The aim of this study was to perform a molecular characterisation of highly polymorphic markers in P. falciparum clinical isolates. In this first study conducted after the introduction of artemisinin-combination therapies, it was observed that P. falciparum field isolates in Cameroon exhibited a high degree of genetic polymorphism in P. falciparum msp-1 and msp-2 markers. These genes encoding for individual functional proteins expressed on the surface of the merozoite appear to play an essential role in the invasion of the red blood cell [59,60] and have captivated a lot of interest as potential vaccine candidates and as drug targets for inhibiting blood-stage replication [61-63]. Our findings revealed a high proportion of multiclonal isolates and MOI and a total of 16 different Pfmsp-1 and 27 different Pfmsp-2 gene types . This figure is not exhaustive since it is obvious that with nested PCR genotyping some of the sub-populations present in mixed infections would not be fully typed by all amplification steps [42,63]. Moreover, more amplified product was found for the msp2 gene than msp1 amplified products, this could be explained by non-synonymous substitutions introduced in template DNA that could jeopardize the proper annealing of the primer at its binding site in msp1 gene or could be explained by the fact that natural selection is more efficient when acting on msp-1 than msp-2 [59,64]suggesting that msp-1 as compared with msp2 proteins are under strong functional constraints in a complex interaction with the host leading to an increased in host's immunological response [65]. The two markers including msp-1 (16 genotypes) and msp-2 (27 genotypes) revealed considerably greater parasite diversity than glurp (5 genotypes). Thus, the observed genetic polymorphism in the these two P. falciparum major merozoite surface proteins Pfmsp1 and pfmsp2 could be explained by balancing selection occurring as a result of different mechanisms of interaction with the host [66]. This distribution of families of msp-1 and msp-2 and their allelic variations were similar to that reported from other countries with meso- to high endemicity of malaria [67, 68]. Genotyping procedure recommended in antimalarial drug trials stipulates consecutive analysis of the three markers starting with msp2 or glurp, and then msp1 [111]. Based on allelic profile of each gene obtained in this study, parasitological outcome assessment could be more accurate when both markers msp-1 and msp-2 are included in the genotyping of recurrent parasitemias in antimalarial drug trials, consistent with other studies [42]. Our findings equally raise concern over the use of glurp genotyping in antimalarial drug trials since this marker revealed limited allelic families as compare to msp-1 and msp-2.
- Genetic diversity for msp-1 and msp-2 Allelic famillies
The present study reported higher numbers of alleles (43 alleles for both MSP-1 and MSP-2) than previously reported to be circulating in the study area and in the central region of Cameroon [22,27] . A peculiarity of in this study is thatRO33 was found to be monomorphic and the most predominant allele type of msp1 compared to the polymorphic K1 and MAD20 allelic families in agreement with previous studies where RO33 was found to be the most predominant allelic family for msp1 locus [39,54,67-70]. However, our finding does not corroborate other studies in which RO33 was identified as the least predominant allelic variant type while MAD20 was the most predominant allelic family of the msp1 locus [22, 34]. This discrepancy could be attributed to the difference in the degree of transmission intensity. Another peculiarity of our study is that, the RO33 family of msp-1 did not show any polymorphism, with only 1 variant (155 bp) detected. This result differs from that of Gabon and West Uganda, where the Ro33 family was polymorphic with three and four allelic variants, respectively [15, 71], but was close to that in Senegal [72] and Brazil [73], showing the monomorphic nature of RO33 family of msp-1 . The allelic variant K1 was the second highly distributed after RO33. This does not corroborate previous finding in the same region of Cameroon [22]. but is consistent with most findings in areas of holoendemic, mesoendemic, and hyperendemic malaria, in which the allelic variant K1 was predominant[39].The predominance of Ro33 and K1 allelic family could be attributed to the balancing selection acting on these two allelic variants[74]. The MAD20 allelic variant was the least predominant among the Pfmsp1 allelic family in agreement with other studies conducted Africa including the Gambia, Nigeria, and Gabon [ 54,70]. This result is does not coroborate other studies reporting the predominance of MAD20 allelic variant over K-1 and RO33 allelic variants [23, 39, 74-76. This low distribution of MAD20 allelic variant could be partly explained by a purifying selection acting on MAD20 allelic type as compared with RO33 and K1 allelic type. This could equally be attributed to a single nucleotide substitution in DNA template that hinders proper annealing with primers designed to amplify the MAD20 allelic type [77].. An association between the distribution of k1 allelic families with severe malaria has been investigated [67,78] while theRO33 allelic family has been frequently reported in asymptomatic malaria cases[21,34]. The predominance of Ro33 seems to be beneficial for the host since the presence of this allele type is related to reduced risk of clinical malaria[67,78].A significant correlation has been established between the genetic diversity of the msp-1 gene of P. falciparum and parasite density. Likewise an association has been observed between MSP-1 allele diversity and age group on one hand, and between msp-1 allele diversity and gender among asymptomatic patients on the other hand[79].
In msp-2 the allelic famillies IC/3D7 and FC27 were almost of equal frequencies consistent with other findings [55], but in contrast with previous reports showing a significant predominance of FC27 over the 3D7/IC allelic family [34]. The frequencies of individual msp2 genotypes were low with 88.9% occurring at a frequency ≤ 10%. However 2 (16.7%) genotypes, belonging to the FC27 and 1 (6.7 %) belonging to the 3D7/IC allelic family were found at frequencies of above 10%. High genetic diversity and low allelic frequencies have been reported previously from other sites including Gabon [15], Uganda [44], Senegal [81], Burkina Faso [82]and Honduras [36].
Diversity, expressed as expected heterozygosity (He), ranged from 0.55 for msp-1 to 0.96 for msp-2 suggesting that the parasite population in Cameroon, exhibits intermediate to higher heterozygosity reflecting intermediate to high transmission pattern [30] consistent with the findings in Uganda, congo and Zimbabwe showing a high heterozygosity ranges between 0.78 and0.8 for msp-2 [34,83]. The correlation between the genetic variation of P. falciparum and malaria endemicity has been described [34-36, 71]. In areas with declining endemicity, it is reported that the number and diversity of alleles (heterozygosity) decrease with decreasing P. falciparum transmission [32,83]. The present study reveals that the msp-2 gene is highly polymorphic as compared with msp 1 gene. This observation is different in low transmission settings where high diversity has been recorded for msp-1 as compared with the msp-2 gene [23,32]. Hence, the high allelic diversity together with the low frequency of individual circulating alleles observed in the present study increase the discriminatory power of msp-1 and msp-2 to differentiate between recrudescence and re-infection. Thus, this study reinforce the importance for the genotyping of P. falciparum based on msp-1 and msp-2 in antimalarial drug efficacy trials, to distinguish between re-infection as recrudescence and emphasize on the importance of implementing msp-1 and msp-2 genotyping in effective malaria management and malaria control strategies in Cameroon and in endemic areas.
- Polyclonal infection expressed as the MOI values
In this study most patients with P. falciparum infections were infected with multiple genetically distinct parasite variants with high level of polyclonal infections observed among msp-2 (88.6%) and msp-1 (33.8%) positive isolates. Such a pattern of parasite structure is typical among P. falciparum populations in areas of high transmission, such as in sub-Saharan Africa where more than 10 variants can be routinely detected in an individual [32,62-64] and selection among these variants in the host is likely to play an important role in parasite diversity. In contrast, in areas of low transmission, such as in Asia or Latin America, patients may have infections with as few as a single variant[37,84]. No data on MOI was available in Cameroon before the introduction of ACT renders it difficult to draw a conclusive statement on the impact of ACTs introduction on MOI. Nevertheless, in this study, the polyclonal infection expressed as the MOI values were heterogeneous across the different loci, and the mean MOI was highest for msp-2 than msp-1 in accordance with previous studies in neighbouring countries with high intensity of malaria transmission including Congo and Gabon [85]. Higher MOI in Cameroon can be the results of multiple infectious mosquito bites or transmission of genetically diverse sporozoite inoculum from a single mosquito bite. Genetically distinct malaria parasites in natural populations have an extremely high rate of genetic recombination during the sexual stages resulting in multiple strains being transmitted simultaneously [86]. Effective recombination of parasites and mutation occurring in several rounds of DNA replication cycles will likely continue to maintain this genetic diversity. During genetic recombination, novel combinations of alleles can be generated offering beneficial features to the parasite, as driven by positive selection enabling spread alleles through the population. Since the evolutionary selection of malaria occurs both within individual hosts and within populations, determining the number of strain in an infection might be an important indicator of transmission intensity [87]. Effective malaria control measures (ACTs deployment, distribution of ITNs) have successfully reduced malaria transmission in many hyperendemic regions of sub-Saharan Africa. After an intervention, the malaria parasite population structure and transmission rate in these regions is expected to become similar to the low transmission rates of the regions of Southeast Asia and South America. The declining malaria transmission, as a result of scaling up interventions, has been shown to affect the genetic diversity pattern and population structure of P. falciparum [38,112]. Therefore the high MOI observed in this study is a reflect of high intensity of malaria transmission in Cameroon, despite several control strategies deployed at the health facilities and in the community such as distribution of free ACTs treatment and mosquito bed nets. This is in agreement with previous findings that observed an increase of MOI with an increase in malaria endemicity [35,36] and a low MOI for msp-1 and msp -2 correlated with a low intensity of malaria transmission [37].
Therefore regular molecular epidemiological surveys need to be performed in order to monitor the genetic diversity of P. falciparum populations in different regions of Cameroon and the world, and then correlate parasite genotypes to the disease phenotypes [79,88-90]. Previous studies reported a reduced risk of clinical malaria associated with polyclonal infections [91] and high rate of severe malaria in individual harbouring mono-infections and very common genotypes [92, 93]. Our findings revealed an increase in the mean MOI according to parasite density but not according to the gender of patients which is consistent with previous studies showing significantly high MOI in patients with moderate to high transmission [55, 94] . In the present study, most of the positive samples were from children aged 0-5 years, and this limited range of age constraint examining the correlation between the MOI and age. However, previous studies showed a pattern of greater MOI in older individuals than younger individuals reflecting more previous exposure to infection [95] while conflicting findings, indicated decreased MOI with age [96,97]. Thus, determining the MOI in endemic area including Cameroon is very important since it can be used to predict clinical outcome and target population with higher attention.
- Genotyping by gel electrophoresis and direct sequencing
In this study, the polymorphic surface antigens msp2, glurp, and msp1 genes for 304 isolates were successfully amplified by Nested PCR and genotyped using agarose gel electrophoresis of which 127 msp-1 and 297 msp-2 gene fragments with single band were selected for further characterisation by direct sequencing. In this study, the lengths of the repeat units, whose number varies between the different allelic variants, was taken into account to set bin width of 20 bp for msp-1 and 50 bp and msp2 then glurp [30]. However, variations in bin width [44,98] and the different fragment sizing methods need to be standardized to facilitate the comparison of data for a particular marker between studies. Although, genotyping method using gel electrophoresis may face some confounding factors including the variability in the electrophoretic migration of a given DNA fragment, this genotyping method could be very important where alternative method are not available and is widely used for P. falciparum genotyping [88,89] beside capillary electrophoresis [40,69,99,100]recognise to have higher resolution power than gel electrophoresis in the ability to distinguish between allelic variants of amplified fragments. However, PCR artefacts are a challenge for this method and MOI determination is often underestimated not only with gel electrophoresis but also with capillary electrophoresis analysis especially when genotyping strategy do not considered separate nPCRs for each allelic family [42]. In contrast to previous studies, a separate nPCR followed by direct sequencing was performed in the present study enabling accurate determination of the fragments sizes for msp1 and msp2 allelic families as well as sequence motifs and nucleotide differences. The development of single-nucleotide polymorphism-based genotyping techniques and next-generation sequencing, might provide highly diverse haplotype markers with sufficient resolution to detect minority population in a mixed infection[101]. Thus, Next generation sequencing technologies and genome-wide characterization could be an alternative strategy to accurately analyse polyclonal infections [103] although complete haplotype characterization of multiclonal infections remains a challenge due to PCR artefacts and sequencing errors [103,104]. Therefore, in order to accurately study the competition and selection between variants in mixed malaria infection, new tools, more sensitive to detect minority populations and quantitative for relative parasite population sizes are required [105]. Another strategy could be the combination of sequencing method and efficient computational tool for an effective characterisation of allelic variants. A good number of software packages are being developed to analyse genome-wide SNP data of field isolates for the estimation of the presence of multiple genotypes, especially minor allele in multiclonal infections [106,107]
- Sequence analysis of genetic polymorphism of msp-1 and msp-2
This study showed highly diverse nature of P. falciparum isolates of Cameroon in respect to length and sequence motifs. Sequencing and gene alignment confirmed the identity Pfmsp1 and Pfmsp2 polymorphisms. Thus, when performing the gene alignment, we found high similarity between the peptides in Pfmsp1 and pfmsp2 of Cameroon and in the other regions in Africa. However, from all peptides analyzed, the region of the alignment corresponding K1 polymorphism had the highest similarity among all the species in the Pfmsp clade included in this study ranging from 93% to 99% homology with previously described polymorphism in isolate from Kenya and Tanzania. The MAD20 peptide sequence polymorphism was the second most conserved with 83% to 100% homology between P. falciparum isolates of Cameroon and other region of Africa and other region of the world. Thus the development of a vaccine based on K1 and MAD20 allelic variant could likely be effective in providing immune protection against malaria in those regions in Africa, although it is not yet known to what extent the high allelic diversity within the K1-like and MAD20-like allelic types is of immunological significance 66]. However, previous analysis indicated more serological variation among the allelic sequences of the K1-like compared to the MAD20-like type [108,109], and more effort has been made to incorporate the repeat sequence variation of the K1-like alleles in recombinant antigens towards design of a multivalent vaccine [64,77,110]and more than 500 different msp1 block 2 allelic sequences has been described, providing a reference for molecular epidemiological studies and potentially for design of a multi-allelic vaccine[66]. Sequencing and immunologic characterization of other allelic variants such as MAD20 for Pfmsp 1, alongside with 3D7/IC for pfmsp 2 should be conducted to obtain more useful information.