Whole genome sequencing of Plasmodium falciparum from dried blood spots using selective whole genome amplification

Sustained clinical benefit of malaria chemoprevention with sulfadoxine-pyrimethamine (SP) in pregnant women in a region with high SP resistance markers

OBJECTIVE

The effectiveness of intermittent preventive treatment of malaria in pregnancy with sulfadoxine-pyrimethamine (IPTp-SP) is threatened by increasing SP-resistance in Africa. We assessed the level of SP-resistance markers, and the clinical and parasitological effectiveness of IPTp-SP in southern Mozambique.

METHODS

P. falciparum infection, antimalarial antibodies and dhfr/dhps SP-resistance mutants were detected by quantitative polymerase chain reaction (qPCR), suspension array technology and targeted deep sequencing, respectively, among 4016 HIV-negative women in Maputo province (2016-2019). Univariate and multivariate regression models were used to assess the association between taking the recommended three or more IPTp-SP doses (IPTp3+) and parasitological and clinical outcomes.

RESULTS

84.3% (3385/4016) women received three or more IPTp-SP doses. The prevalence of quintuple mutants at first antenatal care (ANC) visit was 94.2%. IPTp3+ was associated with a higher clearance rate of qPCR-detected infections from first ANC visit to delivery (adjusted odds ratio [aOR]=5.9, 95% CI: 1.5-33.3; p = 0.012), lower seroprevalence at delivery of antibodies against the pregnancy-specific antigen VAR2CSADBL34 (aOR=0.72, 95% CI: 0.54-0.95; p = 0.022), and lower prevalence of low birth weight deliveries (aOR: 0.61, 95% CI: 0.41-0.90; p = 0.013).

CONCLUSION

A sustained parasitological effect of IPTp-SP contributes to the clinical effectiveness of IPTp3+ in areas with high prevalence of SP-resistance markers.

Sterile protection against P. vivax malaria by repeated blood stage infection in the Aotus monkey model

The malaria parasite Plasmodium vivax remains a major global public health challenge, and no vaccine is approved for use in humans. Here, we assessed whether P. vivax strain-transcendent immunity can be achieved by repeated infection in Aotus monkeys. Sterile immunity was achieved after two homologous infections, whereas subsequent heterologous challenge provided only partial protection. IgG levels based on P. vivax lysate ELISA and protein microarray increased with repeated infections and correlated with the level of homologous protection. Parasite transcriptional profiles provided no evidence of major antigenic switching upon homologous or heterologous challenge. However, we observed significant sequence diversity and transcriptional differences in the P. vivax core gene repertoire between the two strains used in the study, suggesting that partial protection upon heterologous challenge is due to molecular differences between strains rather than immune evasion by antigenic switching. Our study demonstrates that sterile immunity against P. vivax can be achieved by repeated homologous blood stage infection in Aotus monkeys, thus providing a benchmark to test the efficacy of candidate blood stage P. vivax malaria vaccines.

Flexible and cost-effective genomic surveillance of P. falciparum malaria with targeted nanopore sequencing

Genomic surveillance of Plasmodium falciparum malaria can provide policy-relevant information about antimalarial drug resistance, diagnostic test failure, and the evolution of vaccine targets. Yet the large and low complexity genome of P. falciparum complicates the development of genomic methods, while resource constraints in malaria endemic regions can limit their deployment. Here, we demonstrate an approach for targeted nanopore sequencing of P. falciparum from dried blood spots (DBS) that enables cost-effective genomic surveillance of malaria in low-resource settings. We release software that facilitates flexible design of amplicon sequencing panels and use this software to design two target panels for P. falciparum. The panels generate 3-4 kbp reads for eight and sixteen targets respectively, covering key drug-resistance associated genes, diagnostic test antigens, polymorphic markers and the vaccine target csp. We validate our approach on mock and field samples, demonstrating robust sequencing coverage, accurate variant calls within coding sequences, the ability to explore P. falciparum within-sample diversity and to detect deletions underlying rapid diagnostic test failure.

Investigation of an airport-associated cluster of falciparum malaria in Frankfurt, Germany, 2022

Airport malaria is uncommon but increasing in Europe and often difficult to diagnose. We describe the clinical, epidemiological and environmental investigations of a cluster of airport malaria cases and measures taken in response. Three Frankfurt International Airport employees without travel histories to malaria-endemic areas were diagnosed with Plasmodium falciparum malaria in Germany in 2022. Two cases were diagnosed within 1 week, and the third one after 10 weeks. Two cases had severe disease, all three recovered fully. The cases worked in separate areas and no specific location for the transmissions could be identified. No additional cases were detected among airport employees. In June and July, direct flights from Equatorial Guinea, Nigeria and Angola and one parcel originating in Ghana arrived at Frankfurt airport. No vector-competent mosquitoes could be trapped to identify the source of the outbreak. Whole genome sequencing of P. falciparum genomes showed a high genetic relatedness between samples of the three cases and suggested the geographical origin closest to Ghana. A diagnosis of airport malaria should prompt appropriate and comprehensive outbreak investigations to identify the source and to prevent severe forms of falciparum malaria.

Temporal and spatial dynamics of Plasmodium falciparum clonal lineages in Guyana

Plasmodium parasites, the causal agents of malaria, are eukaryotic organisms that obligately undergo sexual recombination within mosquitoes. However, in low transmission settings where most mosquitoes become infected with only a single parasite clone, parasites recombine with themselves, and the clonal lineage is propagated rather than broken up by outcrossing. We investigated whether stochastic/neutral factors drive the persistence and abundance of Plasmodium falciparum clonal lineages in Guyana, a country with relatively low malaria transmission, but the only setting in the Americas in which an important artemisinin resistance mutation (pfk13 C580Y) has been observed. To investigate whether this clonality was potentially associated with the persistence and spatial spread of the mutation, we performed whole genome sequencing on 1,727 Plasmodium falciparum samples collected from infected patients across a five-year period (2016-2021). We characterized the relatedness between each pair of monoclonal infections (n=1,409) through estimation of identity by descent (IBD) and also typed each sample for known or candidate drug resistance mutations. A total of 160 clones (mean IBD ≥ 0.90) were circulating in Guyana during the study period, comprising 13 highly related clusters (mean IBD ≥ 0.40). In the five-year study period, we observed a decrease in frequency of a mutation associated with artemisinin partner drug (piperaquine) resistance (pfcrt C350R) and limited co-occurence of pfcrt C350R with duplications of plasmepsin 2/3, an epistatic interaction associated with piperaquine resistance. We additionally report polymorphisms exhibiting evidence of selection for drug resistance or other phenotypes and reported a novel pfk13 mutation (G718S) as well as 61 nonsynonymous substitutions that increased markedly in frequency. However, P. falciparum clonal dynamics in Guyana appear to be largely driven by stochastic factors, in contrast to other geographic regions. The use of multiple artemisinin combination therapies in Guyana may have contributed to the disappearance of the pfk13 C580Y mutation.

Molecular markers for malaria genetic epidemiology: progress and pitfalls

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

Evidence of artemisinin partial resistance in North-western Tanzania: clinical and drug resistance markers study

Background

Artemisinin-based combination therapies (ACTs) are the recommended antimalarial drugs for the treatment of uncomplicated malaria. The recent emergence of artemisinin partial resistance (ART-R) in Rwanda, Uganda and Eritrea is of great concern. In Tanzania, a nationwide molecular malaria surveillance in 2021 showed a high prevalence of the Kelch13 (K13) 561H mutation in Plasmodium falciparum from the north-western region, close to the border with Rwanda and Uganda. This study was conducted in 2022 to evaluate the efficacy of artemether-lumefantrine (AL) and artesunate-amodiaquine (ASAQ) for the treatment of uncomplicated falciparum malaria and to confirm the presence of ART-R in Tanzania.

Methods

This single-arm study evaluated the efficacy of AL and ASAQ in eligible children aged six months to 10 years at Bukangara Dispensary in Karagwe District, Kagera Region. Clinical and parasitological responses were monitored for 28 days according to standard WHO protocol. Mutations in K13 gene and extended haplotypes with these mutations were analysed using Sanger and whole genome sequencing data, respectively.

Findings

176 children (88 in each AL and ASAQ group) were enrolled and all achieved the defined outcomes. PCR-corrected adequate clinical and parasitological response (ACPR) was 98.3% (95% CI: 90.8-100) and 100.0% (95% CI: 95.8-100) for AL and ASAQ, respectively. Parasitaemia on day 3 was observed in 11/88 (12.5%) and 17/88 (19.3%) in the AL and ASAQ groups, respectively. The half-life of parasitaemia was significantly higher (>6.5 hrs) in patients with parasitaemia on day 3 and/or mutations in K13 gene at enrolment. Most patients with parasitaemia on day 3 (8/11 = 72.7% in the AL group and 10/17 = 58.8% in the ASAQ group) had 561H mutation at enrolment. The parasites with K13 mutations were not similar to those from south-east Asia and Rwanda, but had the same core haplotype of a new 561H haplotype reported in Kagera in 2021.

Interpretation

These findings confirm the presence of ART-R in Tanzania. A context-specific strategy to respond to artemisinin partial resistance is urgently needed. Although both AL and ASAQ showed high efficacy, increased vigilance for reduced efficacy of these ACTs and detection of ART-R in other parts of the country is critical.

Long-Read Sequencing and De Novo Genome Assembly Pipeline of Two Plasmodium falciparum Clones (Pf3D7, PfW2) Using Only the PromethION Sequencer from Oxford Nanopore Technologies without Whole-Genome Amplification

Antimalarial drug resistance has become a real public health problem despite WHO measures. New sequencing technologies make it possible to investigate genomic variations associated with resistant phenotypes at the genome-wide scale. Based on the use of hemisynthetic nanopores, the PromethION technology from Oxford Nanopore Technologies can produce long-read sequences, in contrast to previous short-read technologies used as the gold standard to sequence Plasmodium. Two clones of P. falciparum (Pf3D7 and PfW2) were sequenced in long-read using the PromethION sequencer from Oxford Nanopore Technologies without genomic amplification. This made it possible to create a processing analysis pipeline for human Plasmodium with ONT Fastq only. De novo assembly revealed N50 lengths of 18,488 kb and 17,502 kb for the Pf3D7 and PfW2, respectively. The genome size was estimated at 23,235,407 base pairs for the Pf3D7 clone and 21,712,038 base pairs for the PfW2 clone. The average genome coverage depth was estimated at 787X and 653X for the Pf3D7 and PfW2 clones, respectively. This study proposes an assembly processing pipeline for the human Plasmodium genome using software adapted to large ONT data and the high AT percentage of Plasmodium. This search provides all the parameters which were optimized for use with the software selected in the pipeline.

Selective whole-genome sequencing of Plasmodium parasites directly from blood samples by nanopore adaptive sampling

IMPORTANCE

Malaria is caused by parasites of the genus Plasmodium, and reached a global disease burden of 247 million cases in 2021. To study drug resistance mutations and parasite population dynamics, whole-genome sequencing of patient blood samples is commonly performed. However, the predominance of human DNA in these samples imposes the need for time-consuming laboratory procedures to enrich Plasmodium DNA. We used the Oxford Nanopore Technologies' adaptive sampling feature to circumvent this problem and enrich Plasmodium reads directly during the sequencing run. We demonstrate that adaptive nanopore sequencing efficiently enriches Plasmodium reads, which simplifies and shortens the timeline from blood collection to parasite sequencing. In addition, we show that the obtained data can be used for monitoring genetic markers, or to generate nearly complete genomes. Finally, owing to its inherent mobility, this technology can be easily applied on-site in endemic areas where patients would benefit the most from genomic surveillance.

Evolution and spread of Plasmodium falciparum mutations associated with resistance to sulfadoxine-pyrimethamine in central Africa: a cross-sectional study

BACKGROUND

Efficacy of sulfadoxine-pyrimethamine, the malaria chemoprophylaxis used in pregnant women, and in children when combined with amodiaquine, is threatened by the accumulation of mutations in the Plasmodium falciparum dihydropteroate synthase (pfdhps) and dihydrofolate reductase (pfdhfr) genes. Data on the prevalence of resistant alleles in central Africa and the new pfdhps I431V mutation, particularly associated with other mutations to form the pfdhps vagKgs allele, are scarce. We explored the frequency and geographical distribution of pfdhps and pfdhfr mutations in central Africa in 2014-18, and assessed the evolutionary origin of the vagKgs allele.

METHODS

Samples were collected at 18 health-care centres in seven countries (Angola, Cameroon, Central African Republic, Democratic Republic of the Congo, Gabon, Nigeria, and Republic of the Congo) from patients who showed possible symptoms of malaria between March 1, 2014, and Oct 31, 2018. Samples that were positive for P falciparum were transported to a laboratory in Toulouse, France, and genotyped. The frequency of pfdhfr and pfdhps mutations was studied in 1749 samples. Microsatellites in pfdhps flanking regions and whole-genome analysis compared with parasite genomes from the data-sharing network MalariaGEN were performed on samples carrying the vagKgs allele.

FINDINGS

Mapping of the prevalence of single nucleotide polymorphisms and corresponding alleles of pfdhfr and pfdhps showed a substantial spread of alleles associated with sulfadoxine-pyrimethamine resistance in central Africa during the 2014-18 period, especially an increase going west to east in pfdhps alleles carrying the K540E and A581G mutations. A high prevalence of the pfdhps I431V mutation was observed in Cameroon (exceeding 50% in the northern region) and Nigeria. Genomic analysis showed a recent African emergence and a clonal expansion of the most frequent pfdhps vagKgs allele.

INTERPRETATION

Reduced sulfadoxine-pyrimethamine efficacy due to increased resistance is a worrying situation, especially because the malaria transmission level is high in central Africa. Although the resistance phenotype remains to be confirmed, the emergence and spread of the vagKgs allele in west and central Africa could challenge the use of sulfadoxine-pyrimethamine.

FUNDING

Toulouse Institute for Infectious and Inflammatory Diseases.

Country wide surveillance reveals prevalent artemisinin partial resistance mutations with evidence for multiple origins and expansion of high level sulfadoxine-pyrimethamine resistance mutations in northwest Tanzania

Background

Emergence of artemisinin partial resistance (ART-R) in Plasmodium falciparum is a growing threat to the efficacy of artemisinin combination therapies (ACT) and the efforts for malaria elimination. The emergence of Plasmodium falciparum Kelch13 (K13) R561H in Rwanda raised concern about the impact in neighboring Tanzania. In addition, regional concern over resistance affecting sulfadoxine-pyrimethamine (SP), which is used for chemoprevention strategies, is high.

Methods

To enhance longitudinal monitoring, the Molecular Surveillance of Malaria in Tanzania (MSMT) project was launched in 2020 with the goal of assessing and mapping antimalarial resistance. Community and clinic samples were assessed for resistance polymorphisms using a molecular inversion probe platform.

Findings

Genotyping of 6,278 samples collected countrywide in 2021 revealed a focus of K13 561H mutants in northwestern Tanzania (Kagera) with prevalence of 7.7% (50/649). A small number of 561H mutants (about 1%) were found as far as 800 km away in Tabora, Manyara, and Njombe. Genomic analysis suggests some of these parasites are highly related to isolates collected in Rwanda in 2015, supporting regional spread of 561H. However, a novel haplotype was also observed, likely indicating a second origin in the region. Other validated resistance polymorphisms (622I and 675V) were also identified. A focus of high sulfadoxine-pyrimethamine drug resistance was also identified in Kagera with a prevalence of dihydrofolate reductase 164L of 15% (80/526).

Interpretation

These findings demonstrate the K13 561H mutation is entrenched in the region and that multiple origins of ART-R, similar as to what was seen in Southeast Asia, have occurred. Mutations associated with high levels of SP resistance are increasing. These results raise concerns about the long-term efficacy of artemisinin and chemoprevention antimalarials in the region.

Funding

This study was funded by the Bill and Melinda Gates Foundation and the National Institutes of Health.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

We here analyze data from the first year of an ongoing nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal. The analysis is based on 1097 samples collected at health facilities during passive malaria case detection in 2019; it provides a baseline for analyzing parasite genetic metrics as they vary over time and geographic space. The study's goal was to identify genetic metrics that were informative about transmission intensity and other aspects of transmission dynamics, focusing on measures of genetic relatedness between parasites. We found the best genetic proxy for local malaria incidence to be the proportion of polygenomic infections (those with multiple genetically distinct parasites), although this relationship broke down at low incidence. The proportion of related parasites was less correlated with incidence while local genetic diversity was uninformative. The type of relatedness could discriminate local transmission patterns: two nearby areas had similarly high fractions of relatives, but one was dominated by clones and the other by outcrossed relatives. Throughout Senegal, 58% of related parasites belonged to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci and at one novel locus, reflective of ongoing selection pressure.

Whole-genome surveillance identifies markers of Plasmodium falciparum drug resistance and novel genomic regions under selection in Mozambique

IMPORTANCE

Malaria is a devastating disease caused by Plasmodium parasites. The evolution of parasite drug resistance continues to hamper progress toward malaria elimination, and despite extensive efforts to control malaria, it remains a leading cause of death in Mozambique and other countries in the region. The development of successful vaccines and identification of molecular markers to track drug efficacy are essential for managing the disease burden. We present an analysis of the parasite genome in Mozambique, a country with one of the highest malaria burdens globally and limited available genomic data, revealing current selection pressure. We contribute additional evidence to limited prior studies supporting the effectiveness of SWGA in producing reliable genomic data from complex clinical samples. Our results provide the identity of genomic loci that may be associated with current antimalarial drug use, including artemisinin and lumefantrine, and reveal selection pressure predicted to compromise the efficacy of current vaccine candidates.

Increasing Prevalence of Artemisinin-Resistant HRP2-Negative Malaria in Eritrea

BACKGROUND

Although the clinical efficacy of antimalarial artemisinin-based combination therapies in Africa remains high, the recent emergence of partial resistance to artemisinin in Plasmodium falciparum on the continent is troubling, given the lack of alternative treatments.

METHODS

In this study, we used data from drug-efficacy studies conducted between 2016 and 2019 that evaluated 3-day courses of artemisinin-based combination therapy (artesunate-amodiaquine or artemether-lumefantrine) for uncomplicated malaria in Eritrea to estimate the percentage of patients with day-3 positivity (i.e., persistent P. falciparum parasitemia 3 days after the initiation of therapy). We also assayed parasites for mutations in Pfkelch13 as predictive markers of partial resistance to artemisinin and screened for deletions in hrp2 and hrp3 that result in variable performance of histidine rich protein 2 (HRP2)-based rapid diagnostic tests for malaria.

RESULTS

We noted an increase in the percentage of patients with day-3 positivity from 0.4% (1 of 273) in 2016 to 1.9% (4 of 209) in 2017 and 4.2% (15 of 359) in 2019. An increase was also noted in the prevalence of the Pfkelch13 R622I mutation, which was detected in 109 of 818 isolates before treatment, from 8.6% (24 of 278) in 2016 to 21.0% (69 of 329) in 2019. The odds of day-3 positivity increased by a factor of 6.2 (95% confidence interval, 2.5 to 15.5) among the patients with Pfkelch13 622I variant parasites. Partial resistance to artemisinin, as defined by the World Health Organization, was observed in Eritrea. More than 5% of the patients younger than 15 years of age with day-3 positivity also had parasites that carried Pfkelch13 R622I. In vitro, the R622I mutation conferred a low level of resistance to artemisinin when edited into NF54 and Dd2 parasite lines. Deletions in both hrp2 and hrp3 were identified in 16.9% of the parasites that carried the Pfkelch13 R622I mutation, which made them potentially undetectable by HRP2-based rapid diagnostic tests.

CONCLUSIONS

The emergence and spread of P. falciparum lineages with both Pfkelch13-mediated partial resistance to artemisinin and deletions in hrp2 and hrp3 in Eritrea threaten to compromise regional malaria control and elimination campaigns. (Funded by the Bill and Melinda Gates Foundation and others; Australian New Zealand Clinical Trials Registry numbers, ACTRN12618001223224, ACTRN12618000353291, and ACTRN12619000859189.).

Molecular surveillance of Plasmodium falciparum drug-resistance markers in Vietnam using multiplex amplicon sequencing (2000-2016)

Emergence and spread of Plasmodium falciparum resistance to artemisinin-based combination therapies (ACT) is a major challenge for Greater Mekong Subregion countries in their goal to eliminate malaria by 2030. Tools to efficiently monitor drug resistance beyond resource-demanding therapeutic efficacy studies are necessary. A custom multiplex amplicon sequencing assay based on Illumina technology was designed to target the marker of partial resistance to artemisinin (K13), five candidate modulators of artemisinin resistance, the marker of resistance to chloroquine (crt), and four neutral microsatellite loci. The assay was used to genotype 635 P. falciparum-positive blood samples collected across seven provinces of Vietnam and one of Cambodia between 2000 and 2016. Markers of resistance to artemisinin partner-drugs piperaquine (copy number of plasmepsin-2) and mefloquine (copy number of multidrug-resistance 1) were determined by qPCR. Parasite population structure was further assessed using a 101-SNP barcode. Validated mutations of artemisinin partial resistance in K13 were found in 48.1% of samples, first detection was in 2000, and by 2015 prevalence overcame > 50% in Central Highlands and Binh Phuoc province. K13-C580Y variant became predominant country-wide, quickly replacing an outbreak of K13-I543T in Central Highlands. Mutations in candidate artemisinin resistance modulator genes paralleled the trends of K13 mutants, whereas resistance to piperaquine and mefloquine remained low (≈ 10%) by 2015-2016. Genomic tools applied to malaria surveillance generate comprehensive information on dynamics of drug resistance and population structure and reflect drug efficacy profiles from in vivo studies.

Drug resistance profiling of asymptomatic and low-density Plasmodium falciparum malaria infections on Ngodhe island, Kenya, using custom dual-indexing next-generation sequencing

Malaria control initiatives require rapid and reliable methods for the detection and monitoring of molecular markers associated with antimalarial drug resistance in Plasmodium falciparum parasites. Ngodhe island, Kenya, presents a unique malaria profile, with lower P. falciparum incidence rates than the surrounding region, and a high proportion of sub-microscopic and low-density infections. Here, using custom dual-indexing and Illumina next generation sequencing, we generate resistance profiles on seventy asymptomatic and low-density P. falciparum infections from a mass drug administration program implemented on Ngodhe island between 2015 and 2016. Our assay encompasses established molecular markers on the Pfcrt, Pfmdr1, Pfdhps, Pfdhfr, and Pfk13 genes. Resistance markers for sulfadoxine-pyrimethamine were identified at high frequencies, including a quintuple mutant haplotype (Pfdhfr/Pfdhps: N51I, C59R, S108N/A437G, K540E) identified in 62.2% of isolates. The Pfdhps K540E biomarker, used to inform decision making for intermittent preventative treatment in pregnancy, was identified in 79.2% of isolates. Several variants on Pfmdr1, associated with reduced susceptibility to quinolones and lumefantrine, were also identified (Y184F 47.1%; D1246Y 16.0%; N86 98%). Overall, we have presented a low-cost and extendable approach that can provide timely genetic profiles to inform clinical and surveillance activities, especially in settings with abundant low-density infections, seeking malaria elimination.

Protocol for a cluster randomised placebo-controlled trial of adjunctive ivermectin mass drug administration for malaria control on the Bijagós Archipelago of Guinea-Bissau: the MATAMAL trial

INTRODUCTION

As malaria declines, innovative tools are required to further reduce transmission and achieve elimination. Mass drug administration (MDA) of artemisinin-based combination therapy (ACT) is capable of reducing malaria transmission where coverage of control interventions is already high, though the impact is short-lived. Combining ACT with ivermectin, an oral endectocide shown to reduce vector survival, may increase its impact, while also treating ivermectin-sensitive co-endemic diseases and minimising the potential impact of ACT resistance in this context.

METHODS AND ANALYSIS

MATAMAL is a cluster-randomised placebo-controlled trial. The trial is being conducted in 24 clusters on the Bijagós Archipelago, Guinea-Bissau, where the peak prevalence of Plasmodium falciparum (Pf) parasitaemia is approximately 15%. Clusters have been randomly allocated to receive MDA with dihydroartemisinin-piperaquine and either ivermectin or placebo. The primary objective is to determine whether the addition of ivermectin MDA is more effective than dihydroartemisinin-piperaquine MDA alone in reducing the prevalence of P. falciparum parasitaemia, measured during peak transmission season after 2 years of seasonal MDA. Secondary objectives include assessing prevalence after 1 year of MDA; malaria incidence monitored through active and passive surveillance; age-adjusted prevalence of serological markers indicating exposure to P. falciparum and anopheline mosquitoes; vector parous rates, species composition, population density and sporozoite rates; prevalence of vector pyrethroid resistance; prevalence of artemisinin resistance in P. falciparum using genomic markers; ivermectin's impact on co-endemic diseases; coverage estimates; and the safety of combined MDA.

ETHICS AND DISSEMINATION

The trial has been approved by the London School of Hygiene and Tropical Medicine's Ethics Committee (UK) (19156) and the Comite Nacional de Eticas de Saude (Guinea-Bissau) (084/CNES/INASA/2020). Results will be disseminated in peer-reviewed publications and in discussion with the Bissau-Guinean Ministry of Public Health and participating communities.

TRIAL REGISTRATION NUMBER

NCT04844905.

Genome-wide genetic variation and molecular surveillance of drug resistance in Plasmodium falciparum isolates from asymptomatic individuals in Ouélessébougou, Mali

Sequence analysis of Plasmodium falciparum parasites is informative in ensuring sustained success of malaria control programmes. Whole-genome sequencing technologies provide insights into the epidemiology and genome-wide variation of P. falciparum populations and can characterise geographical as well as temporal changes. This is particularly important to monitor the emergence and spread of drug resistant P. falciparum parasites which is threatening malaria control programmes world-wide. Here, we provide a detailed characterisation of genome-wide genetic variation and drug resistance profiles in asymptomatic individuals in South-Western Mali, where malaria transmission is intense and seasonal, and case numbers have recently increased. Samples collected from Ouélessébougou, Mali (2019-2020; n = 87) were sequenced and placed in the context of older Malian (2007-2017; n = 876) and African-wide (n = 711) P. falciparum isolates. Our analysis revealed high multiclonality and low relatedness between isolates, in addition to increased frequencies of molecular markers for sulfadoxine-pyrimethamine and lumefantrine resistance, compared to older Malian isolates. Furthermore, 21 genes under selective pressure were identified, including a transmission-blocking vaccine candidate (pfCelTOS) and an erythrocyte invasion locus (pfdblmsp2). Overall, our work provides the most recent assessment of P. falciparum genetic diversity in Mali, a country with the second highest burden of malaria in West Africa, thereby informing malaria control activities.

Targeted and whole-genome sequencing reveal a north-south divide in P. falciparum drug resistance markers and genetic structure in Mozambique

Mozambique is one of the four African countries which account for over half of all malaria deaths worldwide, yet little is known about the parasite genetic structure in that country. We performed P. falciparum amplicon and whole genome sequencing on 2251 malaria-infected blood samples collected in 2015 and 2018 in seven provinces of Mozambique to genotype antimalarial resistance markers and interrogate parasite population structure using genome-wide microhaplotyes. Here we show that the only resistance-associated markers observed at frequencies above 5% were pfmdr1-184F (59%), pfdhfr-51I/59 R/108 N (99%) and pfdhps-437G/540E (89%). The frequency of pfdhfr/pfdhps quintuple mutants associated with sulfadoxine-pyrimethamine resistance increased from 80% in 2015 to 89% in 2018 (p < 0.001), with a lower expected heterozygosity and higher relatedness of microhaplotypes surrounding pfdhps mutants than wild-type parasites suggestive of recent selection. pfdhfr/pfdhps quintuple mutants also increased from 72% in the north to 95% in the south (2018; p < 0.001). This resistance gradient was accompanied by a concentration of mutations at pfdhps-436 (17%) in the north, a south-to-north increase in the genetic complexity of P. falciparum infections (p = 0.001) and a microhaplotype signature of regional differentiation. The parasite population structure identified here offers insights to guide antimalarial interventions and epidemiological surveys.

Malaria surveillance reveals parasite relatedness, signatures of selection, and correlates of transmission across Senegal

Parasite genetic surveillance has the potential to play an important role in malaria control. We describe here an analysis of data from the first year of an ongoing, nationwide program of genetic surveillance of Plasmodium falciparum parasites in Senegal, intended to provide actionable information for malaria control efforts. Looking for a good proxy for local malaria incidence, we found that the best predictor was the proportion of polygenomic infections (those with multiple genetically distinct parasites), although that relationship broke down in very low incidence settings (r = 0.77 overall). The proportion of closely related parasites in a site was more weakly correlated ( r = -0.44) with incidence while the local genetic diversity was uninformative. Study of related parasites indicated their potential for discriminating local transmission patterns: two nearby study areas had similarly high fractions of relatives, but one area was dominated by clones and the other by outcrossed relatives. Throughout the country, 58% of related parasites proved to belong to a single network of relatives, within which parasites were enriched for shared haplotypes at known and suspected drug resistance loci as well as at one novel locus, reflective of ongoing selection pressure.

Genomic analysis of Indian isolates of Plasmodium falciparum: Implications for drug resistance and virulence factors

The emergence of drug resistance to frontline treatments such as Artemisinin-based combination therapy (ACT) is a major obstacle to the control and eradication of malaria. This problem is compounded by the inherent genetic variability of the parasites, as many established markers of resistance do not accurately predict the drug-resistant status. There have been reports of declining effectiveness of ACT in the West Bengal and Northeast regions of India, which have traditionally been areas of drug resistance emergence in the country. Monitoring the genetic makeup of a population can help to identify the potential for drug resistance markers associated with it and evaluate the effectiveness of interventions aimed at reducing the spread of malaria. In this study, we performed whole genome sequencing of 53 isolates of Plasmodium falciparum from West Bengal and compared their genetic makeup to isolates from Southeast Asia (SEA) and Africa. We found that the Indian isolates had a distinct genetic makeup compared to those from SEA and Africa, and were more similar to African isolates, with a high prevalence of mutations associated with antigenic variation genes. The Indian isolates also showed a high prevalence of markers of chloroquine resistance (mutations in Pfcrt) and multidrug resistance (mutations in Pfmdr1), but no known mutations associated with artemisinin resistance in the PfKelch13 gene. Interestingly, we observed a novel L152V mutation in PfKelch13 gene and other novel mutations in genes involved in ubiquitination and vesicular transport that have been reported to support artemisinin resistance in the early stages of ACT resistance in the absence of PfKelch13 polymorphisms. Thus, our study highlights the importance of region-specific genomic surveillance for artemisinin resistance and the need for continued monitoring of resistance to artemisinin and its partner drugs.

Genomic variation during culture adaptation of genetically complex Plasmodium falciparum clinical isolates

Experimental studies on the biology of malaria parasites have mostly been based on laboratory-adapted lines, but there is limited understanding of how these may differ from parasites in natural infections. Loss-of-function mutants have previously been shown to emerge during culture of some Plasmodium falciparum clinical isolates in analyses focusing on single-genotype infections. The present study included a broader array of isolates, mostly representing multiple-genotype infections, which are more typical in areas where malaria is highly endemic. Genome sequence data from multiple time points over several months of culture adaptation of 28 West African isolates were analysed, including previously available sequences along with new genome sequences from additional isolates and time points. Some genetically complex isolates eventually became fixed over time to single surviving genotypes in culture, whereas others retained diversity, although proportions of genotypes varied over time. Drug resistance allele frequencies did not show overall directional changes, suggesting that resistance-associated costs are not the main causes of fitness differences among parasites in culture. Loss-of-function mutants emerged during culture in several of the multiple-genotype isolates, affecting genes (including AP2-HS, EPAC and SRPK1) for which loss-of-function mutants were previously seen to emerge in single-genotype isolates. Parasite clones were derived by limiting dilution from six of the isolates, and sequencing identified de novo variants not detected in the bulk isolate sequences. Interestingly, several of these were nonsense mutants and frameshifts disrupting the coding sequence of EPAC, the gene with the largest number of independent nonsense mutants previously identified in laboratory-adapted lines. Analysis of genomic identity by descent to explore relatedness among clones revealed co-occurring non-identical sibling parasites, illustrative of the natural genetic structure within endemic populations.

Population dynamics and drug resistance mutations in Plasmodium falciparum on the Bijagós Archipelago, Guinea-Bissau

Following integrated malaria control interventions, malaria burden on the Bijagós Archipelago has significantly decreased. Understanding the genomic diversity of circulating Plasmodium falciparum malaria parasites can assist infection control, through identifying drug resistance mutations and characterising the complexity of population structure. This study presents the first whole genome sequence data for P. falciparum isolates from the Bijagós Archipelago. Amplified DNA from P. falciparum isolates sourced from dried blood spot samples of 15 asymptomatic malaria cases were sequenced. Using 1.3 million SNPs characterised across 795 African P. falciparum isolates, population structure analyses revealed that isolates from the archipelago cluster with samples from mainland West Africa and appear closely related to mainland populations; without forming a separate phylogenetic cluster. This study characterises SNPs associated with antimalarial drug resistance on the archipelago. We observed fixation of the PfDHFR mutations N51I and S108N, associated with resistance to sulphadoxine-pyrimethamine, and the continued presence of PfCRT K76T, associated with chloroquine resistance. These data have relevance for infection control and drug resistance surveillance; particularly considering expected increases in antimalarial drug use following updated WHO recommendations, and the recent implementation of seasonal malaria chemoprevention and mass drug administration in the region.

Comparison of molecular surveillance methods to assess changes in the population genetics of Plasmodium falciparum in high transmission

A major motivation for developing molecular methods for malaria surveillance is to measure the impact of control interventions on the population genetics of Plasmodium falciparum as a potential marker of progress towards elimination. Here we assess three established methods (i) single nucleotide polymorphism (SNP) barcoding (panel of 24-biallelic loci), (ii) microsatellite genotyping (panel of 12-multiallelic loci), and (iii) varcoding (fingerprinting var gene diversity, akin to microhaplotyping) to identify changes in parasite population genetics in response to a short-term indoor residual spraying (IRS) intervention. Typical of high seasonal transmission in Africa, multiclonal infections were found in 82.3% (median 3; range 1-18) and 57.8% (median 2; range 1-12) of asymptomatic individuals pre- and post-IRS, respectively, in Bongo District, Ghana. Since directly phasing multilocus haplotypes for population genetic analysis is not possible for biallelic SNPs and microsatellites, we chose ~200 low-complexity infections biased to single and double clone infections for analysis. Each genotyping method presented a different pattern of change in diversity and population structure as a consequence of variability in usable data and the relative polymorphism of the molecular markers (i.e., SNPs < microsatellites < var). Varcoding and microsatellite genotyping showed the overall failure of the IRS intervention to significantly change the population structure from pre-IRS characteristics (i.e., many diverse genomes of low genetic similarity). The 24-SNP barcode provided limited information for analysis, largely due to the biallelic nature of SNPs leading to a high proportion of double-allele calls and a view of more isolate relatedness compared to microsatellites and varcoding. Relative performance, suitability, and cost-effectiveness of the methods relevant to sample size and local malaria elimination in high-transmission endemic areas are discussed.

A fast machine-learning-guided primer design pipeline for selective whole genome amplification

Addressing many of the major outstanding questions in the fields of microbial evolution and pathogenesis will require analyses of populations of microbial genomes. Although population genomic studies provide the analytical resolution to investigate evolutionary and mechanistic processes at fine spatial and temporal scales-precisely the scales at which these processes occur-microbial population genomic research is currently hindered by the practicalities of obtaining sufficient quantities of the relatively pure microbial genomic DNA necessary for next-generation sequencing. Here we present swga2.0, an optimized and parallelized pipeline to design selective whole genome amplification (SWGA) primer sets. Unlike previous methods, swga2.0 incorporates active and machine learning methods to evaluate the amplification efficacy of individual primers and primer sets. Additionally, swga2.0 optimizes primer set search and evaluation strategies, including parallelization at each stage of the pipeline, to dramatically decrease program runtime. Here we describe the swga2.0 pipeline, including the empirical data used to identify primer and primer set characteristics, that improve amplification performance. Additionally, we evaluate the novel swga2.0 pipeline by designing primer sets that successfully amplify Prevotella melaninogenica, an important component of the lung microbiome in cystic fibrosis patients, from samples dominated by human DNA.

Selective whole-genome amplification reveals population genetics of Leishmania braziliensis directly from patient skin biopsies

In Brazil, Leishmania braziliensis is the main causative agent of the neglected tropical disease, cutaneous leishmaniasis (CL). CL presents on a spectrum of disease severity with a high rate of treatment failure. Yet the parasite factors that contribute to disease presentation and treatment outcome are not well understood, in part because successfully isolating and culturing parasites from patient lesions remains a major technical challenge. Here we describe the development of selective whole genome amplification (SWGA) for Leishmania and show that this method enables culture-independent analysis of parasite genomes obtained directly from primary patient skin samples, allowing us to circumvent artifacts associated with adaptation to culture. We show that SWGA can be applied to multiple Leishmania species residing in different host species, suggesting that this method is broadly useful in both experimental infection models and clinical studies. SWGA carried out directly on skin biopsies collected from patients in Corte de Pedra, Bahia, Brazil, showed extensive genomic diversity. Finally, as a proof-of-concept, we demonstrated that SWGA data can be integrated with published whole genome data from cultured parasite isolates to identify variants unique to specific geographic regions in Brazil where treatment failure rates are known to be high. SWGA provides a relatively simple method to generate Leishmania genomes directly from patient samples, unlocking the potential to link parasite genetics with host clinical phenotypes.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

Assessment of artemisinin tolerance in Plasmodium falciparum clinical isolates in children with uncomplicated malaria in Ghana

BACKGROUND

Artemisinin-based combination therapy (ACT) is the first-line treatment for uncomplicated malaria in Ghana. Artemisinin (ART) tolerance in Plasmodium falciparum has arisen in Southeast Asia and recently, in parts of East Africa. This is ascribed to the survival of ring-stage parasites post treatment. The present study sought to assess and characterize correlates of potential ART tolerance based on post-treatment parasite clearance, ex vivo and in vitro drug sensitivity, and molecular markers of drug resistance in P. falciparum isolates from children with uncomplicated malaria in Ghana.

METHODS

Six months to fourteen years old children presenting with acute uncomplicated malaria (n = 115) were enrolled in two hospitals and a Health Centre in Ghana's Greater Accra region and treated with artemether-lumefantrine (AL) according to body weight. Pre- and post-treatment parasitaemia (day 0 and day 3) was confirmed by microscopy. The ex vivo ring-stage survival assay (RSA) was used to detect percent ring survival while the 72 h SYBR Green I assay was used to measure the 50% inhibition concentration (IC_50s) of ART and its derivatives and partner drugs. Genetic markers of drug tolerance /resistance were evaluated using selective whole genome sequencing.

RESULTS

Of the total of 115 participants, 85 were successfully followed up on day 3 post-treatment and 2/85 (2.4%) had parasitaemia. The IC₅₀ values of ART, artesunate (AS), artemether (AM), dihydroartemisinin (DHA), amodiaquine (AQ), and lumefantrine (LUM) were not indicative of drug tolerance. However, 7/90 (7.8%) pre-treatment isolates had > 10% ring survival rates against DHA. Of the four isolates (2 RSA positive and 2 RSA negative) with high genomic coverage, P. falciparum (Pf) kelch 13 K188* and Pfcoronin V424I mutations were only present in the two RSA positive isolates with > 10% ring survival rates.

CONCLUSIONS

The observed low proportion of participants with day-3 post-treatment parasitaemia is consistent with rapid ART clearance. However, the increased rates of survival observed in the ex vivo RSA against DHA, maybe a pointer of an early start of ART tolerance. Furthermore, the role of two novel mutations in PfK13 and Pfcoronin genes, harboured by the two RSA positive isolates that had high ring survival in the present study, remains to be elucidated.

Nosocomial Malaria Transmissions Resolved by Genomic Analyses-A Retrospective Case Report Study in France: 2007-2021

BACKGROUND

Exposure of blood to malaria parasites can lead to infection even in the absence of the mosquito vector. During a stay in a healthcare facility, accidental inoculation of the skin with blood from a malaria patient might occur, referred to as nosocomial malaria.

METHODS

Between 2007 and 2021, we identified 6 autochthonous malaria cases that occurred in different French hospitals, originating from nosocomial transmission and imported malaria cases being the infection source. Four cases were observed during the coronavirus disease 2019 pandemic. The genetic relatedness between source and nosocomial infections was evaluated by genome-wide short tandem repeats (STRs) and single-nucleotide polymorphisms (SNPs).

RESULTS

None of the patients with autochthonous malaria had travel history to an endemic area nor had been transfused. For each case, both the source and recipient patients stayed a few hours in the same ward. After diagnosis, autochthonous cases were treated with antimalarials and all recovered except 1. Genetically, each pair of matched source/nosocomial parasite infections showed <1% of different STRs and <6.9% (<1.5% for monoclonal infections) of different SNPs. Similar levels of genetic differences were obtained for parasite DNA samples that were independently sequenced twice as references of identical infections. Parasite phylogenomics were consistent with travel information reported by the source patients.

CONCLUSIONS

Our study demonstrates that genomics analyses may resolve nosocomial malaria transmissions, despite the uncertainty regarding the modes of contamination. Nosocomial transmission of potentially life-threatening parasites should be taken into consideration in settings or occasions where compliance with universal precautions is not rigorous.

Plasmodium malariae: the persisting mysteries of a persistent parasite

Plasmodium malariae is a 'neglected malaria parasite' in as much as the amount of research conducted on it pales into insignificance when compared to that pertaining to Plasmodium falciparum and Plasmodium vivax, its more notorious and pathogenic cousins. There has, however, been an increase in interest in this parasite over the past decade. Principally, this is because of the increasing use of sensitive molecular detection techniques that have revealed a wider than previously recorded prevalence in some regions (particularly in Africa), and high numbers of chronic, asymptomatic infections.

Pf7: an open dataset of Plasmodium falciparum genome variation in 20,000 worldwide samples

We describe the MalariaGEN Pf7 data resource, the seventh release of Plasmodium falciparum genome variation data from the MalariaGEN network.  It comprises over 20,000 samples from 82 partner studies in 33 countries, including several malaria endemic regions that were previously underrepresented.  For the first time we include dried blood spot samples that were sequenced after selective whole genome amplification, necessitating new methods to genotype copy number variations.  We identify a large number of newly emerging crt mutations in parts of Southeast Asia, and show examples of heterogeneities in patterns of drug resistance within Africa and within the Indian subcontinent.  We describe the profile of variations in the C-terminal of the csp gene and relate this to the sequence used in the RTS,S and R21 malaria vaccines.  Pf7 provides high-quality data on genotype calls for 6 million SNPs and short indels, analysis of large deletions that cause failure of rapid diagnostic tests, and systematic characterisation of six major drug resistance loci, all of which can be freely downloaded from the MalariaGEN website.

Resolving drug selection and migration in an inbred South American Plasmodium falciparum population with identity-by-descent analysis

The human malaria parasite Plasmodium falciparum is globally widespread, but its prevalence varies significantly between and even within countries. Most population genetic studies in P. falciparum focus on regions of high transmission where parasite populations are large and genetically diverse, such as sub-Saharan Africa. Understanding population dynamics in low transmission settings, however, is of particular importance as these are often where drug resistance first evolves. Here, we use the Pacific Coast of Colombia and Ecuador as a model for understanding the population structure and evolution of Plasmodium parasites in small populations harboring less genetic diversity. The combination of low transmission and a high proportion of monoclonal infections means there are few outcrossing events and clonal lineages persist for long periods of time. Yet despite this, the population is evolutionarily labile and has successfully adapted to changes in drug regime. Using newly sequenced whole genomes, we measure relatedness between 166 parasites, calculated as identity by descent (IBD), and find 17 distinct but highly related clonal lineages, six of which have persisted in the region for at least a decade. This inbred population structure is captured in more detail with IBD than with other common population structure analyses like PCA, ADMIXTURE, and distance-based trees. We additionally use patterns of intra-chromosomal IBD and an analysis of haplotypic variation to explore past selection events in the region. Two genes associated with chloroquine resistance, crt and aat1, show evidence of hard selective sweeps, while selection appears soft and/or incomplete at three other key resistance loci (dhps, mdr1, and dhfr). Overall, this work highlights the strength of IBD analyses for studying parasite population structure and resistance evolution in regions of low transmission, and emphasizes that drug resistance can evolve and spread in small populations, as will occur in any region nearing malaria elimination.

Genetic mapping of determinants in drug resistance, virulence, disease susceptibility, and interaction of host-rodent malaria parasites

Genetic mapping has been widely employed to search for genes linked to phenotypes/traits of interest. Because of the ease of maintaining rodent malaria parasites in laboratory mice, many genetic crosses of rodent malaria parasites have been performed to map the parasite genes contributing to malaria parasite development, drug resistance, host immune response, and disease pathogenesis. Drs. Richard Carter, David Walliker, and colleagues at the University of Edinburgh, UK, were the pioneers in developing the systems for genetic mapping of malaria parasite traits, including characterization of genetic markers to follow the inheritance and recombination of parasite chromosomes and performing the first genetic cross using rodent malaria parasites. Additionally, many genetic crosses of inbred mice have been performed to link mouse chromosomal loci to the susceptibility to malaria parasite infections. In this chapter, we review and discuss past and recent advances in genetic marker development, performing genetic crosses, and genetic mapping of both parasite and host genes. Genetic mappings using models of rodent malaria parasites and inbred mice have contributed greatly to our understanding of malaria, including parasite development within their hosts, mechanism of drug resistance, and host-parasite interaction.

Gene Coverage Count and Classification (GC3): a locus sequence coverage assessment tool using short-read whole genome sequencing data, and its application to identify and classify histidine-rich protein 2 and 3 deletions in Plasmodium falciparum

Background

The ability of malaria rapid diagnostic tests (RDTs) to effectively detect active infections is being compromised by the presence of malaria strains with genomic deletions at the hrp2 and hrp3 loci, encoding the antigens most commonly targeted in diagnostics for Plasmodium falciparum detection. The presence of such deletions can be determined in publically available P. falciparum whole genome sequencing (WGS) datasets. A computational approach was developed and validated, termed Gene Coverage Count and Classification (GC3), to analyse genome-wide sequence coverage data and provide informative outputs to assess presence and coverage profile of a target locus in WGS data. GC3 was applied to detect deletions at hrp2 and hrp3 (hrp2/3) and flanking genes in different geographic regions and across time points.

Methods

GC3 uses Python and R scripts to extract locus read coverage metrics from mapped WGS data according to user-defined parameters and generates relevant tables and figures. GC3 was tested using WGS data for laboratory reference strains with known hrp2/3 genotypes, and its results compared to those of a hrp2/3-specific qPCR assay. Samples with at least 25% of coding region positions with zero coverage were classified as having a deletion. Publicly available sequence data was analysed and compared with published deletion frequency estimates.

Results

GC3 results matched the expected coverage of known laboratory reference strains. Agreement between GC3 and a hrp2/3-specific qPCR assay reported for 19/19 (100%) hrp2 deletions and 18/19 (94.7%) hrp3 deletions. Among Cambodian (n = 127) and Brazilian (n = 20) WGS datasets, which had not been previously analysed for hrp2/3 deletions, GC3 identified hrp2 deletions in three and four samples, and hrp3 deletions in 10 and 15 samples, respectively. Plots of hrp2/3 coding regions, grouped by year of sample collection, showed a decrease in median standardized coverage among Malawian samples (n = 150) suggesting the importance of a careful, properly controlled follow up to determine if an increase in frequency of deletions has occurred between 2007–2008 and 2014–2015. Among Malian (n = 90) samples, median standardized coverage was lower in 2002 than 2010, indicating widespread deletions present at the gene locus in 2002.

Conclusions

The GC3 tool accurately classified hrp2/3 deletions and provided informative tables and figures to analyse targeted gene coverage. GC3 is an appropriate tool when performing preliminary and exploratory assessment of locus coverage data.

Gravidity and malaria trends interact to modify P. falciparum densities and detectability in pregnancy: a 3-year prospective multi-site observational study

BACKGROUND

Low-density Plasmodium falciparum infections prevail in low transmission settings, where immunity is expected to be minimal, suggesting an immune-independent effect on parasite densities. We aimed to describe parasite densities in pregnancy, and determine how gravidity and antibody-mediated immunity affect these, during a period of declining malaria transmission in southern Mozambique.

METHODS

We documented P. falciparum infections at first antenatal care visits (n = 6471) between November 2016 and October 2019 in Ilha Josina (high-to-moderate transmission area), Manhiça (low transmission area), and Magude (pre-elimination area). Two-way interactions in mixed-effects regression models were used to assess gravidity-dependent differences in quantitative PCR-determined P. falciparum positivity rates (PfPR_qPCR) and densities, in the relative proportion of detectable infections (pDi) with current diagnostic tests (≥ 100 parasites/μL) and in antimalarial antibodies.

RESULTS

PfPR_qPCR declined from 28 to 13% in Ilha Josina and from 5-7 to 2% in Magude and Manhiça. In primigravidae, pDi was highest in Ilha Josina at the first study year (p = 0.048), which declined with falling PfPR_qPCR (relative change/year: 0.41, 95% CI [0.08; 0.73], p = 0.029), with no differences in antibody levels. Higher parasite densities in primigravidae from Ilha Josina during the first year were accompanied by a larger reduction of maternal hemoglobin levels (- 1.60, 95% CI [- 2.49; - 0.72; p < 0.001), than in Magude (- 0.76, 95% CI [- 1.51; - 0.01]; p = 0.047) and Manhiça (- 0.44, 95% CI [- 0.99; 0.10; p = 0.112). In contrast, multigravidae during the transmission peak in Ilha Josina carried the lowest pDi (p = 0.049). As PfPR_qPCR declined, geometric mean of parasite densities increased (4.63, 95% CI [1.28; 16.82], p = 0.020), and antibody levels declined among secundigravidae from Ilha Josina.

CONCLUSIONS

The proportion of detectable and clinically relevant infections is the highest in primigravid women from high-to-moderate transmission settings and decreases with declining malaria. In contrast, the falling malaria trends are accompanied by increased parasite densities and reduced humoral immunity among secundigravidae. Factors other than acquired immunity thus emerge as potentially important for producing less detectable infections among primigravidae during marked declines in malaria transmission.

Development and Optimization of a Selective Whole-Genome Amplification To Study Plasmodium ovale Spp

Since 2010, the human-infecting malaria parasite Plasmodium ovale spp. has been divided into two genetically distinct species, P. ovale wallikeri and P. ovale curtisi. In recent years, application of whole-genome sequencing (WGS) to P. ovale spp. allowed to get a better understanding of its evolutionary history and discover some specific genetic patterns. Nevertheless, WGS data from P. ovale spp. are still scarce due to several drawbacks, including a high level of human DNA contamination in blood samples, infections with commonly low parasite density, and the lack of robust in vitro culture. Here, we developed two selective whole-genome amplification (sWGA) protocols that were tested on six P. ovale wallikeri and five P. ovale curtisi mono-infection clinical samples. Blood leukodepletion by a cellulose-based filtration was used as the gold standard for intraspecies comparative genomics with sWGA. We also demonstrated the importance of genomic DNA preincubation with the endonuclease McrBC to optimize P. ovale spp. sWGA. We obtained high-quality WGS data with more than 80% of the genome covered by ≥5 reads for each sample and identified more than 5,000 unique single-nucleotide polymorphisms (SNPs) per species. We also identified some amino acid changes in pocdhfr and powdhfr for which similar mutations in P. falciparum and P. vivax are associated with pyrimethamine or cycloguanil resistance. In conclusion, we developed two sWGA protocols for P. ovale spp. WGS that will help to design much-needed large-scale P. ovale spp. population studies. IMPORTANCE Plasmodium ovale spp. has the ability to cause relapse, defined as recurring asexual parasitemia originating from liver-dormant forms. Whole-genome sequencing (WGS) data are of importance to identify putative molecular markers associated with relapse or other virulence mechanisms. Due to low parasitemia encountered in P. ovale spp. infections and no in vitro culture available, WGS of P. ovale spp. is challenging. Blood leukodepletion by filtration has been used, but no technique exists yet to increase the quantity of parasite DNA over human DNA when starting from genomic DNA extracted from whole blood. Here, we demonstrated that selective whole-genome amplification (sWGA) is an easy-to-use protocol to obtain high-quality WGS data for both P. ovale spp. species from unprocessed blood samples. The new method will facilitate P. ovale spp. population genomic studies.

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.

A GPU-accelerated compute framework for pathogen genomic variant identification to aid genomic epidemiology of infectious disease: a malaria case study

As recently demonstrated by the COVID-19 pandemic, large-scale pathogen genomic data are crucial to characterize transmission patterns of human infectious diseases. Yet, current methods to process raw sequence data into analysis-ready variants remain slow to scale, hampering rapid surveillance efforts and epidemiological investigations for disease control. Here, we introduce an accelerated, scalable, reproducible, and cost-effective framework for pathogen genomic variant identification and present an evaluation of its performance and accuracy across benchmark datasets of Plasmodium falciparum malaria genomes. We demonstrate superior performance of the GPU framework relative to standard pipelines with mean execution time and computational costs reduced by 27× and 4.6×, respectively, while delivering 99.9% accuracy at enhanced reproducibility.

Whole Genome Sequencing Contributions and Challenges in Disease Reduction Focused on Malaria

Malaria elimination remains an important goal that requires the adoption of sophisticated science and management strategies in the era of the COVID-19 pandemic. The advent of next generation sequencing (NGS) is making whole genome sequencing (WGS) a standard today in the field of life sciences, as PCR genotyping and targeted sequencing provide insufficient information compared to the whole genome. Thus, adapting WGS approaches to malaria parasites is pertinent to studying the epidemiology of the disease, as different regions are at different phases in their malaria elimination agenda. Therefore, this review highlights the applications of WGS in disease management, challenges of WGS in controlling malaria parasites, and in furtherance, provides the roles of WGS in pursuit of malaria reduction and elimination. WGS has invaluable impacts in malaria research and has helped countries to reach elimination phase rapidly by providing required information needed to thwart transmission, pathology, and drug resistance. However, to eliminate malaria in sub-Saharan Africa (SSA), with high malaria transmission, we recommend that WGS machines should be readily available and affordable in the region.

Worms and bugs of the gut: the search for diagnostic signatures using barcoding, and metagenomics-metabolomics

Gastrointestinal (GI) helminth infections cause significant morbidity in both humans and animals worldwide. Specific and sensitive diagnosis is central to the surveillance of such infections and to determine the effectiveness of treatment strategies used to control them. In this article, we: (i) assess the strengths and limitations of existing methods applied to the diagnosis of GI helminth infections of humans and livestock; (ii) examine high-throughput sequencing approaches, such as targeted molecular barcoding and shotgun sequencing, as tools to define the taxonomic composition of helminth infections; and (iii) discuss the current understanding of the interactions between helminths and microbiota in the host gut. Stool-based diagnostics are likely to serve as an important tool well into the future; improved diagnostics of helminths and their environment in the gut may assist the identification of biomarkers with the potential to define the health/disease status of individuals and populations, and to identify existing or emerging anthelmintic resistance.

Tenfold difference in DNA recovery rate: systematic comparison of whole blood vs. dried blood spot sample collection for malaria molecular surveillance

BACKGROUND

Molecular and genomic surveillance is becoming increasingly used to track malaria control and elimination efforts. Blood samples can be collected as whole blood and stored at - 20 °C until DNA extraction, or as dried blood spots (DBS), circumventing the need for a cold chain. Despite the wide use of either method, systematic comparisons of how the method of blood sample preservation affects the limit of detection (LOD) of molecular diagnosis and the proportion of DNA recovered for downstream applications are lacking.

METHODS

Extractions based on spin columns, magnetic beads, Tween-Chelex, and direct PCR without prior extraction were compared for whole blood and dried blood spots (DBS) using dilution series of Plasmodium falciparum culture samples. Extracted DNA was quantified by qPCR and droplet digital PCR (ddPCR).

RESULTS

DNA recovery was 5- to 10-fold higher for whole blood compared to DBS, resulting in a 2- to 3-fold lower LOD for both extraction methods compared to DBS. For whole blood, a magnetic bead-based method resulted in a DNA recovery rate of 88-98% when extracting from whole blood compared to 17-33% for a spin-column based method. For extractions from DBS, the magnetic bead-based method resulted in 8-20% DNA recovery, while the spin-column based method resulted in only 2% DNA recovery. The Tween-Chelex method was superior to other methods with 15-21% DNA recovery, and even more sensitive than extractions from whole blood samples. The direct PCR method was found to have the lowest LOD overall for both, whole blood and DBS.

CONCLUSIONS

Pronounced differences in LOD and DNA yield need to be considered when comparing prevalence estimates based on molecular methods and when selecting sampling protocols for other molecular surveillance applications.

5WBF: a low-cost and straightforward whole blood filtration method suitable for whole-genome sequencing of Plasmodium falciparum clinical isolates

Background

Whole-genome sequencing (WGS) is becoming increasingly helpful to assist malaria control programmes. A major drawback of this approach is the large amount of human DNA compared to parasite DNA extracted from unprocessed whole blood. As red blood cells (RBCs) have a diameter of about 7–8 µm and exhibit some deformability, it was hypothesized that cheap and commercially available 5 µm filters might retain leukocytes but much less of Plasmodium falciparum-infected RBCs. This study aimed to test the hypothesis that such a filtration method, named 5WBF (for 5 µm Whole Blood Filtration), may provide highly enriched parasite material suitable for P. falciparum WGS.

Methods

Whole blood was collected from five patients experiencing a P. falciparum malaria episode (ring-stage parasitaemia range: 0.04–5.5%) and from mock samples obtained by mixing synchronized, ring-stage cultured P. falciparum 3D7 parasites with uninfected human whole blood (final parasitaemia range: 0.02–1.1%). These whole blood samples (50 to 400 µL) were diluted in RPMI 1640 medium or PBS 1× buffer and filtered with a syringe connected to a 5 µm commercial filter. DNA was extracted from 5WBF-treated and unfiltered counterpart blood samples using a commercial kit. The 5WBF method was evaluated on the ratios of parasite:human DNA assessed by qPCR and by sequencing depth and percentages of coverage from WGS data (Illumina NextSeq 500). As a comparison, the popular selective whole-genome amplification (sWGA) method, which does not rely on blood filtration, was applied to the unfiltered counterpart blood samples.

Results

After applying 5WBF, qPCR indicated an average of twofold loss in the amount of parasite template DNA (Pf ARN18S gene) and from 4096- to 65,536-fold loss of human template DNA (human β actin gene). WGS analyses revealed that > 95% of the  parasite nuclear and organellar genomes were all covered at ≥ 10× depth for all samples tested. In sWGA counterparts, the organellar genomes were poorly covered and from 47.7 to 82.1% of the nuclear genome was covered at ≥ 10× depth depending on parasitaemia. Sequence reads were homogeneously distributed across gene sequences for 5WBF-treated samples (n = 5460 genes; mean coverage: 91×; median coverage: 93×; 5th percentile: 70×; 95th percentile: 103×), allowing the identification of gene copy number variations such as for gch1. This later analysis was not possible for sWGA-treated samples, as a much more heterogeneous distribution of reads across gene sequences was observed (mean coverage: 80×; median coverage: 51×; 5th percentile: 7×; 95th percentile: 245×).

Conclusions

The novel 5WBF leucodepletion method is simple to implement and based on commercially available, standardized 5 µm filters which cost from 1.0 to 1.7€ per unit depending on suppliers. 5WBF permits extensive genome-wide analysis of P. falciparum ring-stage isolates from minute amounts of whole blood even with parasitaemias as low as 0.02%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12936-022-04073-1.

Plasmodium falciparum genomic surveillance reveals spatial and temporal trends, association of genetic and physical distance, and household clustering

Molecular epidemiology using genomic data can help identify relationships between malaria parasite population structure, malaria transmission intensity, and ultimately help generate actionable data to assess the effectiveness of malaria control strategies. Genomic data, coupled with geographic information systems data, can further identify clusters or hotspots of malaria transmission, parasite genetic and spatial connectivity, and parasite movement by human or mosquito mobility over time and space. In this study, we performed longitudinal genomic surveillance in a cohort of 70 participants over four years from different neighborhoods and households in Thiès, Senegal—a region of exceptionally low malaria transmission (entomological inoculation rate less than 1). Genetic identity (identity by state, IBS) was established using a 24-single nucleotide polymorphism molecular barcode, identity by descent was calculated from whole genome sequence data, and a hierarchical Bayesian regression model was used to establish genetic and spatial relationships. Our results show clustering of genetically similar parasites within households and a decline in genetic similarity of parasites with increasing distance. One household showed extremely high diversity and warrants further investigation as to the source of these diverse genetic types. This study illustrates the utility of genomic data with traditional epidemiological approaches for surveillance and detection of trends and patterns in malaria transmission not only by neighborhood but also by household. This approach can be implemented regionally and countrywide to strengthen and support malaria control and elimination efforts.

Preservation and Extraction of Malaria Parasite DNA from Dried Blood Spots

Molecular studies related to diagnosis and research rely on collection of blood samples and extraction of high-quality DNA. In Africa, where the populations carried 94% of the total burden of cases and deaths due to malaria in 2019, collection of samples is often challenged by remote study areas and lack of a cold chain to transport and store samples. Collection of blood on filter paper is a technique that is less invasive and has simpler requirements regarding training of staff, storage, and transport of samples than collection of venous blood samples. Dried blood spots (DBS) are therefore commonly used in many research projects. However, DNA quality can be affected by duration and conditions of storage. The quality of the DNA for molecular analyses also depends on a DNA extraction methodology that provides high-quality DNA with high purity and yield. Several protocols for DNA extraction have been described, and many comparative studies have analyzed and optimized the different methodologies to find an alternative to the more costly commercial extraction kits. This chapter describes recommendations for storage and preservation of DBS, and a Chelex-based protocol for extraction of DNA from DBS.

Malaria protection due to sickle haemoglobin depends on parasite genotype

Host genetic factors can confer resistance against malaria¹, raising the question of whether this has led to evolutionary adaptation of parasite populations. Here we searched for association between candidate host and parasite genetic variants in 3,346 Gambian and Kenyan children with severe malaria caused by Plasmodium falciparum. We identified a strong association between sickle haemoglobin (HbS) in the host and three regions of the parasite genome, which is not explained by population structure or other covariates, and which is replicated in additional samples. The HbS-associated alleles include nonsynonymous variants in the gene for the acyl-CoA synthetase family member^2–4PfACS8 on chromosome 2, in a second region of chromosome 2, and in a region containing structural variation on chromosome 11. The alleles are in strong linkage disequilibrium and have frequencies that covary with the frequency of HbS across populations, in particular being much more common in Africa than other parts of the world. The estimated protective effect of HbS against severe malaria, as determined by comparison of cases with population controls, varies greatly according to the parasite genotype at these three loci. These findings open up a new avenue of enquiry into the biological and epidemiological significance of the HbS-associated polymorphisms in the parasite genome and the evolutionary forces that have led to their high frequency and strong linkage disequilibrium in African P. falciparum populations.

A strong association has been found between three regions of the Plasmodium falciparum genome and sickle haemoglobin in children with severe malaria, suggesting parasites have adapted to overcome natural host immunity.

Characterizing the genomic variation and population dynamics of Plasmodium falciparum malaria parasites in and around Lake Victoria, Kenya

Characterising the genomic variation and population dynamics of Plasmodium falciparum parasites in high transmission regions of Sub-Saharan Africa is crucial to the long-term efficacy of regional malaria elimination campaigns and eradication. Whole-genome sequencing (WGS) technologies can contribute towards understanding the epidemiology and structural variation landscape of P. falciparum populations, including those within the Lake Victoria basin, a region of intense transmission. Here we provide a baseline assessment of the genomic diversity of P. falciparum isolates in the Lake region of Kenya, which has sparse genetic data. Lake region isolates are placed within the context of African-wide populations using Illumina WGS data and population genomic analyses. Our analysis revealed that P. falciparum isolates from Lake Victoria form a cluster within the East African parasite population. These isolates also appear to have distinct ancestral origins, containing genome-wide signatures from both Central and East African lineages. Known drug resistance biomarkers were observed at similar frequencies to those of East African parasite populations, including the S160N/T mutation in the pfap2mu gene, which has been associated with delayed clearance by artemisinin-based combination therapy. Overall, our work provides a first assessment of P. falciparum genetic diversity within the Lake Victoria basin, a region targeting malaria elimination.

Genetic surveillance in the Greater Mekong subregion and South Asia to support malaria control and elimination

Background

National Malaria Control Programmes (NMCPs) currently make limited use of parasite genetic data. We have developed GenRe-Mekong, a platform for genetic surveillance of malaria in the Greater Mekong Subregion (GMS) that enables NMCPs to implement large-scale surveillance projects by integrating simple sample collection procedures in routine public health procedures.

Methods

Samples from symptomatic patients are processed by SpotMalaria, a high-throughput system that produces a comprehensive set of genotypes comprising several drug resistance markers, species markers and a genomic barcode. GenRe-Mekong delivers Genetic Report Cards, a compendium of genotypes and phenotype predictions used to map prevalence of resistance to multiple drugs.

Results

GenRe-Mekong has worked with NMCPs and research projects in eight countries, processing 9623 samples from clinical cases. Monitoring resistance markers has been valuable for tracking the rapid spread of parasites resistant to the dihydroartemisinin-piperaquine combination therapy. In Vietnam and Laos, GenRe-Mekong data have provided novel knowledge about the spread of these resistant strains into previously unaffected provinces, informing decision-making by NMCPs.

Conclusions

GenRe-Mekong provides detailed knowledge about drug resistance at a local level, and facilitates data sharing at a regional level, enabling cross-border resistance monitoring and providing the public health community with valuable insights. The project provides a rich open data resource to benefit the entire malaria community.

Funding

The GenRe-Mekong project is funded by the Bill and Melinda Gates Foundation (OPP11188166, OPP1204268). Genotyping and sequencing were funded by the Wellcome Trust (098051, 206194, 203141, 090770, 204911, 106698/B/14/Z) and Medical Research Council (G0600718). A proportion of samples were collected with the support of the UK Department for International Development (201900, M006212), and Intramural Research Program of the National Institute of Allergy and Infectious Diseases.

Advances and opportunities in malaria population genomics

Almost 20 years have passed since the first reference genome assemblies were published for Plasmodium falciparum, the deadliest malaria parasite, and Anopheles gambiae, the most important mosquito vector of malaria in sub-Saharan Africa. Reference genomes now exist for all human malaria parasites and nearly half of the ~40 important vectors around the world. As a foundation for genetic diversity studies, these reference genomes have helped advance our understanding of basic disease biology and drug and insecticide resistance, and have informed vaccine development efforts. Population genomic data are increasingly being used to guide our understanding of malaria epidemiology, for example by assessing connectivity between populations and the efficacy of parasite and vector interventions. The potential value of these applications to malaria control strategies, together with the increasing diversity of genomic data types and contexts in which data are being generated, raise both opportunities and challenges in the field. This Review discusses advances in malaria genomics and explores how population genomic data could be harnessed to further support global disease control efforts.

Failure of rapid diagnostic tests in Plasmodium falciparum malaria cases among travelers to the UK and Ireland: Identification and characterisation of the parasites

OBJECTIVES

Our objective was to systematically investigate false-negative histidine-rich protein 2 rapid diagnostic tests (HRP2-RDT) in imported Plasmodium falciparum malaria cases from travelers to the UK and the Republic of Ireland (RoI).

METHODS

Five imported malaria cases in travellers returning to the UK and RoI from East Africa were reported to the PHE Malaria Reference Laboratory as negative according to histidine-rich protein (HRP2)-RDT. The cases were systematically investigated using microscopic, RDT, molecular, genomic, and in in vitro approaches.

RESULTS

In each case, HRP2-RDT was negative, whereas microscopy confirmed the presence of P. falciparum. Further analysis revealed that the genes encoding HRP2 and HRP3 were deleted in three of the five cases. Whole-genome sequencing in one of these isolates confirmed deletions in P. falciparum chromosomes 8 and 13. Our study produced evidence that the fourth case, which had high parasitemia at clinical presentation, was a rare example of antigen saturation ('prozone-like effect'), leading to a false negative in the HRP2-RDT, while the fifth case was due to low parasitemia.

CONCLUSIONS

False-negative HRP2-RDT results with P. falciparum are concerning. Our findings emphasise the necessity of supporting the interpretation of RDT results with microscopy, in conjunction with clinical observations, and sets out a systematic approach to identifying parasites carrying pfhrp2 and pfhrp3 deletions.

Whole-genome analysis of Malawian Plasmodium falciparum isolates identifies possible targets of allele-specific immunity to clinical malaria

Individuals acquire immunity to clinical malaria after repeated Plasmodium falciparum infections. Immunity to disease is thought to reflect the acquisition of a repertoire of responses to multiple alleles in diverse parasite antigens. In previous studies, we identified polymorphic sites within individual antigens that are associated with parasite immune evasion by examining antigen allele dynamics in individuals followed longitudinally. Here we expand this approach by analyzing genome-wide polymorphisms using whole genome sequence data from 140 parasite isolates representing malaria cases from a longitudinal study in Malawi and identify 25 genes that encode possible targets of naturally acquired immunity that should be validated immunologically and further characterized for their potential as vaccine candidates.