Lineage-informative microhaplotypes for spatio-temporal surveillance of Plasmodium vivax malaria parasites
==========================================================================================================

* Sasha V. Siegel
* Roberto Amato
* Hidayat Trimarsanto
* Edwin Sutanto
* Mariana Kleinecke
* Kathryn Murie
* Georgia Whitton
* Aimee R. Taylor
* James A. Watson
* Mallika Imwong
* Ashenafi Assefa
* Awab Ghulam Rahim
* Nguyen Hoang Chau
* Tran Tinh Hien
* Justin A Green
* Gavin Koh
* Nicholas J. White
* Nicholas Day
* Dominic P. Kwiatkowski
* Julian C. Rayner
* Ric N. Price
* Sarah Auburn

## Abstract

Challenges in understanding the origin of recurrent *Plasmodium vivax* infections constrains the surveillance of antimalarial efficacy and transmission of this neglected parasite. Recurrent infections within an individual may arise from activation of dormant liver stages (relapse), blood-stage treatment failure (recrudescence) or new inoculations (reinfection). Molecular inference of familial relatedness (identity-by-descent or IBD) based on whole genome sequence data, together with analysis of the intervals between parasitaemic episodes (“time-to-event” analysis), can help resolve the probable origin of recurrences. Whole genome sequencing of predominantly low-density *P. vivax* infections is challenging, so an accurate and scalable genotyping method to determine the origins of recurrent parasitaemia would be of significant benefit. We have developed a *P. vivax* genome-wide informatics pipeline to select specific microhaplotype panels that can capture IBD within small, amplifiable segments of the genome. Using a global set of 615 *P. vivax* genomes, we derived a panel of 100 microhaplotypes, each comprising 3-10 high frequency SNPs within <200 bp sequence windows. This panel exhibits high diversity in regions of the Asia-Pacific, Latin America and the horn of Africa (median *H*E = 0.70-0.81) and it captured 89% (273/307) of the polyclonal infections detected with genome-wide datasets. Using data simulations, we demonstrate lower error in estimating pairwise IBD using microhaplotypes, relative to traditional biallelic SNP barcodes. Our panel exhibited high accuracy in predicting the country of origin (median Matthew’s correlation coefficient >0.9 in 90% countries tested) and it also captured local infection outbreak and bottlenecking events. The informatics pipeline is available open-source and yields microhaplotypes that can be readily transferred to high-throughput amplicon sequencing assays for surveillance in malaria-endemic regions.

Key words
*   *Plasmodium vivax*
*   malaria
*   recurrence
*   relapse
*   relatedness
*   genotyping
*   identity-by-descent
*   microhaplotype
*   barcode
*   molecular surveillance

## Introduction

The malaria parasite *Plasmodium vivax* is a major public health threat affecting the poorest and most vulnerable populations in more than 49 endemic countries1. Over the past decade, enhanced malaria control efforts in areas outside of sub-Saharan Africa have achieved a marked decline in *P. falciparum* infections, but a relative rise in the proportion of *P. vivax* cases1. Several biological attributes of *P. vivax* render this species more resilient to interventions than *P. falciparum**2*. *P. vivax* forms dormant liver stage parasites (hypnozoites) that can reactivate weeks to months after initial inoculation causing recurrent episodes of malaria (relapses). A single mosquito inoculation can cause multiple relapses, sustaining transmission over extended periods. Relapses are thought to cause over 60% of clinical cases3. Knowledge of the biology and epidemiology of *P. vivax* relapse, including the underlying reactivation mechanism(s) and host and parasite determinants, are vital for informing public health strategies to combat this parasite, but our understanding of these processes is limited4. A major obstacle to our increased understanding has been the difficulty of classifying recurrent *P. vivax* infections accurately. These recurrences can arise from reinfection (new mosquito inoculations), recrudescence (blood-stage treatment failures), or relapse (reactivation of dormant liver stages). Discriminating between these causes is challenging.

Accurate methods to disentangle relapse from reinfection and recrudescence events are critical to improving our understanding of the therapeutic efficacy of current treatment regimens for *P. vivax*. Accumulating reports suggest that chloroquine, the most widely used drug for treating the *P. vivax* blood stages, is failing in several endemic regions, but recrudescence (drug failure) can be confused with relapse, confounding efficacy studies5. Accurate diagnosis of the cause of recurrent infection is also essential to clinical studies of hypnozoiticidal (anti relapse) regimens (primaquine and tafenoquine) as this relies on distinguishing relapse from reinfection6,7. Reinfection dilutes observed differences in recurrence between hypnozoiticidal interventions and biases treatment effect estimates towards interventions which provide longer post-treatment prophylaxis8.

The ability to discriminate between recurrent *P. falciparum* infections deriving from reinfections (which are likely genetically heterologous to the initial infection given sufficient population diversity) and recrudescences (which contain homologous parasites) using genotyping at a handful of polymorphic markers, represented a crucial advance for *P. falciparum* clinical research9. It allowed therapeutic efficacy studies to be conducted in endemic areas without the need for concomitant detailed epidemiological assessment. However, the situation is more complex for *P. vivax* malaria since relapsing parasites can be genetically homologous or heterologous to the initial infection10,11. However, recent genomic studies have revealed that a proportion of paired *P. vivax* isolates (from acute and recurrent infections) that would be classified as heterologous using traditional genotyping approaches share homology in large segments of the genome, inferring familial relatedness12-14. In mosquito hosts capturing blood meals with mixtures of parasite genomes, the obligate meiotic stage will generate meiotic recombinants, producing sporozoites that share parents. Every natural infection deriving from more than one sporozoite may comprise mixtures of meiotic recombinants. Pairs of infections with evidence of recent identity by descent (IBD) consistent with close relationships such as siblings or half-siblings (as much as ≥50% and ≥25% IBD respectively), are more likely to have derived from the same mosquito inoculation than from different inoculations and are, therefore, more likely to reflect relapse than reinfection events. IBD therefore has significant potential to enable more accurate classification of recurrent *P. vivax* infections. IBD-based measures would also allow finer resolution of the spatial connectivity between *P. vivax* populations than is possible using classical methods such as the fixation index or phylogenetic approaches15. With appropriate genetic data, IBD can be used to characterise major transmission networks (foci) for targeted intervention, or the risks of infection spread between communities and across borders.

While genomic data provides the greatest information to estimate IBD, *P. vivax* patient isolates often have low parasite densities which currently precludes whole genome sequencing at large scale, even when using selective whole genome amplification methods16. Restricting analyses to infections with high parasite densities is not ideal as these are atypical of *P. vivax* infection and hence may not be representative of the true diversity in patient infections. High-throughput genotyping offers a more operationally feasible approach that can be applied to low volume samples such as dried blood spots and would be more readily implemented in surveillance frameworks in malaria-endemic countries17. However, selecting parsimonious marker sets that can capture genome wide IBD effectively is challenging18. To date, genotyping methods for *P. vivax* have relied on either capillary sequencing of microsatellite markers or next generation sequencing (NGS) of Single Nucleotide Polymorphisms (SNPs)2. In a recent *P. falciparum* study, targeted NGS of short regions comprising multiple highly variable SNPs (microhaplotypes) provided a simpler, cheaper, and higher-throughput approach than microsatellite typing, with substantially higher resolution of individual clones than with single SNPs19.

Using 1,816 publicly available *P. vivax* genomes, we identified microhaplotype markers that could estimate IBD across pairs of infections in diverse parasite populations. The goal was to establish a universal *P. vivax* IBD barcode to improve the interpretation of therapeutic evaluations of hypnozoiticidal and schizonticidal antimalarial drugs, elucidate the biology and epidemiology of *P. vivax* relapses, and provide NMCPs and other agencies with actionable information on parasite transmission within and across national borders. *In silico* analyses confirmed that *P. vivax* microhaplotypes can provide much higher resolution to estimate relatedness consistent with close relationships compared to biallelic SNP barcodes that are currently used.

## Results

### Microhaplotype candidate selection

From an initial set of 1,816 isolates from the MalariaGEN Pv4 data release (**Figure 1a**), we identified 615 high quality apparently monoclonal infection samples (>50% genome positions callable, *F*ws ≥0.95)20. This sample set was used to select variants and identify potential microhaplotype candidate regions for targeted sequencing, with the aim of identifying 200 bp regions of the *P. vivax* genome with high diversity in isolates from all geographic regions. In addition to country classifications, the samples were assigned regional population classifications derived from the MalariaGEN Pv4 project and based on genetic clustering patterns20. Although most parasites were from the Asia-Pacific region, the sample set exhibited broad geographic representation, spanning 17 countries and with all regional populations represented by at least 30 isolates (range n=31-151) (**Table 1, Supplementary Figure 1**).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F1)

Figure 1. Microhaplotype discovery pipeline.
Panel a) provides an overview of the marker selection process. Criteria for selecting samples, variants, and windows for potential microhaplotype candidate windows result in a total of 5,460 windows (200 bp). The MalariaGEN Pv4 dataset was filtered to use only high-quality monoclona samples (*F*ws ≥0.95) that had at least 50% of the core genome positions callable. SNP variants from this sample subset were then identified as biallelic, have low genotype missingness (<0.1), had high global minor allele frequencies (MAF ≥0.1), and FILTER = PASS in the MalariaGEN dataset, resulting in 13,090 total SNPs. The co genome was then scanned in coding regions for all 200 bp windows in which > 1 of the identified variants were found and filtered for high diversity (global heterozygosity ≥0.5). Panel b) provides a schematic representation of microhaplotypes with 3 SNPs. Microhaplotypes leverage SNP information content in small-windowed regions of the genome to provide a high-resolution reconstruction of the parasite genome. Three high-diversity SNPs in a single microhaplotype can have as many as 8 distinct combinations of alleles which, when combined with 100 microhaplotypes across the genome, results in high discriminatory power to characterize relatedness. Panel c) provides a map of the *P. vivax* heterozygome. Chromosomal distribution of all windows identified in the global set of high-quality, independent monoclonal infections with at least 1 SNP. Each point is an identified window, with the size increasing as the number of SNPs within the window increases. Potential microhaplotype regions are well distributed across the 14 chromosomes. The microhaplotypes with the highest SNP densities tend to be located at the ends of the chromosomes. Note, microhaplotypes were selected only from the accessible regions of the genome i.e., excluding highly diverse telomeric and sub-telomeric regions where sequence reads could not be mapped accurately. Panel d) illustrates the chromosome distribution of three panels evaluated for their capacity to reconstruct parasite relatedness. Two new microhaplotype panels were created from the microhaplotype discovery pipeline, named “Random mhaps” plotted in orange, and “High-diversity mhaps” plotted in purple. These two panels were selected to have 100 microhaplotype markers using windows that were well-spaced and had between 3-10 SNPS, with even distribution across all 14 chromosomes and a minimum diversity with heterozygosity ≥0.5. The markers for the Random mhaps panel were selected randomly, while the High-diversity mhaps panel were optimised to have the highest heterozygosity possible in each region. These two high-resolution panels were compared to the currently used biallelic 42-SNP panel, named “Broad” in green. Only 38 markers of this panel are considered informative and included in this representation.

View this table:
[Table 1.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T1)

Table 1. Geographic distribution of the sample set.
The countries represented within each regional group are Ethiopia (AF), Cambodia, Thailand and Vietnam (ESEA), Malaysia (MSEA), Indonesia and Papua New Guinea (OCE), Brazil, Colombia, El Salvador, Mexico, Nicaragua and Peru (SAM), Afghanistan, India and Iran (WAS), and Thailand (WSEA). The Thai samples in ESEA and in WSEA derive from provinces east and west respectively of a previously described malaria-free corridor that runs through the centre of Thailand42.

A total of 13,090 candidate high quality biallelic SNPs with global (pooled across regions) high minor allele frequency (MAF >0.1) and low genotype failure rate (missingness < 0.1) were identified in coding regions of the genome (excluding hypervariable regions) to minimise potential changes in primer-binding sites and missingness) (**Figure 1b)**. We then conducted a windowed scan in partially overlapping windows (200 bp windows, 50 bp scanning increments) across all identified variants of interest (39,280 windows). As illustrated in **Figure 1c**, these windows were relatively uniformly distributed across the genome. We then filtered further by heterozygosity to be above the theoretical maximum for a single biallelic SNP of 0.5. Together, these filters identified a total of 5,460 microhaplotype candidate “heterozygome windows” within the nuclear genome. Of these candidates, 100 microhaplotypes with roughly even spacing across all 14 chromosomes were chosen for further panel evaluation that each contained between 3-10 candidate SNPs, optimised for highest heterozygosity (High-diversity SNP panel: mean het = 0.82, min = 0.61, max = 0.96). A second panel of 100 randomly selected microhaplotypes (Random-SNP panel: containing any number of SNPs, heterozygosity > 0.5, evenly spaced) was also included for comparative analysis to investigate the impact of heterozygosity in panel informativeness for several use cases (**Figure 1d**). Details on the coordinates, SNPs and genes covered by the High-diversity and Random-SNP microhaplotype panels are provided in Supplementary **Table 1**. The two hypothetical microhaplotype panels (High-diversity SNP panel, Random-SNP panel) as well as a previously developed biallelic SNP panel (the 42-SNP Broad barcode)21 were then benchmarked for their ability to discriminate relatedness with IBD estimation among real and simulated infection pairs. Four SNPs within the Broad barcode had to be excluded from analyses because of high genotype failure rates (i.e., less than 5 reads covering the SNP in >10% of the available global *P. vivax* sample set), or multi-allelic status (i.e., not bi-allelic). This left 38 SNPs for evaluation; herein we refer to this panel as the 38-SNP Broad barcode.

### Effective estimation of identity by descent with microhaplotype panels

A formal evaluation of the accuracy of the three panels to estimate IBD between monoclonal *P. vivax* infections was first conducted using simulated data generated on pairs of infections across a range of relatedness (*r*) values using the *paneljudge* package22. All data were simulated under the model used to estimate relatedness. In each geographic region, the 95% confidence intervals (CIs) around the estimate ![Graphic][1]</img> were calculated as a measure of the informativeness of a given marker panel (**Figure 2**). In all geographic regions and for all marker panels, the CI intervals were shortest when *r*=1 (identical infections), closely followed by *r*=0 (strangers), and highest when *r*=0.5 (siblings) indicating lowest accuracy in predicting sibling relationships. In all geographic regions, the CI intervals were shortest around the estimate ![Graphic][2]</img> at the High-Diversity SNP microhaplotype panel followed by the Random-SNP microhaplotype panel, and the 38-SNP Broad barcode. At the High-Diversity panel, all (100%) of the estimates of *r*=0.5 had 95% CIs between 0.1 and 0.8. Whilst the Random-SNP microhaplotype panel followed similar trends to the High-diversity SNP panel in each geographic region, the Broad barcode displayed particularly large CI intervals relative to the microhaplotype panels in Maritime Southeast Asia and Oceania. The same trends in marker panel performance were observed using the root mean squared error (RMSE) of the estimates of *r* ![Graphic][3]</img> compared to the data-generating *r* (**Supplementary Figure 2**). In all regions, the High-diversity SNP microhaplotype panel exhibited consistently lower RMSE values than the Random-SNP microhaplotype and the 38-SNP Broad barcode. These results indicate that panels of ∼100 uniformly spaced microhaplotypes provide more accurate estimation of IBD than the current single SNP panel used for *P. vivax* infections; this advantage may be more pronounced in certain geographic regions, possibly reflecting the geographic representation informing SNP selection.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F2)

Figure 2. 
Confidence intervals around relatedness estimates based on data simulated using various data-generating relatedness parameters, r. Data are presented on 3 marker panels: High-diversity microhaplotype panel, Random-SNP microhaplotype panel, and 38 Broad barcode biallelic SNPs. Separate plots are provided for each *r* and geographic region; AF (Africa), ESEA (East Southeast Asia), MSEA (Maritime Southeast Asia), OCE (Oceania), SAM (South America), WAS (West Asia) and WSEA (West Southeast Asia).

In addition to the simulated data, the efficacy of the three panels in estimating IBD between infections was also evaluated using ‘real’ data from sequenced pairs of primary and recurrent *P. vivax* infections. These samples came from a range of clinical trials, as well as from returning travellers23-26. A total of 14 infection pairs satisfied the criteria for selection which required monoclonal infections with high-quality genomic data at both time points. The infection pairs originated from a range of geographic locations in the Asia-Pacific and exhibited a range of durations between the primary and recurrent episode (**Table 2**). Genome-wide IBD measures using *hmmIBD* revealed that 11 infection pairs were clones (IBD ≥0.95) reflecting potential recrudescence or relapse events, whilst 1 pair was distant relatives (0.05≤ IBD <0.25) and 2 pairs were strangers (IBD <0.05). The distant relatives and strangers reflect potential reinfection or relapse events. When the data were restricted to the 100-microhaplotype markers, the resulting *hmmIBD*-based estimates correlated highly with the genome-wide measures of IBD (rho=0.790, Spearman’s rank correlation). One pair of infections (PJ0167-C and PJ0166-C) that were classified as strangers with genomic data (IBD = 0.048) were classified as distant relatives (IBD = 0.156) with the microhaplotypes, but the other 13 pairs had concordant recurrence classifications.

View this table:
[Table 2.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T2)

Table 2. Genome-wide versus microhaplotype-based identity by descent (IBD) estimates in *P. vivax* infection pairs.
List of sample pairs evaluated for relatedness as a fraction of the genome IBD using genome-wide and microhaplotype-based data. IBD fractions were determined using *hmmIBD* software, with microhaplotype variants analysed as biallelic SNPs. Clinical survey treatments; TQ (tafenoquine), PQ (primaquine), AS (artesunate), CQ (chloroquine), AM (artesunate-mefloquine). *Recurrence classification; clone (IBD≥0.95), close relative (0.25≤IBD<0.95), distant relative (0.05≤IBD<0.25), stranger (IBD<0.05). **Returning travellers presenting at the Royal Darwin Hospital, Darwin, Australia.

### Microhaplotype panels can effectively capture diversity and differentiation

Further evaluation of the highest performance panel, the High-diversity SNP microhaplotype panel, was conducted to assess its utility in capturing key population genetic features. The genetic diversity of the High-diversity SNP panel was assessed for isolates from each of the 7 geographic regions using measures of heterozygosity and effective cardinality. Heterozygosities across the microhaplotypes varied, but overall diversity was high in all geographic regions (median > 0.70, i.e., >70% chance on average that two randomly selected isolates will differ at a marker) (**Table 3, Figure 3a**). The lowest diversity was observed in the horn of Africa (median heterozygosity = 0.70) and the highest in Oceania (median heterozygosity = 0.81). Similar trends occurred in effective cardinality (roughly equivalent to the number of alleles, with adjustment for minor allele frequency), with median values ranging from 3.88 in Africa to 5.18 in Oceania (**Table 3, Figure 3b**). The highest effective cardinality observed in the data set was a microhaplotype with a score of 36.9 (roughly 36 different alleles at a single microhaplotype marker) in West Southeast Asia (**Table 4**). A list of the microhaplotypes amongst the top three most diverse markers in each region is presented in Supplementary Table 2. Microhaplotype PvP01_05_v1:1369384, located in PVP01_0532400 (cysteine-rich protective antigen), is amongst the 3 most diverse markers in all 7 regions.

View this table:
[Table 3.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T3)

Table 3. Regional patterns of population diversity using the High-diversity microhaplotype panel.
Measures were conducted on n= 615 high-quality biologically independent monoclonal samples.

View this table:
[Table 4.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T4)

Table 4. Regional patterns of within-host infection diversity.
Measures were conducted on n=922 high-quality biologically independent samples. *Pearson’s Chi-squared test with Yates’ continuity correction. **1-sided Mann-Whitney U testing the hypothesis that the *F*ws in the COI=1 infection group is not larger than in the COI>1 infection group.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F3)

Figure 3. Comparative diversity between the High-diversity microhaplotype panel and the 38-SNP Broad Barcode.
Panel a) presents heterozygosity measures and panel b) presents effective cardinality scores in n=615 high-quality biologically independent monoclonal samples by panel and region. Panel labels; 38-SNP Broad barcode (BR38) and High-diversity microhaplotype panel (MHAP-3_10). Regional labels; Africa (AF), East Southeast Asia (ESEA), Maritime Southeast Asia (MSEA), Oceania (OCE), South America (SAM), West Asia (WAS) and West Southeast Asia (WSEA).

Within-host diversity estimates using the High-diversity SNP microhaplotype panel were generated for each isolate using Complexity of Infection (COI) measures and correlated with the genome-wide estimates derived from the *F**WS* score. Using thresholds of COI = 1 and *F**WS* ≥ 0.95 to define a monoclonal infection, there were no significant differences in the proportion of polyclonal infections defined by the microhaplotype and genome-wide data in any geographic region (**Table 4**). Furthermore, as illustrated in **Figure 4**, the *F**WS* scores of the infections with COI=1 was significantly larger than those with COI>1 in all regions, demonstrating a strong alignment between the microhaplotype and genomic data (**Table 4**).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F4)

Figure 4. Genome-wide *F*WS distribution by microhaplotype-based complexity of infection (COI).
Data from n=922 high-quality biologically independent samples from Africa (AF), East Southeast Asia (ESEA), Maritime Southeast Asia (MSEA), Oceania (OCE), South America (SAM), West Asia (WAS) and West Southeast Asia (WSEA). Panel a) provides boxplots illustrating the distribution of genome wide *F*WS scores in each of the monoclonal and polyclonal infection subsets as determined by THEREALMcCOIL analysis of the SNPs in the 100 microhaplotypes using the proportional function. In all geographic regions, the median genome-wide *F*WS scores are closer to 1 (little to no within-host diversity) in the infections defined as monoclonal. Panel b) illustrates the correlation between genome-wide *F*WS and microhaplotype-based COI estimates; a trend of decreasing COI is observed with increasing *F*WS (i.e., decreasing within-host diversity).

The capacity of the High-diversity SNP microhaplotype panel to capture spatial transmission patterns including differentiation of geographically distinct populations was also assessed, using Principal Coordinate Analysis (PCoA) and IBD-based networks. The microhaplotype-based PCoA trends were consistent with spatial trends observed with genome-wide datasets20. Analysis of all the available *P. vivax* data demonstrated three major clusters with PC1 and PC2 representing South America, Africa and West Asia (group 1), East and West Southeast Asia (group 2), and Oceania and Maritime Southeast Asia (group 3) (**Figure 5a**). Clear differentiation of the regional groups was observed within each of the major clusters. To retain accuracy in MAF estimates, IBD-based analyses were restricted to within each of the 7 regional groupings. In all geographic regions, the microhaplotype-based estimates of pairwise IBD demonstrated a significant positive correlation with the genome-wide estimates (all *P*<0.05, **Supplementary Figure 3**). In Maritime Southeast Asia, a major clonal outbreak in Malaysia (defined K2) that we have described previously with genomic data27, was also distinguished clearly with the microhaplotypes (**Figure 5b**). Furthermore, more subtle sub-structure, including the resolution of previously defined K3 and K4 Malaysian sub-populations was also captured with the microhaplotypes. In other geographic regions, patterns of infection relatedness were largely consistent with genomic patterns described in previous studies (**Supplementary Figure 4**)20. The largest infection networks were observed in South America and Papua New Guinea. Small networks were observed in other regions, but most infections were unrelated at a minimal IBD threshold of 0.25 (consistent with half-siblings or closer). In East Southeast Asia, cross-border networks were observed between Thailand and Vietnam.

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F5)

Figure 5. Spatial trends in *P. vivax* connectivity using microhaplotypes.
Panels a) presents a PCoA plot presenting PC1 (46.4%) against PC2 (25.3%) in n= 615 high-quality biologically independent monoclonal isolates from Pv4.0 at the High-diversity SNP microhaplotype set. The combinations of PC1 and PC2 provides marked separation of all 7 regional groups. Panels b) and c) present whole genome sequencing (WGS) and microhaplotype-based IBD infection networks in Malaysia. The network plots were generated at a set of 224,612 SNPs (WGS) and the High-diversity SNP microhaplotype panel in single clone Malaysian infections (n=57) at a connectivity threshold of minimum IBD 0.5 (siblings or greater relatedness). Infections are colour-coded according to sub-structure definitions based on previously described ADMIXTURE analysis with genomic data27. Infections defined as “New” (purple) were not available in the previous analysis. The clustering patterns of the WGS and microhaplotype-based data are highly consistent; both data sets capture high connectivity amongst the K2 outbreak strains, a distinct K3 sub-population, and divergent K4 infections27. The new infections appear to derive from across the different sub-populations.

Further spatial evaluations of the High-diversity SNP microhaplotype panel were undertaken to assess the accuracy of the panel in detecting imported *P. vivax* cases by determining an infection’s country of origin. Using the data from 21 countries with a minimum sample size of 4 from the available *P. vivax* global dataset, we determined the accuracy of the 494 SNPs within the microhaplotype panel in predicting country of origin using a recently developed Bi-Allele Likelihood (BALK) classifier28. The predictive accuracy of the 494 microhaplotype SNPs was compared to the 38-SNP Broad barcode and three recently identified *P. vivax* geographic barcodes selected specifically for determining country of origin (GEO33, GEO50 and GEO55 comprising 33, 50 and 55 SNPs respectively)28. The median Matthew’s correlation coefficient (MCC, ranging from -1 when prediction is 0% correct to 1 when prediction is 100% correct) of the 494 microhaplotype SNPs were greater than or equal to the 38-SNP Broad barcode, GEO33 and GEO50 in 21 (100%) countries (**Table 5, Figure 6**). The GEO55 panel exhibited higher MCCs than the microhaplotype panel in 3/21 (14%) countries (Cambodia, Indonesia and Vietnam) but lower values in 4/21 (19%) countries (Afghanistan, India, Myanmar and Thailand). The median MCCs of the microhaplotype panel exceeded 0.9 in all countries except Cambodia and Vietnam (19/21, 90%); the median MCCs in these countries were also <0.9 at the 38-SNP Broad barcode and the three GEO panels.

View this table:
[Table 5.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T5)

Table 5. Summary of median MCC scores from comparative evaluations of country prediction performance between the SNP panels.
Comparisons were undertaken between the 494 SNPs in the High-diversity microhaplotype panel (MHAP-3_10), the 38-SNP Broad barcode (BR38), and the 33-, 50- and 55-SNP GEO panels (GEO33, GEO50 and GEO55 respectively). The median Mathews correlation coefficient (MCC) summary statistics are based on 500 repeats of the stratified 10-fold cross-validation using the Bi-Allele Likelihood (BALK) classifier in a set of *n*⍰=⍰799 biologically independent samples from 21 countries (each with n ≥4)28.

![Figure 6.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F6.medium.gif)

[Figure 6.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F6)

Figure 6. Comparison of country prediction performance between SNP panels.
Comparisons were undertaken between the 494 SNPs in the High-diversity microhaplotype panel (MHAP-3_10), the 38-SNP Broad barcode (BR38), and the 33-, 50- and 55-SNP GEO panels (GEO33, GEO50 and GEO55 respectively). The boxplo present the Mathews correlation coefficient (MCC) scores from 500 repeats with stratified 10-fold cross validation for each SNP set using the Bi-Allele Likelihood classifier28. Country labels are provided on the y-axis. Each bar presents the median, interquartile range and minimum and maximum MCC for the given country and model. The analyses were based on *n*⍰=799 biologically independent samples from 21 countries (each with n ≥4 samples).

## Discussion

Our study provides the first description in silico of *P. vivax* microhaplotype panels which can be used to estimate “identity by descent” relatedness between pairs of acute and recurrent infection isolates, and thus help to discriminate between different causes of vivax malaria recurrence. A systematic genome-wide selection process was used to identify panels of 100 globally diverse microhaplotypes using an expansive *P. vivax* genome dataset. The utility of these panels was assessed against a pre-existing *P. vivax* SNP panel using both simulated and whole genome data. These relatively parsimonious panels have significant potential to improve the interpretation of clinical trials and surveillance data. Further details on the integration of microhaplotype data in the proposed use cases, and areas requiring further research and development are described herein.

A key requirement of the *P. vivax* microhaplotype panels was the ability to estimate IBD accurately between paired isolates and thereby help to classify the likely origin of vivax malaria recurrence (i.e., relapse, recrudescence or reinfection) in clinical trials. Our simulations showed that the microhaplotype panels yielded consistently higher accuracy in IBD estimation (across different IBD levels and in different populations) than the 42-SNP Broad barcode that is currently the most widely used SNP barcode for *P. vivax*. This is not surprising as the microhaplotypes have substantially greater genetic information content compared with the 38 evaluable Broad SNPs. We also demonstrated greater efficacy of IBD estimation when microhaplotypes were selected to meet specified diversity criteria (High-diversity SNP panel) relative to microhaplotypes with random SNP selections (any-SNP panel). In accordance with previous predictions18, our simulation-based results demonstrated error rates below 12% in the estimation of pairwise IBD in all populations tested using the high-diversity microhaplotype panel. The error rates were highest for the prediction of infections with 50% IBD. Further work is needed to understand how this might impact clinical predictions of individual treatment response or population-level drug efficacy both in the context of descriptive data analyses and mathematical models29. In descriptive analyses, data summaries (e.g., allelic homology observations or relatedness estimates) and a set of rules (e.g., relatedness values greater than 25% are suggestive of relapse) are used to classify recurrences categorically. Rules-based classification is not the same as estimating the probability of relapse given the data under a statistical model (as in Taylor et al29, which also accounts for the added complexity associated with multiclonal infections). Nevertheless, improved data informativeness for IBD estimation likely translates into improved model-based probabilistic classification performance, at least in the case of the model of Taylor, because that model machinery includes an intermediary that evaluates the probability of the data given IBD partitions compatible with networks of sibling, clonal and stranger parasites29.

It should be acknowledged that panel evaluation based on simulated data was done in a highly idealised manner that captures panel performance in its most favourable light. Data were simulated under the same model used to estimate relatedness, such that the model was perfectly specified in relation to the simulated data. Real data are generated by an ancestral process that the model does not capture, i.e., the model is miss-specified in relation to real data. The impact of miss-specified allele frequencies, miss-classified multiclonal infections, etc. was not evaluated. In addition, it should be acknowledged that genotyping errors and failures (missing data) were not included in our simulations. However, we anticipate that genotyping failures should be infrequent with NGS-based amplicon sequencing methods, where each locus is typically covered by hundreds of reads19.

Data from pairs of *P. vivax* isolates, collected from recurrent infections in clinical trials, confirmed the ability of the high diversity microhaplotype panel to estimate pairwise IBD. High correlations were observed between microhaplotype and genome wide IBD estimates in a set of 14 paired *P. vivax* isolates from the same patient before and after treatment. Using our assigned IBD thresholds, only one of 14 pairs of infections had mismatching classifications of ‘stranger’ versus ‘distant relative’, which are both likely to reflect either reinfections or relapses. However, the available genomic data did not comprise any infection pairs with IBD estimates around 50% (i.e., siblings), which our simulations predicted to be the most difficult relationships to determine accurately. The number of genomic pairs was limited owing to the difficulty in obtaining enough DNA and thus high-quality sequence data from recurrent infections, which typically exhibit low parasite densities. This is a major incentive for using targeted genotyping approaches. Nonetheless, we demonstrated high correlations between the microhaplotype and genome wide IBD estimates in the assessments of day 0 samples, where a range of IBD relationships were observed.

In some scenarios, genetic data alone will not be informative of the likely origin of a recurrence. A single mosquito inoculation may carry clonal parasites or related parasites generated through recombination of heterologous clones. If a human host carries hypnozoites from a single inoculation, the relapsing parasites will therefore either be homologous or related (meiotic siblings) to the incident infection. In situations where the relapses are homologous, they will be genetically indistinguishable from recrudescence. If a host carries hypnozoites from one or more previous mosquito inoculations, relapsing parasites are liable to be unrelated/heterologous to the incident infection, and thus genetically indistinguishable from reinfections. In our study, a high proportion of infection pairs (79%) had homologous genomes, although the majority (86%) of pairs derived from low seasonal transmission settings in Malaysian Borneo, southern Vietnam, Cambodia and western Thailand. A therapeutic efficacy study undertaken in Peru also identified a moderately high proportion of genetically homologous recurrences (52%, 12/23)12. The proportion of relapses amongst the recurrences in our and the Peruvian study is uncertain. In a therapeutic efficacy study undertaken in Cambodia, 20 patients were relocated to a malaria-free area excluding the possibility of reinfection; the authors were able to define confidently 5 recurrences as relapses and demonstrated that 4 (80%) of these were related to the initial infection30. As more comprehensive data become available, a clearer picture of the epidemiology of recurrent *P. vivax* infections will emerge. Even in areas where the prevalence of related recurrences is low, mathematical modelling approaches that combine genetic data with time-to-event data will help to resolve the probable cause of recurrence29.

Another key requirement of the microhaplotype panel was the ability to capture spatial *P. vivax* transmission dynamics. Strategies to contain *P. vivax* effectively will be assisted by a more comprehensive understanding of the major routes of infection spread within and across borders. The SNP-based data from our microhaplotype panel displayed clear geographic trends and high accuracy in predicting the country of origin, suggesting utility in detecting and mapping imported *P. vivax* cases. Although the high-diversity microhaplotype panel was not intentionally selected for country prediction, we captured 494 high diversity bi-allelic SNPs with 100 microhaplotypes and these highly informative genetic data enabled high performance country prediction compared to recently described geographic marker panels (GEO33, GEO50 and GEO55) that were selected specifically to differentiate *P. vivax* populations by country28. The accuracy in pairwise IBD estimation using the microhaplotype data also demonstrates a unique potential for tracking infection spread at micro-epidemiological spatial resolution, to inform the dispersal of infections within and between communities. For example, the microhaplotype data from Malaysia effectively captured a previously described clonal expansion, as well as more subtle population structure reflecting different foci of infection. The spatial analyses conducted here used biallelic SNP data at the microhaplotypes; whilst the high density of SNPs (n=494) provided rich genetic information, even greater information content can be achieved using multiallelic microhaplotypes once new software to deal with these complex datasets becomes available.

Some global regions are not well represented in the Pv4 genomic dataset that were used for marker selection, such as the Indian subcontinent, Central and South America. It is therefore unclear how well the microhaplotypes described here will capture IBD in these regions. However, the informatics pipeline that was established for the microhaplotype selection can be applied readily to update the panel as needed once additional whole genome data become available from new geographical regions. The pipeline can also be used to select country-or region-specific panels where needed.

Information on within-host infection complexity is important to capture epidemiologically relevant transmission dynamics. We observed high concordance in the proportion of polyclonal infections captured by microhaplotype-based COI and genome-wide *F**WS* measures when thresholds of COI=1 and *F**WS* ≥0.95 were applied. Only a few of infections displayed differences in classification of polyclonality between the microhaplotype and genomic datasets. In interpreting these differences, it should be acknowledged that the 0.95 *F**WS* threshold is only a guideline, and that population distributions of within-host diversity generally reflect a continuum, not discrete clusters of monoclonal and polyclonal infections.

Within-host diversity at SNP barcodes also follows a continuum, but this is less marked owing to the lower genetic resolution. When monoclonal thresholds were only applied to the COI data, the *F**WS* demonstrated significantly larger values in the COI=1 (monoclonal) vs COI>1 (polyclonal) infection group in all geographic regions investigated, highlighting consistency between the two data sets. Further work is needed to determine how to phase microhaplotype profiles in highly complex infections where any number of clones may be present in varying proportions, but tools such as the *Dcifer* software provide an important step forward31. Indeed, phasing of polyclonal malaria infections is not a challenge unique to microhaplotypes.

The *P. vivax* genome has an abundance of globally diverse microhaplotype regions that can effectively capture information on infection lineages and spatial connectivity, overcoming the previous requirement to generate genomic data in samples that are often notoriously difficult to sequence. With targeted, deep sequencing platforms, these markers have great potential to inform on the complex diversity within individual infections and associated insights on transmission. The establishment of targeted microhaplotype genotyping tools for *P. vivax* will transform the assessment of clinical surveys in this species, enhance knowledge of relapse biology, and greatly improve surveillance of infections.

## Methods

### Data preparation

The initial dataset comprised 1,816 samples from 17 countries derived from the *P. vivax* community study (Pv4.0) in the vivax and Malaria Genomic Epidemiology Networks (vivaxGEN and MalariaGEN)20, as well as previously published external studies32-34. Of the 1,895 samples described in Pv4.0 data note, the 1,816 samples reflect data for which we had access prior to the open release. All genomic datasets were generated using Illumina short-read sequencing platforms. Sequence alignment, SNP discovery and variant calling, population assignments, and *F*ws (within-sample F statistic) calculations for within-host allele infection complexity were undertaken using previously described methods within the MalariaGEN framework, producing a dataset referred to as *P. vivax* Genome Variation Project release 4.0 (Pv4.0). The *F**WS* estimates the fixation of alleles within each infection relative to the diversity observed in the total population on a scale from 0 to 1 and were provided as part of the Pv4 dataset20,35. Using a threshold of *F*WS ≥0.95 as a proxy to identify a monoclonal infection, all polyclonal infections were excluded from subsequent analysis. High-quality samples were selected using a threshold of ≥50% core genome, notated by the “analysis-set” flag in the MalariaGEN Pv4 dataset. The patient metadata provided by the contributing VivaxGEN partner studies was used to identify independent isolates for marker selection, and recurrent isolates for downstream evaluation of selected marker sites. The MalariaGEN curated metadata was used to define 7 regional-level geographic assignments with the following categories: Africa (AF), South and Central America (SAM), West Southeast Asia (WSEA), East Southeast Asia (ESEA), Maritime Southeast Asia (MSEA), Oceania (OCE) and West Asia (WAS). The regional classifications were based on genomic clustering patterns in the Pv4.0 data notated as “analysis population” flag.

### Marker selection

We tailored our marker selection to the Illumina amplicon-based sequencing workflow as this methodology is widely used in the malaria community and has demonstrated utility and feasibility for malaria molecular surveillance in low- and middle-income countries17. The maximum amplicon size for Illumina amplicon-based sequencing is 200 bp, which set the criteria for the maximum length of each microhaplotype. The target number of microhaplotypes was set to 100 based on a mathematical modelling study, which demonstrated that under idealised settings (i.e., estimating relatedness from data simulated under the model used for estimation), IBD estimates with low root mean squared error (RMSE<0.1) between monoclonal malaria parasites with relatedness equal to 0.5 (relatedness with the highest RMSE) are obtainable using ∼100 polyallelic markers18. Diminishing returns in RMSE reduction were observed above 100 polyallelic markers, highlighting this target number as a pragmatic balance of IBD accuracy against the economic cost of primers18. Panels of microhaplotypes (schematic **Figure 1a**) were subsequently selected to optimise marker selection with the previously identified criteria using a three-phase approach as described in **Figure 1b** with respect to sample, variant and window selection. In the first phase, samples were subset from the MalariaGEN Pv4.0 dataset that consisted of high-quality samples as identified in Pv4 as being independent, 50% of the genome callable, and likely monoclonal with *F*ws ≥0.95. In the second phase, candidate SNPs were identified by filtering the high-quality biallelic SNPs (filter pass), with minor allele frequency (MAF) <0.1 and genotype failure rates >0.1. Finally, candidate microhaplotype windows were identified by scanning the core genome (excluding hypervariable regions) in coding regions with 50 bp overlapping, 200 bp windows to identify amplicon-sized segments containing at least one candidate SNP, producing what we define as the heterozygome map of the accessible regions of the *P. vivax* genome (**Figure 1c**). Once the windows were identified, they were further filtered to a minimum heterozygosity of 0.5 for a total of 5,460 windows. Panels of 100 microhaplotypes were then selected from the *P. vivax* heterozygome by selecting only windows that had between 3-10 SNPs, and approximately evenly distanced spacing across 14 chromosomes. The “Random mhaps” panel had 100 markers chosen irrespectively of diversity (with a minimum of 0.5 heterozygosity) and the “High-diversity” panel had each of its markers chosen with the highest possible heterozygosity in a given genomic region (Figure 1d). Measures of MAF and heterozygosity were calculated using the scikitallel package.

### Population genetic evaluation of the microhaplotype panel

To evaluate the efficacy of the selected microhaplotype panel relative to the Broad *P. vivax* barcode as a global barcode, population diversity was assessed at each panel (microhaplotype and Broad) for each of the 7 geographic regions using measures of effective cardinality and heterozygosity computed with the R-based *paneljudge* package ([https://github.com/aimeertaylor/paneljudge](https://github.com/aimeertaylor/paneljudge)). Analyses were conducted on the high-quality monoclonal samples using one of the two alleles selected at random to reconstruct microhaplotypes at heterozygote genotype positions using the “haploidify” function in Scikit-allel. Measures of within-host diversity were also generated using *THE REAL McCOIL* package on the biallelic SNPs in the High-diversity microhaplotype panel and correlated against the genome-wide *F**WS* estimates of infection complexity35,36.

The potential of the biallelic SNPs within the microhaplotype panel in capturing spatial patterns of transmission was assessed using Principal Coordinate Analysis (PCoA) to illustrate the regional level geographic clustering patterns. PCoA analysis was performed on the same high-quality monoclonal samples (n=615) by conducting PCA measures using *sklearn*.*decomposition PCA* method of scikit-learn package version 1.2.0 on a microhaplotype-based pairwise distance matrix calculated from the proportion of non-identical microphalotype alleles across all positions. Spatial clustering patterns were also assessed using IBD-based measures of differentiation implemented with the hmmIBD package as described below and plotted using the R-based *igraph* package ([https://igraph.org](https://igraph.org)).

A Bi-allele Likelihood (BALK) classifier was used to predict the country of origin of *P. vivax* isolates using SNP data (not multi-allelic microhaplotypes) at biallelic SNPs in the High-diversity SNP microhaplotype panel, Broad *P. vivax* barcode, and three recently identified *P. vivax* geographic barcodes (GEO33, GEO50 and GEO55)21,28. Details of the BALK classifier can be found in the original paper28. Country prediction performance was assessed using data from 799 of the 922 high-quality samples in the available Pv4.0 genomic dataset. The 799 sample-set corresponds with the sample set used in the original paper describing the BALK classifier28. Briefly, the sample set was derived by filtering samples to include a single representative of samples with near-identical genomes (defined as pairs with genetic distance <0.001), subjecting to iterative data quality filtering to obtain the best representative number of samples and informative SNPs without any missing genotype by removing samples with higher missingness iteratively, and then removing samples that appeared to be imported based on genome-wide data clustering patterns. A total of 21 countries had ≥4 samples. The comparative predictive performance of the SNP panels was evaluated using a stratified 10-fold cross-validation with 500 repeats, reporting the Matthews Correlation Coefficient (MCC).

### Evaluation of IBD estimation with the microhaplotype panel

The 2 selected microhaplotype panels (Random microhaplotype panel and High-diversity microhaplotype panel), as well as the 42-SNP Broad barcode panel, were evaluated with *paneljudge* package for their ability to estimate relatedness (IBD) using a range of simulations to compute the uncertainty (error rate) in estimations of relatedness (*r*) ranging from 0 (strangers) to 1 (identical) under varying allele frequency distributions in different geographic regions. Details of the equations, assumptions, and parameter options in the *paneljudge* package can be found on GitHub ([https://github.com/aimeertaylor/paneljudge](https://github.com/aimeertaylor/paneljudge)). Briefly, for each of a range of data-generating *r* values (*r* = 0, 0.25, 0.5, 0.75 and 1.0), data were simulated on a pair of haploid genotypes, generating 100 haploid genotype pairs. Estimates of *r* ![Graphic][4]</img>, switch-rate parameter ![Graphic][5]</img>, the 95% confidence intervals (CIs) around the estimates ![Graphic][6]</img> and ![Graphic][7]</img>, and root mean square errors (RMSE) were then computed. The switch-rate parameter modifies rate of switching between latent IBD and non-IBD states in the hidden Markov model by multiplying the probability of a crossover per base pair (default value 7.4e-7 as per 37) and the inter-marker distance in base pairs. The accuracy of a given marker panel in estimating *r* was measured by the CI width and RMSE, ranging from ∼0 (near perfect informativeness) to 1 (uninformative). The CI width and RMSE of ![Graphic][8]</img> compared to the data-generating *r* never reach 0 as a genome of finite length may have many realised relatednesses compatible with a given probability of IBD. This means that there will always be some variance in realised relatedness around *r*. For comparative value, *paneljudge* simulations were run on a panel of 100 single SNP markers (i.e., biallelic markers) and for a set of 38 assayable biallelic SNPs from the 42-SNP Broad barcode that has been widely used by the *P. vivax* surveillance community21,38,39. The 100-SNP biallelic panel was composed by selecting the highest diversity SNP in each microhaplotype.

In addition to simulation-based data, the efficacy of the selected microhaplotype panel in estimating IBD was assessed relative to ‘gold-standard’ genome-wide data in pairs of primary (day 0) and recurrent *P. vivax* infections with high-quality genomic data (≥50% genome covered by at least 5 reads) which were monoclonal infections (*F**WS* ≥0.95). IBD estimates between infection pairs were generated using the *hmmIBD* version 3.0 software40. As the version 3 software does not enable adjustments in the genotyping error rate for polyallelic markers (which may reduce artificially the IBD estimates between infections) the microhaplotype data were run as biallelic SNP markers rather than polyallelic markers. Sample pairs were run in geographically defined batches using allele frequency estimates based on the regional estimates from the baseline samples (high-quality, monoclonal, independent infection) using default parameters. Genome-wide and microhaplotype-based recurrence classifications were assigned using the following infection pair IBD thresholds; clone (pairwise IBD ≥0.95), close relative (0.25≤ IBD <0.95), distant relative (0.05≤ IBD <0.25), stranger (IBD <0.05). The concordance between the microhaplotype and genome wide IBD estimates was evaluated by the correlation coefficient between the datasets, and the proportion of recurrence classification mismatches.

## Supporting information

Supplemental Table 1 [[supplements/287179_file08.xls]](pending:yes)

Supplemental Table 2 [[supplements/287179_file09.xlsx]](pending:yes)

Supplemental Figure 1 [[supplements/287179_file10.tif]](pending:yes)

Supplemental Figure 2 [[supplements/287179_file11.tif]](pending:yes)

Supplemental Figure 3 [[supplements/287179_file12.tif]](pending:yes)

Supplemental Figure 4 [[supplements/287179_file13.tif]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## Ethics

As detailed in a corresponding data note for the Pv4.0 data set, all isolates were collected in the framework of studies with ethical approval from the local research ethics boards, and with written, informed consent provided by patients or guardians for individuals of age 18 or younger20.

## Data availability

Sequencing data from these samples has been made publicly available in the European Nucleotide Archive (ENA), with details provided in a data note describing the MalariaGEN Pv4.0 data set20.

## Code availability

The custom, in-house scripts used to select the microhaplotype panels are available from [https://github.com/svsiegel/vivax-mhaps](https://github.com/svsiegel/vivax-mhaps).

## Funding

The study was supported by the National Health and Medical Research Council of Australia (APP2001083 supporting SA and SVS), the Wellcome Trust (200909 and ICRG GR071614MA Senior Fellowships in Clinical Science to RNP, 206194/Z17/Z supporting JCR and SVS) the National Institutes of Health (R01AI137154 to JCR) and the Bill & Melinda Gates Foundation (INV-043618 supporting SA and RNP).The whole genome sequencing component of the study was supported by the Medical Research Council and UK Department for International Development (award number M006212 to DK) and the Wellcome Trust (award numbers 206194 and 204911 to DK). The IMPROV clinical trial was supported by the Bill & Melinda Gates Foundation (OPP1054404 awarded to RNP).

## Supplementary Material

**Supplementary Table 1. Microhaplotype panel marker selection information**.

View this table:
[Supplementary Table 2.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/T6)

Supplementary Table 2. Properties of the three most diverse markers in each geographic region for the High-diversity microhaplo type panel.
*The PvP01\_05\_v1:1369384:7,52,99,114,121,126,136,164 microhaplo type was one of the top three most diverse markers in all seven geographic regions.

![Supplementary Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F7.medium.gif)

[Supplementary Figure 1.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F7)

Supplementary Figure 1. *P. vivax* incidence map illustrating regional country groupings.
The baseline *P. vivax* incidence map was derived from the Malaria Atlas Proj (MAP) and presents the number of newly diagnosed *P. vivax* cases per 1,000 population in 202041. The labelled, dashed red lines indicate the boundaries of the geographic regions included in the identity by descent (IBD) analyses: SAM (South America), AF (Africa), WAS (West Asia), WSEA (West Southeast Asia), ESEA (East Southeast Asia), MSEA (Maritime Southeast Asia) and OCE (Oceania). Pinpoints indicate the countries included in each regional grouping.

![Supplementary Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F8.medium.gif)

[Supplementary Figure 2.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F8)

Supplementary Figure 2. Comparative accuracy in relatedness prediction at Random and High-diversity microhaplotype panels, and the 38-SNP Broad barcode.
Root mean square error (RMSE) of relatedness estimates based on data simulated using various data-generating relatedness and switch rate parameters, *r* and *k*, respectively. Data are presented on 3 marker panels: High-diversity SNP microhaplotype panel, Random-SNP microhaplotype panel and 38 Broad barcode biallelic SNPs. Panel comparisons are presented by geographic region; AF (Africa), ESEA (East Southeast Asia), MSEA (Maritime Southeast Asia), OCE (Oceania), SAM (South America), WAS (West Asia) and WSEA (West Southeast Asia).

![Supplementary Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F9.medium.gif)

[Supplementary Figure 3.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F9)

Supplementary Figure 3. Correlations between microhaplotype and genomic estimates of IBD in the regional datasets.
Panels a) to g) represent Africa (AF), East Southeast Asia (ESEA), Maritime Southeast Asia (MSEA), Oceania (OCE), South America (SAM), West Asia (WAS), and West Southeast Asia (WSEA). The microhaplotype and whole genome sequence (WGS) IBD estimates reflect pairwise estimates at the High-diversity SNP microhaplotype panel and a set of 898,448 genome-wide SNPs calculations were performed using *hmmIBD* on the 615 monoclonal sample set. Correlations were assessed with Spearman’s rho statistic (using a paired test) and presented with the associated p-value. At an alpha of 0.05, significantly positive correlations were observed in all regions.

![Supplementary Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/16/2023.03.13.23287179/F10.medium.gif)

[Supplementary Figure 4.](http://medrxiv.org/content/early/2023/03/16/2023.03.13.23287179/F10)

Supplementary Figure 4. Regional microhaplotype-based infection networks.
Panels a) to e) represent South America (SAM), Africa (AF), West Asia (WAS), West Southeast Asia (WSEA), East Southeast Asia (ESEA) and Oceania (OCE). The networks were generated from the High-diversity SNP microhaplotype panel in the 615 monoclonal sample set. The isolates from Maritime Southeast Asia (MSEA) were mostly from Malaysia (57/59) and are represented in Figure 5. Each circle reflects an infection, colour-coded by country, and line lengths reflect relatedness (shorter lines reflect greater relatedness) at a connectivity threshold of minimum identity by descent (IBD) 0.25 (half-siblings or greater relatedness). Each circle reflects an infection, with colour-coding by country.

## Acknowledgements

We thank the patients who contributed their samples to the study, and the health workers and field teams who assisted with the sample collections. Genome sequencing was undertaken by the Wellcome Sanger Institute and we thank the staff of the Wellcome Sanger Institute Sample Logistics, Sequencing, and Informatics facilities for their contribution.

*   Received March 13, 2023.
*   Revision received March 13, 2023.
*   Accepted March 16, 2023.


*   © 2023, 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.  WHO. World Malaria Report 2022. World Health Organization; Geneva 2022. (2022).
    
    

2.  Auburn, S., Cheng, Q., Marfurt, J. & Price, R. N. The changing epidemiology of Plasmodium vivax: Insights from conventional and novel surveillance tools. PLoS medicine 18, e1003560, doi:10.1371/journal.pmed.1003560 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1003560&link_type=DOI) 

3.  Commons, R. J., Simpson, J. A., Watson, J., White, N. J. & Price, R. N. Estimating the Proportion of Plasmodium vivax Recurrences Caused by Relapse: A Systematic Review and Meta-Analysis. The American journal of tropical medicine and hygiene 103, 1094–1099, doi:10.4269/ajtmh.20-0186 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.20-0186&link_type=DOI) 

4.  White, N. J. Determinants of relapse periodicity in Plasmodium vivax malaria. Malaria journal 10, 297, doi:10.1186/1475-2875-10-297 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-10-297&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21989376&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

5.  Price, R. N. et al. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. The Lancet. Infectious diseases 14, 982–991, doi:10.1016/S1473-3099(14)70855-2 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1473-3099(14)70855-2&link_type=DOI) 

6.  von Seidlein, L. et al. Review of key knowledge gaps in glucose-6-phosphate dehydrogenase deficiency detection with regard to the safe clinical deployment of 8-aminoquinoline treatment regimens: a workshop report. Malaria journal 12, 112, doi:10.1186/1475-2875-12-112 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1475-2875-12-112&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23537118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

7.  Thriemer, K. et al. Challenges for achieving safe and effective radical cure of Plasmodium vivax: a round table discussion of the APMEN Vivax Working Group. Malaria journal 16, 141, doi:10.1186/s12936-017-1784-1 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-017-1784-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28381261&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

8.  Cheah, P. Y., Steinkamp, N., von Seidlein, L. & Price, R. N. The ethics of using placebo in randomised controlled trials: a case study of a Plasmodium vivax antirelapse trial. BMC medical ethics 19, 19, doi:10.1186/s12910-018-0259-4 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12910-018-0259-4&link_type=DOI) 

9.  Snounou, G. & Beck, H. P. The use of PCR genotyping in the assessment of recrudescence or reinfection after antimalarial drug treatment. Parasitol Today 14, 462–467, doi:10.1016/s0169-4758(98)01340-4 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0169-4758(98)01340-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17040849&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076800500009&link_type=ISI) 

10. Chen, N., Auliff, A., Rieckmann, K., Gatton, M. & Cheng, Q. Relapses of Plasmodium vivax infection result from clonal hypnozoites activated at predetermined intervals. The Journal of infectious diseases 195, 934–941, doi:10.1086/512242 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/512242&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17330782&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244721700005&link_type=ISI) 

11. Imwong, M. et al. Relapses of Plasmodium vivax infection usually result from activation of heterologous hypnozoites. The Journal of infectious diseases 195, 927–933, doi:10.1086/512241 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/512241&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17330781&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000244721700004&link_type=ISI) 

12. Cowell, A. N., Valdivia, H. O., Bishop, D. K. & Winzeler, E. A. Exploration of Plasmodium vivax transmission dynamics and recurrent infections in the Peruvian Amazon using whole genome sequencing. Genome Med 10, 52, doi:10.1186/s13073-018-0563-0 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-018-0563-0&link_type=DOI) 

13. Bright, A. T. et al. A high resolution case study of a patient with recurrent Plasmodium vivax infections shows that relapses were caused by meiotic siblings. PLoS neglected tropical diseases 8, e2882, doi:10.1371/journal.pntd.0002882 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0002882&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24901334&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

14. Popovici, J. et al. Recrudescence, Reinfection, or Relapse? A More Rigorous Framework to Assess Chloroquine Efficacy for Plasmodium vivax Malaria. The Journal of infectious diseases 219, 315–322, doi:10.1093/infdis/jiy484 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiy484&link_type=DOI) 

15. Taylor, A. R. et al. Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent. PLoS genetics 13, e1007065, doi:10.1371/journal.pgen.1007065 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1007065&link_type=DOI) 

16. Noviyanti, R. et al. Implementing parasite genotyping into national surveillance frameworks: feedback from control programmes and researchers in the Asia-Pacific region. Malaria journal 19, 271, doi:10.1186/s12936-020-03330-5 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-020-03330-5&link_type=DOI) 

17. Jacob, C. G. et al. Genetic surveillance in the Greater Mekong subregion and South Asia to support malaria control and elimination. eLife 10, doi:10.7554/eLife.62997 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.62997&link_type=DOI) 

18. Taylor, A. R., Jacob, P. E., Neafsey, D. E. & Buckee, C. O. Estimating Relatedness Between Malaria Parasites. Genetics 212, 1337–1351, doi:10.1534/genetics.119.302120 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiZ2VuZXRpY3MiO3M6NToicmVzaWQiO3M6MTA6IjIxMi80LzEzMzciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMy8xNi8yMDIzLjAzLjEzLjIzMjg3MTc5LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

19. Tessema, S. K. et al. Sensitive, Highly Multiplexed Sequencing of Microhaplotypes From the Plasmodium falciparum Heterozygome. The Journal of infectious diseases 225, 1227–1237, doi:10.1093/infdis/jiaa527 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiaa527&link_type=DOI) 

20. MalariaGen et al. An open dataset of Plasmodium vivax genome variation in 1,895 worldwide samples. Wellcome open research 7, 136, doi:10.12688/wellcomeopenres.17795.1 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.12688/wellcomeopenres.17795.1&link_type=DOI) 

21. Baniecki, M. L. et al. Development of a single nucleotide polymorphism barcode to genotype Plasmodium vivax infections. PLoS neglected tropical diseases 9, e0003539, doi:10.1371/journal.pntd.0003539 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0003539&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25781890&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

22. LaVerriere, E. et al. Design and implementation of multiplexed amplicon sequencing panels to serve genomic epidemiology of infectious disease: A malaria case study. Mol Ecol Resour 22, 2285–2303, doi:10.1111/1755-0998.13622 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1755-0998.13622&link_type=DOI) 

23. Lacerda, M. V. G. et al. Single-Dose Tafenoquine to Prevent Relapse of Plasmodium vivax Malaria. The New England journal of medicine 380, 215–228, doi:10.1056/NEJMoa1710775 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1710775&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

24. Taylor, W. R. J. et al. Short-course primaquine for the radical cure of Plasmodium vivax malaria: a multicentre, randomised, placebo-controlled non-inferiority trial. Lancet 394, 929–938, doi:10.1016/S0140-6736(19)31285-1 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(19)31285-1&link_type=DOI) 

25. Chu, C. S. et al. Comparison of the Cumulative Efficacy and Safety of Chloroquine, Artesunate, and Chloroquine-Primaquine in Plasmodium vivax Malaria. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America 67, 1543–1549, doi:10.1093/cid/ciy319 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciy319&link_type=DOI) 

26. Grigg, M. J. et al. Efficacy of Artesunate-mefloquine for Chloroquine-resistant Plasmodium vivax Malaria in Malaysia: An Open-label, Randomized, Controlled Trial. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America, doi:10.1093/cid/ciw121 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciw121&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27107287&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

27. Auburn, S. et al. Genomic analysis of a pre-elimination Malaysian Plasmodium vivax population reveals selective pressures and changing transmission dynamics. Nature communications 9, 2585, doi:10.1038/s41467-018-04965-4 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-04965-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29968722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

28. Trimarsanto, H. et al. A molecular barcode and web-based data analysis tool to identify imported Plasmodium vivax malaria. Commun Biol 5, 1411, doi:10.1038/s42003-022-04352-2 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s42003-022-04352-2&link_type=DOI) 

29. Taylor, A. R. et al. Resolving the cause of recurrent Plasmodium vivax malaria probabilistically. Nature communications 10, 5595, doi:10.1038/s41467-019-13412-x (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-13412-x&link_type=DOI) 

30. Popovici, J. et al. Genomic Analyses Reveal the Common Occurrence and Complexity of Plasmodium vivax Relapses in Cambodia. mBio 9, doi:10.1128/mBio.01888-17 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoibWJpbyI7czo1OiJyZXNpZCI7czoxMzoiOS8xL2UwMTg4OC0xNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzE2LzIwMjMuMDMuMTMuMjMyODcxNzkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

31. Gerlovina, I., Gerlovin, B., Rodriguez-Barraquer, I. & Greenhouse, B. Dcifer: an IBD-based method to calculate genetic distance between polyclonal infections. Genetics 222, doi:10.1093/genetics/iyac126 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/genetics/iyac126&link_type=DOI) 

32. Parobek, C. M. et al. Selective sweep suggests transcriptional regulation may underlie Plasmodium vivax resilience to malaria control measures in Cambodia. Proceedings of the National Academy of Sciences of the United States of America 113, E8096–E8105, doi:10.1073/pnas.1608828113 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTEzLzUwL0U4MDk2IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMTYvMjAyMy4wMy4xMy4yMzI4NzE3OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

33. Lo, E. et al. Frequent expansion of Plasmodium vivax Duffy Binding Protein in Ethiopia and its epidemiological significance. PLoS neglected tropical diseases 13, e0007222, doi:10.1371/journal.pntd.0007222 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0007222&link_type=DOI) 

34. Hupalo, D. N. et al. Population genomics studies identify signatures of global dispersal and drug resistance in Plasmodium vivax. Nature genetics 48, 953–958, doi:10.1038/ng.3588 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3588&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27348298&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

35. Auburn, S. et al. Characterization of within-host Plasmodium falciparum diversity using next-generation sequence data. PloS one 7, e32891, doi:10.1371/journal.pone.0032891 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0032891&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22393456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

36. Chang, H. H. et al. THE REAL McCOIL: A method for the concurrent estimation of the complexity of infection and SNP allele frequency for malaria parasites. PLoS Comput Biol 13, e1005348, doi:10.1371/journal.pcbi.1005348 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1005348&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28125584&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

37. Miles, A. et al. Indels, structural variation, and recombination drive genomic diversity in Plasmodium falciparum. Genome research 26, 1288–1299, doi:10.1101/gr.203711.115 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ2Vub21lIjtzOjU6InJlc2lkIjtzOjk6IjI2LzkvMTI4OCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzAzLzE2LzIwMjMuMDMuMTMuMjMyODcxNzkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

38. Dewasurendra, R. L. et al. Use of a Plasmodium vivax genetic barcode for genomic surveillance and parasite tracking in Sri Lanka. Malaria journal 19, 342, doi:10.1186/s12936-020-03386-3 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-020-03386-3&link_type=DOI) 

39. Ba, H. et al. Multi-locus genotyping reveals established endemicity of a geographically distinct Plasmodium vivax population in Mauritania, West Africa. PLoS neglected tropical diseases 14, e0008945, doi:10.1371/journal.pntd.0008945 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0008945&link_type=DOI) 

40. Schaffner, S. F., Taylor, A. R., Wong, W., Wirth, D. F. & Neafsey, D. E. hmmIBD: software to infer pairwise identity by descent between haploid genotypes. Malaria journal 17, 196, doi:10.1186/s12936-018-2349-7 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12936-018-2349-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom) 

41. Battle, K. E. et al. Mapping the global endemicity and clinical burden of Plasmodium vivax, 2000-17: a spatial and temporal modelling study. Lancet 394, 332–343, doi:10.1016/S0140-6736(19)31096-7 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(19)31096-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31229233&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F16%2F2023.03.13.23287179.atom)

 [1]: /embed/inline-graphic-1.gif
 [2]: /embed/inline-graphic-2.gif
 [3]: /embed/inline-graphic-3.gif
 [4]: /embed/inline-graphic-4.gif
 [5]: /embed/inline-graphic-5.gif
 [6]: /embed/inline-graphic-6.gif
 [7]: /embed/inline-graphic-7.gif
 [8]: /embed/inline-graphic-8.gif