Microsatellite markers reveal low levels of population sub-structuring of Plasmodium falciparum in southwestern Nigeria

High genetic complexity but low relatedness in Plasmodium falciparum infections from Western Savannah Highlands and coastal equatorial Lowlands of Cameroon

To determine the diversity and connectivity of infections in Northwestern and Southwestern Cameroon, 232 Plasmodium falciparum infections, collected in 2018 from the Ndop Health District (NHD) in the western savannah highlands in the Northwest and the Limbe Health District (LHD) in the coastal lowland forests in the Southwest of Cameroon were genotyped for nine neutral microsatellite markers. Overall infection complexity and genetic diversity was significantly (p < 0.05) lower in NHD than LHD, (Mean MOI = 2.45 vs. 2.97; Fws = 0.42 vs. 0.47; Mean He = 0.84 vs. 0.89, respectively). Multi-locus linkage disequilibrium was generally low but significantly higher in the NHD than LHD population (mean ISA= 0.376 vs 0.093). Consequently, highly related pairs of isolates were observed in NHD (mean IBS = 0.086) compared to those from the LHD (mean IBS = 0.059). Infections from the two regions were mostly unrelated (mean IBS = 0.059), though the overall genetic differentiation across the geographical range was low. Indices of differentiation between the populations were however significant (overall pairwise Fst = 0.048, Jost's D = 0.133, p < 0.01). Despite the high human migration across the 270km separating the study sites, these results suggest significant restrictions to gene flow against contiguous geospatial transmission of malaria in west Cameroon. Clonal infections in the highland sites could be driven by lower levels of malaria prevalence and seasonal transmission. How these differences in genetic diversity and complexity affect responses to interventions such as drugs will require further investigations from broader community sampling.

Molecular and epidemiological characterization of imported malaria cases in Chile

Background

Chile is one of the South American countries certified as malaria-free since 1945. However, the recent increase of imported malaria cases and the presence of the vector Anopheles pseudopunctipennis in previously endemic areas in Chile require an active malaria surveillance programme.

Methods

Specimens from 268 suspected malaria cases—all imported—collected between 2015 and 2018 at the Public Health Institute of Chile (ISP), were diagnosed by microscopy and positive cases were included for epidemiological analysis. A photo-induced electron transfer fluorogenic primer real-time PCR (PET-PCR) was used to confirm the presence of malaria parasites in available blood samples. Sanger sequencing of drug resistance molecular markers (pfk13, pfcrt and pfmdr1) and microsatellite (MS) analysis were performed in confirmed Plasmodium falciparum samples and results were related to origin of infection.

Results

Out of the 268 suspected cases, 65 were Plasmodium spp. positive by microscopy. A total of 63% of the malaria patients were male and 37% were female; 43/65 of the patients acquired infections in South American endemic countries. Species confirmation of available blood samples by PET-PCR revealed that 15 samples were positive for P. falciparum, 27 for Plasmodium vivax and 4 were mixed infections. The P. falciparum samples sequenced contained four mutant pfcrt genotypes (CVMNT, CVMET, CVIET and SVMNT) and three mutant pfmdr1 genotypes (Y184F/S1034C/N1042D/D1246Y, Y184F/N1042D/D1246Y and Y184F). MS analysis confirmed that all P. falciparum samples presented different haplotypes according to the suspected country of origin. Four patients with P. vivax infection returned to the health facilities due to relapses.

Conclusion

The timely detection of polymorphisms associated with drug resistance will contribute to understanding if current drug policies in the country are appropriate for treatment of imported malaria cases and provide information about the most frequent resistant genotypes entering Chile.

Genome-wide microsatellite characteristics of five human Plasmodium species, focusing on Plasmodium malariae and P. ovale curtisi

Microsatellites can be utilized to explore genotypes, population structure, and other genomic features of eukaryotes. Systematic characterization of microsatellites has not been a focus for several species of Plasmodium, including P. malariae and P. ovale, as the majority of malaria elimination programs are focused on P. falciparum and to a lesser extent P. vivax. Here, five human malaria species (P. falciparum, P. vivax, P. malariae, P. ovale curtisi, and P. knowlesi) were investigated with the aim of conducting in-depth categorization of microsatellites for P. malariae and P. ovale curtisi. Investigation of reference genomes for microsatellites with unit motifs of 1–10 base pairs indicates high diversity among the five Plasmodium species. Plasmodium malariae, with the largest genome size, displays the second highest microsatellite density (1421 No./Mbp; 5% coverage) next to P. falciparum (3634 No./Mbp; 12% coverage). The lowest microsatellite density was observed in P. vivax (773 No./Mbp; 2% coverage). A, AT, and AAT are the most commonly repeated motifs in the Plasmodium species. For P. malariae and P. ovale curtisi, microsatellite-related sequences are observed in approximately 18–29% of coding sequences (CDS). Lysine, asparagine, and glutamic acids are most frequently coded by microsatellite-related CDS. The majority of these CDS could be related to the gene ontology terms “cell parts,” “binding,” “developmental processes,” and “metabolic processes.” The present study provides a comprehensive overview of microsatellite distribution and can assist in the planning and development of potentially useful genetic tools for further investigation of P. malariae and P. ovale curtisi epidemiology.

Diversity and Multiplicity of P. falciparum infections among asymptomatic school children in Mbita, Western Kenya

Multiplicity of infection (MOI) and genetic diversity of P. falciparum infections are important surrogate indicators for assessing malaria transmission intensity in different regions of endemicity. Determination of MOI and diversity of P. falciparum among asymptomatic carriers will enhance our understanding of parasite biology and transmission to mosquito vectors. This study examined the MOI and genetic diversity of P. falciparum parasite populations circulating in Mbita, a region characterized as one of the malaria hotspots in Kenya. The genetic diversity and multiplicity of P. falciparum infections in 95 asymptomatic school children (age 5–15 yrs.) residing in Mbita, western Kenya were assessed using 10 polymorphic microsatellite markers. An average of 79.69% (Range: 54.84–95.74%) of the isolates analysed in this study were polyclonal infections as detected in at least one locus. A high mean MOI of 3.39 (Range: 2.24–4.72) and expected heterozygosity (He) of 0.81 (Range: 0.57–0.95) was reported in the study population. The analysed samples were extensively polyclonal infections leading to circulation of highly genetically diverse parasite populations in the study area. These findings correlated with the expectations of high malaria transmission intensity despite scaling up malaria interventions in the area thereby indicating the need for a robust malaria interventions particularly against asymptomatic carriers in order to attain elimination in the region.

The Diversity, Multiplicity of Infection and Population Structure of P. falciparum Parasites Circulating in Asymptomatic Carriers Living in High and Low Malaria Transmission Settings of Ghana

Background: Diversity in Plasmodium falciparum poses a major threat to malaria control and elimination interventions. This study utilized 12 polymorphic microsatellite (MS) markers and the Msp2 marker to examine diversity, multiplicity of infection (MOI) as well as the population structure of parasites circulating in two sites separated by about 92 km and with varying malaria transmission intensities within the Greater Accra Region of Ghana. Methods: The diversity and MOI of P. falciparum parasites in 160 non-symptomatic volunteers living in Obom (high malaria transmission intensity) and Asutsuare (low malaria transmission intensity) aged between 8 and 60 years was determined using Msp2 genotyping and microsatellite analysis. Results: The prevalence of asymptomatic P. falciparum carriers as well as the parasite density of infections was significantly higher in Obom than in Asutsuare. Samples from Asutsuare and Obom were 100% and 65% clonal, respectively, by Msp2 genotyping but decreased to 50% and 5%, respectively, when determined by MS analysis. The genetic composition of parasites from Obom and Asutsuare were highly distinct, with parasites from Obom being more diverse than those from Asutsuare. Conclusion: Plasmodium falciparum parasites circulating in Obom are genetically more diverse and distinct from those circulating in Asutsuare. The MOI in samples from both Obom and Asutsuare increased when assessed by MS analysis relative to MSP2 genotyping. The TA40 and TA87 loci are useful markers for estimating MOI in high and low parasite prevalence settings.

Genetic diversity and population structure of Plasmodium falciparum in Kenyan-Ugandan border areas

Kenya has, in the last decade, made tremendous progress in the fight against malaria. Nevertheless, continued surveillance of the genetic diversity and population structure of Plasmodium falciparum is required to refine malaria control and to adapt and improve elimination strategies. Twelve neutral microsatellite loci were genotyped in 201 P. falciparum isolates obtained from the Kenyan-Ugandan border (Busia) and from two inland malaria-endemic sites situated in western (Nyando) and coastal (Msambweni) Kenya. Analyses were done to assess the genetic diversity (allelic richness and expected heterozygosity, [H_e ]), multilocus linkage disequilibrium ( I S A ) and population structure. A similarly high degree of genetic diversity was observed among the three parasite populations surveyed (mean H_e  = 0.76; P > 0.05). Except in Msambweni, random association of microsatellite loci was observed, indicating high parasite out-breeding. Low to moderate genetic structure (F_ST   = 0.022-0.076; P < 0.0001) was observed with only 5% variance in allele frequencies observed among the populations. This study shows that the genetic diversity of P. falciparum populations at the Kenyan-Ugandan border is comparable to the parasite populations from inland Kenya. In addition, high genetic diversity, panmixia and weak population structure in this study highlight the fitness of Kenyan P. falciparum populations to successfully withstand malaria control interventions.

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.

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.

Pooled-DNA sequencing identifies genomic regions of selection in Nigerian isolates of Plasmodium falciparum

Background

The burden of falciparum malaria is especially high in sub-Saharan Africa. Differences in pressure from host immunity and antimalarial drugs lead to adaptive changes responsible for high level of genetic variations within and between the parasite populations. Population-specific genetic studies to survey for genes under positive or balancing selection resulting from drug pressure or host immunity will allow for refinement of interventions.

Methods

We performed a pooled sequencing (pool-seq) of the genomes of 100 Plasmodium falciparum isolates from Nigeria. We explored allele-frequency based neutrality test (Tajima’s D) and integrated haplotype score (iHS) to identify genes under selection.

Results

Fourteen shared iHS regions that had at least 2 SNPs with a score > 2.5 were identified. These regions code for genes that were likely to have been under strong directional selection. Two of these genes were the chloroquine resistance transporter (CRT) on chromosome 7 and the multidrug resistance 1 (MDR1) on chromosome 5. There was a weak signature of selection in the dihydrofolate reductase (DHFR) gene on chromosome 4 and MDR5 genes on chromosome 13, with only 2 and 3 SNPs respectively identified within the iHS window. We observed strong selection pressure attributable to continued chloroquine and sulfadoxine-pyrimethamine use despite their official proscription for the treatment of uncomplicated malaria. There was also a major selective sweep on chromosome 6 which had 32 SNPs within the shared iHS region. Tajima’s D of circumsporozoite protein (CSP), erythrocyte-binding antigen (EBA-175), merozoite surface proteins - MSP3 and MSP7, merozoite surface protein duffy binding-like (MSPDBL2) and serine repeat antigen (SERA-5) were 1.38, 1.29, 0.73, 0.84 and 0.21, respectively.

Conclusion

We have demonstrated the use of pool-seq to understand genomic patterns of selection and variability in P. falciparum from Nigeria, which bears the highest burden of infections. This investigation identified known genomic signatures of selection from drug pressure and host immunity. This is evidence that P. falciparum populations explore common adaptive strategies that can be targeted for the development of new interventions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2260-z) contains supplementary material, which is available to authorized users.

Population genetic structure and adaptation of malaria parasites on the edge of endemic distribution

To determine whether the major human malaria parasite Plasmodium falciparum exhibits fragmented population structure or local adaptation at the northern limit of its African distribution where the dry Sahel zone meets the Sahara, samples were collected from diverse locations within Mauritania over a range of ~1000 km. Microsatellite genotypes were obtained for 203 clinical infection samples from eight locations, and Illumina paired‐end sequences were obtained to yield high coverage genomewide single nucleotide polymorphism (SNP) data for 65 clinical infection samples from four locations. Most infections contained single parasite genotypes, reflecting low rates of transmission and superinfection locally, in contrast to the situation seen in population samples from countries further south. A minority of infections shared related or identical genotypes locally, indicating some repeated transmission of parasite clones without recombination. This caused some multilocus linkage disequilibrium and local divergence, but aside from the effect of repeated genotypes there was minimal differentiation between locations. Several chromosomal regions had elevated integrated haplotype scores (|iHS|) indicating recent selection, including those containing drug resistance genes. A genomewide F_ST scan comparison with previous sequence data from an area in West Africa with higher infection endemicity indicates that regional gene flow prevents genetic isolation, but revealed allele frequency differentiation at three drug resistance loci and an erythrocyte invasion ligand gene. Contrast of extended haplotype signatures revealed none to be unique to Mauritania. Discrete foci of infection on the edge of the Sahara are genetically highly connected to the wider continental parasite population, and local elimination would be difficult to achieve without very substantial reduction in malaria throughout the region.

Genetic Diversity and Population Structure of Plasmodium falciparum in Lake Victoria Islands, A Region of Intense Transmission

Understanding the genetic structure and transmission dynamics of Plasmodium falciparum parasites in malaria-endemic regions is crucial before the implementation of interventions. Located in a high-transmission region of western Kenya where P. falciparum is the predominant species, the Lake Victoria islands are ideal for feasibility of malaria elimination studies. We analyzed genetic variation in eight microsatellite loci to examine parasite population structure and gene flow patterns across five sites. High levels of genetic diversity were measured throughout the region (mean heterozygosity index = 0.84). The overall fixation index value between the sites was 0.044, indicating that approximately 5% of the overall allelic variation is due to differences between the populations. Based on these results, we concluded that parasite population structure in the studied islands is shaped by human migration patterns that maintain extensive parasite gene flow between the sites. Consequently, any malaria elimination and interventions strategies in the study area will have to be carried out broadly on all four islands and adjoining mainland region.

Genetic diversity of Plasmodium falciparum in human malaria cases in Mali

Background

In Mali, Plasmodium falciparum malaria is highly endemic and remains stable despite the implementation of various malaria control measures. Understanding P. falciparum population structure variations across the country could provide new insights to guide malaria control programmes. In this study, P. falciparum genetic diversity and population structure in regions of varying patterns of malaria transmission in Mali were analysed.

Methods

A total of 648 blood isolates adsorbed onto filter papers during population surveillance surveys (December 2012–March 2013, October 2013) in four distinct sites of Mali were screened for the presence of P. falciparum via quantitative PCR (qPCR). Multiple loci variable number of tandem repeats analysis (MLVA) using eight microsatellite markers was then performed on positive qPCR samples. Complete genotypes were then analysed for genetic diversity, genetic differentiation and linkage disequilibrium.

Results

Of 156 qPCR-positive samples, complete genotyping of 112 samples was achieved. The parasite populations displayed high genetic diversity (mean He = 0.77), which was consistent with a high level of malaria transmission in Mali. Genetic differentiation was low (F_ST < 0.02), even between sites located approximately 900 km apart, thereby illustrating marked gene flux amongst parasite populations. The lack of linkage disequilibrium further revealed an absence of local clonal expansion, which was corroborated by the genotype relationship results. In contrast to the stable genetic diversity level observed throughout the country, mean multiplicity of infection increased from north to south (from 1.4 to 2.06) and paralleled malaria transmission levels observed locally.

Conclusions

In Mali, the high level of genetic diversity and the pronounced gene flux amongst P. falciparum populations may represent an obstacle to control malaria. Indeed, results suggest that parasite populations are polymorphic enough to adapt to their host and to counteract interventions, such as anti-malarial vaccination. Additionally, the panmictic parasite population structure imply that resistance traits may disseminate freely from one area to another, making control measures performed at a local level ineffective.

Electronic supplementary material

The online version of this article (doi:10.1186/s12936-016-1397-0) contains supplementary material, which is available to authorized users.

Genetic Diversity of Plasmodium falciparum Populations in Malaria Declining Areas of Sabah, East Malaysia

Malaysia has a national goal to eliminate malaria by 2020. Understanding the genetic diversity of malaria parasites in residual transmission foci can provide invaluable information which may inform the intervention strategies used to reach elimination targets. This study was conducted to determine the genetic diversity level of P. falciparum isolates in malaria residual foci areas of Sabah. Malaria active case detection was conducted in Kalabakan and Kota Marudu. All individuals in the study sites were screened for malaria infection by rapid diagnostic test. Blood from P. falciparum-infected individuals were collected on filter paper prior to DNA extraction. Genotyping was performed using merozoite surface protein-1 (MSP-1), merozoite surface protein-2 (MSP-2), glutamate rich protein (GLURP) and 10 neutral microsatellite loci markers. The size of alleles, multiplicity of infection (MOI), mean number of alleles (Na), expected heterozygosity (He), linkage disequilibrium (LD) and genetic differentiation (FST) were determined. In Kalabakan, the MSP-1 and MSP-2 alleles were predominantly K1 and FC27 family types, respectively. The GLURP genotype VI (751-800 bp) was predominant. The MOI for MSP-1 and MSP-2 were 1.65 and 1.20, respectively. The Na per microsatellite locus was 1.70. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.17, 0.37, 0.70 and 0.33, respectively. In Kota Marudu, the MSP-1 and MSP-2 alleles were predominantly MAD20 and 3D7 family types, respectively. The GLURP genotype IV (651-700 bp) was predominant. The MOI for both MSP-1 and MSP-2 was 1.05. The Na per microsatellite locus was 3.60. The He values for MSP-1, MSP-2, GLURP and neutral microsatellites were 0.24, 0.25, 0.69 and 0.30, respectively. A significant LD was observed in Kalabakan (0.495, p<0.01) and Kota Marudu P. falciparum populations (0.601, p<0.01). High genetic differentiation between Kalabakan and Kota Marudu P. falciparum populations was observed (FST = 0.532). The genetic data from the present study highlighted the limited diversity and contrasting genetic pattern of P. falciparum populations in the malaria declining areas of Sabah.

Analyzing Deoxyribose Nucleic Acid from Malaria Rapid Diagnostic Tests to Study Plasmodium falciparum Genetic Diversity in Mali

We evaluated the use of positive malaria rapid diagnostic tests (mRDTs) to determine genetic diversity of Plasmodium falciparum in Mali. Genetic diversity was assessed via multiple loci variable number of tandem repeats analysis (MLVA). We performed DNA extraction from 104 positive and 30 negative used mRDTs that had been stored at ambient temperature for up to 14 months. Extracted DNA was analyzed via quantitative polymerase chain reaction (qPCR), and MLVA genotyping was then assessed on positive qPCR samples. Eighty-three of the positive mRDTs (83/104, 79.8%) and none of the negative mRDTs were confirmed P. falciparum positive via qPCR. We achieved complete genotyping of 90.4% (75/83) of the qPCR-positive samples. Genotyping revealed high genetic diversity among P. falciparum populations in Mali and an absence of population clustering. We show that mRDTs are useful to monitor P. falciparum genetic diversity and thereby can provide essential data to guide malaria control programs.