Population Genomics of Plasmodium malariae from Four African Countries
======================================================================

* Zachary R. Popkin-Hall
* Kelly Carey-Ewend
* Farhang Aghakhanian
* Eniyou C. Oriero
* Misago D. Seth
* Melchior M. Kashamuka
* Billy Ngasala
* Innocent M. Ali
* Eric Sompwe Mukomena
* Celine I. Mandara
* Oksana Kharabora
* Rachel Sendor
* Alfred Simkin
* Alfred Amambua-Ngwa
* Antoinette Tshefu
* Abebe A. Fola
* Deus S. Ishengoma
* Jeffrey A. Bailey
* Jonathan B. Parr
* Jessica T. Lin
* Jonathan J. Juliano

## Abstract

*Plasmodium malariae* is geographically widespread but neglected and may become more prevalent as *P. falciparum* declines. We completed the largest genomic study of African *P. malariae* to-date by performing hybrid capture and sequencing of 77 isolates from Cameroon (n=7), the Democratic Republic of the Congo (n=16), Nigeria (n=4), and Tanzania (n=50) collected between 2015 and 2021. There is no evidence of geographic population structure. Nucleotide diversity was significantly lower than in co-localized *P. falciparum* isolates, while linkage disequilibrium was significantly higher. Genome-wide selection scans identified no erythrocyte invasion ligands or antimalarial resistance orthologs as top hits; however, targeted analyses of these loci revealed evidence of selective sweeps around four erythrocyte invasion ligands and six antimalarial resistance orthologs. Demographic inference modeling suggests that African *P. malariae* is recovering from a bottleneck. Altogether, these results suggest that *P. malariae* is genomically atypical among human *Plasmodium* spp. and panmictic in Africa.

*Plasmodium malariae* is a neglected malaria parasite species with a broad but irregular global distribution1 and the ability to cause persistent infections2. Parasite densities are typically lower than *Plasmodium falciparum*3 and most infections involve multiple species3. While *P. malariae* causes less severe clinical disease than *P. falciparum*, it can still be deadly owing to severe complications such as glomerulonephritis and/or anemia4–6. Prevalence estimates in sub-Saharan Africa vary, but may rival *P. falciparum* in some settings3,7–11, with a particularly high contribution to malaria morbidity during the dry season in areas with seasonal transmission12. Evidence from Tanzania suggests that as control measures reduce *P. falciparum* cases, *P. malariae* may become more prevalent13, as has occurred with *P. knowlesi* in Malaysia14 and *P. vivax* in the Solomon Islands15 and elsewhere16.

*P. malariae* is most closely related to non-human primate malarias found throughout African apes, as well as *Plasmodium brasilianum*, which may represent a recent anthroponosis17. These species form a distinct clade from other human *Plasmodium* spp.17. While the *P. malariae* genome is incomplete, it shares many 1-1 orthologous genes with *P. falciparum*, including genes that are putatively implicated in antimalarial resistance and erythrocyte invasion based on their *P. falciparum* orthologs18. Though genetic studies of *P. malariae* are extremely limited and only one incomplete reference genome is available, previous analyses suggest that it experienced a bottleneck after spilling over from non-human apes to humans17.

To date, genomic analysis of *P. malariae* has been limited. A recent microsatellite study of 75 *P. malariae* isolates from seven African countries identified no geographic structure, high diversity among the microsatellite markers, and strong linkage disequilibrium (LD)19. Microsatellite markers also showed higher diversity in African *P. malariae* than South American and Asian populations19. A single previous whole-genome sequencing study of *P. malariae* leveraged selective whole-genome amplification of 18 isolates from sub-Saharan Africa and Thailand. This study identified that the Thai isolates clustered independently from the African isolates in a maximum-likelihood phylogeny20. The study also identified mutations in the putative antimalarial resistance genes dihydrofolate reductase (*pmdhfr*), dihydropteroate synthase (*pmdhps*), and multridrug resistance protein 1 (*pmmdr1*)20.

To improve understanding of the population genetics and demographic history of *P. malariae*, we conducted the largest genomic analysis of *P. malariae* in Africa to-date, incorporating 77 whole genomes generated using a custom hybrid capture protocol. These 77 isolates span four high-transmission African countries: Cameroon, the Democratic Republic of the Congo, Nigeria, and Tanzania. These genomes were used to characterize the *P. malariae* population in Africa, and compared to *P. falciparum* from similar geographic regions.

## Results

### High sequencing coverage of African *P. malariae* isolates

Of 81 genomic DNA samples that underwent hybrid capture enrichment and sequencing, 77 yielded usable sequences and were used for downstream genomic analysis. These 77 samples include 7 samples from Cameroon, 16 from the Democratic Republic of the Congo, 4 from Nigeria, and 50 from Tanzania. Tanzanian samples spanned 14 regions across Tanzania.Enrichment was overall highly successful, but a lower proportion of *P. malariae* reads were extracted from lower density samples (**Supplemental Figure 1**). In the 77 samples that were successfully enriched and sequenced, sequencing coverage was high, with the mean coverage across all chromosomes surpassing 50-fold (**Supplemental Figure 2**). Forty samples (51.9%) had >50X coverage across all chromosomes, 58 samples (75.3%) had >20X coverage across all chromosomes, and 67 samples (87.0%) had >10X coverage across all chromosomes. A total of 178,179 high-quality SNPs were identified following quality and missingness filtering and repeat masking (see **Methods**).

### Low complexity of infection

Complexity of infection (COI) was estimated using coiaf21 and was compared to that of 662 geographically matched *P. falciparum* isolates from the publicly available MalariaGEN *Pf7* dataset 22 **(Figure 1)**. The majority (92.2%, n = 71) of *P. malariae* isolates were monoclonal (COI= 1). Among the remaining six polyclonal samples, all but one contained two clones, while the remaining sample was estimated to contain three clones. By contrast, polyclonal infections were significantly (ANOVA F = 12.5, p < 0.001, df = 1) more common in *P. falciparum*. Of 662 geographically matched *P. falciparum* isolates, only 384 (58.0%) were monoclonal, while the remaining 278 (42.0%) were polyclonal. Infection by two clones was the most common form of *P. falciparum* polyclonal infection (71.6%, n = 199), similar to *P. malariae*. While country of origin was significantly associated with COI (F = 2.90, p = 0.034, df = 3), the interaction between species and country was not (F = 1.08, p = 0.340, df = 2). A Tukey post-hoc test identified no significant differences in COI by country within species.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/09/2024.09.07.24313132/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/F1)

Figure 1. Complexity of infection in A) *P. malariae* by country and B) geographically matched *P. falciparum* and *P. malariae* isolates overall.
COI values are significantly lower in *P. malariae* than *P. falciparum* (ANOVA F = 12.5, p < 0.001, df = 1), but there is no significant variation by country within species.

### Low nucleotide diversity

Among 3,763 candidate one-to-one orthologs between *P. malariae* and *P. falciparum*, 1,377 ortholog pairs were retained for nucleotide diversity (π) analysis following masking for coverage and repetitive regions (see **Methods**). The average π in *P. malariae* was 1.72 x 10-4, significantly lower (t = -113, p < 0.001, df = 5340) than the average π in *P. falciparum*, which was 6.11 x 10-3 (**Figure 2A**). SNP density was also lower in *P. malariae* than *P. falciparum*(**Supplemental Figure 3**).

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/09/2024.09.07.24313132/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/F2)

Figure 2. A) Nucleotide diversity (π) of orthologous genes among *P. malariae* and *P. falciparum* isolates.
A log-transformed boxplot of π for each gene is shown for each of the 1,377 orthologs retained after masking. 68 orthologs where missing data precluded π calculation for *P. malariae* are not shown. Boxes represent the 25th, 50th, and 75th percentiles, with outliers represented by dots. The difference in π between the species is highly significant (t= -113, p < 0.001, df = 5340). **B) Linkage disequilibrium (LD) decay in *P. falciparum and P. malariae***. R2 values were calculated in PLINK for each distance. LD is higher for *P. malariae* compared to *P. falciparum*, with rapid decay of LD in *P. falciparum* over short distances. The LD decay difference is highly significant (t = -108, p < 0.001, df = 73,920). The insert shows LD from 2 to 10bp, demonstrating that linkage in *P. falciparum* is higher than in *P. malariae* only at very short distances (<8 bp).

### High linkage disequilibrium (LD)

LD was calculated in PLINK23 for both *P. malariae* and the geographically matched *P. falciparum* isolates. Mean R2 values supporting linkage were plotted by distance between base pairs for both species (**Figure 2B**). While the mean R2 is higher in *P. falciparum* for distances < 5 bp (see **Figure 2B inset**), it is consistently higher in *P. malariae* at distances ≥ 5 bp. This difference in LD between the two species is highly significant (t = -108, p < 0.001, df = 73,920).

### No evidence of population structure or geographic differentiation

Principal component analysis (PCA) of the 71 monoclonal *P. malariae* isolates shows no clear structure and no clear relationship between geography and genetic relatedness among isolates, as well as low explanatory value, likely pointing to the lack of genetic variation (**Figure 3**). Discriminant PCA (DAPC) identified seven population clusters, three clusters with multiple samples and four singletons, as the most likely possibility according to the Bayes Information Criterion, but these putative clusters do not correspond to geography (**Supplemental Figure 4**). ADMIXTURE calculated a lower cross-validation error for K=2 (0.244) than K=1 (0.260), but again the putative populations do not correspond to geography (**Supplemental Figure 5**). Finally, a genetic distance matrix identified three distinct clusters, but these clusters also do not correspond to geography and are not borne out in a maximum-likelihood phylogeny(**Supplemental Figure 6**).

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/09/2024.09.07.24313132/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/F3)

Figure 3. Principal component analysis of monoclonal *P. malariae* isolates.
71 monoclonal *P. malariae* isolates and 178,036 biallelic SNPs are included. The first two principal components (percent of total variation explained) are depicted with isolates colored by country of origin (Cameroon n = 6, DRC n = 16, Nigeria n = 3, Tanzania n = 45).

### *P. malariae* is recovering from a genetic bottleneck

In the absence of compelling evidence for population structure within African *P. malariae* populations, we performed demographic inference modeling to identify the most likely history of this population. Of the five models tested (**Figure 4**, the “three epoch” model, indicating recovery from a bottleneck, was the best fit for African *P. malariae* based on both log likelihood (LL) and the composite-likelihood Akaike Information Criterion (CL-AIC) (**Supplemental Table 1**).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/09/2024.09.07.24313132/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/F4)

Figure 4. Schematics of demographic models tested for goodness of fit.
Plots are organized horizontally, with time in generations on the X-axis and width of colored shape corresponding to effective population size (*N*e). *TA* indicates ancient population while *TC* indicates contemporary population. **A)** Standard neutral model for one population, with no change in population size. **B)** Growth model where population growth begins at time *T*. **C) “**Bottlegrowth” model where an instantaneous size change is followed by exponential time growth at time *T*. **D)** Two epoch model where an instantaneous size change occurs at time *T* followed by a constant *Ne*. **E)** Three epoch model where an instantaneous size change occurs prior to time *T*, with *TB* corresponding to the length of the bottleneck, and bottleneck recovery begins after time *T*, with the time since bottleneck recovery represented by *TF*. Among the five models tested, the three epoch model was the best fit (LL = -514, CL-AIC = -59,663).

### Selective sweeps in orthologs of blood-stage vaccine target and putative antimalarial resistance genes

Genome-wide scans to detect selection identified numerous genes across the genome substantially deviating from neutral expectations. Genome-wide Tajima’s D scans across both genes and exons identified a negative skew (average D = -1.03 and -0.99, respectively), consistent with expectations of population expansion following a bottleneck, or potentially indicating directional selection and/or sweep. However, contrary to our expectations, none of the top hits identified by genome-wide *n*SL or Tajima’s D scans were red blood cell (RBC) invasion ligands or putative antimalarial resistance genes, other than *P. malariae myosin A* (*myoA*)24, which has a Tajima’s D value < -2 when scanning across exons, consistent with directional selection (**Supplemental Figure 7, Supplemental Figure 8, Supplemental Table 2, Supplemental Table 3, Supplemental Table 4**).

As such, we adopted a candidate gene approach to scan specific genes of interest due to their putative role in RBC invasion or antimalarial resistance. These genes included apical membrane antigen 1 (*ama1*), chloroquine resistance transporter (*crt*), circumsporozoite protein (*csp*), bifunctional dihydrofolate reductase-thymidylate synthase (*dhfr-ts*), hydroxymethyldihydropterin pyrophosphokinase-dihydropteroate synthase (*pppk-dhps*), Kelch13 (*K13*), liver surface antigen 1 (*lsa1*), multidrug resistance protein 1 (*mdr1*), multidrug resistance protein 2 (*mdr2*), multidrug resistance-associated protein 1 (*mrp1*), multidrug resistance-associated protein 2 (*mrp2*), merozoite surface protein 1 (*msp1*), ookinete surface protein P25 (*P25*), 6-cysteine protein P48/45 (*P48/45*), and thrombospondin-related anonymous protein (*TRAP*). For each of these genes, we performed McDonald-Kreitman (MK) tests to detect evidence of directional selection, as well as all of the tests incorporated in the *DH* software package to identify evidence of sweep.

A significant MK result was found for *lsa1*, but no other genes. *lsa1* also had the second-highest direction of selection, 0.184, which is consistent with weak positive selection (**Table 1**).

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/T1)

Table 1. Output of McDonald-Kreitman Tests on Genes of Interest.
P values <0.05 are bolded.

However, significant D, H, and DH values were identified for six, nine, and ten genes, respectively (**Table 2**). Significant DH p-values consistent with selective sweep were calculated for the RBC invasion genes *ama1, lsa1, msp1*, and *trap*. DH results suggestive of sweep were also found in the putative antimalarial resistance genes *crt, dhfr-ts, pppk-dhps, mdr1, mrp1*, and *mrp2*.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/T2)

Table 2. Output of Tests for Selection on Genes of Interest.
P values <0.05 are bolded.

## Discussion

We present the largest genomic study of *P. malariae* in Africa to-date, enabling a comprehensive analysis of this neglected malaria pathogen in the four high-transmission countries of Cameroon, DRC, Nigeria, and Tanzania. Using hybrid-capture enrichment, we successfully sequenced 77 samples, 71 of which were monoclonal and therefore included in further, more comprehensive analyses. Our study finds that *P. malariae* population genomics differ substantially from those of other human malaria species, including *P. falciparum* isolates from the same countries, with markedly low nucleotide diversity, no geographic population structure, high linkage disequilibrium, and no erythrocyte invasion ligands or antimalarial resistance gene orthologs identified in genome-wide selection scans. This suggests *P. malariae* may rely on markedly different strategies to persist in human populations, perhaps related to its ability to persist in human hosts undetected for months to years.

As has been seen in other non-falciparum species, COI was significantly lower in *P. malariae* than in *P. falciparum* isolates from the same countries. This trend was also observed in African *P. ovale* spp. infections25, and may reflect the lower prevalence and lighter transmission intensity of these non-falciparum species. In *P. malariae*, low COI may also correspond to its low genetic diversity and/or frequent transmission of the same clones. It is also notable that *P. malariae* has most often been detected as part of mixed infections with *P. falciparum*7–10.

In contrast to *P. ovale curtisi*25 and *P. vivax*26, *P. malariae* exhibits significantly lower nucleotide diversity than *P. falciparum* when comparing 1-1 orthologs. While a similar trend is also observed in *P. ovale wallikeri*25, the magnitude of the difference in *P. malariae* is much greater. The extremely low nucleotide diversity within *P. malariae* is likely a relic of the relatively recent bottleneck that occurred following zoonosis from the non-human primates where *P. malariae* originated17 and corresponds to a low effective population size (*Ne*). The lingering impacts of this bottleneck and corresponding low *Ne* are also reflected in the apparent lack of geographic population structure within African *P. malariae*. That lack of geographic population structure is supported by the failure of multiple methods to identify a consistent and parsimonious trend, including PCA, DAPC, ADMIXTURE, and maximum likelihood phylogenetics. Further support comes from the high linkage disequilibrium observed across the *P. malariae* genome; this delayed decay of LD over distances in the *P. malariae* genome suggests fewer recombination events between diverse strains over time.

While our analyses provide strong support for *P. malariae* experiencing a recent bottleneck, they also indicate that it is currently in recovery. The average genome-wide Tajima’s D value is negative, which is consistent with recovery from a bottleneck and/or selective sweep. We find significant evidence of selective sweep in genes of interest, as well as strong evidence from demographic inference modeling that *P. malariae* is in the recovery phase, with the three-epoch model having the best fit. Specifically, within *P. malariae*, that model could imply 1) an ancient population of non-human primate parasites that 2) experienced a severe bottleneck due to spilling over to human hosts and finally 3) began to recover from that bottleneck. While the African *P. malariae* population is still low in diversity and lacking in population structure, the model and our results suggest that *Ne* is on an upward trajectory.

*P. malariae* is highly unusual among human malaria parasites in that genome-wide selection scans identified neither erythrocyte invasion ligands nor antimalarial resistance orthologs as sites of notable balancing or directional selection using Tajima’s D and *n*SL. This difference is particularly notable in contrast with African *P. ovale* spp. from some of the same study regions, where multiple such genes had some of the highest values for both statistics25. However, the candidate gene approach detected signatures of selection in four erythrocyte invasion ligands (*pmama1, pmlsa1, pmmsp1*, and *pmtrap*) and six antimalarial resistance orthologs (*pmcrt, pmdhfr-ts, pmpppk-dhps, pmmdr1, pmmrp1*, and *pmmrp2*). The lack of identification of these genes in genome-wide scans may be due to the bottleneck recovery, along with the remarkably high LD in *P. malariae*.

Of the genes analyzed with the candidate gene approach, only *pmlsa1* shows evidence of positive selection. This trend is unusual given that erythrocyte invasion ligands in other human *Plasmodium* spp. typically show evidence of balancing selection25,27,28. This finding may, however, reflect the susceptibility of the McDonald-Kreitman test to bias from demographic factors and slightly deleterious mutations29,30. Nonetheless, the signatures of selective sweep in *pmlsa1*, as well as signatures of positive selection at four erythrocyte invasion ligands point to a consistent trend within *P. malariae*. This sweep may have arisen as part of the process of adapting to human hosts following the zoonotic event, in which case mutations in either these genes or those with which they are linked would likely have provided a selective advantage.However, it is also possible that the sweep stems from the bottleneck itself, particularly if there has been more than one severe bottleneck31.

Selective sweep in the antimalarial resistance orthologs is more precedented25,32–34 within other *Plasmodium* spp. Given that *P. malariae* is most commonly found in mixed-species infections, it is subject to drug pressure like *P. falciparum*. However, as with the erythrocyte invasion ligands, it is also possible that sweep in these antimalarial resistance orthologs is more related to bottleneck recovery than to drug pressure.

This study is subject to several limitations. As with all *Plasmodium* genomic studies, our ability to perform robust enrichment and high-coverage sequencing was limited by parasitemia, meaning that we could only thoroughly analyze relatively high-density infections. However, the hybrid capture enrichment method is able to enrich much lower density infections than direct sequencing. Second, sampling for the parent studies was not nationally representative. While we have samples spanning 14 regions of mainland Tanzania, the geographic spread of the other three countries in this analysis is restricted to specific regions. While it is possible that we have failed to detect population structure that would only be apparent with the inclusion of samples from other regions, such as West Africa, the Horn of Africa, and Southern Africa, the observed lack of population differentiation between such disparate locations as coastal Tanzania and Nigeria suggests that such structure would be weakly geographically correlated if present at all. In addition, a previous microsatellite study that includes isolates from West Africa found the same lack of geographic separation19. Third, the incomplete nature of the available *P. malariae* genome assembly used both to design our RNA baits and to align sequencing reads, may have led us to miss variants and trends in biologically important regions of the genome which are currently not assembled in chromosomes. Fourth, the study design used for collection of these isolates varied, which may have introduced unexpected biases. Finally, the comparatively low positivity rates for *P. malariae* as compared to *P. falciparum* inherently limit which samples can be used for analysis. Nevertheless, this study remains the largest conducted to-date and enables genomic analyses across multiple geographies and populations.

This study augments our understanding of the population genomics and demographic history of the largely neglected malaria parasite, *P. malariae*, in Africa. As has been suggested by other studies, we find that *P. malariae* is a genomically atypical human malaria parasite, likely owing to the impact of a zoonotic spillover event and associated genetic bottleneck in its past. The lack of geographic population structure within Africa and the extremely low nucleotide diversity present unique challenges for genomic surveillance and molecular epidemiology, as they are likely to complicate efforts to track importation and transmission networks. While *P. malariae* is currently a relatively minor problem in Africa, there is evidence that this may change in the future13,35–38. As such, it may become desirable to tailor national drug strategies and/or design vaccines to target *P. malariae* as well as *P. falciparum*, in addition to monitoring the impact of interventions targeting *P. falciparum* on the *P. malariae* population. Such efforts will need to contend with the selective sweeps identified in many of the relevant loci. While there are limited studies of *P. malariae* genomics within Africa, there are even fewer elsewhere. Evidence of population differentiation between African and Thai isolates suggest that other malaria-endemic regions may not necessarily exhibit the same trends seen in this analysis. As such, epidemiological as well as genomic studies in Asia and the Americas are warranted.

## Methods

DNA was isolated from 16,596 blood samples from six malaria studies conducted across sub-Saharan Africa using standard Chelex methods39. Samples were either dried blood spots or whole blood samples leukodepleted at the time of collection by CF11 filtration40. DNA isolates were screened by real-time polymerase chain reaction for amplification of the *P. malariae* 18S rRNA gene (**Table 3**), as described elsewhere11. Positive samples were considered for sequencing based on geographic diversity and amplification of the 18S rRNA gene before 36 cycles.

View this table:
[Table 3.](http://medrxiv.org/content/early/2024/09/09/2024.09.07.24313132/T3)

Table 3. Details of sample collection for samples included in study.

### Library preparation and sequencing

Selected DNA isolates were fragmented enzymatically and underwent library preparation using the Twist Library Preparation EF 2.0 kit (Twist Bioscience, South San Francisco, CA). Libraries then underwent parasite DNA enrichment using a custom-designed Twist hybridization capture protocol using 280,042 RNA baits specifically designed to amplify *P. malariae* genomic DNA across the full genome without reacting with background human DNA. Enriched samples were then amplified, purified, and submitted for Illumina short-read 150bp sequencing on the NovaSeq 6000 S4-XP system with paired-end chemistry.

### Sequencing data alignment and variant calling

*fastqc* v0.12.1 ([https://www.bioinformatics.babraham.ac.uk/projects/fastqc/](https://www.bioinformatics.babraham.ac.uk/projects/fastqc/)) was used to check the quality of raw sequencing reads before trimming sequencing adapters using *Trim Galore* v0.6.7 ([https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/](https://www.bioinformatics.babraham.ac.uk/projects/trim_galore/)). Trimmed reads were then competitively aligned to the *P. malariae* (PmUG01 strain), *P. falciparum* (Pf3D7), and *Homo sapiens* (Hg38 strain) using BBSplit within BBMap v38.96 ([https://sourceforge.net/projects/bbmap/](https://sourceforge.net/projects/bbmap/)). Reads that best aligned to the *P. falciparum* and *H. sapiens* genomes were discarded, and the remaining reads were then aligned to the PmUG01 reference genome with *bwa-mem2* v2.2.146. *picard* v2.26.11 ([https://broadinstitute.github.io/picard/](https://broadinstitute.github.io/picard/)) was then used to sort and deduplicate reads before variant calling with GATK v4.5.0.047. Coverage statistics were calculated using *samtools* v1.2048,49.

Variants were called across each sample following the *GATK* best practices pipeline47. Variants were called using gVCF mode of the Haplotype Caller function within each chromosome of each individual sample before calling variants across all samples50. SNPs were then filtered with GATK if they met the following thresholds: quality by depth <2.5, Fisher strand bias >10, mapping quality >50, mapping quality rank sum <-2.5, read position rank sum <-2.5. The resulting filtered VCF was then filtered to biallelic SNPs using bcftools v1.2.049 and SNPs with greater than 20% missingness were excluded with vcftools v0.1.1551. SNPs falling within the hypervariable *PIRs* and SNPs within tandem repeats as identified by *Tandem Repeats Finder* v4.09.152 were masked in vcftools. Finally, only SNPs within the 14 chromosomes were retained for further analysis, with all SNPs on the mitochondrion, apicoplast, and extrachromosmal contigs excluded with bcftools.

### Selection of matched *P. falciparum* samples

For comparison with sympatric *Plasmodium falciparum* parasites, whole-genome sequencing data from 384 *P. falciparum* isolates from the same countries of origin as the *P. malariae* isolates (121 from Cameroon, 93 from DRC, and 170 from Tanzania -while Nigerian data were available, the Cameroonian data were more geographically proximate to the Nigerian *P. malariae* isolates) were selected from the publicly-available *Pf7* dataset22.

### Nucleotide diversity

The *OrthoMCL*53 database implemented in *PlasmoDB*54 was used to identify orthologous protein-coding genes between the *P. malariae* and *P. falciparum*, and 3,763 one-to-one orthologs were analyzed to compare nucleotide diversity between the two species. Orthologs were masked based on the masking steps outlined above and also excluded if they were found outside of the *P. falciparum* core genome. In addition, ortholog pairs were excluded unless 60% of samples of each species had ≥ 5X coverage of the gene at least every 10 base pairs. Of the 3,763 candidate orthologs, 1,377 unmasked ortholog pairs were retained for final analyses. Nucleotide diversity (π) was then calculated in vcftools over each ortholog in each species, and the distribution of π by species was compared with a Student’s t-test.

### Linkage disequilibrium and complexity of infection

Linkage disequilibrium across genomic intervals in each species was calculated using *PLINK* v1.90b7.223 with 100 bp windows and a minor allele frequency cutoff of 0.01. Complexity of infection was estimated for each isolate using *coiaf* v0.1.221 with a minor allele frequency cutoff of 0.05. Data visualizations were constructed in *R* v4.2.255.

### Population Structure and Demography

*vcfdo* ([https://github.com/IDEELResearch/vcfdo](https://github.com/IDEELResearch/vcfdo)) was used to calculate within sample allele frequency (WSAF) among monoclonal samples, then filtered to SNPs with WSAF equal to 0 or 1, as heterozygous calls in monoclonal haploid infections are likely to represent sequencing errors or paralogous misalignments. Principal component analysis (PCA) of monoclonal isolates was performed in *PLINK*. Population structure among monoclonal isolates was also assessed using *ADMIXTURE* v1.3.056, discriminant PCA in *adegenet* v2.1.1057,58, genetic distance in *fastreeR* ([https://github.com/gkanogiannis/fastreeR](https://github.com/gkanogiannis/fastreeR)), and maximum likelihood phylogeny in *RAxML Next Generation* v1.2.2 ([https://github.com/amkozlov/raxml-ng](https://github.com/amkozlov/raxml-ng)) following format conversion using *vcf2phylip* v2.8 ([https://github.com/edgardomortiz/vcf2phylip](https://github.com/edgardomortiz/vcf2phylip)).Demographic history was initially inferred using *donni*59 to determine rough estimates for one-population model parameters using the following models implemented in ∂a∂i60: “bottlegrowth\_1d”, “growth”, “snm_1d”, “two_epoch”, and “three_epoch”. The 95% confidence intervals (CI) from the *donni* output were then used as bounds for optimization in *dadi-cli*61.Because of the high linkage in *P. malariae*62, model fit was compared using the composite likelihood Akaike Information Criterion (CL-AIC)63,64, calculated in a custom R script. To perform this calculation, the Godambe Information Matrix62 script in *dadi-cli* was modified to extract the *H* and *J* matrices for the 95% CI with 0.01 step size. All converged optimizations from the best fit file for each model were then used to determine median log likelihood values and median demographic parameter estimates.

### Signatures of Selection

Among monoclonal isolates of each species, *n*SL65 was calculated at each SNP using *selscan* v2.2.066 and the SNPs with absolute values in the top 0.5% were intersected with the PmUG01 genome annotation using the *ape* v5.867, *GenomicRanges* v1.50.268, and *plyranges* v1.18.069 Rpackages. Tajima’s D was calculated in windows of 300 bp with a 10bp step size across the entire genome with *vcftools*; the top 0.5% (absolute value) of windows were analyzed as above.

To detect signs of selection within blood-stage vaccine target orthologs and antimalarial resistance gene orthologs (*ama1, crt, csp, dhfr-ts, Kelch13, lsa1, mdr1, mdr2, mrp1, mrp2, msp1, p25, p48-45, pppk-dhps, trap*), *bcftools* was used to generate fasta sequences for both coding sequences (CDS) and full gene sequences for each gene within each isolate for both *P. malariae* and *P. falciparum*, then sequences were aligned using MAFFT v7.49070. The CDS alignments were imported into DnaSP v6.12.0371, where the McDonald-Kreitman72 test for positive selection was calculated. Because the neutrality index is subject to bias, we calculated and reported direction of selection (DoS) instead73. The full sequence alignments of *P. malariae* sequences with one randomly selected *P. falciparum* outgroup sequence were analyzed with the *Readms* module of *DH* ([https://github.com/drkaizeng/publications-and-software/blob/main/dh/dh.zip](https://github.com/drkaizeng/publications-and-software/blob/main/dh/dh.zip)), which calculates Tajima’s D74, Fay and Wu’s H75, a p-value for the DH test76, and the E test76 over 10,000 coalescent simulations. The DH test is unique in its sensitivity to selective sweeps and insensitivity to demographic forces, whereas Tajima’s D and Fay and Wu’s H are both prone to demographic bias76.

## Supporting information

Supplemental Material [[supplements/313132_file02.docx]](pending:yes)

## Data Availability

Parasite sequence data is available through SRA (BioProject ID PRJNA1157442).

## Data availability

Parasite sequence data is available through SRA (BioProject ID PRJNA1157442).

## Code availability

Code used for analysis is available from GitHub at [https://github.com/IDEELResearch/PmPopGen](https://github.com/IDEELResearch/PmPopGen)

## Ethics declarations

As part of the parent studies from which samples were derived, written informed consent, assent, and/or parental consent was obtained from all participants.Institutional review board approvals were obtained for parent studies as follows: Cameroon, from the Gambia Government/MRCG Joint Ethics Committee (REF SCC1626); DRC, from the University of North Carolina at Chapel Hill (IRB#: 14-0489) and the Kinshasa School of Public Health (ESP/CE/015/014); Nigeria, from the Cross River State Health Research Ethics Committee (REC No CRSMOH/RP/REC/2017/809). For the MSMT samples, the study protocol was submitted to the Tanzanian Medical Research Coordinating Committee (MRCC) of the National Institute for Medical Research (NIMR) for review and ethical approval, which was granted. The protocol was also submitted for review and approval by the ethics committee of WHO in Geneva, Switzerland, which was also granted. All research participants were asked and provided individual consent (or assent for children aged 7–17 years of age) for their participation in the survey and biobanking for future research. For children under the legal age of adulthood in Tanzania (< 18 years), consent was obtained from a parent or guardian. An informed consent form was developed in English and translated in Kiswahili and used to obtain consent both verbally and in writing from all participants. All participants agreed and signed the consent or assent form or provided a thumbprint in conjunction with the signature of an independent witness in case the study participant was illiterate. All experiments were performed in accordance with relevant guidelines and regulations in accordance with the Declaration of Helsinki. For the GANC samples, the study protocol was approved by the National Health Research Ethics Sub-Committee (NatHREC) of the Ministry of Health, Community Development, Gender, Elderly and Children (Dar es Salaam, Tanzania) and the US Centers for Disease Control and Prevention (CDC; Atlanta, GA, USA) institutional review board.

## Funding

This work was supported, in part, by the Bill & Melinda Gates Foundation [grant number 002202]. Under the grant conditions of the Foundation, a Creative Commons Attribution 4.0 Generic License has already been assigned to the Author Accepted Manuscript version that might arise from this submission. Data collection in Geita, Tanzania was funded by USAID/PMI through Jhpiego and CDC. This study was funded in part by the National Institutes of Health (NIH) [T32AI007151 to Z.R.P-H.; T32AI070114 to R.S.; R01AI107949 and R01AI129812 to A.T.; R21 AI148579 to J.B.P. and J.T.L.; R01AI137395 and R21AI152260 to J. T. L.; R01AI132547 and K24AI134990 to J.J.J]. This work was supported through the DELTAS Africa initiative (DELGEME Grant 107740/Z/15/Z). The DELTAS Africa initiative is an independent funding scheme of the African Academy of Sciences (AAS)’s Alliance for Accelerating Excellence in Science in Africa (AESA) and supported by the New Partnership for Africa’s Development Planning and Coordinating Agency (NEPAD) with funding from the Wellcome Trust (DELGEME Grant 107740/Z/15/Z) and the UK government.

## Competing interests

JBP reports research support from Gilead Sciences, non-financial support from Abbott Laboratories, and consulting for Zymeron Corporation, all outside the scope of the current manuscript.

## Acknowledgements

Kyaw Lay Thwai, Claudia Gaither, Meredith Muller, and Srijana Chhetri provided laboratory assistance at UNC. Andrew Guinness wrote the R function to calculate CL-AIC. Professor Martin Meremikwu (University of Calabar Teaching Hospital) and Dr. Dr Tobias Apinjoh (University of Buea) provided invaluable support for field collections in Nigeria and Cameroon, respectively. The authors wish to thank participants and parents or guardians of all children who took part in the studies which were undertaken. We appreciate the support of the in-country teams, and other colleagues including community and health facility staff who took part in data collection or laboratory processing of samples.

## Footnotes

*   † Co-first authors

*   Received September 7, 2024.
*   Revision received September 7, 2024.
*   Accepted September 9, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Fuehrer, H.-P., Campino, S. & Sutherland, C. J. The primate malaria parasites Plasmodium malariae, Plasmodium brasilianum and Plasmodium ovale spp.: genomic insights into distribution, dispersal and host transitions. Malar J 21, 138 (2022).
    
    

2.  2.Sutherland, C. J. Persistent Parasitism: The Adaptive Biology of Malariae and Ovale Malaria. Trends Parasitol 32, 808–819 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2016.07.001&link_type=DOI) 

3.  3.Oriero, E. C., Amenga-Etego, L., Ishengoma, D. S. & Amambua-Ngwa, A. Plasmodium malariae, current knowledge and future research opportunities on a neglected malaria parasite species. Crit Rev Microbiol , 1–13 (2021).
    
    

4.  4.Ehrich, J. H. H. & Eke, F. U. Malaria-induced renal damage: Facts and myths. Pediatric Nephrology 22, 626–637 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00467-006-0332-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17205283&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246077100002&link_type=ISI) 

5.  5.Hedelius, R., Fletcher, J. J., Glass, W. F., Susanti, A. I. & Maguire, J. D. Nephrotic Syndrome and Unrecognized Plasmodium malariae Infection in a US Navy Sailor 14 Years After Departing Nigeria. J Travel Med 18, 288–291 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1708-8305.2011.00526.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21722243&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

6.  6.Langford, S. et al. Plasmodium malariae Infection Associated with a High Burden of Anemia: A Hospital-Based Surveillance Study. PLoS Negl Trop Dis 9, e0004195 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0004195&link_type=DOI) 

7.  7.Sendor, R. et al. Epidemiology of Plasmodium malariae and Plasmodium ovale spp. in Kinshasa Province, Democratic Republic of Congo. Nature Communications 2023 14:1 14, 1–13 (2023).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-023-37194-5&link_type=DOI) 

8.  8.Sendor, R. et al. Similar Prevalence of Plasmodium falciparum and Non-P. falciparum Malaria Infections among Schoolchildren, Tanzania. Emerg Infect Dis 29, 1143–1153 (2023).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2906.221016&link_type=DOI) 

9.  9.Popkin-Hall, Z. R. et al. Malaria Species Positivity Rates Among Symptomatic Individuals Across Regions of Differing Transmission Intensities in Mainland Tanzania. J Infect Dis (2023) doi:10.1093/infdis/jiad522.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiad522&link_type=DOI) 

10. 10.Popkin-Hall, Z. R. et al. Prevalence of non-falciparum malaria infections among asymptomatic individuals in four regions of Mainland Tanzania. Parasit Vectors 17, 1–6 (2024).
    
    

11. 11.Gumbo, A. et al. Occurrence and Distribution of Nonfalciparum Malaria Parasite Species Among Adolescents and Adults in Malawi. J Infect Dis 225, 257–268 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiab353&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34244739&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

12. 12.Greenwood, B. M. et al. Mortality and morbidity from malaria among children in a rural area of The Gambia, West Africa. Trans R Soc Trop Med Hyg 81, 478–486 (1987).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0035-9203(87)90170-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3318021&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

13. 13.Yman, V. et al. Persistent transmission of Plasmodium malariae and Plasmodium ovale species in an area of declining Plasmodium falciparum transmission in eastern Tanzania. PLoS Negl Trop Dis 13, e0007414 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0007414&link_type=DOI) 

14. 14.Cooper, D. J. et al. Plasmodium knowlesi Malaria in Sabah, Malaysia, 2015–2017: Ongoing Increase in Incidence Despite Near-elimination of the Human-only Plasmodium Species. Clinical Infectious Diseases 70, 361–367 (2020).
    
    

15. 15.Quah, Y. W. et al. Molecular epidemiology of residual Plasmodium vivax transmission in a paediatric cohort in Solomon Islands. Malar J 18, 106 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-019-2727-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30922304&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

16. 16.Price, R. N., Commons, R. J., Battle, K. E., Thriemer, K. & Mendis, K. Plasmodium vivax in the Era of the Shrinking P. falciparum Map. Trends Parasitol 36, 560–570 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pt.2020.03.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32407682&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

17. 17.Plenderleith, L. J. et al. Zoonotic origin of the human malaria parasite Plasmodium malariae from African apes. Nat Commun 13, 1868 (2022).
    
    

18. 18.Rutledge, G. G. et al. Plasmodium malariae and P. ovale genomes provide insights into malaria parasite evolution. Nature 542, 101–104 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature21038&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28117441&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

19. 19.Oriero, E. C. et al. Plasmodium malariae structure and genetic diversity in sub-Saharan Africa determined from microsatellite variants and linked SNPs in orthologues of antimalarial resistance genes. Sci Rep 12, 21881 (2022).
    
    

20. 20.Ibrahim, A. et al. Selective whole genome amplification of Plasmodium malariae DNA from clinical samples reveals insights into population structure. Sci Rep 10, 10832 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-67568-4&link_type=DOI) 

21. 21.Paschalidis, A., Watson, O. J., Aydemir, O., Verity, R. & Bailey, J. A. coiaf: Directly estimating complexity of infection with allele frequencies. PLoS Comput Biol 19, e1010247 (2023).
    
    

22. 22.Malaria GEN et al. An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples. Wellcome Open Research 2021 6:42 6, 42 (2021).
    
    

23. 23.Purcell, S. et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81, 559–575 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

24. 24.Thomas, J. C. et al. Interaction and dynamics of the Plasmodium falciparumMTIP–MyoA complex, a key component of the invasion motor in the malaria parasite. Mol Biosyst 6, 494–498 (2010).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1039/b922093c&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20174678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

25. 25.Carey-Ewend, K. et al. Population genomics of Plasmodium ovale species in sub-Saharan Africa. bioRxiv 2024.04.10.588912 (2024) doi:10.1101/2024.04.10.588912.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyNC4wNC4xMC41ODg5MTJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

26. 26.Neafsey, D. E. et al. The malaria parasite Plasmodium vivax exhibits greater genetic diversity than Plasmodium falciparum. Nat Genet 44, 1046–1050 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2373&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22863733&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

27. 27.Tapaopong, P. et al. Genetic diversity and molecular evolution of Plasmodium vivax Duffy Binding Protein and Merozoite Surface Protein-1 in northwestern Thailand. Infection, Genetics and Evolution 113, 105467 (2023).
    
    

28. 28.Naung, M. T. et al. Global diversity and balancing selection of 23 leading Plasmodium falciparum candidate vaccine antigens. PLoS Comput Biol 18, e1009801 (2022).
    
    

29. 29.Charlesworth, J. & Eyre-Walker, A. The McDonald-Kreitman test and slightly deleterious mutations. Mol Biol Evol 25, 1007–1015 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msn005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18195052&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000255758200002&link_type=ISI) 

30. 30.Eyre-Walker, A. Changing Effective Population Size and the McDonald-Kreitman Test. Genetics 162, 2017–2024 (2002).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjE2Mi80LzIwMTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8wOS8yMDI0LjA5LjA3LjI0MzEzMTMyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

31. 31.Wilson, B. A., Petrov, D. A. & Messer, P. W. Soft Selective Sweeps in Complex Demographic Scenarios. Genetics 198, 669–684 (2014).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMTk4LzIvNjY5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMDkvMjAyNC4wOS4wNy4yNDMxMzEzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. 32.Parobek, C. M. et al. Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia. Proc Natl Acad Sci U S A 113, E8096–E8105 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEzLzUwL0U4MDk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMDkvMjAyNC4wOS4wNy4yNDMxMzEzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

33. 33.Shaukat, A. et al. Selective sweep and phylogenetic models for the emergence and spread of pyrimethamine resistance mutations in Plasmodium vivax. Infection, Genetics and Evolution 68, 221–230 (2019).
    
    

34. 34.Alam, M. T. et al. Selective sweeps and genetic lineages of Plasmodium falciparum drug - resistant alleles in Ghana. Journal of Infectious Diseases 203, 220–227 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiq038&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21288822&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

35. 35.Nguiffo-Nguete, D. et al. Plasmodium malariae contributes to high levels of malaria transmission in a forest–savannah transition area in Cameroon. Parasites & Vectors 2023 16:1 16, 1–10 (2023).
    
    

36. 36.Molineaux, L., Storey, J., Cohen, J. E. & Thomas, A. A Longitudinal Study of Human Malaria in the West African Savanna in the Absence of Control Measures: Relationships between Different Plasmodium Species, in Particular P. Falciparum and P. Malariae. Am J Trop Med Hyg 29, 725–737 (1980).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiIyOS81LzcyNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

37. 37.Roucher, C., Rogier, C., Sokhna, C., Tall, A. & Trape, J. F. A 20-Year Longitudinal Study of Plasmodium ovale and Plasmodium malariae Prevalence and Morbidity in a West African Population. PLoS One 9, e87169 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0087169&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24520325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

38. 38.Hawadak, J., Dongang Nana, R. R. & Singh, V. Global trend of Plasmodium malariae and Plasmodium ovale spp. malaria infections in the last two decades (2000–2020): a systematic review and meta-analysis. Parasit Vectors 14, 1–14 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13071-021-04992-z&link_type=DOI) 

39. 39.Teyssier, N. B. et al. Optimization of whole-genome sequencing of Plasmodium falciparum from low-density dried blood spot samples. Malar J 20, 1–8 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-021-03646-w&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

40. 40.Venkatesan, M. et al. Using CF11 cellulose columns to inexpensively and effectively remove human DNA from Plasmodium falciparum-infected whole blood samples. Malar J 11, 1–7 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-11-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22212246&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

41. 41.Emerson, C. et al. Women attending antenatal care as a sentinel surveillance population for malaria in Geita region, Tanzania: feasibility and acceptability to women and providers. Malar J 22, 1–10 (2023).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-023-04445-1&link_type=DOI) 

42. 42.Gutman, J. R. et al. Using antenatal care as a platform for malaria surveillance data collection: study protocol. Malar J 22, 1–10 (2023).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-023-04445-1&link_type=DOI) 

43. 43.Markwalter, C. F. et al. Direct Comparison of Standard and Ultrasensitive PCR for the Detection of Plasmodium falciparum from Dried Blood Spots in Bagamoyo, Tanzania. Am J Trop Med Hyg 104, 1371–1374 (2021).
    
    

44. 44.Rapp, T. et al. Microheterogeneity of Transmission Shapes Submicroscopic Malaria Carriage in Coastal Tanzania. J Infect Dis 230, 485–496 (2024).
    
    

45. 45.Oriero, E. C. et al. Seroprevalence and Parasite Rates of Plasmodium malariae in a High Malaria Transmission Setting of Southern Nigeria. Am J Trop Med Hyg 103, 2208–2216 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.20-0593&link_type=DOI) 

46. 46.Md, V., Misra, S., Li, H. & Aluru, S. Efficient architecture-aware acceleration of BWA-MEM for multicore systems. Proceedings - 2019 IEEE 33rd International Parallel and Distributed Processing Symposium, IPDPS 2019 314–324 (2019) doi:10.1109/IPDPS.2019.00041.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/IPDPS.2019.00041&link_type=DOI) 

47. 47.Van der Auwera, G. A. & O’Connor, B. D. Genomics in the Cloud: Using Docker, GATK, and WDL in Terra. (O’Reilly Media, 2020).
    
    

48. 48.Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078– 2079 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btp352&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19505943&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000268808600014&link_type=ISI) 

49. 49.Danecek, P. et al. Twelve years of SAMtools and BCFtools. Gigascience 10, 1–4 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giab016&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

50. 50.Poplin, R. et al. Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv 201178 (2018) doi:10.1101/201178.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czo4OiIyMDExNzh2MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

51. 51.Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr330&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21653522&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000292778700023&link_type=ISI) 

52. 52.Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res 27, 573–580 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/27.2.573&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9862982&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078181500028&link_type=ISI) 

53. 53.Li, L., Stoeckert, C. J. & Roos, D. S. OrthoMCL: Identification of Ortholog Groups for Eukaryotic Genomes. Genome Res 13, 2178–2189 (2003).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjEzLzkvMjE3OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

54. 54.Alvarez-Jarreta, J. et al. VEuPathDB: the eukaryotic pathogen, vector and host bioinformatics resource center in 2023. Nucleic Acids Res 52, D808–D816 (2024).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkad1003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37953350&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

55. 55.R Core Team. R: A Language and Environment for Statistical Computing. Preprint at [https://www.r-project.org/](https://www.r-project.org/) (2020).
    
    

56. 56.Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res 19, 1655–1664 (2009).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjE5LzkvMTY1NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

57. 57.Jombart, T. & Bateman, A. adegenet: a R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405 (2008).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btn129&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18397895&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256169300015&link_type=ISI) 

58. 58.Jombart, T. & Ahmed, I. adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics 27, 3070–3071 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr521&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21926124&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296099300022&link_type=ISI) 

59. 59.Tran, L. N., Sun, C. K., Struck, T. J., Sajan, M. & Gutenkunst, R. N. Computationally Efficient Demographic History Inference from Allele Frequencies with Supervised Machine Learning. Mol Biol Evol 41, (2024).
    
    

60. 60.Gutenkunst, R. N., Hernandez, R. D., Williamson, S. H. & Bustamante, C. D. Inferring the Joint Demographic History of Multiple Populations from Multidimensional SNP Frequency Data. PLoS Genet 5, e1000695 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1000695&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19851460&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

61. 61.Huang, X., Struck, T. J., Davey, S. W. & Gutenkunst, R. N. dadi-cli: Automated and distributed population genetic model inference from allele frequency spectra. bioRxiv 2023.06.15.545182 (2023) doi:10.1101/2023.06.15.545182.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiYmlvcnhpdiI7czo1OiJyZXNpZCI7czoxOToiMjAyMy4wNi4xNS41NDUxODJ2MSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzA5LzIwMjQuMDkuMDcuMjQzMTMxMzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

62. 62.Coffman, A. J., Hsieh, P. H., Gravel, S. & Gutenkunst, R. N. Computationally Efficient Composite Likelihood Statistics for Demographic Inference. Mol Biol Evol 33, 591–593 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msv255&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26545922&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

63. 63.Varin, C. & Vidoni, P. A Note on Composite Likelihood Inference and Model Selection. 92, 519–528 (2005).
    
    

64. 64.Noskova, E., Ulyantsev, V., Koepfli, K. P., O’brien, S. J. & Dobrynin, P. GADMA: Genetic algorithm for inferring demographic history of multiple populations from allele frequency spectrum data. Gigascience 9, 1–18 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/gigascience/giaa056&link_type=DOI) 

65. 65.Ferrer-Admetlla, A., Liang, M., Korneliussen, T. & Nielsen, R. On Detecting Incomplete Soft or Hard Selective Sweeps Using Haplotype Structure. Mol Biol Evol 31, 1275–1291 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msu077&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24554778&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000335914400020&link_type=ISI) 

66. 66.Szpiech, Z. A. selscan 2.0: scanning for sweeps in unphased data. Bioinformatics 40, (2024).
    
    

67. 67.Paradis, E. & Schliep, K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty633&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30016406&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

68. 68.Lawrence, M. et al. Software for Computing and Annotating Genomic Ranges. PLoS Comput Biol 9, e1003118 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1003118&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23950696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

69. 69.Lee, S., Cook, D. & Lawrence, M. Plyranges: A grammar of genomic data transformation. Genome Biol 20, (2019).
    
    

70. 70.Katoh, K. & Standley, D. M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol Biol Evol 30, 772–780 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/mst010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23329690&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000317002300004&link_type=ISI) 

71. 71.Rozas, J. et al. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol Biol Evol 34, 3299–3302 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msx248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29029172&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 

72. 72.McDonald, J. H. & Kreitman, M. Adaptive protein evolution at the Adh locus in Drosophila. Nature 351, 652–654 (1991).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/351652a0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1904993&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991FT11200066&link_type=ISI) 

73. 73.Stoletzki, N. & Eyre-Walker, A. Estimation of the Neutrality Index. Mol Biol Evol 28, 63–70 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msq249&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20837603&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F09%2F2024.09.07.24313132.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285418600014&link_type=ISI) 

74. 74.Tajima, F. Statistical Method for Testing the Neutral Mutation Hypothesis by DNA Polymorphism. Genetics 123, 585–595 (1989).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6OToiMTIzLzMvNTg1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMDkvMjAyNC4wOS4wNy4yNDMxMzEzMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

75. 75.Fay, J. C. & Wu, C. I. Hitchhiking Under Positive Darwinian Selection. Genetics 155, 1405–1413 (2000).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjE1NS8zLzE0MDUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8wOS8yMDI0LjA5LjA3LjI0MzEzMTMyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

76. 76.Zeng, K., Fu, Y. X., Shi, S. & Wu, C. I. Statistical tests for detecting positive selection by utilizing high-frequency variants. Genetics 174, 1431–1439 (2006).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjE3NC8zLzE0MzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8wOS8yMDI0LjA5LjA3LjI0MzEzMTMyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==)