Genetic analysis reveals unique characteristics of Plasmodium falciparum parasite populations in Haiti

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.

Plasmodium falciparum pfhrp2 and pfhrp3 Gene Deletions and Relatedness to Other Global Isolates, Djibouti, 2019-2020

Deletions of pfhrp2 and paralogue pfhrp3 (pfhrp2/3) genes threaten Plasmodium falciparum diagnosis by rapid diagnostic test. We examined 1,002 samples from suspected malaria patients in Djibouti City, Djibouti, to investigate pfhrp2/3 deletions. We performed assays for Plasmodium antigen carriage, pfhrp2/3 genotyping, and sequencing for 7 neutral microsatellites to assess relatedness. By PCR assay, 311 (31.0%) samples tested positive for P. falciparum infection, and 296 (95.2%) were successfully genotyped; 37 (12.5%) samples were pfhrp2+/pfhrp3+, 51 (17.2%) were pfhrp2+/pfhrp3-, 5 (1.7%) were pfhrp2-/pfhrp3+, and 203 (68.6%) were pfhrp2-/pfhrp3-. Histidine-rich protein 2/3 antigen concentrations were reduced with corresponding gene deletions. Djibouti P. falciparum is closely related to Ethiopia and Eritrea parasites (pairwise G_ST 0.68 [Ethiopia] and 0.77 [Eritrea]). P. falciparum with deletions in pfhrp2/3 genes were highly prevalent in Djibouti City in 2019-2020; they appear to have arisen de novo within the Horn of Africa and have not been imported.

Evaluation of a parasite-density based pooled targeted amplicon deep sequencing (TADS) method for molecular surveillance of Plasmodium falciparum drug resistance genes in Haiti

Sequencing large numbers of individual samples is often needed for countrywide antimalarial drug resistance surveillance. Pooling DNA from several individual samples is an alternative cost and time saving approach for providing allele frequency (AF) estimates at a population level. Using 100 individual patient DNA samples of dried blood spots from a 2017 nationwide drug resistance surveillance study in Haiti, we compared codon coverage of drug resistance-conferring mutations in four Plasmodium falciparum genes (crt, dhps, dhfr, and mdr1), for the same deep sequenced samples run individually and pooled. Samples with similar real-time PCR cycle threshold (Ct) values (+/- 1.0 Ct value) were combined with ten samples per pool. The sequencing success for samples in pools were higher at a lower parasite density than the individual samples sequence method. The median codon coverage for drug resistance-associated mutations in all four genes were greater than 3-fold higher in the pooled samples than in individual samples. The overall codon coverage distribution for pooled samples was wider than the individual samples. The sample pools with < 40 parasites/μL blood showed more discordance in AF calls for dhfr and mdr1 between the individual and pooled samples. This discordance in AF estimation may be due to low amounts of parasite DNA, which could lead to variable PCR amplification efficiencies. Grouping samples with an estimated ≥ 40 parasites/μL blood prior to pooling and deep sequencing yielded the expected population level AF. Pooling DNA samples based on estimates of > 40 parasites/μL prior to deep sequencing can be used for rapid genotyping of a large number of samples for these four genes and possibly other drug resistant markers in population-based studies. As Haiti is a low malaria transmission country with very few mixed infections and continued chloroquine sensitivity, the pooled sequencing approach can be used for routine national molecular surveillance of resistant parasites.