Genetic Polymorphisms in Plasmodium vivax Dihydrofolate Reductase and Dihydropteroate Synthase in Isolates from the Philippines,Bangladesh, and Nepal

Genetic polymorphisms of pvdhfr and pvdhps genes of Plasmodium vivax were investigated in 83 blood samples collected from patients in the Philippines, Bangladesh, and Nepal. The SNP-haplotypes of the pvdhfr gene at the amino acid positions 13, 33, 57, 58, 61, 117, and 173, and that of the pvdhps gene at the positions 383 and 553 were analyzed by nested PCR-RFLP. Results suggest diverse polymorphic patterns of pvdhfr alone as well as the combination patterns with pvdhps mutant alleles in P. vivax isolates collected from the 3 endemic countries in Asia. All samples carried mutant combination alleles of pvdhfr and pvdhps. The most prevalent combination alleles found in samples from the Philippines and Bangladesh were triple mutant pvdhfr combined with single mutant pvdhps allele and triple mutant pvdhfr combined with double wild-type pvdhps alleles, respectively. Those collected from Nepal were quadruple mutant pvdhfr combined with double wild-type pvdhps alleles. New alternative antifolate drugs which are effective against sulfadoxine-pyrimethamine (SP)-resistant P. vivax are required.     

Plasmodium vivax Drug Resistance Genes; Pvmdr1 and Pvcrt-o Polymorphisms in Relation to Chloroquine Sensitivity from a Malaria Endemic Area of Thailand

The aim of the study was to explore the possible molecular markers of chloroquine resistance in Plasmodium vivax isolates in Thailand. A total of 30 P. vivax isolates were collected from a malaria endemic area along the Thai-Myanmar border in Mae Sot district of Thailand. Dried blood spot samples were collected for analysis of Pvmdr1 and Pvcrt-o polymorphisms. Blood samples (100 μl) were collected by finger-prick for in vitro chloroquine susceptibility testing by schizont maturation inhibition assay. Based on the cut-off IC₅₀ of 100 nM, 19 (63.3%) isolates were classified as chloroquine resistant P. vivax isolates. Seven non-synonymous mutations and 2 synonymous were identified in Pvmdr1 gene. Y976F and F1076L mutations were detected in 7 (23.3%) and 16 isolates (53.3%), respectively. Analysis of Pvcrt-o gene revealed that all isolates were wild-type. Our results suggest that chloroquine resistance gene is now spreading in this area. Monitoring of chloroquine resistant molecular markers provide a useful tool for future control of P. vivax malaria.     

Sequence signatures involved in targeting the male-specific lethal complex to X-chromosomal genes in Drosophila melanogaster

Background

Drug resistance in the malaria parasite Plasmodium falciparum severely compromises the treatment and control of malaria. A knowledge of the critical mutations conferring resistance to particular drugs is important in understanding modes of drug action and mechanisms of resistances. They are required to design better therapies and limit drug resistance. A mutation in the gene (pfcrt) encoding a membrane transporter has been identified as a principal determinant of chloroquine resistance in P. falciparum, but we lack a full account of higher level chloroquine resistance. Furthermore, the determinants of resistance in the other major human malaria parasite, P. vivax, are not known. To address these questions, we investigated the genetic basis of chloroquine resistance in an isogenic lineage of rodent malaria parasite P. chabaudi in which high level resistance to chloroquine has been progressively selected under laboratory conditions.

Results

Loci containing the critical genes were mapped by Linkage Group Selection, using a genetic cross between the high-level chloroquine-resistant mutant and a genetically distinct sensitive strain. A novel high-resolution quantitative whole-genome re-sequencing approach was used to reveal three regions of selection on chr11, chr03 and chr02 that appear progressively at increasing drug doses on three chromosomes. Whole-genome sequencing of the chloroquine-resistant parent identified just four point mutations in different genes on these chromosomes. Three mutations are located at the foci of the selection valleys and are therefore predicted to confer different levels of chloroquine resistance. The critical mutation conferring the first level of chloroquine resistance is found in aat1, a putative aminoacid transporter.

Conclusions

Quantitative trait loci conferring selectable phenotypes, such as drug resistance, can be mapped directly using progressive genome-wide linkage group selection. Quantitative genome-wide short-read genome resequencing can be used to reveal these signatures of drug selection at high resolution. The identities of three genes (and mutations within them) conferring different levels of chloroquine resistance generate insights regarding the genetic architecture and mechanisms of resistance to chloroquine and other drugs. Importantly, their orthologues may now be evaluated for critical or accessory roles in chloroquine resistance in human malarias P. vivax and P. falciparum.     

Molecular Evidence of Plasmodium vivax Mono and Mixed Malaria Parasite Infections in Duffy-Negative Native Cameroonians

The malaria parasite Plasmodium vivax is known to be majorly endemic to Asian and Latin American countries with no or very few reports of Africans infected with this parasite. Since the human Duffy antigens act as receptors for P. vivax to invade human RBCs and Africans are generally Duffy-negative, non-endemicity of P. vivax in Africa has been attributed to this fact. However, recent reports describing P. vivax infections in Duffy-negative Africans from West and Central parts of Africa have been surfaced including a recent report on P. vivax infection in native Cameroonians. In order to know if Cameroonians living in the southern regions are also susceptible to P. vivax infection, we collected finger-prick blood samples from 485 malarial symptomatic patients in five locations and followed PCR diagnostic assays with DNA sequencing of the 18S ribosomal RNA gene. Out of the 201 malaria positive cases detected, 193 were pure P. falciparum, six pure P. vivax and two mixed parasite infections (P. falciparum + P. vivax). The eight P. vivax infected samples (six single + two mixed) were further subjected to DNA sequencing of the P. vivax multidrug resistance 1 (pvmdr1) and the P.vivax circumsporozoite (pvcsp) genes. Alignment of the eight Cameroonian pvmdr1 sequences with the reference sequence showed high sequence similarities, reconfirming P. vivax infection in all the eight patients. DNA sequencing of the pvcsp gene indicated all the eight P. vivax to be of VK247 type. Interestingly, DNA sequencing of a part of the human Duffy gene covering the promoter region in the eight P. vivax-infected Cameroonians to identify the T-33C mutation revealed all these patients as Duffy-negative. The results provide evidence of single P. vivax as well as mixed malaria parasite infection in native Cameroonians and add knowledge to the growing evidences of P. vivax infection in Duffy-negative Africans.     

Plasmodium vivax: Prevalence of mutations associated with sulfadoxine–pyrimethamine resistance in Plasmodium vivax clinical isolates from Pakistan

The main aim of the present study was to investigate the frequency of SNPs-haplotypes of dhfr and dhps genes associated to sulfadoxine–pyrimethamine (SP) resistance in Plasmodium vivax clinical isolates circulating in a malaria endemic area, Pakistan. All 164 collected isolates were analyzed for SNPs-haplotypes at positions 13, 33, 57, 58, 61, 117 and 173 of pvdhfr and 383 and 553 of pvdhps genes using PCR–RFLP methods. All examined isolates were found to carry wild-type amino acids at positions 13, 33, 57, 61 and 173, while 58R and 117N mutations were detected among 15.2% and 53.6% of isolates, respectively. Based on the size polymorphism of pvdhfr genes at repeat region, type B (79.3%) was the most prevalent variant. The combination of pvdhfr and pvdhps haplotypes demonstrated nine distinct haplotypes. The three most prevalent haplotypes were I₁₃P₃₃F₅₇S₅₈T₆₁S₁₁₇I₁₇₃/A₃₈₃A₅₅₃ (43.9%), I₁₃P₃₃F₅₇S₅₈T₆₁N₁₁₇I₁₇₃/A₃₈₃A₅₅₃ (33.6%) and I₁₃P₃₃F₅₇R₅₈T₆₁N₁₁₇I₁₇₃/A₃₈₃A₅₅₃ (12.2%). The presence of mutant haplotypes is worrying and indicates the emergence of drug tolerant/resistant P. vivax isolates in Pakistan in near future.     

Prevalence and Clinical Manifestations of Malaria in Aligarh,India

Malaria is one of the most widespread infectious diseases of tropical countries with an estimated 207 million cases globally. In India, there are endemic pockets of this disease, including Aligarh. Hundreds of Plasmodium falciparum and P. vivax cases with severe pathological conditions are recorded every year in this district. The aim of this study is to find out changes in liver enzymes and kidney markers. Specific diagnosis for P. falciparum and P. vivax was made by microscopic examination of Giemsa stained slides. Clinical symptoms were observed in both of these infections. Liver enzymes, such as AST, ALT, and ALP, and kidney function markers, such as creatinine and urea, were estimated by standard biochemical techniques. In Aligarh district, P. vivax, P. falciparum, and mixed infections were 64%, 34%, and 2%, respectively. In case of P. falciparum infection, the incidences of anemia, splenomegaly, renal failure, jaundice, and neurological sequelae were higher compared to those in P. vivax infection. Recrudescence and relapse rates were 18% and 20% in P. falciparum and P. vivax infections, respectively. Liver dysfunctions and renal failures were more common in P. falciparum patients, particularly in elderly patients. Artesunate derivatives must, therefore, be introduced for the treatment of P. falciparum as they resist to chloroquine as well as sulfadoxine-pyrimethamine combinations.     

Genetic Polymorphisms in VIR Genes among Indian Plasmodium vivax Populations

The vir genes are antigenic genes and are considered to be possible vaccine targets. Since India is highly endemic to Plasmodium vivax, we sequenced 5 different vir genes and investigated DNA sequence variations in 93 single-clonal P. vivax isolates. High variability was observed in all the 5 vir genes; the vir 1/9 gene was highly diverged across Indian populations. The patterns of genetic diversity do not follow geographical locations, as geographically distant populations were found to be genetically similar. The results in general present complex genetic diversity patterns in India, requiring further in-depth population genetic and functional studies.     

A High Resolution Case Study of a Patient with Recurrent Plasmodium vivax Infections Shows That Relapses Were Caused by Meiotic Siblings

Plasmodium vivax infects a hundred million people annually and endangers 40% of the world''s population. Unlike Plasmodium falciparum, P. vivax parasites can persist as a dormant stage in the liver, known as the hypnozoite, and these dormant forms can cause malaria relapses months or years after the initial mosquito bite. Here we analyze whole genome sequencing data from parasites in the blood of a patient who experienced consecutive P. vivax relapses over 33 months in a non-endemic country. By analyzing patterns of identity, read coverage, and the presence or absence of minor alleles in the initial polyclonal and subsequent monoclonal infections, we show that the parasites in the three infections are likely meiotic siblings. We infer that these siblings are descended from a single tetrad-like form that developed in the infecting mosquito midgut shortly after fertilization. In this natural cross we find the recombination rate for P. vivax to be 10 kb per centimorgan and we further observe areas of disequilibrium surrounding major drug resistance genes. Our data provide new strategies for studying multiclonal infections, which are common in all types of infectious diseases, and for distinguishing P. vivax relapses from reinfections in malaria endemic regions. This work provides a theoretical foundation for studies that aim to determine if new or existing drugs can provide a radical cure of P. vivax malaria.     

FastD: Fast detection of insecticide target‐site mutations and overexpressed detoxification genes in insect populations from RNA‐Seq data

Target‐site mutations and detoxification gene overexpression are two major mechanisms conferring insecticide resistance. Molecular assays applied to detect these resistance genetic markers are time‐consuming and with high false‐positive rates. RNA‐Seq data contains information on the variations within expressed genomic regions and expression of detoxification genes. However, there is no corresponding method to detect resistance markers at present. Here, we collected 66 reported resistance mutations of four insecticide targets (AChE, VGSC, RyR, and nAChR) from 82 insect species. Next, we obtained 403 sequences of the four target genes and 12,665 sequences of three kinds of detoxification genes including P450s, GSTs, and CCEs. Then, we developed a Perl program, FastD, to detect target‐site mutations and overexpressed detoxification genes from RNA‐Seq data and constructed a web server for FastD (http://www.insect-genome.com/fastd). The estimation of FastD on simulated RNA‐Seq data showed high sensitivity and specificity. We applied FastD to detect resistant markers in 15 populations of six insects, Plutella xylostella, Aphis gossypii, Anopheles arabiensis, Musca domestica, Leptinotarsa decemlineata and Apis mellifera. Results showed that 11 RyR mutations in P. xylostella, one nAChR mutation in A. gossypii, one VGSC mutation in A. arabiensis and five VGSC mutations in M. domestica were found to be with frequency difference >40% between resistant and susceptible populations including previously confirmed mutations G4946E in RyR, R81T in nAChR and L1014F in VGSC. And 49 detoxification genes were found to be overexpressed in resistant populations compared with susceptible populations including previously confirmed detoxification genes CYP6BG1, CYP6CY22, CYP6CY13, CYP6P3, CYP6M2, CYP6P4 and CYP4G16. The candidate target‐site mutations and detoxification genes were worth further validation. Resistance estimates according to confirmed markers were consistent with population phenotypes, confirming the reliability of this program in predicting population resistance at omics‐level.     

Primaquine Administration after Falciparum Malaria Treatment in Malaria Hypoendemic Areas with High Incidence of Falciparum and Vivax Mixed Infection: Pros and Cons

Mixed infections of Plasmodium falciparum and Plasmodium vivax is high (~30%) in some malaria hypoendemic areas where the patients present with P. falciparum malaria diagnosed by microscopy. Conventional treatment of P. falciparum with concurrent chloroquine and 14 days of primaquine for all falciparum malaria patients may be useful in areas where mixed falciparum and vivax infections are high and common and also with mild or moderate G6PD deficiency in the population even with or without subpatent vivax mixed infection. It will be possibly cost-effective to reduce subsequent vivax illness if the patients have mixed vivax infection. Further study to prove this hypothesis may be warranted.     

A Case of Plasmodium ovale wallikeri Infection in a Chinese Worker Returning from West Africa

In contrast to the gradual reduction in the number of locally transmitted malaria cases in China, the number of imported malaria cases has been increasing since 2008. Here, we report a case of a 39-year-old Chinese man who acquired Plasmodium ovale wallikeri infection while staying in Ghana, West Africa for 6 months in 2012. Microscopic examinations of Giemsa-stained thin and thick blood smears indicated Plasmodium vivax infection. However, the results of rapid diagnostic tests, which were conducted 3 times, were not in agreement with P. vivax. To further check the diagnosis, standard PCR analysis of the small-subunit rRNA gene was conducted, based on which a phylogeny tree was constructed. The results of gene sequencing indicated that this malaria is a variant of P. ovale (P. ovale wallikeri). The infection in this patient was not a new infection, but a relapse of the infection from the one that he had contracted in West Africa.     

Geographical spread and structural basis of sulfadoxine-pyrimethamine drug-resistant malaria parasites

The global spread of sulfadoxine (Sdx, S) and pyrimethamine (Pyr, P) resistance is attributed to increasing number of mutations in DHPS and DHFR enzymes encoded by malaria parasites. The association between drug resistance mutations and SP efficacy is complex. Here we provide an overview of the geographical spread of SP resistance mutations in Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) encoded dhps and dhfr genes. In addition, we have collated the mutation data and mapped it on to the three-dimensional structures of DHPS and DHFR which have become available. Data from genomic databases and 286 studies were collated to provide a comprehensive landscape of mutational data from 2005 to 2019. Our analyses show that the Pyr-resistant double mutations are widespread in Pf/PvDHFR (P. falciparum ～61% in Asia and the Middle East, and in the Indian sub-continent; in P. vivax ～33% globally) with triple mutations prevailing in Africa (～66%) and South America (～33%). For PfDHPS, triple mutations dominate South America (～44%), Asia and the Middle East (～34%) and the Indian sub-continent (～27%), while single mutations are widespread in Africa (～45%). Contrary to the status for P. falciparum, Sdx-resistant single point mutations in PvDHPS dominate globally. Alarmingly, highly resistant quintuple and sextuple mutations are rising in Africa (PfDHFR-DHPS) and Asia (Pf/PvDHFR-DHPS). Structural analyses of DHFR and DHPS proteins in complexes with substrates/drugs have revealed that resistance mutations map proximal to Sdx and Pyr binding sites. Thus new studies can focus on discovery of novel inhibitors that target the non-substrate binding grooves in these two validated malaria parasite drug targets.     

Evaluation of rapid diagnostics for Plasmodium falciparum and P. vivax in Mae Sot Malaria endemic area, Thailand

Prompt and accurate diagnosis of malaria is the key to prevent disease morbidity and mortality. This study was carried out to evaluate diagnostic performance of 3 commercial rapid detection tests (RDTs), i.e., Malaria Antigen Pf/Pan™, Malaria Ag-Pf™, and Malaria Ag-Pv™ tests, in comparison with the microscopic and PCR methods. A total of 460 blood samples microscopically positive for Plasmodium falciparum (211 samples), P. vivax (218), mixed with P. falciparum and P. vivax (30), or P. ovale (1), and 124 samples of healthy subjects or patients with other fever-related infections, were collected. The sensitivities of Malaria Ag-Pf™ and Malaria Antigen Pf/Pan™ compared with the microscopic method for P. falciparum or P. vivax detection were 97.6% and 99.0%, or 98.6% and 99.0%, respectively. The specificities of Malaria Ag-Pf™, Malaria Ag-Pv™, and Malaria Antigen Pf/Pan™ were 93.3%, 98.8%, and 94.4%, respectively. The sensitivities of Malaria Ag-Pf™, Malaria Antigen Pf/Pan™, and microscopic method, when PCR was used as a reference method for P. falciparum or P. vivax detection were 91.8%, 100%, and 96.7%, or 91.9%, 92.6%, and 97.3%, respectively. The specificities of Malaria Ag-Pf™, Malaria Ag-Pv™, Malaria Antigen Pf/Pan™, and microscopic method were 66.2%, 92.7%, 73.9%, and 78.2%, respectively. Results indicated that the diagnostic performances of all the commercial RDTs are satisfactory for application to malaria diagnosis.     

Coexistence of Malaria and Thalassemia in Malaria Endemic Areas of Thailand

Hemoglobinopathy and malaria are commonly found worldwide particularly in malaria endemic areas. Thalassemia, the alteration of globin chain synthesis, has been reported to confer resistance against malaria. The prevalence of thalassemia was investigated in 101 malaria patients with Plasmodium falciparum and Plasmodium vivax along the Thai-Myanmar border to examine protective effect of thalassemia against severe malaria. Hemoglobin typing was performed using low pressure liquid chromatography (LPLC) and α-thalassemia was confirmed by multiplex PCR. Five types of thalassemia were observed in malaria patients. The 2 major types of thalassemia were Hb E (18.8%) and α-thalassemia-2 (11.9%). There was no association between thalassemia hemoglobinopathy and malaria parasitemia, an indicator of malaria disease severity. Thalassemia had no significant association with P. vivax infection, but the parasitemia in patients with coexistence of P. vivax and thalassemia was about 2-3 times lower than those with coexistence of P. falciparum and thalassemia and malaria without thalassemia. Furthermore, the parasitemia of P. vivax in patients with coexistence of Hb E showed lower value than coexistence with other types of thalassemia and malaria without coexistence. Parasitemia, hemoglobin, and hematocrit values in patients with coexistence of thalassemia other than Hb E were significantly lower than those without coexistence of thalassemia. Furthermore, parasitemia with coexistence of Hb E were 2 times lower than those with coexistence of thalassemia other than Hb E. In conclusion, the results may, at least in part, support the protective effect of thalassemia on the development of hyperparasitemia and severe anemia in malaria patients.     

Resistome analysis of Mycobacterium tuberculosis: Identification of aminoglycoside 2'-Nacetyltransferase (AAC) as co-target for drug desigining

The emergence of multidrug resistant tuberculosis (MDRTB) highlights the urgent need to understand the mechanisms of resistance to the drugs and to develop a new arena of therapeutics to treat the disease. Ethambutol, isonazid, pyrazinamide, rifampicin are first line of drugs against TB, whereas aminoglycoside, polypeptides, fluoroquinolone, ethionamide are important second line of bactericidal drugs used to treat MDRTB, and resistance to one or both of these drugs are defining characteristic of extensively drug resistant TB. We retrieved 1,221 resistant genes from Antibiotic Resistance Gene Database (ARDB), which are responsible for resistance against first and second line antibiotics used in treatment of Mycobacterium tuberculosis infection. From network analysis of these resistance genes, 53 genes were found to be common. Phylogenetic analysis shows that more than 60% of these genes code for acetyltransferase. Acetyltransferases detoxify antibiotics by acetylation, this mechanism plays central role in antibiotic resistance. Seven acetyltransferase (AT-1 to AT-7) were selected from phylogenetic analysis. Structural alignment shows that these acetyltransferases share common ancestral core, which can be used as a template for structure based drug designing. From STRING analysis it is found that acetyltransferase interact with 10 different proteins and it shows that, all these interaction were specific to M. tuberculosis. These results have important implications in designing new therapeutic strategies with acetyltransferase as lead co-target to combat against MDR as well as Extreme drug resistant (XDR) tuberculosis.

Abbreviations

AA - amino acid, AT - Acetyltransferase, AAC - Aminoglycoside 2''-N-acetyltransferase, XDR - Extreme drug-resistant, MDR - Multidrug-resistant, Mtb - Mycobacterium tuberculosis, TB - Tuberculosis.     

Mutational and acquired carbapenem resistance mechanisms in multidrug resistant Pseudomonas aeruginosa clinical isolates from Recife, Brazil

An investigation was carried out into the genetic mechanisms responsible for multidrug resistance in nine carbapenem-resistant Pseudomonas aeruginosaisolates from different hospitals in Recife, Brazil. Susceptibility to antimicrobial agents was determined by broth microdilution. Polymerase chain reaction (PCR) was employed to detect the presence of genes encoding β-lactamases, aminoglycoside-modifying enzymes (AMEs), 16S rRNA methylases, integron-related genes and OprD. Expression of genes coding for efflux pumps and AmpC cephalosporinase were assessed by quantitative PCR. The outer membrane proteins were separated by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. The blaSPM-1, blaKPC-2 and blaGES-1 genes were detected in P. aeruginosaisolates in addition to different AME genes. The loss of OprD in nine isolates was mainly due to frameshift mutations, premature stop codons and point mutations. An association of loss of OprD with the overexpression of MexAB-OprM and MexXY-OprM was observed in most isolates. Hyper-production of AmpC was also observed in three isolates. Clonal relationship of the isolates was determined by repetitive element palindromic-PCR and multilocus sequence typing. Our results show that the loss of OprD along with overexpression of efflux pumps and β-lactamase production were responsible for the multidrug resistance in the isolates analysed.     

Positional cloning of rp2 QTL associates the P450 genes CYP6Z1, CYP6Z3 and CYP6M7 with pyrethroid resistance in the malaria vector Anopheles funestus

Pyrethroid resistance in Anopheles funestus is threatening malaria control in Africa. Elucidation of underlying resistance mechanisms is crucial to improve the success of future control programs. A positional cloning approach was used to identify genes conferring resistance in the uncharacterised rp2 quantitative trait locus (QTL) previously detected in this vector using F6 advanced intercross lines (AIL). A 113 kb BAC clone spanning rp2 was identified and sequenced revealing a cluster of 15 P450 genes and one salivary protein gene (SG7-2). Contrary to A. gambiae, AfCYP6M1 is triplicated in A. funestus, while AgCYP6Z2 orthologue is absent. Five hundred and sixty-five new single nucleotide polymorphisms (SNPs) were identified for genetic mapping from rp2 P450s and other genes revealing high genetic polymorphisms with one SNP every 36 bp. A significant genotype/phenotype association was detected for rp2 P450s but not for a cluster of cuticular protein genes previously associated with resistance in A. gambiae. QTL mapping using F6 AIL confirms the rp2 QTL with an increase logarithm of odds score of 5. Multiplex gene expression profiling of 15 P450s and other genes around rp2 followed by individual validation using qRT–PCR indicated a significant overexpression in the resistant FUMOZ-R strain of the P450s AfCYP6Z1, AfCYP6Z3, AfCYP6M7 and the glutathione-s-transferase GSTe2 with respective fold change of 11.2, 6.3, 5.5 and 2.8. Polymorphisms analysis of AfCYP6Z1 and AfCYP6Z3 identified amino acid changes potentially associated with resistance further indicating that these genes are controlling the pyrethroid resistance explained by the rp2 QTL. The characterisation of this rp2 QTL significantly improves our understanding of resistance mechanisms in A. funestus.     

Identification of Immunodominant B-cell Epitope Regions of Reticulocyte Binding Proteins in Plasmodium vivax by Protein Microarray Based Immunoscreening

Plasmodium falciparum can invade all stages of red blood cells, while Plasmodium vivax can invade only reticulocytes. Although many P. vivax proteins have been discovered, their functions are largely unknown. Among them, P. vivax reticulocyte binding proteins (PvRBP1 and PvRBP2) recognize and bind to reticulocytes. Both proteins possess a C-terminal hydrophobic transmembrane domain, which drives adhesion to reticulocytes. PvRBP1 and PvRBP2 are large (> 326 kDa), which hinders identification of the functional domains. In this study, the complete genome information of the P. vivax RBP family was thoroughly analyzed using a prediction server with bioinformatics data to predict B-cell epitope domains. Eleven pvrbp family genes that included 2 pseudogenes and 9 full or partial length genes were selected and used to express recombinant proteins in a wheat germ cell-free system. The expressed proteins were used to evaluate the humoral immune response with vivax malaria patients and healthy individual serum samples by protein microarray. The recombinant fragments of 9 PvRBP proteins were successfully expressed; the soluble proteins ranged in molecular weight from 16 to 34 kDa. Evaluation of the humoral immune response to each recombinant PvRBP protein indicated a high antigenicity, with 38-88% sensitivity and 100% specificity. Of them, N-terminal parts of PvRBP2c (PVX_090325-1) and PvRBP2 like partial A (PVX_090330-1) elicited high antigenicity. In addition, the PvRBP2-like homologue B (PVX_116930) fragment was newly identified as high antigenicity and may be exploited as a potential antigenic candidate among the PvRBP family. The functional activity of the PvRBP family on merozoite invasion remains unknown.