Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2017 Oct 27;13(10):e1007065.
doi: 10.1371/journal.pgen.1007065. eCollection 2017 Oct.

Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent

Affiliations

Quantifying connectivity between local Plasmodium falciparum malaria parasite populations using identity by descent

Aimee R Taylor et al. PLoS Genet. .

Abstract

With the rapidly increasing abundance and accessibility of genomic data, there is a growing interest in using population genetic approaches to characterize fine-scale dispersal of organisms, providing insight into biological processes across a broad range of fields including ecology, evolution and epidemiology. For sexually recombining haploid organisms such as the human malaria parasite P. falciparum, however, there have been no systematic assessments of the type of data and methods required to resolve fine scale connectivity. This analytical gap hinders the use of genomics for understanding local transmission patterns, a crucial goal for policy makers charged with eliminating this important human pathogen. Here we use data collected from four clinics with a catchment area spanning approximately 120 km of the Thai-Myanmar border to compare the ability of divergence (FST) and relatedness based on identity by descent (IBD) to resolve spatial connectivity between malaria parasites collected from proximal clinics. We found no relationship between inter-clinic distance and FST, likely due to sampling of highly related parasites within clinics, but a significant decline in IBD-based relatedness with increasing inter-clinic distance. This association was contingent upon the data set type and size. We estimated that approximately 147 single-infection whole genome sequenced parasite samples or 222 single-infection parasite samples genotyped at 93 single nucleotide polymorphisms (SNPs) were sufficient to recover a robust spatial trend estimate at this scale. In summary, surveillance efforts cannot rely on classical measures of genetic divergence to measure P. falciparum transmission on a local scale. Given adequate sampling, IBD-based relatedness provides a useful alternative, and robust trends can be obtained from parasite samples genotyped at approximately 100 SNPs.

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. SMRU clinics on the Thai-Myanmar border.
The border between Myanmar and Thailand (also the Moei river) is depicted in blue. Grey edges are proportional to inter-clinic proportions of highly related barcode parasite sample pairs (pairs with π^IBD>0.5). Latitudes and longitudes, respectively, are: 17.128107, 98.382152 (Maela); 16.83014, 98.53737 (Wang Pha); 16.5781479, 98.5846176 (Mae Kon Ken); 16.3258896, 98.670166 (Mawker Thai).
Fig 2
Fig 2. FST estimates based on 2001–2014 barcode data plotted with respect to inter-clinic distance.
Annotations refer to site comparisons using abbreviated clinic names (MLA for Maela, 212 parasite samples; WPA for Wang Pha, 457 parasite samples; MKK for Mae Kon Ken, 116 parasite samples; and MKT for Mawker Thai, 388 parasite samples). All parasite samples were single-infection. For a clinic pair, A and B say, the FST estimate was based on nA + nB parasite samples, where n denotes the number of parasite samples per clinic. Error bars represent 95% confidence intervals based on bootstrapping over SNPs.
Fig 3
Fig 3. FST estimates based on 2001–2014 WGS data plotted with respect to inter-clinic distance.
Annotations refer to site comparisons using abbreviated clinic names (MLA for Maela, 55 parasite samples; WPA for Wang Pha, 103 parasite samples; MKK for Mae Kon Ken, 4 parasite samples; and MKT for Mawker Thai, 16 parasite samples). All parasite samples were single-infection. For clinic pair, A and B say, the FST estimate was based on nA + nB parasite samples, where n denotes the number of parasite samples per clinic. Error bars represent 95% confidence intervals, based on boostrapping over SNPs.
Fig 4
Fig 4. Logit-transformed proportions of highly related 2001–2014 barcode parasite sample pairs with respect to inter-clinic distance.
Annotations refer to site comparisons using abbreviated clinic names (MLA for Maela, 212 parasite samples; WPA for Wang Pha, 457 parasite samples; MKK for Mae Kon Ken, 116 parasite samples; and MKT for Mawker Thai, 388 parasite samples). All parasite samples were single-infection. For inter-clinic pair A and B say, the proportion was based on nA × nB parasite sample pairs, where n denotes the number of parasite samples per clinic. Error bars represent 95% confidence intervals based on bootstrapping over highly related parasite sample pair labels (equal to 1 if π^IBD>0.5 and 0 otherwise), and are therefore zero where there are no π^IBD>0.5.
Fig 5
Fig 5. Logit-transformed proportions of highly related 2001–2014 WGS parasite sample pairs plotted with respect to inter-clinic distance.
Annotations refer to site comparisons using abbreviated clinic names (MLA for Maela, 55 parasite samples; WPA for Wang Pha, 103 parasite samples; MKK for Mae Kon Ken, 4 parasite samples; and MKT for Mawker Thai, 16 parasite samples). All parasite samples were single-infection. For inter-clinic pair A and B say, the proportion was based on nA × nB parasite sample pairs, where n denotes the number of parasite samples per clinic. Error bars represent 95% confidence intervals based on bootstrapping over highly related parasite sample pair labels (equal to 1 if π^IBD>0.5 and 0 otherwise), and are therefore zero where there are no π^IBD>0.5.
Fig 6
Fig 6. Proportions of significant negative spatial trend estimates with respect to data subset sample size.
Spatial trend estimates were based on regression of highly related parasite sample pair labels (equal to 1 if π^IBD>0.5  and 0 otherwise) onto distance, within temporally adjusted models (βadjusted ΔDistance (km) for barcode data, and βadjusted 2014 ΔDistance (km) for WGS data).

Similar articles

Cited by

References

    1. Gardy J, Loman NJ, Rambaut A. Real-time digital pathogen surveillance—the time is now. Genome Biol. 2015;16(1):155 doi: 10.1186/s13059-015-0726-x - DOI - PMC - PubMed
    1. Neher RA, Bedford T. nextflu: Real-time tracking of seasonal influenza virus evolution in humans. Bioinformatics. 2015;31(21):3546–8. doi: 10.1093/bioinformatics/btv381 - DOI - PMC - PubMed
    1. Bielejec F, Baele G, Vrancken B, Suchard MA, Rambaut A, Lemey P. SpreaD3: Interactive Visualization of Spatiotemporal History and Trait Evolutionary Processes. Mol Biol Evol. 2016;33(8):2167–9. doi: 10.1093/molbev/msw082 - DOI - PMC - PubMed
    1. Hall M, Woolhouse M, Rambaut A. Epidemic Reconstruction in a Phylogenetics Framework: Transmission Trees as Partitions of the Node Set. PLoS Comput Biol. 2015;11(12):1–36. - PMC - PubMed
    1. De Maio N, Wu CH, Wilson DJ. SCOTTI: Efficient Reconstruction of Transmission within Outbreaks with the Structured Coalescent. PLoS Comput Biol. 2016;12(9):1–23. - PMC - PubMed

LinkOut - more resources