Genetic polymorphism of the N-terminal region in circumsporozoite surface protein of Plasmodium falciparum field isolates from Sudan.

BACKGROUND: Malaria caused by Plasmodium falciparum parasite is still known to be one of the most significant public health problems in sub-Saharan Africa. Genetic diversity of the Sudanese P. falciparum based on the diversity in the circumsporozoite surface protein (PfCSP) has not been previously studied. Therefore, this study aimed to investigate the genetic diversity of the N-terminal region of the pfcsp gene. METHODS: A cross-sectional molecular study was conducted; 50 blood samples have been analysed from different regions in Sudan. Patients were recruited from the health facilities of Khartoum, New Halfa, Red Sea, White Nile, Al Qadarif, Gezira, River Nile, and Ad Damazin during malaria transmission seasons between June to October and December to February 2017-2018. Microscopic and nested PCR was performed for detection of P. falciparum. Merozoite surface protein-1 was performed to differentiate single and multiple clonal infections. The N-terminal of the pfcsp gene has been sequenced using PCR-Sanger dideoxy method and analysed to sequences polymorphism including the numbers of haplotypes (H), segregating sites (S), haplotypes diversity (Hd) and the average number of nucleotide differences between two sequences (Pi) were obtained using the software DnaSP v5.10. As well as neutrality testing, Tajima's D test, Fu and Li's D and F statistics. RESULTS: PCR amplification resulted in 1200 bp of the pfcsp gene. Only 21 PCR products were successfully sequenced while 29 were presenting multiple clonal P. falciparum parasite were not sequenced. The analysis of the N-terminal region of the PfCSP amino acids sequence compared to the reference strains showed five different haplotypes. H1 consisted of 3D7, NF54, HB3 and 13 isolates of the Sudanese pfcsp. H2 comprised of 7G8, Dd2, MAD20, RO33, Wellcome strain, and 5 isolates of the Sudanese pfcsp. H3, H4, and H5 were found in 3 distinct isolates. Hd was 0.594 ± 0.065, and S was 12. The most common polymorphic site was A98G; other sites were D82Y, N83H, N83M, K85L, L86F, R87L, R87F, and A98S. Fu and Li's D* test value was - 2.70818, Fu and Li's F* test value was - 2.83907, indicating a role of negative balancing selection in the pfcsp N-terminal region. Analysis with the global pfcsp N-terminal regions showed the presence of 13 haplotypes. Haplotypes frequencies were 79.4%, 17.0%, 1.6% and 1.0% for H1, H2, H3 and H4, respectively. Remaining haplotypes frequency was 0.1% for each. Hd was 0.340 ± 0.017 with a Pi of 0.00485, S was 18 sites, and Pi was 0.00030. Amino acid polymorphisms identified in the N-terminal region of global pfcsp were present at eight positions (D82Y, N83H/M, K85L/T/N, L86F, R87L/F, A98G/V/S, D99G, and G100D). CONCLUSIONS: Sudanese pfcsp N-terminal region was well-conserved with only a few polymorphic sites. Geographical distribution of genetic diversity showed high similarity to the African isolates, and this will help and contribute in the deployment of RTS,S, a PfCSP-based vaccine, in Sudan.

Diversity pattern of Duffy binding protein sequence among Duffy-negatives and Duffy-positives in Sudan.

BACKGROUND: Vivax malaria is a leading public health concern worldwide. Due to the high prevalence of Duffy-negative blood group population, Plasmodium vivax in Africa historically is less attributable and remains a neglected disease. The interaction between Duffy binding protein and its cognate receptor, Duffy antigen receptor for chemokine plays a key role in the invasion of red blood cells and serves as a novel vaccine candidate against P. vivax. However, the polymorphic nature of P. vivax Duffy binding protein (DBP), particularly N-terminal cysteine-rich region (PvDBPII), represents a major obstacle for the successful design of a DBP-based vaccine to enable global protection. In this study, the level of pvdbpII sequence variations, Duffy blood group genotypes, number of haplotypes circulating, and the natural selection at pvdbpII in Sudan isolates were analysed and the implication in terms of DBP-based vaccine design was discussed. METHODS: Forty-two P. vivax-infected blood samples were collected from patients from different areas of Sudan during 2014-2016. For Duffy blood group genotyping, the fragment that indicates GATA-1 transcription factor binding site of the FY gene (- 33T > C) was amplified by PCR and sequenced by direct sequencing. The region II flanking pvdbpII was PCR amplified and sequenced by direct sequencing. The genetic diversity and natural selection of pvdbpII were done using DnaSP ver 5.0 and MEGA ver 5.0 programs. Based on predominant, non-synonymous, single nucleotide polymorphisms (SNPs), prevalence of Sudanese haplotypes was assessed in global isolates. RESULTS: Twenty SNPs (14 non-synonymous and 6 synonymous) were identified in pvdbpII among the 42 Sudan P. vivax isolates. Sequence analysis revealed that 11 different PvDBP haplotypes exist in Sudan P. vivax isolates and the region has evolved under positive selection. Among the identified PvDBP haplotypes five PvDBP haplotypes were shared among Duffy-negative as well as Duffy-positive individuals. The high selective pressure was mainly found on the known B cell epitopes (H3) of pvdbpII. Comparison of Sudanese haplotypes, based on 10 predominant non-synonymous SNPs with 10 malaria-endemic countries, demonstrated that Sudanese haplotypes were prevalent in most endemic countries. CONCLUSION: This is the first pvdbp genetic diversity study from an African country. Sudanese isolates display high haplotype diversity and the gene is under selective pressure. Haplotype analysis indicated that Sudanese haplotypes are a representative sample of the global population. However, studies with a large number of samples are needed. These findings would be valuable for the development of PvDBP-based malaria vaccine.

Genetic diversity of Plasmodium vivax metacaspase 1 and Plasmodium vivax multi-drug resistance 1 genes of field isolates from Mauritania, Sudan and Oman.

BACKGROUND: Plasmodium vivax is the second most important human malaria parasite, widely spread across the world. This parasite is associated with important issues in the process toward malaria elimination, including potential for relapse and increased resistance to chloroquine. Plasmodium vivax multi-drug resistant (pvmdr1) is suspected to be a marker of resistance although definitive evidence is lacking. Progress has been made in knowledge of biological factors affecting parasite growth, including mechanisms of regulated cell death and the suspected role of metacaspase. Plasmodium vivax metacaspase1 (PvMCA1-cd) has been described with a catalytic domain composed of histidine (H372) and cysteine (C428) residues. The aim of this study was to test for a link between the conserved histidine and cysteine residues in PvMCA1-cd, and the polymorphism of the P. vivax multi-drug resistant gene (pvmdr1). RESULTS: Thirty P. vivax isolates were collected from Mauritania, Sudan, and Oman. Among the 28 P. vivax isolates successfully sequenced, only 4 samples showed the conserved His (372)-Cys (428) residues in PvMCA1-cd. Single nucleotide polymorphisms observed were H372T (46.4%), H372D (39.3%), and C428R (85.7%). A new polymorphic catalytic domain was observed at His (282)-Cys (305) residues. Sequences alignment analysis of pvmdr1 showed SNP in the three codons 958, 976 and 1076. A single SNP was identified at the codon M958Y (60%), 2 SNPs were found at the position 976: Y976F (13%) and Y976V (57%), and 3 SNPs were identified at the position 1076: F1076L (40%), F1076T (53%) and F1076I (3%). Only one isolate was wildtype in all three codons (MYF), 27% were single MYL mutants, and 10% were double MFL mutants. Three new haplotypes were also identified: the triple mutant YVT was most prevalent (53.3%) distributed in the three countries, while triple YFL and YVI mutants (3%), were only found in samples from Sudan and Mauritania. CONCLUSIONS: Triple or quadruple mutants for metacaspase genes and double or triple mutants for Pvmdr1 were observed in 24/28 and 19/28 samples. There was no difference in the frequency of mutations between PvMCA1-cd and Pvmdr1 (P > 0.2). Histidine and cysteine residues in PvMCA1-cd are highly polymorphic and linkage disequilibrium with SNPs of Pvmdr1 gene may be expected from these three areas with different patterns of P. vivax transmission.

Extensive genetic diversity in Plasmodium vivax from Sudan and its genetic relationships with other geographical isolates.

Plasmodium vivax, traditionally overlooked has experienced a notable increase in cases in East Africa. This study investigated the geographical origin and genetic diversity of P. vivax in Sudan using 14 microsatellite markers. A total of 113 clinical P. vivax samples were collected from two different ecogeographical zones, New Halfa and Khartoum, in Sudan. Additionally, 841 geographical samples from the database were incorporated for a global genetic analysis to discern genetic relationships among P. vivax isolates on regional and worldwide scales. On the regional scale, our findings revealed 91 unique and 8 shared haplotypes among the Sudan samples, showcasing a remarkable genetic diversity compared to other geographical isolates and supporting the hypothesis that P. vivax originated from Africa. On a global scale, distinct genetic clustering of P. vivax isolates from Africa, South America, and Asia (including Papua New Guinea and Solomon Island) was observed, with limited admixture among the three clusters. Principal component analysis emphasized the substantial contribution of African isolates to the observed global genetic variation. The Sudanese populations displayed extensive genetic diversity, marked by significant multi-locus linkage disequilibrium, suggesting an ancestral source of P. vivax variation globally and frequent recombination among the isolates. Notably, the East African P. vivax exhibited similarity with some Asian isolates, indicating potential recent introductions. Overall, our results underscore the effectiveness of utilizing microsatellite markers for implementing robust control measures, given their ability to capture extensive genetic diversity and linkage disequilibrium patterns.

Prevalence of glucose-6-phosphate dehydrogenase deficiency (G6PDd), CareStart qualitative rapid diagnostic test performance, and genetic variants in two malaria-endemic areas in Sudan.

Glucose-6-phosphate dehydrogenase deficiency (G6PDd) is the most common enzymopathy globally, and deficient individuals may experience severe hemolysis following treatment with 8-aminoquinolines. With increasing evidence of Plasmodium vivax infections throughout sub-Saharan Africa, there is a pressing need for population-level data at on the prevalence of G6PDd. Such evidence-based data will guide the expansion of primaquine and potentially tafenoquine for radical cure of P. vivax infections. This study aimed to quantify G6PDd prevalence in two geographically distinct areas in Sudan, and evaluating the performance of a qualitative CareStart rapid diagnostic test as a point-of-care test. Blood samples were analyzed from 491 unrelated healthy persons in two malaria-endemic sites in eastern and central Sudan. A pre-structured questionnaire was used which included demographic data, risk factors and treatment history. G6PD levels were measured using spectrophotometry (SPINREACT) and first-generation qualitative CareStart rapid tests. G6PD variants (202 G>A; 376 A>G) were determined by PCR/RFLP, with a subset confirmed by Sanger sequencing. The prevalence of G6PDd by spectrophotometry was 5.5% (27/491; at 30% of adjusted male median, AMM); 27.3% (134/491; at 70% of AMM); and 13.1% (64/490) by qualitative CareStart rapid diagnostic test. The first-generation CareStart rapid diagnostic test had an overall sensitivity of 81.5% (95%CI: 61.9 to 93.7) and negative predictive value of 98.8% (97.3 to 99.6). All persons genotyped across both study sites were wild type for the G6PD G202 variant. For G6PD A376G all participants in New Halfa had wild type AA (100%), while in Khartoum the AA polymorphism was found in 90.7%; AG in 2.5%; and GG in 6.8%. Phenotypic G6PD B was detected in 100% of tested participants in New Halfa while in Khartoum, the phenotypes observed were B (96.2%), A (2.8%), and AB (1%). The African A- phenotype was not detected in this study population. Overall, G6PDd prevalence in Sudan is low-to-moderate but highly heterogeneous. Point-of-care testing with the qualitative CareStart rapid diagnostic test demonstrated moderate performance with moderate sensitivity and specificity but high negative predicative value. The two sites harbored primarily the African B phenotype. A country-wide survey is recommended to understand GP6PD deficiencies more comprehensively in Sudan.

The Duffy T-33C is an insightful marker of human history and admixture.

A contrasting genotype and allele frequency pattern between Africans and non-Africans in the Duffy (T-33C) locus is reported. Its near fixation in various populations suggest is no longer under natural selection, and that current distribution is possibly a relic of distant extreme selection combined with genetic drift during the out of Africa. We put this difference into the utility to infer the ancestral state of ambiguous loci in different populations.