Skip to main content

Genetic polymorphisms in malaria vaccine candidate Plasmodium falciparum reticulocyte-binding protein homologue-5 among populations in Lagos, Nigeria



Vaccines are the most reliable alternative to elicit sterile immunity against malaria but their development has been hindered by polymorphisms and strain-specificity in previously studied antigens. New vaccine candidates are therefore urgently needed. Highly conserved Plasmodium falciparum reticulocyte-binding protein homologue-5 (PfRH5) has been identified as a potential candidate for anti-disease vaccine development. PfRH5 is essential for erythrocyte invasion by merozoites and crucial for parasite survival. However, there is paucity of data on the extent of genetic variations on PfRH5 in field isolates of Plasmodium falciparum. This study described genetic polymorphisms at the high affinity binding polypeptides (HABPs) 36718, 36727, 36728 of PfRH5 in Nigerian isolates of P. falciparum. This study tested the hypothesis that only specific conserved B and T cell epitopes on PfRH5 HABPs are crucial for vaccine development.


One hundred and ninety-five microscopically confirmed P. falciparum samples collected in a prospective cross-sectional study of three different populations in Lagos, Nigeria. Genetic diversity and haplotype construct of Pfrh5 gene were determined using bi-directional sequencing approach. Tajima’s D and the ratio of nonsynonymous vs synonymous mutations were utilized to estimate the extent of balancing and directional selection in the pfrh5 gene.


Sequence analysis revealed three haplotypes of PfRH5 with negative Tajima’s D and dN/dS value of − 1.717 and 0.011 ± 0.020, respectively. A single nucleotide polymorphism, SNP (G → A) at position 608 was observed, which resulted in a change of the amino acid cysteine at position 203 to tyrosine. Haplotype and nucleotide diversities were 0.318 ± 0.016 and 0.0046 ± 0.0001 while inter-population genetic differentiation ranged from 0.007 to 0.037. Five polypeptide variants were identified, the most frequent being KTKYH with a frequency of 51.3%. One B-cell epitope, 151 major histocompatibility complex (MHC) class II T-cell epitopes, four intrinsically unstructured regions (IURs) and six MHC class I T-cell epitopes were observed in the study. Phylogenetic analysis of the sequences showed clustering and evidence of evolutionary relationship with 3D7, PAS-2 and FCB-2 RH5 sequences.


This study has revealed low level of genetic polymorphisms in PfRH5 antigen with B- and T-cell epitopes in intrinsically unstructured regions along the PfRH5 gene in Lagos, Nigeria. A broader investigation is however required in other parts of the country to support the possible inclusion of PfRH5 in a cross-protective multi-component vaccine.


The pathology associated with Plasmodium falciparum malaria is elicited by erythrocyte invasion during the asexual stage of parasite development, which is an important target for vaccine development [1]. Owing to the widespread of multi-drug resistant strains of P. falciparum in malaria endemic regions, vaccines may be the reliable intervention strategy to mount sterile immunity against malaria through both humoral and T-cell responses.

The newly introduced malaria vaccine, the RTS, S/ASO1 (RTS,S), recently recommended by the World Health Organization (WHO) elicits some protective effect against malaria. However, its field efficacy is limited to 39% among children with uncomplicated malaria and 29% in severe malaria cases and only protective against P. falciparum [2]. The field efficacies of vaccine against P. falciparum, the deadliest species of malaria parasite, have been significantly hampered by genetic polymorphisms in candidate antigens. These polymorphisms alter the epitope expression leading to loss of vaccine efficacy [3], hence vaccines designed with sequence compositions identical to the infecting parasites can offer protection against infections by such parasites [4, 5].

Population genetic analysis of vaccine candidate antigens offers an insight into the extent and dynamics of genetic polymorphisms and guides vaccine design [4]. This often involves investigation of location and distribution of these genetic variations. However, genetically distinct P. falciparum strains involved in the development of strain-specific immunity are most likely carried by individuals living in malaria endemic areas [4, 6, 7].

Chemotherapy, host immunity and intensity of competitions within hosts impose selection pressure on vaccine candidate antigens [8]. Epidemiological survey of genetic diversity associated with P. falciparum at study sites for vaccine trials is utilized in molecular identification of prevalent genetic variants of vaccine candidate proteins [4]. Plasmodium falciparum reticulocyte binding-like protein family (PfRH) has been in the focus for vaccine development because it is a key determinant of erythrocyte invasion [9]. Plasmodium falciparum Rh comprises of five members: PfRH1, PfRH2a, PfRH2, PfRH4 and PfRH5. PfRH5 is currently a promising vaccine candidate antigens fast approaching clinical development [9]. Most of the clinical manifestations in malaria is as a result of erythrocyte invasion, which is a complex process involving merozoite attachment to red blood cell (RBC), re-orientation, tight junction formation and parasite entry into the RBC. The receptor PfRH5 on the RBC is Basigin (CD 147).

PfRH5 is a 63-kD protein, encoded by the rh5 gene located on chromosome 4. The protein is located at the tight junction together with its associated binding proteins, AMA1 and RON2. PfRH5 lack trans-membrane domain in its sequence, its binding with AMA1 and RON2 may be important for its membrane anchorage. Antibodies against PfRH5 from different parasite strains, inhibited the growth of different strains of P. falciparum, out-performing other antigens on the vaccine development platforms [9].

Nonetheless, limited genetic polymorphisms have been reported on pfrh5 gene, however, these SNPs significantly modify the protein affecting its binding to basigin and in the active complex. Certain peptides in PfRH5 have been identified as having high red blood cell binding affinity (HABPs) [10]. They include the HABPs 36718, 36727 and 36728 in which genetic mutations affect the protein structure of the whole PfRH5 protein [10]. Currently, a total of five informative non-synonymous SNPs have been reported in PfRH5 [10] but there is a paucity of information on PfRH5 polymorphisms among P. falciparum isolates circulating in most African countries.

This cross-sectional study aimed to evaluate the status and extent of polymorphisms in PfRH5 HABPs 36718, 36727 and 36728 among P. falciparum isolates from Lagos, Nigeria. The study also investigated the presence of B- and T-cell epitopes that are immunologically important with respect to effective malaria vaccine design.


Sampling and study design

In this cross-sectional study, participants were recruited at the peak of malaria transmission from March to September 2016 at three health facilities in Lagos State, South-West Nigeria. The study areas included (a) Ijede General Hospital, Ijede, (b) Ajeromi General Hospital, Ajegunle, (c) St. Kizito’s Primary Healthcare Centre, Lekki (Fig. 1a). Malaria transmission is high in these areas and occurs all-year round, the selected facilities have over 500 outpatients per week. Informed consent (parental consent and/or assent from under-16 participants) were obtained from the patients presenting at these centres with symptoms suggestive of malaria (axillary temperature ≥ 37.5 or history of fever in the last 48 h). Participants with P. falciparum mono-infection (infection by P. falciparum alone as confirmed by microscopy and pan-specific SD Bioline malaria rapid diagnostic test kit, Wagle, India) and a minimum of 2000 P. falciparum parasites/µL estimated by thick film examination were enrolled into the study, after which they were treated with artemether-lumefantrine according to national guidelines [11]. Finger-pricked blood samples collected from the participants (Fig. 1b) were used to prepare smears according to WHO standard protocol [12].

Fig. 1
figure 1

a Geographical location of the three study sites. b Participant sampling

The slides were independently read by two microscopists for the estimation of parasite densities and species determination using the thick and thin blood films, respectively, assuming an average white blood cell count of 8000/μL. gDNA extraction, parasite genotyping and polymerase chain reaction (PCR) Amplification of pfrh5 gene were done at the Malaria Research and Immunology Laboratory of the Nigerian Institute of Medical Research, Lagos, Nigeria. DNA re-amplification, cleaning and sequencing of pfrh5 were done at the Department of Biomedical Chemistry, Graduate school of Medicine, University of Tokyo, Tokyo, Japan.

DNA extraction and parasite clonal differentiation

Extraction of whole-blood genomic DNA (gDNA) was carried out using QIAGEN QIAamp® DNA Extraction mini-kit (Qiagen, UK, CAT NO: 51306) from the individual dried blood spots. The extracted DNA samples were stored at − 20 °C until PCR. The quality as well as quantity of each DNA sample was determined using Nanodrop 1000 (Thermo Fisher Scientific, USA). To differentiate the parasites based on clonal origin, oligonucleotide primers (MSP2-OF:5′-ATGAAGGTAATTAAAACATTGTCTATTATA-3′ and MSP2-OR: 5′-CTTTGTTACC-ATCGGTACATTCTT-3′ for the primary reaction. The secondary reactions involved: FC27-OF:5′-AATACTAAGAGTGTAGGTGCARATGCTCCA-3′ and FC27-OR: 5′-TTTTATTTGGT GCATTGCCAGAACTTGAAC-3′, IC/3D7M2-F: 5′ AGAAGTATGG-CAGAAAGTAAK CCTYCTACT-3′and IC/3D7M2-R:5′-GATTGTAATTCGGGGGATTTGTTCG-3′ from published sequences in the GenBank [13], were used to amplify the block 3 of MSP-2 in nested-PCR. The conditions for the primary reactions were an initial denaturation step at 95 °C for 5 min, followed by 30 cycles consisting of 1 min of denaturation step at 94 °C, 2 min of annealing at 55 °C, and 2 min of extension at 72 °C, followed by a final extension at 72 °C for 5 min. The same cycling conditions were used for the secondary PCR (Nest 2) reaction except that the extension in the cycle was for 1 min. The DNA regions were amplified in 50 µL reaction mixtures containing: 1.25 U PfuUltra II Fusion HS DNA polymerase (Agilent Technologies, Santa Clara, USA), 1× polymerase reaction buffer, 2.5 mM MgCl2, 0.25 mM dNTPs and 1 µM of each primer. Each amplification with conserved or family-specific primer pair being done separately as previously described [14], using 2 μL of 5 ng/μL DNA sample. PCR products were electrophoresed on 1.5% agarose gels, and DNA bands visualized with the gel documentation system (Infinity 3026, Vilber Lourmat, Marne-la-Vallée, France) after ethidium bromide staining. The sizes of the fragments obtained were estimated by comparison to the 1.3 Kb (100 bp) DNA ladder (Jena Bioscience GmbH, Germany) as marker in comparison with standard purified genomic DNA from 3D7 laboratory strains similarly electrophoresed.

Amplification of pfrh5 gene

Plasmodium falciparum reticulocyte binding protein homologue 5 (pfrh5) gene was amplified from the whole-blood genomic DNA of samples with single parasite line using the primer set sequences namely: RH5-F1: 5′-ACCCATGAGGAATTGAGTC-3′ and RH5-R1: 5′-CGGTTTCATCATCTGTCTC-3′ amplifying the exon region encoding HABPs 36718, 36727 and 36728 [10]. The cycling parameters for the PCR were: 10 min initial denaturation at 94 °C followed by 35 cycles of 1 min denaturation at 94 °C, 1 min annealing at 55 °C, 2 min extension at 72 °C and a 10 min final extension at 72 °C. All PCR reactions were carried out with 2 μL of DNA in a total volume of 50 μL containing 0.2 mM dNTPs, 2 mM MgCl2, 1 μM of each primer, and 1 unit of PfuUltra II Fusion HS DNA polymerase (Agilent Technologies, Santa Clara, USA). The gel slices containing the DNA fragment of interest were excised and purified with the MagExtractor- PCR & Gel Clean up kit (TOYOBO Co. LTD., Osaka, Japan).

The purified DNA was sequenced using both forward and reverse primers after precipitation using an ABI3777 automatic DNA sequencer (Applied Biosystems). After sequencing, the gene sequences were edited using the Bioedit software and aligned amongst themselves using CLUSTAL W in MEGA 6 [15]. Blast searches were done to compare the RH5 haplotypes with the corresponding sequences available in the Genbank database and also with P. falciparum 3D7 reference RH5 sequence (GenBank accession number XM_001351508.1). The sequences were deposited in the GenBank database of NCBI (National Centre for Biotechnological Information) (GenBank: KX418644–KX418646).

Genetic analyses

The number of segregating sites (S), observed nucleotide diversity per site between any two sequences assuming that the samples are random was estimated as (∏), number of haplotype (H: uninterrupted sequence of nucleotides or amino acid sequence variants using the PfRH5 locus of the P. falciparum 3D7 strain (GenBank accession number: XM_001351508.1) as the reference sequence for all the comparisons), number of mutations, haplotype diversity (Hd), Tajima D (positive values show excess of alleles and balancing selection), the ratio of the rate of non-synonymous to synonymous mutations (dN/dS) (1 values are indication of negative selection), Linkage disequilibrium, haplotype prevalence, Estimate of recombination effect ‘R’ between adjacent sites (to assess the rate of recombination which may influence Tajima’s D) were determined using DnaSP version 5.10.1 while phylogenetic analysis was done by neighbour-joining (NJ) method with Kimura-2-parameter distance matrix in MEGA 6.

Immunological epitopes (B- and T-cell)

The ABC predserver [16] using a threshold of 0.81 and the Human Leukocyte antigens (HLA) binding prediction algorithm (Tmhc Prediction server) were used to predict the PfRH5 linear B-cell epitopes and screen the PfRH5 sequences for potential CD4 +/CD8 + epitopes with respect to related HLA. Major histocompatibility complexes (MHC-I and MHC-II) display specificity for binding with their respective epitopes. In humans, these MHC molecules are known by HLA alleles. The HLA binding peptides were predicted (threshold setting = 3.000) using the Tmhc Pred server. This prediction is based on the virtual and quantitative matrices based on 97 MHC alleles using position specific scoring matrices (PSSMs) and utilizes supervised learning method called support vector machine [16]. The relative binding ability of the different peptides to specific MHC molecule were assessed. The peptide concentrations that leads to 50% inhibition, IC50, were taken. The top binders with the cut-off value of IC50 < 100 nM were shortlisted and those with IC50 < 50 nM selected to check the potential impact of the SNPs with respect to alteration of the core sequences of lead epitopes.

The intrinsically unstructured regions (IURs) were screened using the Biocomputing package for Accurate Detection of Protein Disorder (Cspritz) server [17].


Demographic profile of patients

Three hundred and eighty-four febrile subjects comprising 185 males and 199 females were positive for P. falciparum malaria out of the 1883 examined in the three study areas. Of the 384 microscopically confirmed malaria patients, 41.93% were < 5 years, 300 (100 per site) were included in this study. Overall, median age was 10 years and prevalence of P. falciparum malaria among the febrile participants was 20.39% ordered as Ijede > Ajeromi > Lekki (Table 1). Approximately 36.25% of the uninfected healthy controls included were < 5 years (Table 1). Three hundred individuals (≥ 1 year old) were randomly selected for genetic analysis.

Table 1 Demographic profile of study population

Parasite clonal distribution

Of the 300 isolates included in this study, 195 were PCR positive for MSP2 (showing bands for one or both families). Quality assessment of the extracted DNA showed an average yield of 80 ng/µL and purity 1.86 (OD 260/280). Nested PCR analysis of the MSP-2-confirmed samples revealed that both 3D7 and FC27 families of P. falciparum MSP-2 were present in the study sites (Fig. 2).

Fig. 2
figure 2

Distribution of MSP2 alleles in isolates from the sites studied

Polymorphisms and haplotypic variations of PfRH5

Analysis of the HABPs 36718, 36727 and 36728 of PfRH5 revealed 810 bp sequence (encoding amino acid residues 23–241 of PFRH5) of which 795 sites were invariable or monomorphic and 15 (1.85%) were polymorphic (6 singleton and 9 parsimony informative) sites. A total of 5 non-synonymous SNPs were found in the populations studied described as: K62R, T81Q, P197S, C203Y and H240R with the mutation C203Y being the most prevalent with 74.87%. Of these resultant 5 amino acid changes, 7 peculiar polypeptide variants were observed, and defined based on the amino acid composition at the 5 variable residue positions with the most frequent variant (variant 1: KTKYH) having a frequency of 51.3% and the least (variant 7: ENKCR) with a frequency of 0.5% (Table 2). Three different RH5 haplotypes (H1, H2 and H3) were found among the P. falciparum populations studied based on these variants (Table 4). A single nucleotide change (G608A) was observed in the major polymorphic nucleotide of PfRH5, HABP 36727, (according to the 3D7 reference strain numbering) corresponding to amino acid change at residue 203 from Cysteine to Tyrosine (Additional file 1: Table S1). The haplotype diversity (Hd), nucleotide diversity (π) and average number of pairwise nucleotide differences within population showed no clear evidence of regional variation (Table 3). The sequences reported in this study have been deposited in the database (GenBank: KX418644–KX418646). Phylogenetic analysis showed evolutionary relationship of the haplotypes with reported PfRH5 around the world (Fig. 3) using Plasmodium reichenowi as the out-group.

Table 2 Frequency of 7 PfRH5 variants of 5 SNPs in Nigerian isolates
Table 3 Measures of DNA sequence polymorphisms at the HABP 36727 of RH5 among Nigerian P. falciparum populations (N = 195)
Fig. 3
figure 3

Molecular phylogenetic analysis of RH5 haplotypes by Maximum Likelihood method. A neighbour-joining (NJ) tree depicting the relationship between different RH5 haplotypes observed among Nigerian P. falciparum populations and with other RH5 sequences including PAS-2, FCB-2, 7G8 and 3D reference strain RH5. The evolutionary history was inferred by using the Maximum Likelihood method. The analysis involved 78 nucleotide sequences. Numbers below the lines indicate percentage bootstrap values for 1000 replications. All positions containing gaps and missing data were eliminated. Evolutionary analyses were conducted in MEGA6 [15]. The RH5 sequence of Plasmodium reichenowi was used as an out-group

Plasmodium falciparum genetic structure, evidence of selection and recombination

The LD (Fig. 4) and dN/dS values for PfRH5 showed negative selection with negative values for Tajima’s D and values < 1 for dN/dS (Table 3). Based on the three PfRH5 haplotypes (H1, H2 and H3), inter-population nucleotide differences (Kxy) varied from 10.662 (between Ijede and Ajeromi) to 11.939 (between Lekki and Ajeromi). Similarly, the average number of nucleotide (Dxy) and net nucleotide substitution per site (Da) between populations ranged from 0.0233 to 0.0261 and 0.00015 to 0.00096, respectively (Table 4).

Fig. 4
figure 4

Linkage disequilibrium (LD) plot showing non-random association between nucleotide variants at different polymorphic sites on RH5. The R2 values are plotted against the nucleotide distances with two-tailed Fisher’s exact test of significance using DnaSP. The value of LD index (range from − 1 to +1) increased with increasing distance, indicating no recombination events taking place. All polymorphic sites were considered in the analysis. a Ajeromi, b Lekki, c Ijede, d total

Table 4 Inter-population genetic differentiation of P. falciparum

Variability of Immunological (B- and T-cell) epitopes in PfRH5

To investigate the representativeness of specific predicted epitopes on the basis of their immunological significance, in comparison with the complete active PfRH5 protein, B- and T-cell epitopes and the IURs were screened for genetic mutations. This analysed the complexity of genetic diversity of predicted epitopes. One B-cell epitope, 151 MHC class II T-cell epitopes, four IURs (highly disordered) and six MHC class I T-cell epitopes were observed (Fig. 5, Table 5). Two out of the 5 genetic mutations in the pfrh5 genes were found located in the epitope region, however, there was no commitment to populating genetic mutations in epitope regions along the gene sequences mapping the pfrh5 HABPs (x2 = 6.72, P = 0.062, df = 2). A large part of the gene sequences was conserved, of which some of the epitopes were located in the disordered region of the protein (Fig. 5).

Fig. 5
figure 5

Distribution of SNPs, B- and T-cell epitopes on PfRH5 from Nigerian P. falciparum

Table 5 Predicted CD8 + T-cell, CD4 + T-cell, and B-cell epitopes in PfRH5 (IC50 < 50 nM)


In natural P. falciparum populations, one major hindrance to effective vaccine design and development is the sequence variation on the surface of vaccine candidate proteins. Therefore, a vaccine targeting a single allelic form of parasite proteins cannot elicit required protective immunity [18].

The difficulty of developing an effective malaria vaccine is largely due to polymorphisms in candidate proteins resulting strain-specificity and failed field trials. However, PfRH5 is highly susceptible to cross-strain neutralizing vaccine-induced antibodies, and out-performs all other antigens delivered by vaccine platforms [19, 20]. Moreover, only limited polymorphisms have been reported in PfRH5. Hence, there is need for population-specific investigation of sequence diversity of PfRH5 in Nigeria to inform design of effective vaccine.

This study investigated the extent of polymorphisms in P. falciparum RH5 in Lagos, Nigeria to improve understanding of the population-specific processes of acquired immunity and vaccine failure.

HABPs 36718, 36727 and 36728 of PfRH5 were genetically analysed from Nigerian isolates. Among the genetic mutations and amino acid changes found, five non-synonymous SNPs were observed giving rise to 5 polypeptide variants with frequencies ranging from 51.3 to 84.9%.

Genetic analysis of the 810 bp PfRH5 sequences revealed that Tajima’s D and dN/dS values for PfRH5 were negative suggesting a purifying selection on this gene. Consequent upon the very low polymorphisms in this gene, there is a high tendency for susceptibility of Nigerian isolates of P. falciparum to universally effective vaccines based on PfRH5 [21] reported a rationale to advance PfRH5-based vaccine to the Phase I/IIa human trials, which if successful, could be combined with vaccines targeting other stages of the P. falciparum life cycle, such as the sporozoite stage in the RTS, S vaccine. This promises a good efficacy of such vaccines going by the results from this study.

Phylogenetic analysis revealed clustering of the haplotypes and evolutionary relationship of the haplotypes with the reference 3D7 and other reported strains including the PAS-2, FCB-2 RH5 sequences [21]. It has been clearly demonstrated that evolutionary analysis may help explain changes in tropism when information on the location and nature of functional genetic variation are provided [22, 23].

The low Fst values suggest high gene flow between two populations in the study sites. This means that any allele introduced spreads fast among the P. falciparum populations with a negative impact on control measures if it is undesirable.

However, this also implies that interbreeding occurs more freely between 2 populations, as studied, which is consistent with the finding of three RH5 haplotypes. Thus, HABP 36727, with no clear evidence of regional variation indicates that most of the malaria infections were actually originating from strains circulating inside Lagos, suggesting limited chances of finding indications of PfRH5 polymorphisms. This has been observed as a potential limitation of this study.

Apart from immunity, occurrence of ancillary causes of selection which may operate on the parasite at other stages of development, is an important limitation of applying genetic technique [21]. These may include rapidly increased frequency of antigenic alleles at a locus which occasionally occurs per chance and may be indistinguishable from selection events. Alleles existing at a locus for a period of time can acquire individual sets of neutral mutations that are unique unless recombination events occur causing linkage between the functionally different alleles.

An important factor to consider in the selection of PfRH5 for involvement in the multi-component malaria vaccine design is the location of epitopes, in conserved regions or/and intrinsically unstructured regions (IURs) of the protein [22, 24]. IURs are important targets of host immune system, including antibody recognition. Disordered epitopes are usually smaller, but more efficient and sensitive to epitope sequence variations [24]. It is therefore imperative to identify population-specific polymorphic immunological hotspots within the disordered antigen–antibody interface from our study population. Although C203Y mutation was most frequent, only one of the codons (197) within the PfRH5–basigin interface was observed to be polymorphic. These findings show a high level of conservation of PfRH5 in comparison with the expectations for a protein with such critical functional role. However, the presence of polymorphisms in the predicted immunologically important peptide requires further attention to understand the impact on the effectiveness of PfRH5-based vaccine.

This study revealed epitopes that are conserved and located in disordered regions on the RH5 sequences from Lagos Nigeria. The inhibition activity of these epitopes or fragments should be assessed by functional assays to ascertain their biological and immunological relevance as candidate epitopes of a PfRH5-based vaccine.

Recognition of a peptide by T cells requires a number of steps, most notably the binding of the peptide to a host-encoded HLA molecule. It was also observed that not significant number of SNPs were found on the epitope regions of PfRH5, which could be inferred to reveal no categorical commitment to populating polymorphisms along these regions. However, larger study populations may be required to confirm this finding.


Globally efficacious RH5-based vaccine may be considerably protective in Nigeria in spite of the complex genetic diversity of Nigerian isolates of P. falciparum. Identification of variants herein reported, is suggestive of the necessity of a wider coverage of studies on the potential efficacy of PfRH5-based vaccine in the study communities in relation to the spatial distribution of immunologically significant genetic mutations in isolates from Nigeria.

Availability of data and materials

All data generated or analysed during this study are included in this published article. The nucleotide sequences of P. falciparum RH5 gene obtained from this study have been submitted in GenBank database under the Accession Numbers: KX418644–KX418646.



reticulocyte-binding protein homolog 5


genomic deoxyribonucleic acid


high affinity binding polypeptides


linkage disequilibrium


major histocompatibility complex


human leukocyte antigen


internally unstructured regions


Molecular Evolutionary Genetics Analysis


  1. Miller LH, Baruch DI, Marsh K, Doumbo OK. The pathogenic basis of malaria. Nature. 2002;415:673–9.

    Article  CAS  Google Scholar 

  2. WHO. Malaria vaccine: WHO position paper. Vaccine. 2018;36:3576–7.

    Article  Google Scholar 

  3. Srisutham S, Saralamba N, Sriprawat K, Mayxay M, Smithuis F, Nosten F, et al. Genetic diversity of three surface protein genes in Plasmodium malariae from three Asian countries. Malar J. 2018;17:24.

    Article  Google Scholar 

  4. Takala SL, Plowe CV. Genetic diversity and malaria vaccine design, testing and efficacy: preventing and overcoming ‘vaccine resistant malaria’. Parasite Immunol. 2009;31:560–73.

    Article  CAS  Google Scholar 

  5. Neafsey DE, Juraska M, Bedford T, Benkeser D, Valim C, Griggs A, et al. Genetic diversity and protective efficacy of the RTS, S/AS01 malaria vaccine. N Engl J Med. 2015;373:2025–37.

    Article  CAS  Google Scholar 

  6. Kyabayinze DJ, Karamagi C, Kiggundu M, Kamya MR, Wabwire-Mangen F, Kironde F, et al. Multiplicity of Plasmodium falciparum infection predicts antimalarial treatment outcome in Ugandan children. Afr Health Sci. 2008;8:200–5.

    PubMed  PubMed Central  Google Scholar 

  7. Sawyer SA, Parsch J, Zhang Z, Hartl DL. Prevalence of positive selection among nearly neutral amino acid replacements in Drosophila. Proc Natl Acad Sci USA. 2007;104:6504–10.

    Article  CAS  Google Scholar 

  8. Mobegi VA, Duffy CW, Amambua-Ngwa A, Loua KM, Laman E, Nwakanma DC, et al. Genome-wide analysis of selection on the malaria parasite Plasmodium falciparum in West African populations of differing infection endemicity. Mol Biol Evol. 2014;31:1490–9.

    Article  CAS  Google Scholar 

  9. Douglas AD, Williams AR, Illingworth JJ, Kamuyu G, Biswas S, Goodman AL, et al. The blood-stage malaria antigen PfRH5 is susceptible to vaccine-inducible cross-strain neutralizing antibody. Nat Commun. 2011;20:601.

    Article  CAS  Google Scholar 

  10. Arévalo-Pinzón G, Curtidor H, Muñoz M, Patarroyo MA, Bermudez A, Patarroyo ME. A single amino acid change in the Plasmodium falciparum RH5 (PfRH5) human RBC binding sequence modifies its structure and determines species-specific binding activity. Vaccine. 2012;30:637–46.

    Article  CAS  Google Scholar 

  11. Federal Ministry of Health. National malaria control programme, 2005 Annual Report. Abuja: Federal Ministry of Health; 2008.

  12. WHO. Assessment and monitoring of antimalarial drug efficacy for the treatment of uncomplicated falciparum malaria. World Health Organization: Geneva, WHO/HTM/RBM/2003.5, 2003.

  13. Snounou G, Viriyakosol S, Jarra W, Thaithong S, Brown KN. Identification of the four human malaria parasite species in field samples by the polymerase chain reaction and detection of a high prevalence of mixed infections. Mol Biochem Parasitol. 1993;58:283–92.

    Article  CAS  Google Scholar 

  14. Mayengue PI, Ndounga M, Malonga FV, Bitemo M, Ntoumi F. Genetic polymorphism of merozoite surface protein-1 and merozoite surface protein-2 in Plasmodium falciparum isolates from Brazzaville, Republic of Congo. Malar J. 2011;10:276.

    Article  CAS  Google Scholar 

  15. Tamura K, Battistuzzi FU, Billing P, Murillo O, Filipski A, Kumar S. Estimating divergence times in large molecular phylogenies. Proc Natl Acad Sci USA. 2012;109:19333–8.

    Article  Google Scholar 

  16. Saha S, Raghava GPS. Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins. 2006;65:40–8.

    Article  CAS  Google Scholar 

  17. Biocomputing package for Accurate Detection of Protein Disorder (Cspritz) ( Accessed 23 Dec 2018.

  18. Wright GJ, Rayner JC. Plasmodium falciparum erythrocyte invasion: combining function with immune evasion. PLoS Pathog. 2014;10:e1003943.

    Article  CAS  Google Scholar 

  19. Patel SD, Ahouidi AD, Bei AK, Dieye TN, Mboup S, Harrison SC, et al. Plasmodium falciparum merozoite surface antigen, PfRH5, elicits detectable levels of invasion-inhibiting antibodies in humans. J Infect Dis. 2013;208:1679–87.

    Article  CAS  Google Scholar 

  20. Ouattara A, Tran TM, Doumbo S, Matthew A, Agrawal S, Niangaly A, et al. Extent and dynamics of polymorphism in the malaria vaccine candidate Plasmodium falciparum reticulocyte-binding protein homologue-5 in Kalifabougou, Mali. Am J Trop Med Hyg. 2018;99:43–50.

    Article  CAS  Google Scholar 

  21. Bustamante LY, Bartholdson SJ, Crosnier C, Campos MG, Wanaguru M, Nguon C. A full-length recombinant Plasmodium falciparum PfRH5 protein induces inhibitory antibodies that are effective across common PfRH5 genetic variants. Vaccine. 2013;31:373–9.

    Article  CAS  Google Scholar 

  22. Sironi M, Cagliani R, Forni D, Clerici M. Evolutionary insights into host-pathogen interactions from mammalian sequence data. Nat Rev Genet. 2015;16:224–36.

    Article  CAS  Google Scholar 

  23. Otsuk H, Kaneko O, Thongkukiatkul A, Tachibana M, Iriko H, Takeo S, et al. Single amino acid substitution in Plasmodium yoelii erythrocyte ligand determines its localization and controls parasite virulence. Proc Natl Acad Sci USA. 2009;106:7167–72.

    Article  Google Scholar 

  24. Ramamurthy M, Sankar S, Abraham AM, Nandagopal B, Sridharan G. B cell epitopes in the intrinsically disordered regions of neuraminidase and hemagglutinin proteins of H5N1 and H9N2 avian influenza viruses for peptide-based vaccine development. J Cell Biochem. 2019;120(10):17534–44.

    PubMed  CAS  Google Scholar 

Download references


The authors are grateful to the patients who contributed samples for this study, and to the staff, nurses and doctors of the hospitals visited who facilitated sample collection. Our sincere gratitude and appreciations go to Department of Biomedical Chemistry, The University of Tokyo for the immense contribution to the gene sequencing aspect of this study. We are grateful to the patients who contributed samples for this study and Authors appreciate the Department of Microbiology of the Nigerian Institute of Medical Research, Nigeria for their support. EOB is a FLAIR Research Fellow of The Royal Society, UK, and supported by a Global Challenge Research Fund (GCRF) Grant (No. FLR\R1\190353) through a partnership between the African Academy of Sciences and the Royal Society, UK.


This work was funded personally.

Author information

Authors and Affiliations



OA conceived the study and carried out the laboratory and molecular analysis. AAO and EOB supervised the work and revised the manuscript, BAI, KH and OA handled data analysis and interpretation. OA wrote the manuscript. KMO, YW, GUE, KK and AA revised the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Olusola Ajibaye.

Ethics declarations

Ethics approval and consent to participate

Ethical approval (Ref No: IRB/10/136) to conduct this study was obtained from the Institutional Review Board of the Nigerian Institute of Medical Research, Lagos. Permission was obtained from the different health facilities where the study was carried out. All the research participants gave written informed consent, guardian/parental consent and/or accent for participants less than 16 years (children). The adults provided written informed consent and the parents/guardians gave informed consent on behalf of their children/wards. Participants who declined participation were not denied treatment. Generally, the study was conducted in line with the principle of Good Clinical Laboratory Practice (GCLP) and the 1994 Helsinki Declaration.

Consent for Publication

Not applicable.

Competing interests

We hereby state that there is no conflict of interest regarding this work. There is no financial or non-financial competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Mutations in PfRH5 sequences from Lagos, Nigeria.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ajibaye, O., Osuntoki, A.A., Balogun, E.O. et al. Genetic polymorphisms in malaria vaccine candidate Plasmodium falciparum reticulocyte-binding protein homologue-5 among populations in Lagos, Nigeria. Malar J 19, 6 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Polymorphisms
  • Reticulocyte-binding protein homologue-5
  • Haplotypes
  • Histocompatibility
  • Gene flow
  • Linkage
  • Unstructured regions