Investigating the drivers of the spatio-temporal patterns of genetic differences between Plasmodium falciparum malaria infections in Kilifi County, Kenya

Knowledge of how malaria infections spread locally is important both for the design of targeted interventions aiming to interrupt malaria transmission and the design of trials to assess the interventions. A previous analysis of 1602 genotyped Plasmodium falciparum parasites in Kilifi, Kenya collected over 12 years found an interaction between time and geographic distance: the mean number of single nucleotide polymorphism (SNP) differences was lower for pairs of infections which were both a shorter time interval and shorter geographic distance apart. We determine whether the empiric pattern could be reproduced by a simple model, and what mean geographic distances between parent and offspring infections and hypotheses about genotype-specific immunity or a limit on the number of infections would be consistent with the data. We developed an individual-based stochastic simulation model of households, people and infections. We parameterized the model for the total number of infections, and population and household density observed in Kilifi. The acquisition of new infections, mutation, recombination, geographic location and clearance were included. We fit the model to the observed numbers of SNP differences between pairs of parasite genotypes. The patterns observed in the empiric data could be reproduced. Although we cannot rule out genotype-specific immunity or a limit on the number of infections per individual, they are not necessary to account for the observed patterns. The mean geographic distance between parent and offspring malaria infections for the base model was 0.5 km (95% CI 0.3-1.5), for a distribution with 68% of distances shorter than the mean. Very short mean distances did not fit well, but mixtures of distributions were also consistent with the data. For a pathogen which undergoes meiosis in a setting with moderate transmission and a low coverage of infections, analytic methods are limited but an individual-based model can be used with genotyping data to estimate parameter values and investigate hypotheses about underlying processes

Evaluating the performance of malaria genetics for inferring changes in transmission intensity using transmission modelling

Substantial progress has been made globally to control malaria, however there is a growing need for innovative new tools to ensure continued progress. One approach is to harness genetic sequencing and accompanying methodological approaches as have been used in the control of other infectious diseases. However, to utilise these methodologies for malaria we first need to extend the methods to capture the complex interactions between parasites, human and vector hosts, and environment, which all impact the level of genetic diversity and relatedness of malaria parasites. We develop an individual-based transmission model to simulate malaria parasite genetics parameterised using estimated relationships between complexity of infection and age from 5 regions in Uganda and Kenya. We predict that cotransmission and superinfection contribute equally to within-host parasite genetic diversity at 11.5% PCR prevalence, above which superinfections dominate. Finally, we characterise the predictive power of six metrics of parasite genetics for detecting changes in transmission intensity, before grouping them in an ensemble statistical model. The model predicted malaria prevalence with a mean absolute error of 0.055. Different assumptions about the availability of sample metadata were considered, with the most accurate predictions of malaria prevalence made when the clinical status and age of sampled individuals is known. Parasite genetics may provide a novel surveillance tool for estimating the prevalence of malaria in areas in which prevalence surveys are not feasible. However, the findings presented here reinforce the need for patient metadata to be recorded and made available within all future attempts to use parasite genetics for surveillance

Evaluating the performance of malaria genomics for inferring changes in transmission intensity using transmission modelling

AbstractAdvances in genetic sequencing and accompanying methodological approaches have resulted in pathogen genetics being used in the control of infectious diseases. To utilise these methodologies for malaria we first need to extend the methods to capture the complex interactions between parasites, human and vector hosts, and environment. Here we develop an individual-based transmission model to simulate malaria parasite genetics parameterised using estimated relationships between complexity of infection and age from 5 regions in Uganda and Kenya. We predict that cotransmission and superinfection contribute equally to within-host parasite genetic diversity at 11.5% PCR prevalence, above which superinfections dominate. Finally, we characterise the predictive power of six metrics of parasite genetics for detecting changes in transmission intensity, before grouping them in an ensemble statistical model. The best performing model successfully predicted malaria prevalence with mean absolute error of 0.055, suggesting genetic tools could be used for monitoring the impact of malaria interventions.</jats:p

Micro-epidemiological structuring of Plasmodium falciparum parasite populations in regions with varying transmission intensities in Africa

Background: The first models of malaria transmission assumed a completely mixed and homogeneous population of parasites.  Recent models include spatial heterogeneity and variably mixed populations. However, there are few empiric estimates of parasite mixing with which to parametize such models. Methods: Here we genotype 276 single nucleotide polymorphisms (SNPs) in 5199 P. falciparum isolates from two Kenyan sites (Kilifi county and Rachuonyo South district) and one Gambian site (Kombo coastal districts) to determine the spatio-temporal extent of parasite mixing, and use Principal Component Analysis (PCA) and linear regression to examine the relationship between genetic relatedness and distance in space and time for parasite pairs. Results: Using 107, 177 and 82 SNPs that were successfully genotyped in 133, 1602, and 1034 parasite isolates from The Gambia, Kilifi and Rachuonyo South district, respectively, we show that there are no discrete geographically restricted parasite sub-populations, but instead we see a diffuse spatio-temporal structure to parasite genotypes.  Genetic relatedness of sample pairs is predicted by relatedness in space and time. Conclusions: Our findings suggest that targeted malaria control will benefit the surrounding community, but unfortunately also that emerging drug resistance will spread rapidly through the population

Evaluating the Performance of Malaria Genetics for Inferring Changes in Transmission Intensity Using Transmission Modeling.

Substantial progress has been made globally to control malaria, however there is a growing need for innovative new tools to ensure continued progress. One approach is to harness genetic sequencing and accompanying methodological approaches as have been used in the control of other infectious diseases. However, to utilize these methodologies for malaria, we first need to extend the methods to capture the complex interactions between parasites, human and vector hosts, and environment, which all impact the level of genetic diversity and relatedness of malaria parasites. We develop an individual-based transmission model to simulate malaria parasite genetics parameterized using estimated relationships between complexity of infection and age from five regions in Uganda and Kenya. We predict that cotransmission and superinfection contribute equally to within-host parasite genetic diversity at 11.5% PCR prevalence, above which superinfections dominate. Finally, we characterize the predictive power of six metrics of parasite genetics for detecting changes in transmission intensity, before grouping them in an ensemble statistical model. The model predicted malaria prevalence with a mean absolute error of 0.055. Different assumptions about the availability of sample metadata were considered, with the most accurate predictions of malaria prevalence made when the clinical status and age of sampled individuals is known. Parasite genetics may provide a novel surveillance tool for estimating the prevalence of malaria in areas in which prevalence surveys are not feasible. However, the findings presented here reinforce the need for patient metadata to be recorded and made available within all future attempts to use parasite genetics for surveillance

An open dataset of Plasmodium falciparum genome variation in 7,000 worldwide samples.

MalariaGEN is a data-sharing network that enables groups around the world to work together on the genomic epidemiology of malaria. Here we describe a new release of curated genome variation data on 7,000 Plasmodium falciparum samples from MalariaGEN partner studies in 28 malaria-endemic countries. High-quality genotype calls on 3 million single nucleotide polymorphisms (SNPs) and short indels were produced using a standardised analysis pipeline. Copy number variants associated with drug resistance and structural variants that cause failure of rapid diagnostic tests were also analysed.  Almost all samples showed genetic evidence of resistance to at least one antimalarial drug, and some samples from Southeast Asia carried markers of resistance to six commonly-used drugs. Genes expressed during the mosquito stage of the parasite life-cycle are prominent among loci that show strong geographic differentiation. By continuing to enlarge this open data resource we aim to facilitate research into the evolutionary processes affecting malaria control and to accelerate development of the surveillance toolkit required for malaria elimination

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

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

Effect of transmission intensity on hotspots and micro-epidemiology of malaria in sub-Saharan Africa

Background: Malaria transmission intensity is heterogeneous, complicating the implementation of malaria control interventions. We provide a description of the spatial micro-epidemiology of symptomatic malaria and asymptomatic parasitaemia in multiple sites. Methods: We assembled data from 19 studies conducted between 1996 and 2015 in 7 countries of sub-Saharan Africa with homestead-level geospatial data. Data from each site were used to quantify spatial autocorrelation and examine the temporal stability of hotspots. Parameters from these analyses were examined to identify trends over varying transmission intensity. Results: Significant hotspots of malaria transmission were observed in most years and sites. The risk ratios of malaria within hotspots were highest at low malaria positive fractions (MPFs) and decreased with increasing MPF (p&lt;0·001). However, statistical significance of hotspots was lowest at extremely low and extremely high MPFs, with a peak in statistical significance at an MPF of ~0·3. In four sites with longitudinal data we noted temporal instability and variable negative correlations between MPF and average age of symptomatic malaria across all sites, suggesting varying degrees of temporal stability. Conclusions: We observed geographical micro-variation in malaria transmission at sites with a variety of transmission intensities across sub-Saharan Africa. Hotspots are marked at lower transmission intensity, but it becomes difficult to show statistical significance when cases are sparse at very low transmission intensity. Given the predictability with which hotspots occur as transmission intensity falls, malaria control programmes should have a low threshold for responding to apparent clustering of cases

Variation in the effectiveness of insecticide treated nets against malaria and outdoor biting by vectors in Kilifi, Kenya

Background: Insecticide treated nets (ITNs) protect humans against bites from the Anopheles mosquito vectors that transmit malaria, thereby reducing malaria morbidity and mortality. It has been noted that ITN use leads to a switch from indoor to outdoor feeding among these vectors. It might be expected that outdoor feeding would undermine the effectiveness of ITNs that target indoors vectors, but data are limited. Methods: We linked homestead level geospatial data to clinical surveillance data at a primary healthcare facility in Kilifi County in order to map geographical heterogeneity in ITN effectiveness and observed vector feeding behaviour using landing catches and CDC light traps in six selected areas of varying ITN effectiveness. We quantified the interaction between mosquitoes and humans to evaluate whether outdoor vector biting is a potential explanation for the variation in ITN effectiveness. Results: We observed 37% and 46% visits associated with positive malaria slides among ITN users and non-ITN-users, respectively; ITN use was associated with 32% protection from malaria (crude OR = 0.68, 95% CI: 0.64, 0.73). We obtained modification of ITN effectiveness by geographical area (p=0.016), and identified 6 hotspots using the spatial scan statistic. Majority of mosquitoes were caught outdoor (60%) and were of the An. funestus group (75%). The overall propensity to feed at times when most people were asleep was high; the vast majority of the Anopheles mosquitoes were caught at times when most people are indoors asleep. Estimates for the proportion of human-mosquito contact between the first and last hour when most humans were asleep was consistently high across all locations, ranging from 0.83 to 1.00. Conclusion: Our data do not provide evidence of an epidemiological association between microgeographical variations in ITN effectiveness and variations in the microgeographical distribution of outdoor biting.</p