Genotyping and assessment of P. falciparum multiple clonal infections
The study included a total of 261 febrile patients consisted of 143 (54.8%) male and 118 (45.2%) females. All study participants have reported fever > 37.5 C. Based on the WHO classification, all patients were classified as mild uncomplicated malaria infection. The microscopic examinations revealed that a total 196 (75.1%) patients were positive for P. falciparum infection whiles the remaining 65 (24.9%) were harboring co-infection of P. falciparum/ P. vivax. Molecular analysis using the semi-nested PCR confirmed the presence of 23.9% (47/196) P. falciparum/ P. vivax co-infections. Confirmed cases of mono-infections of P. falciparum were further amplified for the msp1 gene to exclude multiple allelic infections. Out of the 149 P. falciparum mono-infections, the msp1 revealed the presence of 83 (55.7%) multiple allelic infection which were excluded from the amplification process of the PfCSP gene. The selected 66 (44.3%) P. falciparum single allelic, mono-infections isolates were amplified for the PfCSP gene. PCR amplification produced a product sizes that vary in length from 1 kb to 1.2 kb.
Phylogenetic network of the Sudanese PfCSP and PfCSP reference sequences
In comparison with the different PfCSP reference strains, the Sudanese isolates showed high similarity with the African isolates. The phylogenetic analysis showed that most of the sequences clustered together forming a major 2 different clades consisting of the majority of the study isolates. the first clade consisted of 24 isolates of them 18 isolates were clustered together with a bootstrap value of 100% similarity and the other 6 isolates clustered together with a bootstrap value of 63. For the second major clade, 18 isolates clustered together with a bootstrap value of 100% and showing slight similarity with the Wellcome and the HB3 reference strains. The remaining isolates were distributed forming 2 different smaller clades with a high similarity with the 3D7, MAD20, RO33, Dd2, 7G8, and the NF54 reference strains (Fig. 1).
Genetic diversity in the N-terminal region of the Sudanese and global PfCSP
The analysis of the N-terminal region showed the presence of 2 different haplotypes; Hap01 included the 3D7 reference strain plus 19 isolates, and Hap02 included 47 isolates. The haplotype diversity score was 0.425 ± 0.00727. The low haplotype diversity among the Sudanese isolates indicates that the N-terminal region is highly conserved. The 2 haplotypes were observed due to the amino acid polymorphism at A98G. The amino acid insertion “61NNNGDNGREGKDEDKRDGNN80” was detected in a total of 47 isolates representing Hap02. Whereas the KLKQP motif was conserved in all the isolates (Fig. 2).
Genetic diversity in the central repeat region of the Sudanese PfCSP
Based on the NANP and NVDP repeats of the central repeat region, a total of 11 haplotypes were seen with a Hapd of 0.779 ± 0.00097. Hap01 and Hap02 were the most frequent haplotypes; 24 isolates and 18 isolates, respectively, followed by haplotypes Hap03, Hap05, and Hap06 which consisted of 9, 7, and 2 isolates, correspondingly. Whereas the remaining haplotypes consisted of a single isolate per each. The numbers of NANP repeats at each haplotype were ranging from 33 to 37 repeats whereas NVDP repeats were ranging from 3 to 6 repeats. Accordingly, the total number of repeats per haplotype was ranging from 37 to 40 repeats (Fig. 3).
Diversity patterns in the C‑terminal region of the Sudanese PfCSP
The analysis of the genetic diversity in the C-terminal region of the Sudanese PfCSP gene showed the presence of 10 different haplotypes with a haplotype diversity of 0.457 ± 0.073. Hap01 was the most frequent haplotype; 49 isolates. Hap02 was reported among 5 isolates, Hap07 was reported among 6 isolates. the remaining haplotypes were only represented by one isolate each (Fig. 4). Non-synonymous substitutions were seen in different amino acids in the C-terminal region. D294N, K298N, and N301S, were seen in Hap03, Hap07, and all the haplotypes except for Hap01 and Hap06, respectively. Several non-synonymous amino acids changes were also seen at the Th2R and the Th3R T-cell epitope regions including T317K, E317K, Q318E, K321N, I322K, T322K, R322K, K324Q, I327L, G352N, S354P, R355K, N356D, Q357E, and E361A (Fig. 4). The analysis of the C-terminal nucleotide diversity was 0.00929 ± 0.00150, and an average number of nucleotide differences of 3.51199 was observed. Neutrality testing indices to investigate the null hypothesis of natural selection using Fu and Li's D*, Fu and Li's F*, and Tajima's D tests were − 0.40167, -0.45855, and − 0.35629. All tests were statistically not significant, P > 0.10.
Haplotype network analysis and diversity patterns in the C‑terminal region of the Sudanese and global PfCSP
Based on the C-terminal region, the analysis of the haplotype network of the Sudanese isolates and the global published C-terminal region of the PfCSP the number of haplotypes detected was 69 different haplotypes with a haplotype diversity of 0.8651 ± 0.008. Hap11 was the most frequent haplotype, followed by Hap43, Hap12. The haplotype that includes the 3D7 reference strain was considered Hap01. Hap01 was including most of the studied Sudanese isolates; 48 isolates out of 66 isolates. The remaining isolates were distributed into Hap02 (5 isolates), Hap03 (1), Hap04 (1), Hap05 (1), Hap06 (1), Hap07 (6), Hap08 (1) Hap09 (1), and Hap10 (1) (Fig. 5).