Systematic review of Plasmodium falciparum and Plasmodium vivax polyclonal infections: Impact of prevalence, study population characteristics, and laboratory procedures

Genetic analysis of the circumsporozoite gene in Plasmodium falciparum isolates from Cameroon: Implications for efficacy and deployment of RTS,S/AS01 vaccine

Current understanding of genetic polymorphisms and natural selection in Plasmodium falciparum circumsporozoite (PfCSP), the leading malaria vaccine, is crucial for the development of next-generation vaccines, and such data is lacking in Africa. Blood samples were collected among Plasmodium-infected individuals living in four Cameroonian areas (Douala, Maroua, Mayo-Oulo, Pette). DNA samples were amplified using nested PCR protocols, sequenced, and BLASTed. Single nucleotide polymorphisms (SNPs) were analysed in each PfCSP region, and their impact on PfCSP function/structure was predicted in silico. The N-terminal region showed a limited polymorphism with four haplotypes, and three novel SNPs (N68Y, R87W, K93E) were found. Thirty-five haplotypes were identified in the central region, with several variants (e.g., NVNP and KANP). The C-terminal region was also highly diverse, with 25 haplotypes and eight novel SNPs (N290D, N308I, S312G, K317A, V344I, D356E, E357L, D359Y). Most polymorphic codon sites were mainly observed in the Th2R subregion in isolates from Douala and Pette. The codon site 321 was under episodic positive selection. One novel (E357L) and three known (K322I, G349D, D359Y) SNPs show an impact on function/structure. This study showed extensive genetic diversity with geographical patterns and evidence of the selection of Cameroonian PfCSP central and C-terminal regions.

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.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

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.

Genetic Profiling of Plasmodium ovale wallikeri Relapses With Microsatellite Markers and Whole-Genome Sequencing

Like Plasmodium vivax, both Plasmodium ovale curtisi and Plasmodium ovale wallikeri have the ability to cause relapse in humans, defined as recurring asexual parasitemia originating from liver-dormant forms subsequent to a primary infection. Here, we investigated relapse patterns in P ovale wallikeri infections from a cohort of travelers who were exposed to the parasite in sub-Saharan Africa and then experienced relapses after their return to France. Using a novel set of 8 highly polymorphic microsatellite markers, we genotyped 15 P ovale wallikeri relapses. For most relapses, the paired primary and relapse infections were highly genetically related (with 12 being homologous), an observation that was confirmed by whole-genome sequencing for the 4 relapses we further studied. This is, to our knowledge, the first genetic evidence of relapses in P ovale spp.

Multiplexed ddPCR-amplicon sequencing reveals isolated Plasmodium falciparum populations amenable to local elimination in Zanzibar, Tanzania

Zanzibar has made significant progress toward malaria elimination, but recent stagnation requires novel approaches. We developed a highly multiplexed droplet digital PCR (ddPCR)-based amplicon sequencing method targeting 35 microhaplotypes and drug-resistance loci, and successfully sequenced 290 samples from five districts covering both main islands. Here, we elucidate fine-scale Plasmodium falciparum population structure and infer relatedness and connectivity of infections using an identity-by-descent (IBD) approach. Despite high genetic diversity, we observe pronounced fine-scale spatial and temporal parasite genetic structure. Clusters of near-clonal infections on Pemba indicate persistent local transmission with limited parasite importation, presenting an opportunity for local elimination efforts. Furthermore, we observe an admixed parasite population on Unguja and detect a substantial fraction (2.9%) of significantly related infection pairs between Zanzibar and the mainland, suggesting recent importation. Our study provides a high-resolution view of parasite genetic structure across the Zanzibar archipelago and provides actionable insights for prioritizing malaria elimination efforts.

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.

Exploring how space, time, and sampling impact our ability to measure genetic structure across Plasmodium falciparum populations

A primary use of malaria parasite genomics is identifying highly related infections to quantify epidemiological, spatial, or temporal factors associated with patterns of transmission. For example, spatial clustering of highly related parasites can indicate foci of transmission and temporal differences in relatedness can serve as evidence for changes in transmission over time. However, for infections in settings of moderate to high endemicity, understanding patterns of relatedness is compromised by complex infections, overall high forces of infection, and a highly diverse parasite population. It is not clear how much these factors limit the utility of using genomic data to better understand transmission in these settings. In particular, further investigation is required to determine which patterns of relatedness we expect to see with high quality, densely sampled genomic data in a high transmission setting and how these observations change under different study designs, missingness, and biases in sample collection. Here we investigate two identity-by-state measures of relatedness and apply them to amplicon deep sequencing data collected as part of a longitudinal cohort in Western Kenya that has previously been analysed to identify individual-factors associated with sharing parasites with infected mosquitoes. With these data we use permutation tests, to evaluate several hypotheses about spatiotemporal patterns of relatedness compared to a null distribution. We observe evidence of temporal structure, but not of fine-scale spatial structure in the cohort data. To explore factors associated with the lack of spatial structure in these data, we construct a series of simplified simulation scenarios using an agent based model calibrated to entomological, epidemiological and genomic data from this cohort study to investigate whether the lack of spatial structure observed in the cohort could be due to inherent power limitations of this analytical method. We further investigate how our hypothesis testing behaves under different sampling schemes, levels of completely random and systematic missingness, and different transmission intensities.

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.

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.

High-throughput Plasmodium falciparum hrp2 and hrp3 gene deletion typing by digital PCR to monitor malaria rapid diagnostic test efficacy

Most rapid diagnostic tests for Plasmodium falciparum malaria target the Histidine-Rich Proteins 2 and 3 (HRP2 and HRP3). Deletions of the hrp2 and hrp3 genes result in false-negative tests and are a threat for malaria control. A novel assay for molecular surveillance of hrp2/hrp3 deletions was developed based on droplet digital PCR (ddPCR). The assay quantifies hrp2, hrp3, and a control gene with very high accuracy. The theoretical limit of detection was 0.33 parasites/µl. The deletion was reliably detected in mixed infections with wild-type and hrp2-deleted parasites at a density of >100 parasites/reaction. For a side-by-side comparison with the conventional nested PCR (nPCR) assay, 248 samples were screened in triplicate by ddPCR and nPCR. No deletions were observed by ddPCR, while by nPCR hrp2 deletion was observed in 8% of samples. The ddPCR assay was applied to screen 830 samples from Kenya, Zanzibar/Tanzania, Ghana, Ethiopia, Brazil, and Ecuador. Pronounced differences in the prevalence of deletions were observed among sites, with more hrp3 than hrp2 deletions. In conclusion, the novel ddPCR assay minimizes the risk of false-negative results (i.e., hrp2 deletion observed when the sample is wild type), increases sensitivity, and greatly reduces the number of reactions that need to be run.