Genetic variability of Plasmodium falciparum histidine-rich proteins 2 and 3 in Central America
Malaria Journal volume 18, Article number: 31 (2019)
Malaria is an important disease in many tropical countries. Rapid diagnostic tests (RDTs) are valuable tools for diagnosing malaria in remote areas. The majority of RDTs used for the diagnosis of Plasmodium falciparum are based on the detection of the specific histidine-rich proteins (PfHRP2 and PfHRP3). During the last decade, the threat posed by the lack of expression of these antigens and the variability of the proteins on the diagnosis of malaria has been widely discussed. The aim of this study was to evaluate the genetic diversity of pfhrp2 and pfhrp3 of P. falciparum isolates collected in three Central American countries.
DNA samples were amplified and sequenced to assess the diversity of nucleotides and amino acids. A search for known epitopes within the amino acid sequence was carried out, and the sensitivity of the sequences was evaluated according to a predictive model. A phylogenetic analysis was carried out including homologous sequences from different regions of the world. Protein structures were predicted in silico.
Five different patterns for PfHRP2 and one pattern for PfHRP3 were identified. Isolates from Central America show a high level of genetic diversity in pfhrp2; however, the amino acid sequences seem to contain enough motifs to be detected by the RDTs currently available.
It is unlikely that the variability of the pfhrp2 and pfhrp3 genes has a significant impact on the ability of the RDTs to detect the PfHRP antigens in Central America.
Malaria is a tropical disease responsible for a vast burden on public health in many developing countries. Although the incidence of malaria cases is still considerably high in Africa, Southeast Asia and the Amazon Basin, the Central American sub-region has experienced a notable decrease in the number of cases in recent decades . These achievements have led to propose the elimination of malaria by the year 2030 by the countries of Central America and the island of Hispaniola .
Microscopy is the gold standard method for malaria diagnosis and has shown excellent performance when adequate infrastructure conditions and trained personnel are available . However, due to the decrease in the number of cases in a geographic region, it is necessary to have diagnostic alternatives that do not require highly qualified personnel and equipment.
Rapid diagnostic tests (RDTs) were developed in the 1990s  and are valuable tools in the timely diagnosis of malarial infections in geographically remote areas. Currently available RDTs for the diagnosis of malaria are based on the detection of one or a combination of parasite antigens circulating in the blood of infected patients: the Plasmodium falciparum histidine-rich proteins PfHRP2 and PfHRP3, Plasmodium lactate dehydrogenase (pLDH), and Plasmodium aldolase . PfHRP2 and PfHRP3 are the antigens most commonly used to detect Plasmodium falciparum infections and most of the RDT products (> 90%) available in the market use specific antibodies against both proteins [5,6,7]. PfHRP proteins are common target antigens for RDT manufacturing companies because of their abundance in blood during blood-stage infections , its specificity, and the presence of repetitive epitopes that enable their detection by multiple antibodies, increasing the sensitivity of the technique [9, 10]. PfHRP2 and its homologous PfHRP3 are soluble proteins encoded in genes with two exons interrupted by one intron, with a peptide signal encoded in exon 1 and a histidine-alanine rich repeat region in exon 2 [11, 12].
A recent study reported 21% of P. falciparum isolates from three Central American countries showing deletions of the exon 1-intron 1 segment of pfhrp2 and pfhrp3, and 6.3% of the isolates lacked the exon 2 region of both genes .
Parasite failure to express the antigen is a well-demonstrated factor that affects the sensitivity of the PfHRP-based RDTs , but it has also been suggested that variations in the amino acid repeats within exon 2 of PfHRP2 could affect the sensitivity [15,16,17], which could have negative repercussions for malaria control and elimination programs. For this reason, a collection of parasite isolates from Central America was analysed to assess the genetic diversity within exon 2 of pfhrp2 and pfhrp3, in order to provide useful information to better understand the possible effect of the diversity of parasites in the Central American region on the sensitivity of RDTs based on the detection of PfHRP antigens.
The aim of the study was to evaluate the genetic diversity within exon 2 of pfhrp2 and pfhrp3 genes of P. falciparum isolates collected in three Central American countries. A convenience sample was collected from a previous study that assessed the deletions of both genes and its flanking regions . In this study, a subset of DNA samples that previously amplified the coding region of one or both genes were selected for further analysis. Thirty-five samples of pfhrp2 (Honduras = 16, Nicaragua = 16, Guatemala = 3) and 5 samples of pfhrp3 (Honduras = 3, Nicaragua = 2) were included. Samples showed heterogeneous parasitic densities. All the samples showed a previous positive result by a PfHRP2-based RDT. Since each of the three countries uses a different RDT brand for the routine diagnosis of malaria, the samples were analysed by the following RDTs: CareStart™ Malaria HRP2/pLDH(Pf/Pv) Combo (Honduras), SD Bioline MALARIA Ag P.f/P.v (Nicaragua) and CareStart™Malaria RDT Single Kit (Guatemala).
Blood samples of patients infected with malaria falciparum were collected on Whatman™ filter paper (GE Healthcare Bio-Sciences Corp, NJ, USA) for routine chloroquine resistance surveillance [18, 19]. DNA was extracted through a Chelex-100 based method . A PCR method was used to confirm parasite species .
DNA amplification, sequencing and in silico translation
Exons 2 of the pfhrp2 and pfhrp3 genes were amplified by a semi-nested PCR approach  (Fig. 1). All amplification reactions were carried out in a volume of 50 µl, with 25 µl of 2X Master mix (Promega Corp.), 2.0 μl of each primer (Additional file 1) at a concentration of 10 μM and 1.0 μl of template DNA.
Reactions of both the primary and secondary PCR were carried out by an initial denaturation temperature at 94 °C for 10 min, followed by 35 cycles (37 cycles for the secondary PCR) of 94 °C for 50 s, an annealing step of 55 °C for 50 s, 72 °C for 1 min; and a final extension step at 72 °C for 10 min. The second amplification product was observed by agarose electrophoresis. Expected PCR products ranged between 720 and 830 bp for pfhrp2 and < 500 bp for pfhrp3. Sequencing was performed with the forward and reverse primers of the secondary PCR at the Macrogen facilities (macrogenusa.com). Sequences were deposited at NCBI and accession numbers were assigned.
Geneious v.9.1.7 (Biomatters Ltd., Auckland, New Zealand) and Mega7  v7.0.26 were used for sequence editing and analysis. The sequence length, nucleotide composition, and degree of polymorphism (π) were calculated. π was calculated as the average number of nucleotide differences per site between DNA sequences in all possible pairs in the sample population. Nucleotide sequences were translated in silico to corresponding amino acids (aa) using the correct open reading frame. Amino acid repeats were given a numeric code (1–20) according to previous reports [16, 22]. The number and frequency of each aa repeat were estimated. The number of type 2 repeats was multiplied by the number of type 7 repeats (2 × 7) in each sequence in order to predict the PfHRP2 detection sensitivity according to Baker’s regression model [15, 16]. Sequences were classified as “sensitive” if the score produced by the number of type 2 times type 7 repeats was greater than 43 and “non-sensitive” if it was less than 43 repeats. Minimal epitopes recognized by 11 PfHRP-specific monoclonal antibodies (MAbs) were searched and quantified in each PfHRP sequence .
The ClustalW tool was used to align the pfhrp2nucleotide sequences. The alignment was used to build a cladogram using the UPGMA (Unweighted Pair Group Method with Arithmetic Mean) method, with a bootstrap of 1000 replicates. A cladogram was constructed including only the sequences obtained in this study. A second set of cladograms was constructed including homologous sequences downloaded from the GenBank database at the National Center for Biotechnology Information (NCBI).
Sequences downloaded from NCBI were classified according to the geographical origin of the parasites based on four WHO Health Regions: Africa, Americas, Western Pacific, and South East Asia. Evolutionary distances were computed using the Maximum Composite Likelihood Method. Non-homologous sequences were selected as outgroups for the cladogram construction.
Protein structure prediction
The secondary structure of PfHRP2 and PfHRP3 aa sequences was predicted through the template-based software RaptorX (http://raptorx.uchicago.edu). Sizes of partial aa sequences are shown in Table 1. The p value and the overall unnormalized global distance test (uGDT) were calculated for each model . Two software were used to predict 3D models for the polypeptide sequences: (a) Protein Homology/analogy Recognition Engine v2.0 (PHYRE2)  (http://www.sbg.bio.ic.ac.uk/phyre2/html/page.cgi?id=index), and (b) I-TASSER (https://zhanglab.ccmb.med.umich.edu/I-TASSER/) .
Nucleotide composition and diversity of pfhrp2 and pfhrp3 genes
Amplification products of 35 and 5 isolates were sequenced for pfhrp2 and pfhrp3 genes respectively. Five different pfhrp2 patterns were identified (Accession Nos MH996688-95), while only one pattern was obtained for pfhrp3 (MH996696-97) (Table 1). No correlation was observed between the nucleotide sequence and the geographical origin of the isolates. Pattern I of pfhrp2 was the most frequently observed in the three countries (overall 62.9%). Only one parasite isolate had a unique pfhrp2 sequence (pattern IV) and the remaining 4 sequence patterns were present in more than 3 isolates. As shown in Table 1, the size of the repeat sequence of exon 2 ranged from 720 to 828 nucleotides (240 to 276 aa) in pfhrp2, and it was shorter in pfhrp3 with only 484 base pairs (161 aa). The most frequent aa within PfHRP2 sequences were Ala (42.4–44.4%), His (39.7–41%), and Asp (10.8–11.6%) (Additional file 2).
Although the five different patterns of pfhrp2 revealed a similar nucleotide composition (Additional file 2), the nucleotide diversity (π) within the 35 isolates indicates this is a highly polymorphic gene (π = 0.044601). π was also calculated for 65 homologous pfhrp2 sequences isolated from 4 WHO Health Regions with a result of 0.049603.
Variation in amino acid repeat patterns
Eleven and eight amino acid repeats were identified from pfhrp2 and pfhrp3 respectively (Table 2 and Fig. 2). Repeats 2, 3, 5—8, 10, and 12 were observed in all the patterns of pfhrp2. Repeat type 1 was absent from pattern II and repeats 4 and 18 were present only in pattern IV. Repeats 2 (36%) and 7 (25%) were the most frequent among the PfHRP2 sequences. All patterns, with the exception of pattern II, started with type 1 repeat and all patterns ended with the 10–12 repeat combination.
The repeats’ pattern of the PfHRP3 gene showed a very different structure compared to PfHRP2. The most common repeat within PfHRP3 was type 16 (40%) followed by type 17 (25%). Types 2, 3, 5, 6, 8, 10, and 12 repeats were absent from PfHRP3, and type 7 represented only 5% of repeats. Motifs 15, 16, 17, and 20 were unique of PfHRP3 and absent from PfHRP2 sequences. A non-repetitive sequence of 28 amino acids interspersed between two blocks of repeats was present in all PfHRP3 sequences (ANHGFHFNLHDNNSHTLHHAKANACFDD) (Fig. 2).
Prediction of sensitivity through the Bakers’ model
A predictive model proposed to assess the detection sensitivity of an RDT under parasite densities > 250 μl−1 was used for the 5 PfHRP2 patterns and the PfHRP3 sequence . Baker’s model is based on a binary logistic regression as a function of the number of type 2× type 7 repeats and classifies a sequence as “sensitive” if the score or product of multiplication is higher than 43. The five PfHRP2 patterns produced scores ranging from 78 to 110. According to Baker’s predictive model, all the sequences could be detected with parasitaemia > 250 μl−1. The PfHRP3 sequence lacks type 2 repeats and reveals only one type 7 repetition.
Epitopes targeted by monoclonal antibodies in RDTs
Linear epitopes (motifs) recognized by eleven MAbs  were searched in the sequences obtained in this study. Each MAb mainly recognize a dominant epitope (listed in Additional file 3) although they may be able to bind to other epitopes as well. The most frequent motif among the five patterns of PfHRP2 was AHHAADAHH, recognized by the MAb C1-13. This motif was present from 19 to 21 times in each pattern. The second and third most frequent motifs were AHHAHHA (14 to 16 times) and DAHHAADAHHA (7 to 10 times). Motifs HAHHAHHAADAHH and AYAHHAHHAAY were absent in all the patterns. Only two motifs (AHHAHHA and AHHAHHV) were present in the PfHRP3 pattern (Additional file 3).
BLAST and cluster analyses
The five pfhrp2 sequence patterns obtained in this study were compared with 65 similar sequences available in NCBI using the BLAST tool. The pattern I revealed 99% identity with the Honduran reference strain HB3 (Acc. No AY816261). Patterns II and III shared 97–99% similarity with Acc. No FJ871227, obtained from a strain isolated in Honduras. Patterns IV and V were 98% similar to the sequence KU723614 of Colombian origin. In contrast, the unique pattern of pfhrp3 showed the highest similarity (99%) with a homologous sequence from Madagascar (Acc. No EU589824).
The putative amino acid sequences corresponding to the five patterns of the pfhrp2 gene from Central America were aligned in order to construct a cladogram that graphically revealed the similarities between them. As shown in Fig. 3, patterns IV and V revealed a high degree of similarity and their only difference was the insertion of a block of three repetitions (4-18-7) in the middle part of the pattern IV (Fig. 2).
For a better understanding of the similarities between pfhrp2 and pfhrp3 sequences obtained in this study and those from parasites from other regions of the world, 65 pfhrp2 and 28 pfhrp3 homologous sequences were downloaded from NCBI. Sequences were classified according to their geographical origin and according to WHO Health geographic Regions (for pfhrp2: African Region n = 27, Americas n = 11, Western Pacific n = 8, South East Asia n = 17, unknown n = 2; for pfhrp3: African Region n = 13, Americas n = 4, Western Pacific n = 7, South East Asia n = 4). A non-homologous sequence was included in each cladogram as an outgroup (Fig. 4). All sequences were aligned through ClustalW and a translation alignment method. Translation alignments were performed to nucleotide sequences that were translated into protein. Both alignments were used to construct separate cladograms for pfhrp2 and pfhrp3.
No clear geographic clustering is observed with any of the two cladograms of both genes. However, sequences from Central America and sequences from some American countries cluster together for pfhrp2.
Protein structure prediction
Based on their sequence and homology, the protein secondary structure was predicted for five partial PfHRP2 patterns and one PfHRP3 partial sequence, through an online template-based software. Partial amino acid sequences of exon 2 downloaded from NCBI were also analysed (Accession numbers AY816261 and KC558597). Complete coding sequences (CDS) including 305 amino acids for PfHRP2 (Acc. No XM_002808697) and 264 amino acids for PfHRP3 (Acc. No U69552) were also predicted for secondary structure using the RaptorX software.
As shown in the additional file 5, the predicted secondary structure for all the polypeptides matches with helices and coils. p-value and uGDT were calculated for each model. p-value is the likelihood of a predicted model being worse than the best of a set of randomly-generated models for a protein. Values < 10−3 are a good indicator of relative quality of a model. uGDT measures the absolute model quality. Values > 50 are a good indicator of quality. The p-values of all models were slightly higher than 10−3, and the uGDT values ranged between 28 and 75. Five models obtained an overall uGDT > 50, therefore the robustness of these predictions is partially reliable (Additional file 4).
Three dimensional models of protein sequences could not be meaningfully predicted with the Phyre2 software since 0% of residues modelled at > 90% confidence, and ~ 89% of each sequence was predicted disordered.
The I-TASSER software was used to predict protein models for complete CDS of PfHRP2 (XM_002808697) and PfHRP3 (U69552). Five models were predicted for each polypeptide (Fig. 5). A C-score was calculated for each model. C-score is a confidence score for estimating the quality of predicted models. It is calculated based on the significance of threading template alignments and the convergence parameters of the structure assembly simulations. C-score is typically in the range of − 5 to 2, where a C-score of higher value signifies a model with a high confidence and vice versa.
The five PfHRP2 predicted models revealed C-scores ranged between − 2.78 to − 3.69 (Model 1: − 2.88, Model 2: − 2.87, Model 3: − 3.06, Model 4: − 3.13, and Model 5: − 3.69). For PfHRP3, the following C-scores were calculated for each model (Model 1: − 3.08, Model 2: − 3.32, Model 3: − 3.14, Model 4: − 2.86, and Model 5: − 2.41).
Four factors have been described as causes of failure for PfHRP-based RDTs for the diagnosis of Plasmodium falciparum infections: (a) Low parasitaemia, (b) gene deletions, (c) lack of antigen expression, and (d) polymorphisms in the amino acid sequences . Previous studies reported the existence of parasite isolates with deletions in the pfhrp2 and pfhrp3 genes for the first time in Central America [13, 27]. This finding could hinder malaria control strategies relying on the detection of PfHRP2 and PfHRP3 antigens. In order to assess whether variations in the amino acid sequences of PfHRP2 and PfHRP3 could cause an added risk for the efficacy of the RDTs in Central America, the exon 2 of both genes was sequenced.
Five different pfhrp2 sequences and a single pfhrp3 sequence were identified from 35 and 5 P. falciparum isolates, respectively. According to this result, the isolates from Central America show a high level of genetic diversity in pfhrp2 even though substantially lower than that reported in parasites from Africa and Asia. Sequence variations in both genes from isolates worldwide have been recently reviewed [28, 29]. When studying 458 isolates collected from 38 countries, African parasites showed greater diversity according to the proportion of unique PfHRP2 sequences per country. American isolates, on the other hand, were the least diverse. A similar trend was observed when pfhrp3 sequences were studied from 80 isolates . The highest levels of genetic variation (87%–98%) [15, 17, 29] were also reported in African countries in a review that included isolates from America, Africa and the Asia–Pacific region . The nucleotide diversity (π) within sequences from Central America and worldwide confirmed that pfhrp2 is highly polymorphic regardless of the geographic origin of the parasites .
The pfhrp2 and pfhrp3 genes codify for tandem repeats of amino acids. Tandem repeats are unstable genetic elements within genomes associated with rapid evolutionary changes . pfhrp2 and pfhrp3 are located in the sub-telomeric/telomeric regions of chromosomes 8 and 13, respectively. Telomeric regions of Plasmodium are highly susceptible to changes due to homologous recombination [31, 32]. This could explain such a high gene diversity regardless of the geographic region . In order to better understand the mechanisms of geographical dispersion of P. falciparum, an alignment of sequences obtained in this study was performed and compared with sequences of parasites from other regions of the world. A set of cladograms was constructed but no clear geographic clustering was observed. This apparent lack of clustering in the populations of the parasite suggests that pfhrp2 and pfhrp3 sequences are not under strong evolutionary selective pressure nor offer survival advantages for the species .
Some authors have hypothesized that variations in the sequence of PfHRP2 may affect the sensitivity of RDTs because of the polymorphisms of amino acid repeats [15, 29, 33]. Parasite isolates from Central America revealed eleven amino acid repeats in PfHRP2, all of them previously described in strains from other geographical regions [15, 22, 33,34,35,36]. The two most common motifs were type 2 (AHHAHHAAD) and type 7 (AHHAAD). The abundance of these motifs has been associated with an increase in the sensitivity of RDTs [15, 16, 36]. Both motifs have been proposed as potential epitopes targeted by monoclonal antibodies [16, 17, 23, 29] assuming that a higher epitope frequency within the sequence may result in greater sensitivity . Bakers’ regression model  is based on the number of type 2 multiplied by the number of type 7 (2 × 7) repeats and it was applied in a case study from Uganda that reported a false negative RDT due to a low value of 2 × 7 repeats . If this model has a good predictive capacity, the five PfHRP2 patterns described for Central American parasites could be detected even at low parasitaemia levels. However, there is evidence that suggests that Bakers’ model does not explain all the failures of non-detection of RDTs based exclusively on the score produced by the number of 2 × 7 repeats [17, 35, 38] and that the sequence diversity of repeats is not likely to negatively affect the performance of RDTs [29, 38]. Moreover, the amino acid sequence found in PfHRP3 is much less diverse than PfHRP2 which is consistent with reports elsewhere [15,16,17, 29, 36].
Specific antibodies against PfHRP2 cross react with PfHRP3 [12, 39, 40] but in circulating parasites from Central America this sequence lacks repeat 2 whilst repeat 7 is underrepresented. Furthermore, more than 90% of the isolates have deleted the loci associated to pfhrp3 , therefore the expression of PfHRP3 does not seem to contribute significantly to the sensitivity of the RDTs in Central America.
A second approach was applied to attempt to correlate the amino acid sequences with the potential capacity of an RDT to detect the PfHRP2 antigen . Eleven epitopes ranging from 8 to 15 amino acids which are recognized by commercially available monoclonal antibodies (MAbs)  were searched within the five PfHRP2 patterns and the unique PfHRP3 sequence. The five PfHRP2 sequence patterns revealed four motifs present at a high frequency, and five other motifs present in a lower proportion in almost all sequences (Additional file 3). The two most frequent epitopes were recognized by the MAbs C1-13 and 3A4. Similar results were reported in a study performed with 120 samples from Papua New Guinea . Given the absence of two motifs (HAHHAHHAADAHH, AYAHHAHHAAY) within the sequences reported here, RDTs using MAbs C2-3 and Genway might be less sensitive, as suggested by Lee et al. . In any case, it is important to remember that MAbs recognize dominant epitopes more avidly but in turn are able to recognize other epitopes with a certain degree of similarity to their dominant epitope. These results are further evidence that, despite the high genetic diversity, the RDTs based on the detection of PfHRP2 show an adequate sensitivity whenever the parasites express the antigen. The PfHRP3 sequence includes only 2 out of 11 motifs recognized by MAbs, consequently, this antigen could contribute less to the sensitivity of the RDT.
The three-dimensional structure of PfHRP2 and PfHRP3 have not been elucidated. Because sensitivity of an RDT depends on the binding of MAbs to specific epitopes, which can be linear or conformational, two softwares were applied to predict in silico the secondary and tertiary structure of the PfHRP proteins. The predicted structures could be useful to better understand the topographic availability of epitopes for MAbs.
According to the RaptorX software, uGDT values > 50 are good indicators of absolute model quality. The highest value was obtained for the complete coding sequence (CDS) of PfHRP2 (uGDT = 75). Four partial amino acid sequences encoded in exon 2 showed values > 50. These values suggest that these proteins could be organized mainly as helices as shown in Additional file 5. Additionally, folded protein models predicted by I-TASSER showed low C-scores limiting the robustness of the models.
Bioinformatic analysis are preliminary tools for predicting protein structures, however, the predicted models indicate that although some epitopes could be hidden inside the protein, there are enough epitopes exposed on the surface available for the binding of specific MAbs.
One limitation of the study is that all the samples were diagnosed as positive through a PfHRP2-based RDT, and different results could have been observed if samples with a negative RDT and a positive PCR or microscopy test would have been included. A second limitation of the study is the small number of samples analysed per country. According to the World Malaria Report 2018 , the percentage of samples analysed in relation to the number of cases of P. falciparum infections during the year in which the samples were collected was of 12.5% (16/128) for Honduras, 5.9% (3/51) for Guatemala, and 4.7% (16/342) for Nicaragua. This limitation could not be overcome due to logistical obstacles, and consequently the objective of the study did not have epidemiological connotations.
Despite these limitations this study provides insights into the genetic variability of pfhrp2 and pfhrp3 in isolates from three countries in Central America and shows that PfHRP2 is highly polymorphic. It is unlikely that this variability has a significant impact on the sensitivity of currently available PfHRP-based RDTs.
rapid diagnostic test
Plasmodium falciparum histidine-rich protein 2 gene
Plasmodium falciparum histidine-rich protein 3 gene
Plasmodium falciparum histidine rich protein
Plasmodium lactate dehydrogenase
polymerase chain reaction
unweighted pair group method with arithmetic mean
National Center for Biotechnology Information
World Health Organization
unnormalized global distance test
- Acc. No:
Basic Local Alignment Search Tool
complete coding sequences
WHO. World malaria report 2017. Geneva: World Health Organization; 2017.
WHO. Global TECHNICAL STRATEGY FOR MALARIA 2016–2030. Geneva: World Health Organization; 2015.
Fontecha GA, Mendoza M, Banegas E, Poorak M, De Oliveira AM, Mancero T, et al. Comparison of molecular tests for the diagnosis of malaria in Honduras. Malar J. 2012;11:119.
Dietze R, Perkins M, Boulos M, Luz F, Reller B, Corey GR. The diagnosis of Plasmodium falciparum infection using a new antigen detection system. Am J Trop Med Hyg. 1995;52:45–9.
WHO. Malaria rapid diagnostic test performance: summary results of WHO product testing of malaria RDTs: rounds 1–6 (2008–2015). Geneva: World Health Organization; 2016.
Wongsrichanalai C, Barcus MJ, Muth S, Sutamihardja A, Wernsdorfer WH. A review of malaria diagnostic tools: microscopy and rapid diagnostic test (RDT). Am J Trop Med Hyg. 2007;77:119–27.
Mouatcho JC, Goldring JP. Malaria rapid diagnostic tests: challenges and prospects. J Med Microbiol. 2013;62:1491–505.
Moody A. Rapid diagnostic tests for malaria parasites. Clin Microbiol Rev. 2002;15:66–78.
Wellems TE, Howard RJ. Homologous genes encode two distinct histidine-rich proteins in a cloned isolate of Plasmodium falciparum. Proc Natl Acad Sci USA. 1986;83:6065–9.
Panton LJ, McPhie P, Maloy WL, Wellems TE, Taylor DW, Howard RJ. Purification and partial characterization of an unusual protein of Plasmodium falciparum: histidine-rich protein II. Mol Biochem Parasitol. 1989;35:149–60.
Rock EP, Marsh K, Saul AJ, Wellems TE, Taylor DW, Maloy WL, et al. Comparative analysis of the Plasmodium falciparum histidine-rich proteins HRP-I, HRP-II and HRP-III in malaria parasites of diverse origin. Parasitology. 1987;95:209–27.
Howard RJ, Uni S, Aikawa M, Aley SB, Leech JH, Lew AM, et al. Secretion of a malarial histidine-rich protein (Pf HRP II) from Plasmodium falciparum-infected erythrocytes. J Cell Biol. 1986;103:1269–77.
Fontecha G, Mejia RE, Banegas E, Ade MP, Mendoza L, Ortiz B, et al. Deletions of pfhrp2 and pfhrp3 genes of Plasmodium falciparum from Honduras, Guatemala and Nicaragua. Malar J. 2018;17:320.
Gamboa D, Ho MF, Bendezu J, Torres K, Chiodini PL, Barnwell JW, et al. A large proportion of P. falciparum isolates in the Amazon region of Peru lack pfhrp2 and pfhrp3: implications for malaria rapid diagnostic tests. PLoS One. 2010;5:e8091.
Mariette N, Barnadas C, Bouchier C, Tichit M, Menard D. Country-wide assessment of the genetic polymorphism in Plasmodium falciparum and Plasmodium vivax antigens detected with rapid diagnostic tests for malaria. Malar J. 2008;7:219.
Baker J, McCarthy J, Gatton M, Kyle DE, Belizario V, Luchavez J, et al. Genetic diversity of Plasmodium falciparum histidine-rich protein 2 (PfHRP2) and its effect on the performance of PfHRP2-based rapid diagnostic tests. J Infect Dis. 2005;192:870–7.
Wurtz N, Fall B, Bui K, Pascual A, Fall M, Camara C, et al. Pfhrp2 and pfhrp3 polymorphisms in Plasmodium falciparum isolates from Dakar, Senegal: impact on rapid malaria diagnostic tests. Malar J. 2013;12:34.
Fontecha GA, Sanchez AL, Mendoza M, Banegas E, Mejia-Torres RE. A four-year surveillance program for detection of Plasmodium falciparum chloroquine resistance in Honduras. Mem Inst Oswaldo Cruz. 2014;109:492–3.
Mejia Torres RE, Banegas EI, Mendoza M, Diaz C, Bucheli ST, Fontecha GA, et al. Efficacy of chloroquine for the treatment of uncomplicated Plasmodium falciparum malaria in Honduras. Am J Trop Med Hyg. 2013;88:850–4.
de Lamballerie X, Zandotti C, Vignoli C, Bollet C, de Micco P. A one-step microbial DNA extraction method using “Chelex 100” suitable for gene amplification. Res Microbiol. 1992;143:785–90.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis Version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.
Bharti PK, Chandel HS, Ahmad A, Krishna S, Udhayakumar V, Singh N. Prevalence of pfhrp2 and/or pfhrp3 Gene deletion in Plasmodium falciparum population in eight highly endemic states in India. PLoS ONE. 2016;11:e0157949.
Lee N, Gatton ML, Pelecanos A, Bubb M, Gonzalez I, Bell D, et al. Identification of optimal epitopes for Plasmodium falciparum rapid diagnostic tests that target histidine-rich proteins 2 and 3. J Clin Microbiol. 2012;50:1397–405.
Kallberg M, Wang H, Wang S, Peng J, Wang Z, Lu H, et al. Template-based protein structure modeling using the RaptorX web server. Nat Protoc. 2012;7:1511–22.
Kelley LA, Mezulis S, Yates CM, Wass MN, Sternberg MJ. The Phyre2 web portal for protein modeling, prediction and analysis. Nat Protoc. 2015;10:845–58.
Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER Suite: protein structure and function prediction. Nat Methods. 2015;12:7–8.
Abdallah JF, Okoth SA, Fontecha GA, Torres RE, Banegas EI, Matute ML, et al. Prevalence of pfhrp2 and pfhrp3 gene deletions in Puerto Lempira, Honduras. Malar J. 2015;14:19.
Gendrot M, Fawaz R, Dormoi J, Madamet M, Pradines B. Genetic diversity and deletion of Plasmodium falciparum histidine-rich protein 2 and 3: a threat to diagnosis of P. falciparum malaria. Clin Microbiol Infect. 2018;1:1. https://doi.org/10.1016/j.cmi.2018.09.009 (Epub ahead of print).
Baker J, Ho MF, Pelecanos A, Gatton M, Chen N, Abdullah S, et al. Global sequence variation in the histidine-rich proteins 2 and 3 of Plasmodium falciparum: implications for the performance of malaria rapid diagnostic tests. Malar J. 2010;9:129.
Deme AB, Park DJ, Bei AK, Sarr O, Badiane AS, Gueye Pel H, et al. Analysis of pfhrp2 genetic diversity in Senegal and implications for use of rapid diagnostic tests. Malar J. 2014;13:34.
Scherf A, Mattei D. Cloning and characterization of chromosome breakpoints of Plasmodium falciparum: breakage and new telomere formation occurs frequently and randomly in subtelomeric genes. Nucleic Acids Res. 1992;20:1491–6.
Anders RF, Brown GV, Coppel RL, Kemp DJ. Repeat structures in malaria antigens. P N G Med J. 1986;29:87–93.
Kumar N, Singh JP, Pande V, Mishra N, Srivastava B, Kapoor R, et al. Genetic variation in histidine rich proteins among Indian Plasmodium falciparum population: possible cause of variable sensitivity of malaria rapid diagnostic tests. Malar J. 2012;11:298.
Atroosh WM, Al-Mekhlafi HM, Al-Jasari A, Sady H, Al-Delaimy AK, Nasr NA, et al. Genetic variation of pfhrp2 in Plasmodium falciparum isolates from Yemen and the performance of HRP2-based malaria rapid diagnostic test. Parasit Vectors. 2015;8:388.
Li P, Xing H, Zhao Z, Yang Z, Cao Y, Li W, et al. Genetic diversity of Plasmodium falciparum histidine-rich protein 2 in the China–Myanmar border area. Acta Trop. 2015;152:26–31.
Nderu D, Kimani F, Thiong’o K, Akinyi M, Karanja E, Meyer CG, et al. PfHRP2-PfHRP3 diversity among Kenyan isolates and comparative evaluation of PfHRP2/pLDH malaria RDT with microscopy and nested PCR methodologies. Parasitol Int. 2018;67:793–9.
Wurtz N, Briolant S, Lemarie D, PommierdeSanti V, Pascual A, Roodt T, et al. Delayed diagnosis of Plasmodium falciparum in a soldier in Uganda: false-positive rapid diagnostic test associated with reduced repeats in pfhrp2. Med Sante Trop. 2013;23:181–4 (in French).
Kumar Bharti P, Singh Chandel H, Krishna S, Nema S, Ahmad A, Udhayakumar V, et al. Sequence variation in Plasmodium falciparum histidine rich proteins 2 and 3 in Indian isolates: implications for malaria rapid diagnostic test performance. Sci Rep. 2017;7:1308.
Lee N, Baker J, Andrews KT, Gatton ML, Bell D, Cheng Q, et al. Effect of sequence variation in Plasmodium falciparum histidine-rich protein 2 on binding of specific monoclonal antibodies: implications for rapid diagnostic tests for malaria. J Clin Microbiol. 2006;44:2773–8.
Sharma YD. Genomic organization, structure and possible function of histidine-rich proteins of malaria parasites. Int J Biochem. 1988;20:471–7.
Willie N, Zimmerman PA, Mehlotra RK. Plasmodium falciparum histidine-rich protein 2 gene variation in a malaria-endemic area of Papua New Guinea. Am J Trop Med Hyg. 2018;99:697–703.
WHO. World malaria report 2018. Geneva: World Health Organization; 2018.
GF, AP, BO, and GM contributed with the study design; AP, GM, and DE performed the experiments. All authors read and approved the final manuscript.
This study was possible thanks to the financial support of DICU-UNAH. Samples were provided by national Health officials of Honduras, Guatemala and Nicaragua.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The need for ethics approval was waived by the CEI-MEIZ ethics committee of the UNAH because the samples used in this study were obtained from the routine diagnosis of malaria carried out by the health ministries of Guatemala, Nicaragua and Honduras. Nevertheless, this study considered all the ethical postulates of the Geneva Declaration.
Funding for this study was provided by the Scientific Research Office of the UNAH (DICYP-UNAH), Honduras, and the office of the Pan American Health Organization in Honduras. The funding institutions did not participate in the design of the study and the collection, analysis and interpretation of the data or in the drafting of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of the primer’s sequences used to amplify the exon 2 of pfhrp2 and pfhrp3.
Nucleotide and amino acid composition of five pfhrp2 and 1 pfhrp3 sequences.
Frequency of epitope motifs present in the sequences of PfHRP2 and PfHRP3.
Template-based secondary structure prediction of PfHRP2 (a–f), and PfHRP3 (g–h) partial proteins using the RaptorX web server. Letters (a–e) show the five patterns of PfHRP2 described in this study; (f) partial sequence of Plasmodium falciparum reference strain 3D7 (Accession No XM_002808697); (g) PfHRP3 sequence from Central American isolates; (h) Plasmodium falciparum strain N569 (Accession No KC558597); (i) Complete CDS of PfHRP2 from Plasmodium falciparum strain 3D7 (Accession No XM_002808697); (j) Complete CDS of PfHRP3 from Plasmodium falciparum accession No U69552. The boxes on the right indicate the prediction of the secondary structure. The red boxes show the probability of helix structures and the grey boxes indicate the probability of coil formation for each amino acid position.
p-value and uGDT of secondary structures of models predicted for PfHRP2 and PfHRP3.
About this article
Cite this article
Fontecha, G., Pinto, A., Escobar, D. et al. Genetic variability of Plasmodium falciparum histidine-rich proteins 2 and 3 in Central America. Malar J 18, 31 (2019). https://doi.org/10.1186/s12936-019-2668-3
- Plasmodium falciparum
- Rapid diagnostic test