Molecular approaches to determine the multiplicity of Plasmodium infections

From multiplicity of infection to force of infection for sparsely sampled Plasmodium falciparum populations at high transmission

High multiplicity of infection or MOI, the number of genetically distinct parasite strains co-infecting a single human host, characterizes infectious diseases including falciparum malaria at high transmission. It accompanies high asymptomatic Plasmodium falciparum prevalence despite high exposure, creating a large transmission reservoir challenging intervention. High MOI and asymptomatic prevalence are enabled by immune evasion of the parasite achieved via vast antigenic diversity. Force of infection or FOI, the number of new infections acquired by an individual host over a given time interval, is the dynamic sister quantity of MOI, and a key epidemiological parameter for monitoring the impact of antimalarial interventions and assessing vaccine or drug efficacy in clinical trials. FOI remains difficult, expensive, and labor-intensive to accurately measure, especially in high-transmission regions, whether directly via cohort studies or indirectly via the fitting of epidemiological models to repeated cross-sectional surveys. We propose here the application of queuing theory to obtain FOI on the basis of MOI, in the form of either a two-moment approximation method or Little's law. We illustrate these methods with MOI estimates obtained under sparse sampling schemes with the recently proposed " v a r coding" method, based on sequences of the v a r multigene family encoding for the major variant surface antigen of the blood stage of malaria infection. The methods are evaluated with simulation output from a stochastic agent-based model, and are applied to an interrupted time-series study from Bongo District in northern Ghana before and immediately after a three-round transient indoor residual spraying (IRS) intervention. We incorporate into the sampling of the simulation output, limitations representative of those encountered in the collection of field data, including under-sampling of v a r genes, missing data, and usage of antimalarial drug treatment. We address these limitations in MOI estimates with a Bayesian framework and an imputation bootstrap approach. We demonstrate that both proposed methods give good and consistent FOI estimates across various simulated scenarios. Their application to the field surveys shows a pronounced reduction in annual FOI during intervention, of more than 70%. The proposed approach should be applicable to the many geographical locations where cohort or cross-sectional studies with regular and frequent sampling are lacking but single-time-point surveys under sparse sampling schemes are available, and for MOI estimates obtained in different ways. They should also be relevant to other pathogens of humans, wildlife and livestock whose immune evasion strategies are based on large antigenic variation resulting in high multiplicity of infection.

Plasmodium falciparum genetic diversity and multiplicity of infection based on msp-1, msp-2, glurp and microsatellite genetic markers in sub-Saharan Africa: a systematic review and meta-analysis

BACKGROUND

In sub-Saharan Africa (SSA), Plasmodium falciparum causes most of the malaria cases. Despite its crucial roles in disease severity and drug resistance, comprehensive data on Plasmodium falciparum genetic diversity and multiplicity of infection (MOI) are sparse in SSA. This study summarizes available information on genetic diversity and MOI, focusing on key markers (msp-1, msp-2, glurp, and microsatellites). The systematic review aimed to evaluate their influence on malaria transmission dynamics and offer insights for enhancing malaria control measures in SSA.

METHODS

The review was conducted following the Preferred Reporting Items for Systematic Review and Meta-Analysis (PRISMA) guidelines. Two reviewers conducted article screening, assessed the risk of bias (RoB), and performed data abstraction. Meta-analysis was performed using the random-effects model in STATA version 17.

RESULTS

The review included 52 articles: 39 cross-sectional studies and 13 Randomized Controlled Trial (RCT)/cohort studies, involving 11,640 genotyped parasite isolates from 23 SSA countries. The overall pooled mean expected heterozygosity was 0.65 (95% CI: 0.51-0.78). Regionally, values varied: East (0.58), Central (0.84), Southern (0.74), and West Africa (0.69). Overall pooled allele frequencies of msp-1 alleles K1, MAD20, and RO33 were 61%, 44%, and 40%, respectively, while msp-2 I/C 3D7 and FC27 alleles were 61% and 55%. Central Africa reported higher frequencies (K1: 74%, MAD20: 51%, RO33: 48%) than East Africa (K1: 46%, MAD20: 42%, RO33: 31%). For msp-2, East Africa had 60% and 55% for I/C 3D7 and FC27 alleles, while West Africa had 62% and 50%, respectively. The pooled allele frequency for glurp was 66%. The overall pooled mean MOI was 2.09 (95% CI: 1.88-2.30), with regional variations: East (2.05), Central (2.37), Southern (2.16), and West Africa (1.96). The overall prevalence of polyclonal Plasmodium falciparum infections was 63% (95% CI: 56-70), with regional prevalences as follows: East (62%), West (61%), Central (65%), and South Africa (71%).

CONCLUSION

The study shows substantial regional variation in Plasmodium falciparum parasite genetic diversity and MOI in SSA. These findings suggest a need for malaria control strategies and surveillance efforts considering regional-specific factors underlying Plasmodium falciparum infection.

Genotyping Plasmodium falciparum gametocytes using amplicon deep sequencing

BACKGROUND

Understanding the dynamics of gametocyte production in polyclonal Plasmodium falciparum infections requires a genotyping method that detects distinct gametocyte clones and estimates their relative frequencies. Here, a marker was identified and evaluated to genotype P. falciparum mature gametocytes using amplicon deep sequencing.

METHODS

A data set of polymorphic regions of the P. falciparum genome was mined to identify a gametocyte genotyping marker. To assess marker resolution, the number of unique haplotypes in the marker region was estimated from 95 Malawian P. falciparum whole genome sequences. Specificity of the marker for detection of mature gametocytes was evaluated using reverse transcription-polymerase chain reaction of RNA extracted from NF54 mature gametocytes and rings from a non-gametocyte-producing strain of P. falciparum. Amplicon deep sequencing was performed on experimental mixtures of mature gametocytes from two distinct parasite clones, as well as gametocyte-positive P. falciparum field isolates to evaluate the quantitative ability and determine the limit of detection of the genotyping approach.

RESULTS

A 400 bp region of the pfs230 gene was identified as a gametocyte genotyping marker. A larger number of unique haplotypes was observed at the pfs230 marker (34) compared to the sera-2 (18) and ama-1 (14) markers in field isolates from Malawi. RNA and DNA genotyping accurately estimated gametocyte and total parasite clone frequencies when evaluating agreement between expected and observed haplotype frequencies in gametocyte mixtures, with concordance correlation coefficients of 0.97 [95% CI: 0.92-0.99] and 0.92 [95% CI: 0.83-0.97], respectively. The detection limit of the genotyping method for male gametocytes was 0.41 pfmget transcripts/µl [95% CI: 0.28-0.72] and for female gametocytes was 1.98 ccp4 transcripts/µl [95% CI: 1.35-3.68].

CONCLUSIONS

A region of the pfs230 gene was identified as a marker to genotype P. falciparum gametocytes. Amplicon deep sequencing of this marker can be used to estimate the number and relative frequency of parasite clones among mature gametocytes within P. falciparum infections. This gametocyte genotyping marker will be an important tool for studies aimed at understanding dynamics of gametocyte production in polyclonal P. falciparum infections.

Molecular speciation of Plasmodium and multiplicity of P. falciparum infection in the Central region of Ghana

Malaria is endemic in the Central region of Ghana, however, the ecological and the seasonal variations of Plasmodium population structure and the intensity of malaria transmission in multiple sites in the region have not been explored. In this cross-sectional study, five districts in the region were involved. The districts were Agona Swedru, Assin Central and Gomoa East (representing the forest zone) and Abura-Asebu-Kwamankese and Cape Coast representing the coastal zone. Systematically, blood samples were collected from patients with malaria. The malaria status was screened with a rapid diagnostic test (RDT) kit (CareStart manufactured by Access Bio in Somerset, USA) and the positive ones confirmed microscopically. Approximately, 200 μL of blood was used to prepare four dried blood spots of 50μL from each microscopy positive sample. The Plasmodium genome was sequenced at the Malaria Genome Laboratory (MGL) of Wellcome Sanger Institute (WSI), Hinxton, UK. The single nucleotide polymorphisms (SNPs) in the parasite mitochondria (PfMIT:270) core genome aided the species identification of Plasmodium. Subsequently, the complexity of infection (COI) was determined using the complexity of infection likelihood (COIL) computational analysis. In all, 566 microscopy positive samples were sequenced. Of this number, Plasmodium genome was detected in 522 (92.2%). However, whole genome sequencing was successful in 409/522 (72.3%) samples. In total, 516/522 (98.8%) of the samples contained P. falciparum mono-infection while the rest (1.2%) were either P. falciparum/P. ovale (Pf/Po) (n = 4, 0.8%) or P. falciparum/P. malariae/P. vivax (Pf/Pm/Pv) mixed-infection (n = 2, 0.4%). All the four Pf/Po infections were identified in samples from the Assin Central municipality whilst the two Pf/Pm/Pv triple infections were identified in Abura-Asebu-Kwamankese district and Cape Coast metropolis. Analysis of the 409 successfully sequenced genome yielded between 1-6 P. falciparum clones per individual infection. The overall mean COI was 1.78±0.92 (95% CI: 1.55-2.00). Among the study districts, the differences in the mean COI between ecological zones (p = 0.0681) and seasons (p = 0.8034) were not significant. However, regression analysis indicated that the transmission of malaria was more than twice among study participants aged 15-19 years (OR = 2.16, p = 0.017) and almost twice among participants aged over 60 years (OR = 1.91, p = 0.021) compared to participants between 20-59 years. Between genders, mean COI was similar except in Gomoa East where females recorded higher values. In conclusion, the study reported, for the first time, P. vivax in Ghana. Additionally, intense malaria transmission was found to be higher in the 15-19 and > 60 years, compared to other age groups. Therefore, active surveillance for P. vivax in Ghana and enhanced malaria control measures in the 15-19 year group years and those over 60 years are recommended.

Sickle cell allele HBB-rs334(T) is associated with decreased risk of childhood Burkitt lymphoma in East Africa

Burkitt lymphoma (BL) is an aggressive B-cell lymphoma that significantly contributes to childhood cancer burden in sub-Saharan Africa. Plasmodium falciparum, which causes malaria, is geographically associated with BL, but the evidence remains insufficient for causal inference. Inference could be strengthened by demonstrating that mendelian genes known to protect against malaria-such as the sickle cell trait variant, HBB-rs334(T)-also protect against BL. We investigated this hypothesis among 800 BL cases and 3845 controls in four East African countries using genome-scan data to detect polymorphisms in 22 genes known to affect malaria risk. We fit generalized linear mixed models to estimate odds ratios (OR) and 95% confidence intervals (95% CI), controlling for age, sex, country, and ancestry. The ORs of the loci with BL and P. falciparum infection among controls were correlated (Spearman's ρ = 0.37, p = .039). HBB-rs334(T) was associated with lower P. falciparum infection risk among controls (OR = 0.752, 95% CI 0.628-0.9; p = .00189) and BL risk (OR = 0.687, 95% CI 0.533-0.885; p = .0037). ABO-rs8176703(T) was associated with decreased risk of BL (OR = 0.591, 95% CI 0.379-0.992; p = .00271), but not of P. falciparum infection. Our results increase support for the etiological correlation between P. falciparum and BL risk.

Amplicon deep sequencing of five highly polymorphic markers of Plasmodium falciparum reveals high parasite genetic diversity and moderate population structure in Ethiopia

BACKGROUND

Plasmodium falciparum genetic diversity can add information on transmission intensity and can be used to track control and elimination interventions.

METHODS

Dried blood spots (DBS) were collected from patients who were recruited for a P. falciparum malaria therapeutic efficacy trial in three malaria endemic sites in Ethiopia from October to December 2015, and November to December 2019. qPCR-confirmed infections were subject to amplicon sequencing of polymorphic markers ama1-D3, csp, cpp, cpmp, msp7. Genetic diversity, the proportion of multiclonal infections, multiplicity of infection, and population structure were analysed.

RESULTS

Among 198 samples selected for sequencing, data was obtained for 181 samples. Mean MOI was 1.38 (95% CI 1.24-1.53) and 17% (31/181) of infections were polyclonal. Mean He across all markers was 0.730. Population structure was moderate; populations from Metema and Metehara 2015 were very similar to each other, but distinct from Wondogent 2015 and Metehara 2019.

CONCLUSION

The high level of parasite genetic diversity and moderate population structure in this study suggests frequent gene flow of parasites among sites. The results obtained can be used as a baseline for additional parasite genetic diversity and structure studies, aiding in the formulation of appropriate control strategies in Ethiopia.

A genomic platform for surveillance and antigen discovery in Plasmodium spp. using long-read amplicon sequencing

Many vaccine candidate proteins in the malaria parasite Plasmodium falciparum are under strong immunological pressure and confer antigenic diversity. We present a sequencing and data analysis platform for the genomic surveillance of the insertion or deletion (indel)-rich antigens merozoite surface protein 1 (MSP1), MSP2, glutamate-rich protein (GLURP), and CSP from P. falciparum using long-read circular consensus sequencing (CCS) in multiclonal malaria isolates. Our platform uses 40 PCR primers per gene to asymmetrically barcode and identify multiclonal infections in pools of up to 384 samples. With msp2, we validated the method using 235 mock infections combining 10 synthetic variants at different concentrations and infection complexities. We applied this strategy to P. falciparum isolates from a longitudinal cohort in Tanzania. Finally, we constructed an analysis pipeline that streamlines the processing and interpretation of epidemiological and antigenic diversity data from demultiplexed FASTQ files. This platform can be easily adapted to other polymorphic antigens of interest in Plasmodium or any other human pathogen.

Complexity of Plasmodium falciparum infection and genetic variations associated with differences in parasite clearance time in two Malian villages

Background

Effective approaches to fight against malaria include disease prevention, an early diagnosis of malaria cases, and rapid management of confirmed cases by treatment with effective antimalarials. Artemisinin-based combination therapies are first-line treatments for uncomplicated malaria in endemic areas. However, cases of resistance to artemisinin have already been described in South-East Asia resulting in prolonged parasite clearance time after treatment. In Mali, though mutations in the K13 gene associated with delayed clearance in Asia are absent, a significant difference in parasite clearance time following treatment with artesunate was observed between two malaria endemic sites, Bougoula-Hameau and Faladje. Hypothetically, differences in complexity of Plasmodium falciparum infections may be accounted for this difference. Hence, the aims of this study were to assess the complexity of infection (COI) and genetic diversity of P. falciparum parasites during malaria treatment in Bougoula-Hameau and Faladje in Mali.

Methods

Thirty (30) patients per village were randomly selected from 221 patients enrolled in a prospective artesunate monotherapy study conducted in Faladje and Bougoula-Hameau in 2016. All parasitemic blood samples of patients from enrollment to last positive slide were retained to assess malaria parasite COI and polymorphisms. DNA were extracted with a Qiagen kit and Pfcsp and Pfama1 encoding gene were amplified by nested PCR and sequenced using the Illumina platform. The parasite clearance time (PCT) was determined using the parasite clearance estimator of Worldwide Antimarial Resistance Network (WWARN). Data were analyzed with R®.

Results

The median number of genetically distinct parasite clones was similar at enrollment, 7 (IQR of 5-9) in Faladje and 6 (IQR of 4-10) in Bougoula-Hameau (p-value = 0.1). On the first day after treatment initiation, the COI was higher in Faladje (6; CI:4-8) than in Bougoula-Hameau (4; CI:4-6) with a p-value =0. 02. Overall, COI was high with higher PCT. Finally, there was a low genetic diversity between Faladje and Bougoula-Hameau.

Conclusion

This study demonstrated that the difference in PCT observed between the two villages could be due to differences in the complexity of infection of these two villages.

Amplicon Deep Sequencing Reveals Multiple Genetic Events Lead to Treatment Failure with Atovaquone-Proguanil in Plasmodium falciparum

Atovaquone-proguanil (AP) is used as treatment for uncomplicated malaria, and as a chemoprophylactic agent against Plasmodium falciparum. Imported malaria remains one of the top causes of fever in Canadian returning travelers. Twelve sequential whole-blood samples before and after AP treatment failure were obtained from a patient diagnosed with P. falciparum malaria upon their return from Uganda and Sudan. Ultradeep sequencing was performed on the cytb, dhfr, and dhps markers of treatment resistance before and during the episode of recrudescence. Haplotyping profiles were generated using three different approaches: msp2-3D7 agarose and capillary electrophoresis, and cpmp using amplicon deep sequencing (ADS). A complexity of infection (COI) analysis was conducted. De novo cytb Y268C mutants strains were observed during an episode of recrudescence 17 days and 16 h after the initial malaria diagnosis and AP treatment initiation. No Y268C mutant reads were observed in any of the samples prior to the recrudescence. SNPs in the dhfr and dhps genes were observed upon initial presentation. The haplotyping profiles suggest multiple clones mutating under AP selection pressure (COI > 3). Significant differences in COI were observed by capillary electrophoresis and ADS compared to the agarose gel results. ADS using cpmp revealed the lowest haplotype variation across the longitudinal analysis. Our findings highlight the value of ultra-deep sequencing methods in the understanding of P. falciparum haplotype infection dynamics. Longitudinal samples should be analyzed in genotyping studies to increase the analytical sensitivity.

coiaf: Directly estimating complexity of infection with allele frequencies

In malaria, individuals are often infected with different parasite strains. The complexity of infection (COI) is defined as the number of genetically distinct parasite strains in an individual. Changes in the mean COI in a population have been shown to be informative of changes in transmission intensity with a number of probabilistic likelihood and Bayesian models now developed to estimate the COI. However, rapid, direct measures based on heterozygosity or FwS do not properly represent the COI. In this work, we present two new methods that use easily calculated measures to directly estimate the COI from allele frequency data. Using a simulation framework, we show that our methods are computationally efficient and comparably accurate to current approaches in the literature. Through a sensitivity analysis, we characterize how the distribution of parasite densities, the assumed sequencing depth, and the number of sampled loci impact the bias and accuracy of our two methods. Using our developed methods, we further estimate the COI globally from Plasmodium falciparum sequencing data and compare the results against the literature. We show significant differences in the estimated COI globally between continents and a weak relationship between malaria prevalence and COI.

Genetic diversity of Plasmodium vivax populations from the China-Myanmar border identified by genotyping merozoite surface protein markers

BACKGROUND

Parasite diversity and population structure influence malaria control measures. Malaria transmission at international borders affects indigenous residents and migrants, defying management efforts and resulting in malaria re-introduction. Here we aimed to determine the extent and distribution of genetic variations in Plasmodium vivax populations and the complexity of infections along the China-Myanmar border.

METHODS

We collected clinical P. vivax samples from local and migrant malaria patients from Laiza and Myitsone, Kachin State, Myanmar, respectively. We characterized the polymorphisms in two P. vivax merozoite surface protein markers, Pvmsp-3α and Pvmsp-3β, by PCR-restriction fragment length polymorphism (PCR-RFLP) analysis. We sought to determine whether these genetic markers could differentiate these two neighboring parasite populations.

RESULTS

PCR revealed three major size variants for Pvmsp-3α and four for Pvmsp-3β among the 370 and 378 samples, respectively. PCR-RFLP resolved 26 fragment-size alleles by digesting Pvmsp-3α with Alu I and Hha I and 28 alleles by digesting Pvmsp-3β with Pst I. PCR-RFLP analysis of Pvmsp-3α found that infections in migrant laborers from Myitsone bore more alleles than did infections in residents of Laiza, while such difference was not evident from genotyping Pvmsp-3β. Infections originating from these two places contained distinct but overlapping subpopulations of P. vivax. Infections from Myitsone had a higher multiplicity of infection as judged by the size of the Pvmsp-3α amplicons and alleles after Alu I/Hha I digestion.

CONCLUSIONS

Migrant laborers from Myitsone and indigenous residents from Laiza harbored overlapping but genetically distinct P. vivax parasite populations. The results suggested a more diverse P. vivax population in Myitsone than in the border town of Laiza. PCR-RFLP of Pvmsp-3α offers a convenient method to determine the complexity of P. vivax infections and differentiate parasite populations.

Low genetic diversity of Plasmodium falciparum merozoite surface protein 1 and 2 and multiplicity of infections in western Ethiopia following effective malaria interventions

BACKGROUND

Genetic diversity of malaria parasites can inform the intensity of transmission and poses a major threat to malaria control and elimination interventions. Characterization of the genetic diversity would provide essential information about the ongoing control efforts. This study aimed to explore allelic polymorphism of merozoite surface protein 1 (msp1) and merozoite surface protein 2 (msp2) to determine the genetic diversity and multiplicity of Plasmodium falciparum infections circulating in high and low transmission sites in western Ethiopia.

METHODS

Parasite genomic DNA was extracted from a total of 225 dried blood spots collected from confirmed uncomplicated P. falciparum malaria-infected patients in western Ethiopia. Of these, 72.4% (163/225) and 27.6% (62/225) of the samples were collected in high and low transmission areas, respectively. Polymorphic msp1 and msp2 genes were used to explore the genetic diversity and multiplicity of falciparum malaria infections. Genotyping of msp1 was successful in 86.5% (141/163) and 88.7% (55/62) samples collected from high and low transmission areas, respectively. Genotyping of msp2 was carried out among 85.3% (139/163) and 96.8% (60/62) of the samples collected in high and low transmission sites, respectively. Plasmodium falciparum msp1 and msp2 genes were amplified by nested PCR and the PCR products were analysed by QIAxcel ScreenGel Software. A P-value of less or equal to 0.05 was considered significant.

RESULTS

High prevalence of falciparum malaria was identified in children less than 15 years as compared with those  ≥ 15 years old (AOR = 2.438, P = 0.005). The three allelic families of msp1 (K1, MAD20, and RO33) and the two allelic families of msp2 (FC27 and 3D7), were observed in samples collected in high and low transmission areas. However, MAD 20 and FC 27 alleles were the predominant allelic families in both settings. Plasmodium falciparum isolates circulating in western Ethiopia had low genetic diversity and mean MOI. No difference in mean MOI between high transmission sites (mean MOI 1.104) compared with low transmission area (mean MOI 1.08) (p > 0.05). The expected heterozygosity of msp1 was slightly higher in isolates collected from high transmission sites (He = 0.17) than in those isolates from low transmission (He = 0.12). However, the heterozygosity of msp2 was not different in both settings (Pfmsp2: 0.04 in high transmission; pfmsp2: 0.03 in low transmission).

CONCLUSION

Plasmodium falciparum from clinical malaria cases in western Ethiopia has low genetic diversity and multiplicity of infection irrespective of the intensity of transmission at the site of sampling. These may be signaling the effectiveness of malaria control strategies in Ethiopia; although further studies are required to determine how specific intervention strategies and other parameters that drive the pattern.

Plasmodium vivax: the potential obstacles it presents to malaria elimination and eradication

Initiatives to eradicate malaria have a good impact on P. falciparum malaria worldwide. P. vivax, however, still presents significant difficulties. This is due to its unique biological traits, which, in comparison to P. falciparum, pose serious challenges for malaria elimination approaches. P. vivax's numerous distinctive characteristics and its ability to live for weeks to years in liver cells in its hypnozoite form, which may elude the human immune system and blood-stage therapy and offer protection during mosquito-free seasons. Many malaria patients are not fully treated because of contraindications to primaquine use in pregnant and nursing women and are still vulnerable to P. vivax relapses, although there are medications that could radical cure P. vivax. Additionally, due to CYP2D6's highly variable genetic polymorphism, the pharmacokinetics of primaquine may be impacted. Due to their inability to metabolize PQ, some CYP2D6 polymorphism alleles can cause patients to not respond to treatment. Tafenoquine offers a radical treatment in a single dose that overcomes the potentially serious problem of poor adherence to daily primaquine. Despite this benefit, hemolysis of the early erythrocytes continues in individuals with G6PD deficiency until all susceptible cells have been eliminated. Field techniques such as microscopy or rapid diagnostic tests (RDTs) miss the large number of submicroscopic and/or asymptomatic infections brought on by reticulocyte tropism and the low parasitemia levels that accompany it. Moreover, P. vivax gametocytes grow more quickly and are much more prevalent in the bloodstream. P. vivax populations also have a great deal of genetic variation throughout their genome, which ensures evolutionary fitness and boosts adaptation potential. Furthermore, P. vivax fully develops in the mosquito faster than P. falciparum. These characteristics contribute to parasite reservoirs in the human population and facilitate faster transmission. Overall, no genuine chance of eradication is predicted in the next few years unless new tools for lowering malaria transmission are developed (i.e., malaria elimination and eradication). The challenging characteristics of P. vivax that impede the elimination and eradication of malaria are thus discussed in this article.

A maximum-likelihood method to estimate haplotype frequencies and prevalence alongside multiplicity of infection from SNP data

The introduction of genomic methods facilitated standardized molecular disease surveillance. For instance, SNP barcodes in Plasmodium vivax and Plasmodium falciparum malaria allows the characterization of haplotypes, their frequencies and prevalence to reveal temporal and spatial transmission patterns. A confounding factor is the presence of multiple genetically distinct pathogen variants within the same infection, known as multiplicity of infection (MOI). Disregarding ambiguous information, as usually done in ad-hoc approaches, leads to less confident and biased estimates. We introduce a statistical framework to obtain maximum-likelihood estimates (MLE) of haplotype frequencies and prevalence alongside MOI from malaria SNP data, i.e., multiple biallelic marker loci. The number of model parameters increases geometrically with the number of genetic markers considered and no closed-form solution exists for the MLE. Therefore, the MLE needs to be derived numerically. We use the Expectation-Maximization (EM) algorithm to derive the maximum-likelihood estimates, an efficient and easy-to-implement algorithm that yields a numerically stable solution. We also derive expressions for haplotype prevalence based on either all or just the unambiguous genetic information and compare both approaches. The latter corresponds to a biased ad-hoc estimate of prevalence. We assess the performance of our estimator by systematic numerical simulations assuming realistic sample sizes and various scenarios of transmission intensity. For reasonable sample sizes, and number of loci, the method has little bias. As an example, we apply the method to a dataset from Cameroon on sulfadoxine-pyrimethamine resistance in P. falciparum malaria. The method is not confined to malaria and can be applied to any infectious disease with similar transmission behavior. An easy-to-use implementation of the method as an R-script is provided.

Temporal dynamics of Plasmodium falciparum population in Metehara, east-central Ethiopia

Background

Plasmodium falciparum is the most serious, genetically most complex and fastest-evolving malaria parasite. Information on genetic diversity of this parasite would guide policy decision and malaria elimination endeavors. This study explored the temporal dynamics of P. falciparum population in two time points in Metehara, east-central Ethiopia.

Methods

The participants were quantitative real-time polymerase chain reaction-confirmed patients who were recruited for uncomplicated falciparum malaria therapeutic efficacy test in 2015 and 2019. Dry blood spot samples were analysed by the nested PCR to genotype P. falciparum merozoite surface protein (msp1, msp2) and glutamate-rich protein (glurp) genes.

Results

While msp1, msp2 and glurp genotypes were successfully detected in 26(89.7%), 24(82.8%) and 14(48.3%) of 2015 samples (n = 29); the respective figures for 2019 (n = 41) were 31(68.3%), 39(95.1%), 25(61.0%). In 2015, the frequencies of K1, MAD20 and RO33 allelic families of msp1, and FC27 and IC/3D7 of msp2 were 19(73.1%), 8(30.6%), 14(53.8%), 21(87.5%), 12(50.5%); and in 2019 it was 15(48.4%), 19(61.3%), 15(48.4%), 30(76.9%), 27(69.2%) respectively. MAD20 has shown dominance over both K1 and RO33 in 2019 compared to the proportion in 2015. Similarly, although FC27 remained dominant, there was shifting trend in the frequency of IC/3D7 from 50.5% in 2015 to 69.2% in 2019. The multiplicity of infection (MOI) and expected heterozygosity index (He) in 2015 and 2019 were respectively [1.43 ± 0.84] and [1.15 ± 0.91], 0.3 and 0.03 for msp1. However, there was no significant association between MOI and age or parasitaemia in both time points.

Conclusion

The lower genetic diversity in P. falciparum population in the two time points and overall declining trend as demonstrated by the lower MOI and He may suggest better progress in malaria control in Metehara. But, the driving force and selective advantage of switching to MAD20 dominance over the other two msp1 allelic families, and the dynamics within msp2 alleles needs further investigation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04277-5.

Population genomics in neglected malaria parasites

Malaria elimination includes neglected human malaria parasites Plasmodium vivax, Plasmodium ovale spp., and Plasmodium malariae. Biological features such as association with low-density infection and the formation of hypnozoites responsible for relapse make their elimination challenging. Studies on these parasites rely primarily on clinical samples due to the lack of long-term culture techniques. With improved methods to enrich parasite DNA from clinical samples, whole-genome sequencing of the neglected malaria parasites has gained increasing popularity. Population genomics of more than 2200 P. vivax global isolates has improved our knowledge of parasite biology and host-parasite interactions, identified vaccine targets and potential drug resistance markers, and provided a new way to track parasite migration and introduction and monitor the evolutionary response of local populations to elimination efforts. Here, we review advances in population genomics for neglected malaria parasites, discuss how the rich genomic information is being used to understand parasite biology and epidemiology, and explore opportunities for the applications of malaria genomic data in malaria elimination practice.

Genetic polymorphism of Plasmodium falciparum msp-1, msp-2 and glurp vaccine candidate genes in pre-artemisinin era clinical isolates from Lakhimpur district in Assam, Northeast India

Background

Northeast India shares its international border with Southeast Asia and has a number of malaria endemic zones. Monitoring genetic diversity of malaria parasites is important in this area as drug resistance and increasing genetic diversity form a vicious cycle in which one favours the development of the other. This retrospective study was done to evaluate the genetic diversity patterns in Plasmodium falciparum strains circulating in North Lakhimpur area of Assam in the pre-artemisinin era and compare the findings with current diversity patterns.

Methods

Genomic DNA extraction was done from archived blood spot samples collected in 2006 from malaria-positive cases in Lakhimpur district of Assam, Northeast India. Three antigenic markers of genetic diversity were studied – msp-1 (block-2), msp-2 (block-3) and the glurp RII region of P. falciparum using nested PCR.

Results

Allelic diversity was examined in 71 isolates and high polymorphism was observed. In msp-1, eight genotypes were detected; K1 (single allele), MAD20 (six different alleles) and RO33 (single allele) allelic families were noted. Among msp-2 genotypes, 22 distinct alleles were observed out of which FC27 had six alleles and IC/3D7 had 16 alleles. In RII region of glurp, nine genotypes were obtained. Expected heterozygosity (H_E) values of the three antigenic markers were 0.72, 0.81 and 0.88, respectively. Multiplicity of infection (MOI) values noted were 1.28, 1.84 and 1.04 for msp-1, msp-2 and glurp, respectively.

Conclusion

Results suggest a high level of genetic diversity in P. falciparum msp (block-2 of msp-1 and block-3 of msp-2) and the glurp RII region in Northeast India in the pre-artemisinin era when chloroqunine was the primary drug used for uncomplicated falciparum malaria. Comparison with current studies have revealed that the genetic diversity in these genes is still high in this region, complicating malaria vaccine research.

Infection age as a predictor of epidemiological metrics for malaria

BACKGROUND

Accurate estimation of the burden of Plasmodium falciparum is essential for strategic planning for control and elimination. Due in part to the extreme heterogeneity in malaria exposure, immunity, other causes of disease, direct measurements of fever and disease attributable to malaria can be difficult. This can make a comparison of epidemiological metrics both within and between populations hard to interpret. An essential part of untangling this is an understanding of the complex time-course of malaria infections.

METHODS

Historic data from malariatherapy infections, in which individuals were intentionally infected with malaria parasites, were reexamined in aggregate. In this analysis, the age of each infection was examined as a potential predictor describing aggregate patterns across all infections. A series of piecewise linear and generalized linear regressions were performed to highlight the infection age-dependent patterns in both parasitaemia and gametocytaemia, and from parasitaemia and gametocytaemia to fever and transmission probabilities, respectively.

RESULTS

The observed duration of untreated patent infection was 130 days. As infections progressed, the fraction of infections subpatent by microscopy was seen to increase steadily. The time-averaged malaria infections had three distinct phases in parasitaemia: a growth phase for the first 6 days of patency, a rapid decline from day 6 to day 18, and a slowly declining chronic phase for the remaining duration of the infection. During the growth phase, parasite densities increased sharply to a peak. Densities sharply decline for a short period of time after the peak. During the chronic phase, infections declined steadily as infections age. gametocytaemia was strongly correlated with lagged asexual parasitaemia. Fever rates and transmission efficiency were strongly correlated with parasitaemia and gametocytaemia. The comparison between raw data and prediction from the age of infection has good qualitative agreement across all quantities of interest for predicting averaged effects.

CONCLUSION

The age of infection was established as a potentially useful covariate for malaria epidemiology. Infection age can be estimated given a history of exposure, and accounting for exposure history may potentially provide a new way to estimate malaria-attributable fever rates, transmission efficiency, and patent fraction in immunologically naïve individuals such as children and people in low-transmission regions. These data were collected from American adults with neurosyphilis, so there are reasons to be cautious about extending the quantitative results reported here to general populations in malaria-endemic regions. Understanding how immune responses modify these statistical relationships given past exposure is key for being able to apply these results more broadly.

Comparative analysis of targeted next-generation sequencing for Plasmodium falciparum drug resistance markers

Well-defined molecular resistance markers are available for a range of antimalarial drugs, and molecular surveillance is increasingly important for monitoring antimalarial drug resistance. Different genotyping platforms are available, but these have not been compared in detail. We compared Targeted Amplicon Deep sequencing (TADs) using Ion Torrent PGM with Illumina MiSeq for the typing of antimalarial drug resistance genes. We developed and validated protocols to type the molecular resistance markers pfcrt, pfdhfr, pfdhps, pfmdr1, pfkelch, and pfcytochrome b, in Plasmodium falciparum for the Ion Torrent PGM and Illumina MiSeq sequencing platforms. With P. falciparum 3D7 and K1 as reference strains, whole blood samples (N = 20) and blood spots from Rapid Diagnostic Test (RDT) samples (N = 5) from patients with uncomplicated falciparum malaria from Ubon Ratchathani were assessed on both platforms and compared for coverage (average reads per amplicon), sequencing accuracy, variant accuracy, false positive rate, false negative rate, and alternative allele detection, with conventional Sanger sequencing as the reference method for SNP calling. Both whole blood and RDT samples could be successfully sequenced using the Ion Torrent PGM and Illumina MiSeq platforms. Coverage of reads per amplicon was higher with Illumina MiSeq (28,886 reads) than with Ion Torrent PGM (1754 reads). In laboratory generated artificial mixed infections, the two platforms could detect the minor allele down to 1% density at 500X coverage. SNPs calls from both platforms were in complete agreement with conventional Sanger sequencing. The methods can be multiplexed with up to 96 samples per run, which reduces cost by 86% compared to conventional Sanger sequencing. Both platforms, using the developed TAD protocols, provide an accurate method for molecular surveillance of drug resistance markers in P. falciparum, but Illumina MiSeq provides higher coverage than Ion Torrent PGM.

Genotyping of Plasmodium falciparum to Assess Clone Composition in Parasite Cultures

Defining clone composition in Plasmodium falciparum cultures is key to verify that in vitro experiments are performed on the parasite line of interest. Genotyping of the highly polymorphic merozoite surface protein 2 gene (msp2) is a widely established method to define P. falciparum clones. Specific size variants from the two msp2 families (IC and FC27) can be used as "fingerprints" to identify individual clones in parasite mixtures. Size variant genotyping of msp2 using fluorescent nested PCR followed by fragment analysis by capillary electrophoresis (CE) provides accurate information about the presence of one or multiple parasite clones. Here, we describe an adaptation of this approach to assess the integrity and purity of P. falciparum lines kept in in vitro culture. In addition, we describe the use of synthetic mock parasite mixtures with the msp2 sequences from the parasite lines kept in culture that can provide a good estimate of the assay sensitivity, specificity, and reproducibility. We suggest that genotyping of P. falciparum lines should be performed on a regular basis as part of the standard procedures of in vitro parasite culture, as a way to secure that the parasite lines of interest are cultivated, and to monitor any cross-contamination and/or recombination events.

Distribution of clones among hosts for the lizard malaria parasite Plasmodium mexicanum

Background

Malaria parasites reproduce asexually, leading to the production of large numbers of genetically identical parasites, here termed a clonal line or clone. Infected hosts may harbor one or more clones, and the number of clones in a host is termed multiplicity of infection (MOI). Understanding the distribution of parasite clones among hosts can shed light on the processes shaping this distribution and is important for modeling MOI. Here, I determine whether the distribution of clones of the lizard malaria parasite Plasmodium mexicanum differ significantly from statistical distributions commonly used to model MOI and logical extensions of these models.

Methods

The number of clones per infection was assessed using four microsatellite loci with the maximum number of alleles at any one locus used as a simple estimate of MOI for each infection. I fit statistical models (Poisson, negative binomial, zero-inflated models) to data from four individual sites to determine a best fit model. I also simulated the number of alleles per locus using an unbiased estimate of MOI to determine whether the simple (but potentially biased) method I used to estimate MOI influenced model fit.

Results

The distribution of clones among hosts at individual sites differed significantly from traditional Poisson and negative binomial distributions, but not from zero-inflated modifications of these distributions. A consistent excess of two-clone infections and shortage of one-clone infections relative to all fit distributions was also observed. Any bias introduced by the simple method for estimating of MOI did not appear to qualitatively alter the results.

Conclusions

The statistical distributions used to model MOI are typically zero-truncated; truncating the Poisson or zero-inflated Poisson yield the same distribution, so the reasonable fit of the zero-inflated Poisson to the data suggests that the use of the zero-truncated Poisson in modeling is adequate. The improved fit of zero-inflated distributions relative to standard distributions may suggest that only a portion of the host population is located in areas suitable for transmission even at small sites (<1 ha). Collective transmission of clones and premunition may also contribute to deviations from standard distributions.

Markov chain Monte Carlo Gibbs sampler approach for estimating haplotype frequencies among multiple malaria infected human blood samples

Background

Malaria patients can have two or more haplotypes in their blood sample making it challenging to identify which haplotypes they carry. In addition, there are challenges in measuring the type and frequency of resistant haplotypes in populations. This study presents a novel statistical method Gibbs sampler algorithm to investigate this issue.

Results

The performance of the algorithm is evaluated on simulated datasets consisting of patient blood samples characterized by their multiplicity of infection (MOI) and malaria genotype. The simulation used different resistance allele frequencies (RAF) at each Single Nucleotide Polymorphisms (SNPs) and different limit of detection (LoD) of the SNPs and the MOI. The Gibbs sampler algorithm presents higher accuracy among high LoD of the SNPs or the MOI, validated, and deals with missing MOI compared to previous related statistical approaches.

Conclusions

The Gibbs sampler algorithm provided robust results when faced with genotyping errors caused by LoDs and functioned well even in the absence of MOI data on individual patients.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-021-03841-9.

Genetic diversity and genetic relatedness in Plasmodium falciparum parasite population in individuals with uncomplicated malaria based on microsatellite typing in Eastern and Western regions of Uganda, 2019-2020

BACKGROUND

Genetic diversity and parasite relatedness are essential parameters for assessing impact of interventions and understanding transmission dynamics of malaria parasites, however data on its status in Plasmodium falciparum populations in Uganda is limited. Microsatellite markers and DNA sequencing were used to determine diversity and molecular characterization of P. falciparum parasite populations in Uganda.

METHODS

A total of 147 P. falciparum genomic DNA samples collected from cross-sectional surveys in symptomatic individuals of 2-10 years were characterized by genotyping of seven highly polymorphic neutral microsatellite markers (n = 85) and genetic sequencing of the Histidine Rich Protein 2 (pfhrp2) gene (n = 62). ArcGIS was used to map the geographical distribution of isolates while statistical testing was done using Student's t-test or Wilcoxon's rank-sum test and Fisher's exact test as appropriate at P ≤ 0.05.

RESULTS

Overall, 75.5% (95% CI 61.1-85.8) and 24.5% (95% CI14.2-38.9) of parasites examined were of multiclonal (mixed genotype) and single clone infections, respectively. Multiclonal infections occurred more frequently in the Eastern region 73.7% (95% CI 48.8-89.1), P < 0.05. Overall, multiplicity of infection (MOI) was 1.9 (95% CI 1.7-2.1), P = 0.01 that was similar between age groups (1.8 vs 1.9), P = 0.60 and regions (1.9 vs 1.8), P = 0.43 for the < 5 and ≥ 5 years and Eastern and Western regions, respectively. Genomic sequencing of the pfhrp2 exon2 revealed a high level of genetic diversity reflected in 96.8% (60/62) unique sequence types. Repeat type AHHAAAHHATD and HRP2 sequence Type C were more frequent in RDT-/PCR + samples (1.9% vs 1.5%) and (13% vs 8%), P < 0.05 respectively. Genetic relatedness analysis revealed small clusters of gene deleted parasites in Uganda, but no clustering with Eritrean parasites.

CONCLUSION

High level of genetic diversity of P. falciparum parasites reflected in the frequency of multiclonal infections, multiplicity of infection and variability of the pfhrp2 gene observed in this study is consistent with the high malaria transmission intensity in these settings. Parasite genetic analysis suggested spontaneous emergence and clonal expansion of pfhrp2 deleted parasites that require close monitoring to inform national malaria diagnosis and case management policies.

Population structure and diversity of Plasmodium falciparum in children with asymptomatic malaria living in different ecological zones of Ghana

BACKGROUND

Genetic diversity in Plasmodium falciparum populations can be used to describe the resilience and spatial distribution of the parasite in the midst of intensified intervention efforts. This study used microsatellite analysis to evaluate the genetic diversity and population dynamics of P. falciparum parasites circulating in three ecological zones of Ghana.

METHODS

A total of 1168 afebrile children aged between 3 to 13 years were recruited from five (5) Primary schools in 3 different ecological zones (Sahel (Tamale and Kumbungu), Forest (Konongo) and Coastal (Ada and Dodowa)) of Ghana. Asymptomatic malaria parasite carriage was determined using microscopy and PCR, whilst fragment analysis of 6 microsatellite loci was used to determine the diversity and population structure of P. falciparum parasites.

RESULTS

Out of the 1168 samples examined, 16.1 and 39.5% tested positive for P. falciparum by microscopy and nested PCR respectively. The genetic diversity of parasites in the 3 ecological zones was generally high, with an average heterozygosity (He) of 0.804, 0.787 and 0.608 the rainy (peak) season for the Sahel, Forest and Coastal zones respectively. The mean He for the dry (off-peak) season were 0.562, 0.693 and 0.610 for the Sahel, Forest and Coastal zones respectively. Parasites from the Forest zone were more closely related to those from the Sahel than from the Coastal zone, despite the Coastal zone being closer in physical distance to the Forest zone. The fixation indexes among study sites ranged from 0.049 to 0.112 during the rainy season and 0.112 to 0.348 during the dry season.

CONCLUSION

A large asymptomatic parasite reservoir was found in the school children during both rainy and dry seasons, especially those in the Forest and Sahel savannah zones where parasites were also found to be related compared to those from the Coastal zone. Further studies are recommended to understand why despite the roll out of several malaria interventions in Ghana, high transmission still persist.

Development of a new barcode-based, multiplex-PCR, next-generation-sequencing assay and data processing and analytical pipeline for multiplicity of infection detection of Plasmodium falciparum

BACKGROUND

Simultaneous infection with multiple malaria parasite strains is common in high transmission areas. Quantifying the number of strains per host, or the multiplicity of infection (MOI), provides additional parasite indices for assessing transmission levels but it is challenging to measure accurately with current tools. This paper presents new laboratory and analytical methods for estimating the MOI of Plasmodium falciparum.

METHODS

Based on 24 single nucleotide polymorphisms (SNPs) previously identified as stable, unlinked targets across 12 of the 14 chromosomes within P. falciparum genome, three multiplex PCRs of short target regions and subsequent next generation sequencing (NGS) of the amplicons were developed. A bioinformatics pipeline including B4Screening pathway removed spurious amplicons to ensure consistent frequency calls at each SNP location, compiled amplicons by SNP site diversity, and performed algorithmic haplotype and strain reconstruction. The pipeline was validated by 108 samples generated from cultured-laboratory strain mixtures in different proportions and concentrations, with and without pre-amplification, and using whole blood and dried blood spots (DBS). The pipeline was applied to 273 smear-positive samples from surveys conducted in western Kenya, then providing results into StrainRecon Thresholding for Infection Multiplicity (STIM), a novel MOI estimator.

RESULTS

The 24 barcode SNPs were successfully identified uniformly across the 12 chromosomes of P. falciparum in a sample using the pipeline. Pre-amplification and parasite concentration, while non-linearly associated with SNP read depth, did not influence the SNP frequency calls. Based on consistent SNP frequency calls at targeted locations, the algorithmic strain reconstruction for each laboratory-mixed sample had 98.5% accuracy in dominant strains. STIM detected up to 5 strains in field samples from western Kenya and showed declining MOI over time (q < 0.02), from 4.32 strains per infected person in 1996 to 4.01, 3.56 and 3.35 in 2001, 2007 and 2012, and a reduction in the proportion of samples with 5 strains from 57% in 1996 to 18% in 2012.

CONCLUSION

The combined approach of new multiplex PCRs and NGS, the unique bioinformatics pipeline and STIM could identify 24 barcode SNPs of P. falciparum correctly and consistently. The methodology could be applied to field samples to reliably measure temporal changes in MOI.

Prevalence of Asymptomatic Malaria Parasitemia in Odisha, India: A Challenge to Malaria Elimination

The prevalence of malaria in India is decreasing, but it remains a major concern for public health administration. The role of submicroscopic malaria and asymptomatic malaria parasitemia and their persistence is being explored. A cross-sectional survey was conducted in the Kandhamal district of Odisha (India) during May-June 2017. Blood samples were collected from 1897 individuals for screening of asymptomatic parasitemia. Samples were screened using rapid diagnostic tests (RDTs) and examined microscopically for Plasmodium species. Approximately 30% of randomly selected samples (n = 586) were analyzed using real-time PCR (qPCR), and the genetic diversity of Plasmodium falciparum was analyzed. The prevalence of Plasmodium species among asymptomatic individuals detected using qPCR was 18%, which was significantly higher than that detected by microscopy examination (5.5%) or RDT (7.3%). Of these, 37% had submicroscopic malaria. The species-specific prevalence among asymptomatic malaria-positive cases for P. falciparum, Plasmodium vivax, and mixed infection (P. falciparum and P. vivax) by qPCR was 57%, 29%, and 14%, respectively. The multiplicity of infection was 1.6 and 1.2 for the merozoite surface protein-1 gene (msp1) and (msp2), respectively. Expected heterozygosity was 0.64 and 0.47 for msp1 and msp2, respectively. A significant proportion of the study population, 105/586 (18%), was found to be a reservoir for malaria infection, and identification of this group will help in the development of elimination strategies.

Use of a Plasmodium vivax genetic barcode for genomic surveillance and parasite tracking in Sri Lanka

Background

Sri Lanka was certified as a malaria-free nation in 2016; however, imported malaria cases continue to be reported. Evidence-based information on the genetic structure/diversity of the parasite populations is useful to understand the population history, assess the trends in transmission patterns, as well as to predict threatening phenotypes that may be introduced and spread in parasite populations disrupting elimination programmes. This study used a previously developed Plasmodium vivax single nucleotide polymorphism (SNP) barcode to evaluate the population dynamics of P. vivax parasite isolates from Sri Lanka and to assess the ability of the SNP barcode for tracking the parasites to its origin.

Methods

A total of 51 P. vivax samples collected during 2005–2011, mainly from three provinces of the country, were genotyped for 40 previously identified P. vivax SNPs using a high-resolution melting (HRM), single-nucleotide barcode method. Minor allele frequencies, linkage disequilibrium, pair-wise F_ST values, and complexity of infection (COI) were evaluated to determine the genetic diversity. Structure analysis was carried out using STRUCTURE software (Version 2.3.4) and SNP barcode was used to identify the genetic diversity of the local parasite populations collected from different years. Principal component analysis (PCA) was used to determine the clustering according to global geographic regions.

Results

The proportion of multi-clone infections was significantly higher in isolates collected during an infection outbreak in year 2007. The minor allele frequencies of the SNPs changed dramatically from year to year. Significant linkage was observed in sample sub-sets from years 2005 and 2007. The majority of the isolates from 2007 consisted of at least two genetically distinct parasite strains. The overall percentage of multi-clone infections for the entire parasite sample was 39.21%. Analysis using STRUCTURE software (Version 2.3.4) revealed the high genetic diversity of the sample sub-set from year 2007. In-silico analysis of these data with those available from other global geographical regions using PCA showed distinct clustering of parasite isolates according to geography, demonstrating the usefulness of the barcode in determining an isolate to be indigenous.

Conclusions

Plasmodium vivax parasite isolates collected during a disease outbreak in year 2007 were more genetically diverse compared to those collected from other years. In-silico analysis using the 40 SNP barcode is a useful tool to track the origin of an isolate of uncertain origin, especially to differentiate indigenous from imported cases. However, an extended barcode with more SNPs may be needed to distinguish highly clonal populations within the country.

Microgeographic Epidemiology of Malaria Parasites in an Irrigated Area of Western Kenya by Deep Amplicon Sequencing

To improve food security, investments in irrigated agriculture are anticipated to increase throughout Africa. However, the extent to which environmental changes from water resource development will impact malaria epidemiology remains unclear. This study was designed to compare the sensitivity of molecular markers used in deep amplicon sequencing for evaluating malaria transmission intensities and to assess malaria transmission intensity at various proximities to an irrigation scheme. Compared to ama1, csp, and msp1 amplicons, cpmp required the smallest sample size to detect differences in infection complexity between transmission risk zones. Transmission intensity was highest within 5 km of the irrigation scheme by polymerase chain reaction positivity rate, infection complexity, and linkage disequilibrium. The irrigated area provided a source of parasite infections for the surrounding 2- to 10-km area. This study highlights the suitability of the cpmp amplicon as a measure for transmission intensities and the impact of irrigation on microgeographic epidemiology of malaria parasites.

Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region

The Asia-Pacific region faces formidable challenges in achieving malaria elimination by the proposed target in 2030. Molecular surveillance of Plasmodium parasites can provide important information on malaria transmission and adaptation, which can inform national malaria control programmes (NMCPs) in decision-making processes. In November 2019 a parasite genotyping workshop was held in Jakarta, Indonesia, to review molecular approaches for parasite surveillance and explore ways in which these tools can be integrated into public health systems and inform policy. The meeting was attended by 70 participants from 8 malaria-endemic countries and partners of the Asia Pacific Malaria Elimination Network. The participants acknowledged the utility of multiple use cases for parasite genotyping including: quantifying the prevalence of drug resistant parasites, predicting risks of treatment failure, identifying major routes and reservoirs of infection, monitoring imported malaria and its contribution to local transmission, characterizing the origins and dynamics of malaria outbreaks, and estimating the frequency of Plasmodium vivax relapses. However, the priority of each use case varies with different endemic settings. Although a one-size-fits-all approach to molecular surveillance is unlikely to be applicable across the Asia-Pacific region, consensus on the spectrum of added-value activities will help support data sharing across national boundaries. Knowledge exchange is needed to establish local expertise in different laboratory-based methodologies and bioinformatics processes. Collaborative research involving local and international teams will help maximize the impact of analytical outputs on the operational needs of NMCPs. Research is also needed to explore the cost-effectiveness of genetic epidemiology for different use cases to help to leverage funding for wide-scale implementation. Engagement between NMCPs and local researchers will be critical throughout this process.

A hybrid model for the effects of treatment and demography on malaria superinfection

•

Treatment, demography impact the distribution of multiplicity of infection (MOI).

•

The MOI can be modeled with an alternative hyper-Poisson distribution.

•

The distribution of MOI determines the average rate of recovery.

•

The average rate of recovery is used to map between exposure and prevalence.

As standard mathematical models for the transmission of vector-borne pathogens with weak or no apparent sterilizing immunity, Susceptible-Infected-Susceptible (SIS) systems such as the Ross-Macdonald equations are a useful starting point for modeling the impacts of interventions on prevalence for diseases that cannot superinfect their hosts. In particular, they are parameterizable from quantities we can estimate such as the force of infection (FOI), the rate of natural recovery from a single infection, the treatment rate, and the rate of demographic turnover. However, malaria parasites can superinfect their host which has the effect of increasing the duration of infection before total recovery. Queueing theory has been applied to capture this behavior, but a problem with current queueing models is the exclusion of factors such as demographic turnover and treatment. These factors in particular can affect the entire shape of the distribution of the multiplicity of infection (MOI) generated by the superinfection process, its transient dynamics, and the population mean recovery rate. Here we show the distribution of MOI can be described by an alternative hyper-Poisson distribution. We then couple our resulting equations to a simple vector transmission model, extending previous Ross-Macdonald theory.

Genetic diversity of Plasmodium falciparum and genetic profile in children affected by uncomplicated malaria in Cameroon

BACKGROUND

Malaria is a major public health problem in Cameroon. The study of the genetic diversity within parasite population is essential for understanding the mechanism underlying malaria pathology and to determine parasite clones profile in an infection, for proper malaria control strategies. The objective of this study was to perform a molecular characterization of highly polymorphic genetic markers of Plasmodium falciparum, and to determine allelic distribution with their influencing factors valuable to investigate malaria transmission dynamics in Cameroon.

METHODS

A total of 350 P. falciparum clinical isolates were characterized by genotyping block 2 of msp-1, block 3 of msp-2, and region II of glurp gene using nested PCR and DNA sequencing between 2012 and 2013.

RESULTS

A total of 5 different genotypes with fragment sizes ranging from 597 to 817 bp were recorded for GLURP. Overall, 16 MSP-1 genotypes, including K1, MAD20 and RO33 were identified, ranging from 153 to 335 bp. A peculiarity about this study is the RO33 monomorphic pattern revealed among the Pfmsp-1 allelic type. Again, this study identified 27 different Pfmsp-2 genotypes, ranging from 140 to 568 bp in size, including 15 belonging to the 3D7-type and 12 to the FC27 allelic families. The analysis of the MSP-1 and MSP-2 peptides indicates that the region of the alignment corresponding K1 polymorphism had the highest similarity in the MSP1and MSP2 clade followed by MAD20 with 93% to 100% homology. Therefore, population structure of P. falciparum isolates is identical to that of other areas in Africa, suggesting that vaccine developed with K1 and MAD20 of Pfmsp1 allelic variant could be protective for Africa children but these findings requires further genetic and immunological investigations. The multiplicity of infection (MOI) was significantly higher (P < 0.05) for Pfmsp-2 loci (3.82), as compare with Pfmsp-1 (2.51) and heterozygotes ranged from 0.55 for Pfmsp-1 to 0.96 for Pfmsp-2.

CONCLUSION

High genetic diversity and allelic frequencies in P. falciparum isolates indicate a persisting high level of transmission. This study advocate for an intensification of the malaria control strategies in Cameroon. Trial registration This study was approved by Cameroon National Ethics Committee. It is a randomized controlled trial retrospectively registered in NIH U.S. National Library of Medicine, ClinicalTrials.gov on the 28/11/2016 at https://clinicaltrials.gov/ct2/show/NCT02974348 with the registration number NCT02974348.

Malaria Molecular Epidemiology: An Evolutionary Genetics Perspective

Malaria is a vector-borne disease that involves multiple parasite species in a variety of ecological settings. However, the parasite species causing the disease, the prevalence of subclinical infections, the emergence of drug resistance, the scale-up of interventions, and the ecological factors affecting malaria transmission, among others, are aspects that vary across areas where malaria is endemic. Such complexities have propelled the study of parasite genetic diversity patterns in the context of epidemiologic investigations. Importantly, molecular studies indicate that the time and spatial distribution of malaria cases reflect epidemiologic processes that cannot be fully understood without characterizing the evolutionary forces shaping parasite population genetic patterns. Although broad in scope, this review in the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology highlights the need for understanding population genetic concepts when interpreting parasite molecular data. First, we discuss malaria complexity in terms of the parasite species involved. Second, we describe how molecular data are changing our understanding of malaria incidence and infectiousness. Third, we compare different approaches to generate parasite genetic information in the context of epidemiologically relevant questions related to malaria control. Finally, we describe a few Plasmodium genomic studies as evidence of how these approaches will provide new insights into the malaria disease dynamics. *This article is part of a curated collection.

Strengthening Surveillance Systems for Malaria Elimination by Integrating Molecular and Genomic Data

Unprecedented efforts in malaria control over the last 15 years have led to a substantial decrease in both morbidity and mortality in most endemic settings. However, these progresses have stalled over recent years, and resurgence may cause dramatic impact on both morbidity and mortality. Nevertheless, elimination efforts are currently going on with the objective of reducing malaria morbidity and mortality by 90% and malaria elimination in at least 35 countries by 2030. Strengthening surveillance systems is of paramount importance to reach those targets, and the integration of molecular and genomic techniques into routine surveillance could substantially improve the quality and robustness of data. Techniques such as polymerase chain reaction (PCR) and quantitative PCR (qPCR) are increasingly available in malaria endemic countries, whereas others such as sequencing are already available in a few laboratories. However, sequencing, especially next-generation sequencing (NGS), requires sophisticated infrastructure with adequate computing power and highly trained personnel for data analysis that require substantial investment. Different techniques will be required for different applications, and cost-effective planning must ensure the appropriate use of available resources. The development of national and sub-regional reference laboratories could help in minimizing the resources required in terms of equipment and trained staff. Concerted efforts from different stakeholders at national, sub-regional, and global level are needed to develop the required framework to establish and maintain these reference laboratories.

Asymptomatic malaria in the clinical and public health context

Introduction: Historically, the global community has focused on the control of symptomatic malaria. However, interest in asymptomatic malaria has been growing, particularly in the context of malaria elimination.Areas covered: We undertook a comprehensive PubMed literature review on asymptomatic malaria as it relates to detection and elimination with emphasis between 2014 and 2019. Diagnostic tools with a low limit of detection (LOD) have allowed us to develop a more detailed understanding of asymptomatic malaria and its impact. These highly sensitive diagnostics have demonstrated that the prevalence of asymptomatic malaria is greater than previously thought. In addition, it is now possible to detect the malaria reservoir in the community, something that was previously not feasible. Asymptomatic malaria has previously not been treated, but research has begun to examine whether treating individuals with asymptomatic malaria may lead to health benefits. Finally, we have begun to understand the importance of asymptomatic malaria in ongoing transmission.Expert opinion: Therefore, with malaria elimination back on the agenda, asymptomatic malaria can no longer be ignored, especially in light of new ultra-sensitive diagnostic tools.

Characterization of drug resistance and genetic diversity of Plasmodium falciparum parasites from Tripura, Northeast India

Monitoring of anti-malarial drug resistance is vital in Northeast India as this region shares its international border with Southeast Asia. Genetic diversity of Plasmodium parasites regulates transmission dynamics, disease severity and vaccine efficacy. P. falciparum chloroquine resistance transporter (Pfcrt), multidrug resistance-1 (Pfmdr-1) and kelch 13 propeller (PfK-13) genes which govern antimalarial drug resistance and three genetic diversity markers, merozoite surface protein 1 and 2 (Pfmsp-1, Pfmsp-2) and glutamate rich protein (Pfglurp) were evaluated from Tripura, Northeast India using molecular tools. In the Pfcrt gene, 87% isolates showed triple mutations at codons M74I, N75E and K76T. 12.5% isolates in Pfmdr-1 gene showed mutation at N86Y. No polymorphism in PfK-13 propeller was found. Polyclonal infections were observed in 53.85% isolates and more commonly in adults (p = 0.0494). In the Pfmsp-1 locus, the K1 allelic family was predominant (71.2%) followed by the 3D7/IC family (69.2%) in the Pfmsp-2 locus. RII region of Pfglurp exhibited nine alleles with expected heterozygosity of 0.85. The multiplicity of infection for Pfmsp-1, Pfmsp-2 and Pfglurp were 1.56, 1.31 and 1.06 respectively. Overall, the study demonstrated a high level of chloroquine resistance and extensive parasite diversity in the region, necessitating regular surveillance in this population group.

Applying next-generation sequencing to track falciparum malaria in sub-Saharan Africa

Next-generation sequencing (NGS) technologies are increasingly being used to address a diverse range of biological and epidemiological questions. The current understanding of malaria transmission dynamics and parasite movement mainly relies on the analyses of epidemiologic data, e.g. case counts and self-reported travel history data. However, travel history data are often not routinely collected or are incomplete, lacking the necessary level of accuracy. Although genetic data from routinely collected field samples provides an unprecedented opportunity to track the spread of malaria parasites, it remains an underutilized resource for surveillance due to lack of local awareness and capacity, limited access to sensitive laboratory methods and associated computational tools and difficulty in interpreting genetic epidemiology data. In this review, the potential roles of NGS in better understanding of transmission patterns, accurately tracking parasite movement and addressing the emerging challenges of imported malaria in low transmission settings of sub-Saharan Africa are discussed. Furthermore, this review highlights the insights gained from malaria genomic research and challenges associated with integrating malaria genomics into existing surveillance tools to inform control and elimination strategies.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing

BACKGROUND

Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/μl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/μl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated.

RESULTS

Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/μl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/μl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/μl, but only 61% at low Plasmodium concentrations (< 5 genomes/μl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases.

CONCLUSIONS

Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

The preparation and clinical application of diagnostic DNA microarray for the detection of pathogens in intracranial bacterial and fungal infections

The present study prepared 2 types of DNA diagnostic chips based on 16S ribosomal DNA (rDNA) and 18S-28S rDNA, and evaluated their values in the detection of pathogens in intracranial bacterial/fungal infections. A total of 14 probes of bacteria (Klebsiella pneumonia, Acinetobacter baumannii, Pseudomonas aeruginosa, Escherichia coli, Haemophilus influenza, Stenotrophomonas maltophilia, Neisseria meningitidis, Enterobacter spp., Enterococcus faecalis, Enterococcus faecium, Listeria monocytogenes, Staphylococcus aureus, Streptococcus pneumonia and coagulase negative staphylococcus) and 4 probes of fungi (Candida albicans, Candida tropicalis, Candida glabrata and Cryptococcus neoformans), determined frequently in cerebrospinal fluid (CSF), were designed and used for preparation of microarrays. CSF samples from 88 patients with clinically suspected intracranial infection and standard strains were used to evaluate the chips. The same samples were also analyzed by culture and sequencing. The results demonstrated that the sensitivity, specificity and false-positive rate of the microarray assay compared with culture method were 100 vs. 68.3% (P<0.05), 97.1 vs. 100%, and 2.9 vs. 0%, respectively. The minimum concentration of detection with the chips was 10 cfu ml^-1 for bacteria and 100 cfu ml^-1 for fungi. The specificity of the probes was confirmed, and no cross-reaction was detected in the present study. Furthermore, 13 cases were positive in the microarray and negative in culture. However, 4 cases were not identified as clear pathogens and only positive in the 16S probe sites. The diagnostic DNA microarray for intracranial infections has proven to be more rapid and sensitive, and it may be a better option for clinical application than culture methods.