- Open Access
Ribosomal and non-ribosomal PCR targets for the detection of low-density and mixed malaria infections
Malaria Journal volume 18, Article number: 154 (2019)
The unexpected high proportion of submicroscopic malaria infections in areas with low transmission intensity challenges the control and elimination of malaria in the Americas. The current PCR-based assays present limitations as most protocols still rely on amplification of few-copies target gene. Here, the hypothesis was that amplification of different plasmodial targets—ribosomal (18S rRNA) and non-ribosomal multi-copy sequences (Pvr47 for Plasmodium vivax and Pfr364 for Plasmodium falciparum)—could increase the chances of detecting submicroscopic malaria infection.
A non-ribosomal real-time PCR assay targeting Pvr47/Pfr364 (NR-qPCR) was established and compared with three additional PCR protocols, two of them based on 18S rRNA gene amplification (Nested-PCR and R-qPCR) and one based on Pvr47/Pfr364 targets (NR-cPCR). The limit of detection of each PCR protocol, at single and artificial mixed P. vivax/P. falciparum infections, was determined by end-point titration curves. Field samples from clinical (n = 110) and subclinical (n = 324) malaria infections were used to evaluate the impact of using multiple molecular targets to detect malaria infections.
The results demonstrated that an association of ribosomal and non-ribosomal targets did not increase sensitivity to detect submicroscopic malaria infections. Despite of that, artificial mixed-malaria infections demonstrated that the NR-qPCR was the most sensitive protocol to detect low-levels of P. vivax/P. falciparum co-infections. Field studies confirmed that submicroscopic malaria represented a large proportion (up to 77%) of infections among asymptomatic Amazonian residents, with a high proportion of infections (~ 20%) identified only by the NR-qPCR.
This study presents a new species-specific non-ribosomal PCR assay with potential to identify low-density P. vivax and P. falciparum infections. As the majority of subclinical infections was caused by P. vivax, the commonest form of malaria in the Amazon area, future studies should investigate the potential of Pvr47/Pfr364 to detect mixed-malaria infections in the field.
Malaria is a worldwide public health concern that is present in roughly 90 countries, mainly in tropical and subtropical regions . While Plasmodium vivax is the most widely distributed parasite causing malaria, Plasmodium falciparum accounts for the most severe forms of the disease . Although malaria incidence rate is estimated to have decreased by 18% globally between 2010 and 2016, a recent increase in case incidence occurred in the Americas, particularly in the Amazon rainforest [1, 3].
In order to progress towards malaria control and elimination, it is critical to understand the sources of transmission (the infectious reservoir) and those at risk of infection at the population level . In this context, the molecular detection of Plasmodium infections in endemic areas have confirmed previous finding of high frequencies of malaria infections at densities below the limit of conventional microscopic diagnostics [5,6,7,8,9]. It is particularly relevant as data from systematic reviews have shown that across different geographic areas low-density infections may represent the majority of malaria infections [10, 11]. Accordingly, a substantial proportion of asymptomatic and submicroscopic malarial infections has been described in peri-urban areas of the Brazilian Amazon . Many are the implications of these findings as submicroscopic malaria carriers may be able to transmit the Plasmodium parasites, acting as reservoirs for malaria [13, 14]. Beyond the practical value of using molecular tools to identify submicroscopic carriers and mixed-species infections, malaria infections at any density may have significant health and socioeconomic consequences .
Historically, the small subunit 18S of the ribosomal RNA gene (18S rRNA) has been the most common target used for molecular diagnosis of malaria [16,17,18,19,20]. As this gene is present in few copies (5 to 8) in the genome of Plasmodium parasites , low sensitivity and reproducibility of standard PCR protocols based on 18S rRNA gene amplification have been described [22, 23]. In the last decade, the genomic data mining of Plasmodium parasites has allowed the discovery of new species-specific multi-copy targets which show potential for molecular diagnosis of P. vivax and P. falciparum malaria [24,25,26]. Among the promising targets include the non-coding subtelomeric repeat sequences Pvr47 and Pfr364 that are present in 14 and 41 copies in the genomes of P. vivax and P. falciparum, respectively . While there is evidence for their location and distribution, the biological functions of Pvr47 and Pfr364 remains to be established. By using a single-step PCR assay to amplify Pvr47/Pfr364 targets, it was possible to demonstrate the relatively higher sensitivity of these targets as compared to the amplification of 18S rRNA gene by the conventional nested-PCR assay .
Since most malaria PCR-based protocols still relies on amplification of 18S rRNA gene, which has low sensitivity to detect low-density infections, we evaluated here how useful Pvr47/Pfr364 targets are to detect single and mixed P. vivax and P. falciparum infections in clinical and subclinical malaria. As the original PCR protocol to amplify Pvr47/Pfr364 involved DNA visualization on gel electrophoresis , here a new qPCR protocol targeting these high-copy non-ribosomal sequences was established. The experimental approach evaluated whether amplification of different plasmodial targets (Pvr47/Pfr364 and 18S rRNA gene) could increase the chances of detecting submicroscopic malaria infections. For that, field samples (clinical and subclinical malaria) were amplified by four different PCR assays, two of them targeting Pvr47/Pfr364 sequences  and two targeting the 18S rRNA gene [16, 17].
Study population and participants
Ethical and methodological aspects of this study were approved by the Ethical Committee of Research on Human Beings from the René Rachou Institute/Oswaldo Cruz Foundation (protocols No 24/2008, and No 1.821.955/2016), according to the Brazilian National Council of Health (Resolutions 196/96 and 466/12). All participants were informed about the objectives and procedures of the study, with voluntary participation through written informed consent.
This group included individuals with clinical suspicion of malaria who sought care at Brazilian malaria reference healthcare facilities located in both endemic (States of Rondônia and Mato Grosso) and non-endemic areas (Minas Gerais). After short-trips to malaria transmission areas, these individuals presented symptoms suggestive of uncomplicated malaria infection, such as fever, myalgia, chills, and headaches. Non-inclusion criteria included: (i) refusal or inability to sign the informed consent; (ii) age below 5 years-old; (iii) pregnant women; and (iv) any other co-morbidity that could be traced. One-hundred-and-ten patients were enrolled in the study, which gives 90% statistical power at 5% significance level assuming 30% of prevalence by light microscopy and an estimative of 50% increase by molecular assays; the majority of study population were adults, with a median age of 40 years (IQR 30.75–48.25), and a proportion female:male of 1:5. For each individual, blood sample was collected at a single time point by venipuncture or finger prick. The period of blood collection varied from 2008 to 2017, and since there, all samples from participants have been maintained in the biorepository of malaria research group at FIOCRUZ-MINAS (Belo Horizonte, MG), Brazilian Ministry of Health, according to the local legislation.
This group included malaria-exposed individuals who had participated of cross-sectional surveys carried-out between 2008 and 2015 in a rural community of the Brazilian Amazon rainforest, Rio Pardo (1°46′S—1°54′S, 60°22′W—60°10′W), Presidente Figueiredo municipality, State of Amazonas. The study site and malaria transmission patterns have been described in details elsewhere [27,28,29]. In this area, malaria transmission is considered hypo to mesoendemic, and the majority of residents were natives from the Amazon region . For the current study, the non-inclusion criteria were: (i) refusal or inability to sign the informed consent; (ii) any signs and/or symptoms that could be related to symptomatic malaria such as fever, myalgia, chills, and headaches; the absence of symptoms was self-reported and obtained during a personal interview conducted through structured questionnaire, as previously described ; (iii) age below 5 years-old, as subclinical malaria infection is not prevalent in Amazon children ; (iv) pregnant women; and (v) any other morbidity that could be traced. This group was composed by 324 participants; considering 7% of prevalence by light microscopy  and assuming that in an area of greater chance of P. vivax infection molecular assays is expected to increase malaria prevalence by a factor of 2 , sample size will give more than 90% of statistical power at 5% significance level. The median age of studied individuals were 40 years (IQR 24–56), with female:male ratio of 1:1.5, and an average of 35 years (IQR 21–52) living in the endemic area. For each individual, blood sample was collected at a single time-point by venipuncture or finger prick. As in the study area the number of malaria cases fluctuated during the last 8 years, reflecting periods of high and low malaria transmission , samples covered the temporal-variation in the profile of malaria transmission; thus, the blood samples were collected from November 2008 to June 2015, and all biological specimens have been maintained in the biorepository of malaria research group at FIOCRUZ-MINAS (Belo Horizonte, MG), Brazilian Ministry of Health, according to the local legislation.
Conventional light microscopy
At the time of blood collection, all individuals (clinical and subclinical groups) were submitted to a finger-prick for malaria diagnosis by light microscopy. The Giemsa-stained thick blood smears were prepared and examined by experienced local microscopists, according to the malaria diagnosis guidelines of the Brazilian Ministry of Health . Parasite density was estimated as the number of parasites per microlitre of blood (parasites/µL), and all microscopically positive cases were treated immediately in local health services, following the national malaria treatment protocols recommended by the Brazilian Ministry of Health .
Extraction of genomic DNA
The experimental approach to the molecular detection of clinical and subclinical malaria infections was included in Fig. 1. Genomic DNA (gDNA) was extracted from either whole blood samples collected in EDTA, or from dried blood spots on filter paper using the Gentra Puregene Blood Kit (Qiagen) and the QIAamp DNA Mini Kit (Qiagen), respectively, according to manufacturer’s instructions. As an internal control of the DNA extractions, 10% of the samples extracted were randomly submitted to a PCR assay for the amplification of a human gene (ABO blood group), according to the protocol previously described . All samples tested amplified the target gene. The extracted DNA was stored at − 20 °C until use.
Panel of Plasmodium reference samples
The following parasites were used as reference in molecular assays: (i) P. falciparum (3D7 strain) from in vitro continuous blood-stage cultures maintained in the routine of the laboratory, according to protocol previously described ; (ii) P. vivax, P. falciparum and Plasmodium malariae from peripheral blood of malaria infected individuals whose species-specific diagnosis was confirmed by PCR protocols (Malaria Biorepository, FIOCRUZ-Minas, MG, Brazil); (iii) Plasmodium brasilianum/P. malariae (Peruvian III strain, MR4-349) kindly provided by the Malaria Research and Reference Reagent Resource Center, MR4 (Biodefense and Emerging Infections Research Resources Repository, BEI Resources, NIAID/NIH, ATCC, USA).
End-point Plasmodium DNA titrations of single and mixed Plasmodium vivax and Plasmodium falciparum infections
DNA samples from high-density P. vivax (12,900 parasites/µL) and P. falciparum (13,400 parasites/µL) infections were used to determine the end-point detection for ribosomal and non-ribosomal PCR assays. For each plasmodial DNA, three-fold dilutions were prepared range from 2000 to 0.3 parasites/µL. Similarly, artificial mixed malaria infections were produced by mixed DNA samples from P. vivax and P. falciparum in different proportions (1:1 until 1:2000; threefold dilutions). Each sample was assayed in triplicate and submitted to all PCR protocols (Fig. 1).
Primers design and real-time PCR conditions to amplify the non-ribosomal Pvr47/Pfr364 targets (NR-qPCR)
Primers and MGB (minor groove binder) probes were designed for Pvr47 and Pfr364 targets using Primer-Blast (NCBI) and OLIGO (version 4.0, 1999) softwares, considering conserved regions in each species alignments. The alignments were based on the sequences previously described , composed of 14 copies of Pvr47 (P. vivax Sal-I strain) and 41 copies of Pfr364 (subfamilies 1 and 2 of P. falciparum 3D7 strain). The set of oligonucleotides obtained for Pvr47 was 5′TCCGCAGCTCACAAATGTTC3′ (forward), 5′ACATGGGGATTCTAAGCCAATTTA3′ (reverse), and 5′HEX-TCCGCGAGGGCTGCAA3′ (probe), which binds respectively to positions 142–161, 245–222 and 189–204 of P. vivax Sal-I (GenBank accession number AAKM01000578), resulting in a 104 bp amplicon. For Pfr364, the set used was 5′ACTCGCAATAACGCTGCAT3′ (forward), 5′TTCCCTGCCCAAAAACGG3′ (reverse), and 5′FAM-GGTGCCGGGGGTTTCTACGC3′ (probe), which binds respectively to positions 313–331, 400–383 and 335–354 of P. falciparum 3D7 (Pf3D7_04_12625_14152, in ), resulting in an 88 bp amplicon. All NR-qPCR reactions were performed in 10 μL volumes containing 2 μL of DNA (~ 3 to 6 μL of whole blood) and 5 μL of TaqMan Universal PCR Master Mix (Thermo Fisher Scientific). For Pvr47 amplification was used 50 nM of forward primer, 900 nM of reverse primer, and 250 nM of probe; for Pfr364, 900 nM of forward primer, 300 nM of reverse primer, and 150 nM of probe were used. The PCR assays were performed using the automatic thermocycler ViiA7 Real-Time PCR System (Thermo Fisher Scientific) and the following cycling parameter: a pre-incubation and initial denaturation, respectively, at 50 °C for 2 min and 95 °C for 10 min, followed by 40 cycles of denaturation at 95 °C for 15 s, primers annealing at 52 °C for 1 min, and extension at 60 °C for 1 min. The fluorescence acquisition was performed at the end of each extension step. Analytical sensitivity and specificity of Pvr47 and Pfr364 NR-qPCR assays were included as an Additional file 1. Specificity was investigated using template DNA from either malaria unexposed individuals (n = 30) or other Plasmodium species (Additional file 2). For the limit of detection, standard curves were prepared from serial dilution of plasmid DNA carrying the Pvr47 or Pfr364 target (Additional file 3). The cycle threshold (Ct) values of 37 and 38 (Ct≤ 37 or Ct≤ 38) were used to define positivity to P. vivax and P. falciparum, respectively.
Amplification of Pvr47/Pfr364 targets by conventional PCR assay (NR-cPCR)
The amplification of the targets Pvr47 (P. vivax) and Pfr364 (P. falciparum) were conducted by using primers previously described . Details about primers and cycling conditions were included in Additional file 4.
18S rRNA gene amplification by nested-PCR (Nested-PCR) and real-time PCR assay (R-qPCR)
Nested-PCR assay for amplification of 18S rRNA gene was performed as described , using genus and species-specific primers previously described . The real-time PCR (R-qPCR) method uses a single pair of genus-specific primers for the detection of the 18S rRNA gene, with two internal species-specific hydrolysis probes for P. falciparum and P. vivax . Primers and cycling conditions for 18SrRNA gene amplification (Nested-PCR and R-qPCR) were included as Additional file 4.
Differences in proportions were evaluated using the Chi square (χ2) test or Fisher’s exact test, as appropriate. The Probit Regression analysis was used to define the limit of detection (LOD) of NR-qPCR assay, with different input concentrations of plasmid DNA used to calculate the predicted proportion of positive results (MedCalc Statistical Software program, Ostend, Belgium). Heat map of malaria prevalence was constructed using the heatmap.2 function available in the R package gplots. The assessment of sensitivity and specificity of PCR assays was determined as described , with the reference standard established by combining the detection of any PCR assay excluding the protocol under evaluation. The analyses of sensitivity and specificity were performed in the GraphPad InStat, version 3.0 (GraphPad Software, San Diego, CA, USA). The significance at the 5% level was considered to all analysis.
Limit of detection of Pvr47 and Pfr364 NR-qPCR
To determine the LOD of NR-qPCR protocols, standard curves were constructed for each DNA plasmid carrying the target gene (Pvr47 or Pfr364) through serial dilution ranging from 20,000 to 0.05 copies/μL (Additional file 3). By using probit regression analysis was possible to demonstrate that NR-qPCR developed here presented a 95% probability of detecting levels as low as 0.66 copies/µL for P. vivax (Additional file 5A) and 3.27 copies/µL for P. falciparum (Additional file 5B). No amplification was observed by using template DNA from either malaria unexposed individuals or other Plasmodium species (data not shown, Additional file 2).
Ribosomal and non-ribosomal detection of monoinfections caused by either P. vivax or P. falciparum
The ability of different targets (18S rRNA vs. Pvr47/Pfr364) to detect P. vivax and P. falciparum monoinfections at low parasite densities was investigated through serial dilutions of field samples containing known amounts of either P. vivax or P. falciparum gDNA. In case of P. vivax monoinfection, no significant difference was observed between parasite targets amplified by PCR (Table 1). Despite of that, the NR-qPCR was the only protocol able to consistently detect the lowest levels of parasite densities (all replicates amplified until 1 parasite/µL). Variability between PCRs assays using the same target resulted in a difference in positivity between conventional and real-time PCR assays targeting Pvr47 (67% vs. 96% for NR-cPCR and NR-qPCR, respectively). Considering P. falciparum titration (Table 2), regardless the PCR assay, the amplification of Pfr364 sequence was more precise than 18S rRNA gene (80% vs. 59%, p = 0.0359, Fisher’s exact test). While the amplification of Pfr364 by NR-qPCR was reproducible until 3 parasites/μL, the amplification of 18S rRNA was inconsistent, alternating between positive and negative results at similar levels of parasite density (both nested-PCR and R-qPCR).
Ribosomal and non-ribosomal detection of artificial mixed P. vivax and P. falciparum infections
Next, the ability of ribosomal and non-ribosomal targets to amplify artificial mixed malaria infections was evaluated. By fixing the amount of P. vivax DNA (1433 parasites/μL) and varying the amount of P. falciparum (1489 to 0.7 parasites/μL), it was possible to demonstrate that the protocols based on non-ribosomal targets were much more precise to identify both parasite species, even when P. falciparum was present at very low densities (Table 3). Taken together, non-ribosomal protocols identified 85% (41 out of 48) artificial mixed infections, while ribosomal protocols identified only 31% (15 out of 48). By comparing the variation intra-target, the amplification of 18S rRNA gene by R-qPCR demonstrated a trend to amplify P. vivax in detriment of P. falciparum, even when P. falciparum DNA was present at relatively high concentrations (Table 3). By fixing the concentration of P. falciparum (1489 parasites/μL) and varying P. vivax densities (1433 to 0.7 parasites/μL), a good performance of both targets was observed until 6 parasites/μL (Table 4). Although there was no significant difference between the amplification of ribosomal and non-ribosomal targets (69% vs. 83%), the NR-qPCR protocol developed here seems to present a better performance as compared to 18S rRNA amplification, detecting almost all replicates in all P. vivax/P. falciparum dilution points. In fact, NR-qPCR was the only PCR assay able to consistently detect mixed infection when one of the species was present in a ratio of about 700-fold lower than the other species (1489 vs. 2 parasites/μL to all replicates) (Table 4).
Field evaluation of ribosomal and non-ribosomal PCR targets in clinical and subclinical malaria infections
Initially, the potential use of both ribosomal and non-ribosomal PCR targets for application in malaria field studies involved samples from 110 symptomatic individuals whose light microscopy (LM) confirmed 35 (32%) malaria infections (Fig. 2). As compared with LM, 18S rRNA gene PCR-based assays (nested-PCR and/or R-qPCR) identified a similar proportion of malaria infections (35%, n = 39) (Fig. 2a); of interest, the majority of positive samples were amplified by both 18S rRNA protocols. Basically, the same proportion of positives was obtained with the non-ribosomal protocols (NR-cPCR and/or NR-qPCR). In addition, the use of different PCR targets (ribosomal and non-ribosomal) did not increase malaria positivity (Fig. 2a, the right-side bar chart, in lilac; p = 1.0 for ribosomal vs. both PCR-targets and p = 0.89 for non-ribosomal vs. both targets). To confirm that an association of ribosomal and non-ribosomal targets did not increase sensitivity to detect submicroscopic malaria infections, the values of sensitivity and specificity were determined for molecular PCR-assays. In accordance, the performance of each PCR assay was similar with almost no detection of false positives or negatives (Table 5). Plasmodium species-specific identification (Fig. 2b) showed the same pattern of positivity, with no significant difference obtained between methods (LM vs. PCR based-protocols) or target amplified (ribosomal vs. non-ribosomal). In this symptomatic malaria patients, P. vivax and P. falciparum were found in similar proportions. Despite of that, a couple of mixed-infections (5 out of 6) identified by LM (Fig. 2b, first panel) could not be confirmed by any PCR-based assay.
Screening for subclinical malaria infections among individuals from cross-sectional surveys identified a positivity of 7% (21 out of 324) by LM (Fig. 3A). In this population, molecular diagnosis of malaria increased three to four times the detection of subclinical malaria carries. However, the amplification of either 18S rRNA gene (nested-PCR/R-qPCR) or non-ribosomal targets (NR-cPCR/NR-qPCR) detected a similar proportion of subclinical infections, being 22% (n = 72) and 27% (n = 89), respectively. As the predicted sensitivity of individual PCR assays varied significantly in this group of low-parasite densities, especially with high false negative rate for nested-PCR assay (Table 5), a significant amount of submicroscopic infections was identified only by the real-time PCR assays, i.e., R-qPCR for 18S rRNA and NR-qPCR for Pvr47/Pfr364 sequences. In addition, the majority of submicroscopic infections were identified by both ribosomal and non-ribosomal targets (Fig. 3A, the right-side bar chart, in lilac; p = 0.09 ribosomal vs. both PCR-targets, and p = 0.86 for non-ribosomal vs. both targets). Overall, while PCR protocols identified 92 (28%) subclinical infections, only 21 (7%) were identified by LM (p < 0.0001, Fig. 3), indicating that 71 out of 92 (77%) infections were missed by the routine LM. Species-specific identification demonstrated a high proportion of P. vivax as compared with P. falciparum and mixed infections (Fig. 3B), and confirmed differences in the sensitivity between protocols. Although real-time PCR assays (R-qPCR and NR-qPCR) presented a better performance than conventional PCR assays (NRc-PCR and Nested-PCR), the results confirmed the gains achieved by using NR-qPCR to detect subclinical malaria infection (Fig. 3B). According, considering only the submicroscopic infections identified by any PCR-based assay (73 out of 92), Venn diagram demonstrated a high proportion of infections (18%, 13 out of 73) identified only by the NR-qPCR (Fig. 4). In terms of amplified target region, 96% (n = 70) of submicroscopic infections were detected by non-ribosomal sequences, while 78% (n = 57) were detected by ribosomal target (p = 0.0024, Fisher’s exact test).
Although major advances have been reached for the molecular detection of malaria parasites [25, 26, 36, 37], most sensitive PCR-based assays require high-volume of venous blood and complex sample processing [8, 23, 38, 39], being not feasible in the context of malaria routine surveillance. The current study involved investigate the hypothesis that the amplification of both ribosomal and non-ribosomal multi-copy PCR targets could increase the chances of detecting low parasite density and mixed P. falciparum and P. vivax infection. For that, a non-ribosomal (NR) qPCR targeting the multi-copy Pvr47/Pfr364 sequences was developed and this new protocol was compared with the original non-ribosomal gel-stained PCR-based protocol  as well as with two species-specific PCR assays based on the 18S rRNA gene amplification.
The end-point titration assays of field samples revealed that the NR-qPCR protocol was able to accurately detect both P. vivax and P. falciparum infections—in single and artificial mixed infections—producing reproducible results at the lowest parasite densities (1–3 parasites/µL). Although there was considerable variation between PCR protocols assayed, the non-ribosomal protocols (NR-cPCR and NR-qPCR) were more accurate than ribosomal (nested-PCR and R-qPCR) to detect mixed-species infections. Of interest, only NR-qPCR assay developed here were able to detect P. falciparum when this species was present in a proportion of 240-fold lower than P. vivax. As the sensitivity of any PCR protocol depends largely on the molecular target used , the high copy number of Pfr364 (around 20 copies of “subfamily 1” targeted by specific primers) probably facilitated the detection of low levels of P. falciparum in co-infections as compared to 18S rRNA (around 5–8 copies). Although different multi-copy targets have been described as sensitive for molecular diagnosis of malaria [23, 25, 36], those studies did not investigate the reliability of these targets in mixed-malaria infections, which precludes any potential comparison with results described here. In addition, most of the studies have been carried-out in endemic areas, such as Papua New Guinea, that currently does not represent an unstable and low-transmission endemic area . More work needs to be done in this field of investigation. An apparent inability of 18S rRNA qPCR to detect low P. falciparum densities in situation where P. vivax was present in much higher densities was observed. The use of a single pair of primers to detect both species may have been a determinant factor in causing failure of R-qPCR to identify mixed infections. A similar phenomenon of primer competition was described in the original protocol , straightening that species-specific primers should be used in field studies in which malaria co-infection is expected to be relevant.
In clinical malaria suspects, the overall prevalence for P. vivax and P. falciparum detected by amplification of non-ribosomal Pvr47/Pfr364 targets was not significantly different than that of either conventional microscopy or 18S rRNA gene amplification. The predicted sensitivity and specificity of each PCR protocol assayed here were also similar, and it was independent of the parasite target. Although the clinical sample size limited the statistical power to detect small differences between protocols, these results were not completely unexpected as symptomatic patients usually present high parasite densities in the peripheral blood; consequently, it may facilitate the confirmation of malaria infection by less sensitive protocols such as microscopy and rapid diagnostic tests (RDTs) . These findings reinforce that submicroscopic malaria infections are not prevalent among symptomatic patients, and LM and RDTs are adequate tools for case management [10, 42]. Nevertheless, the limited sensitivity of microscopy in correct identification of mixed-species malaria should be considered in areas where more than one Plasmodium species is circulating [43, 44], a result that was confirmed here.
While low frequencies of submicroscopic infections were observed in the group of clinical malaria cases (3–5%), screening for malaria in cross-sectional surveys demonstrated a large proportion (> 70%) of malaria cases in the study area that was not detected by conventional microscopy. The majority of subclinical infections were caused by P. vivax, the commonest malaria parasite in the Amazon basin, and frequently associated with low-density infections [27, 28, 45,46,47]. These findings are in accordance with recent reports showing high proportions of submicroscopic P. vivax infections across different endemic settings, particularly areas with relatively low transmission intensity [10,11,12, 26]. Although the reason for this high rate of asymptomatic P. vivax infections is unknown, it is probably associated with the unique biology of P. vivax that includes a fast acquisition of clinical immunity as compared with P. falciparum . It is particularly relevant because in different epidemiological settings there are perspectives on treating asymptomatic infections for malaria elimination . In the study area, the results demonstrated that, in general, multiple molecular targets (i.e., ribosomal plus non-ribosomal) did not increase sensitivity to detect subclinical malaria infections. Despite of that, the NR-qPCR developed here was the most sensitive protocol to detect submicroscopic asymptomatic malaria infections, which resulted in a significantly higher prevalence of submicroscopic infections (70 out of 73, 96%) when compared to that detected by ribosomal PCR assays (57 out of 73, 78%). While more sensitive amplification of 18S rRNA gene has been described [39, 49], the likelihood of amplify 18S rRNA gene was dependent on (i) large blood volume (2 mL); (ii) careful removal of plasma and buffy coat as prerequisite to avoid interference during PCR processing; (iii) concentration of purified DNA dehydrated in a centrifugal vacuum concentrator; additionally, these “high-volume” 18S rRNA PCR strategy did not allow the detection by species (only Plasmodium spp.) .
The apparent ability of Pvr47/Pfr364 NR-qPCR to increase sensitivity to investigate the true prevalence of malaria infection is relevant as an unexpectedly large reservoir of infections may hinder control and elimination of malaria in the Americas [3, 50]. These findings are critical as both subclinical and submicroscopic malaria carriers remain untreated in the Brazilian Amazon region and therefore might remain infective over long periods of time . As parasite densities cannot be assumed as a static parameter and thus may fluctuate over time falling below the detection threshold of the assay , future studies should approach longitudinal PCR-malaria surveys. Although the NR-qPCR developed here may constitute powerful additive tools to identify endemic sites where relevant control measures have to be settled and monitored , the costs of PCR-based assays limited such type of study. In general, nucleic acid amplification tests (NAATs) require expensive equipment available, well-equipped laboratories, qualified personnel, and large quantities of disposable supplies that need to be frozen or refrigerated, which is sometimes difficult in the country [41, 53]. Currently, WHO recommends that the use of NAATs be considered only for epidemiological research or surveys mapping submicroscopic infections in low transmission areas . Innovative and cost-effective strategies that identify the real burden of malaria infections (those detected by qPCR) are required to reach malaria elimination goals, but remain a challenge .
Assuming that the NR-qPCR developed here seems to be the most sensitive method—as it was positive in a number of samples not detected by other PCR protocols—the results suggested that NR-qPCR has a lower detection threshold. Despite of that, it is important to clarify the technical limitations that apply for the definition of “reference standard” for PCR-detection of submicroscopic malaria infections. In general, the estimative of test accuracy are based on the assumptions that the reference standard is 100% sensitive and that specific disagreements between the reference standard and the diagnostic test being evaluated (index test) result from incorrect classification by the index test . However, this statement cannot be applied for the detection of low-density malaria infections because there is no “gold-standard”. While the conventional microscopy diagnostic present high number of false negatives at low parasite density , there is no consensus about a PCR assay able to detect all malaria infections . Due to these inherent limitations, the “reference standard” for each molecular diagnostic method was defined as a combination of positive detections by any PCR assay, excluding the method under evaluation, as described before . Consequently, “false positive” in this type of analysis is considered in the sense that no other PCR method found these infections. In this scenario, the findings unlikely represent a tendency to false positive by NR-qPCR as it was established by (i) end-point titration of well-characterized field samples, including mono and artificial mixed-infections; (ii) reproducibility of replicates at low levels of parasitaemia; (iii) no amplification with gDNA samples from malaria-free volunteers; (iv) no cross-reactivity with other Plasmodium species. Furthermore, considering the rules for quality assessment of diagnostic accuracy studies (QUADAS-2) , the risk of bias of the present study was reduced as methodological design involved: (i) structured sample size calculations and random selection of malaria-exposed individuals, with explicit exclusion criteria defined in methods; (ii) in the estimative of sensitivity and/or specificity, the diagnostic test being evaluated was clearly interpreted before the reference standard was known; (iii) the execution of the PCR-based assays and the definition of reference standard were described in sufficient detail to permit replication of the test. Consequently, valuable malaria information can be retrieved from the current study.
Finally, relatively low frequencies of mixed-malaria infections were detected here, which precluded a more detailed evaluation of the potential of Pvr47/Pfr364 to detect mixed-malaria infections in the field. In the Amazon area, besides P. vivax being the predominant malaria parasite [27, 29], the progress achieved in malaria control has decreased the number of P. falciparum cases in recent years . Notwithstanding this study limitation, it is highly relevant the results from the end-point titration experiments showing the ability of Pvr47/Pfr364 to consistently detect P. vivax/P. falciparum co-infection, as the accurate detection of malaria mixed-infections seems to be critical for control and management of malaria [43, 44]. In fact, disease burden due to mixed species infections remains largely unknown, and this limitation have the potential to influence decisions on testing vaccines and new antimalarial drugs . As malaria has been re-emerging in areas where it was previously controlled, dealing with mixed malaria infection cannot be bypassed, as recent evidence suggest that the frequency of these infections may be much higher than previously expected , including in the Amazon region . Due to the risk of P. falciparum reemergence from Amazonian neighboring countries with high transmission rates, a cross-border malaria study to evaluate the relevance of NR-qPCR in mixed-malaria infections are on progress.
Although the simultaneous use of ribosomal and non-ribosomal PCR-targets did not impact the molecular diagnosis of malaria, the amplification of Pfr364 and Pvr47 multi-copy targets by the NR-qPCR seems to be a valuable tool in detection of subclinical and mixed P. vivax/P. falciparum infections, even though one of the species was present in a ratio of hundred-times lower than the other species.
- 18S rRNA:
small subunit 18S of the ribosomal RNA gene
polymerase chain reaction
non-ribosomal real-time PCR
non-ribosomal conventional PCR
ribosomal real-time PCR
limit of detection
WHO. World Malaria Report 2018. Geneva: World Health Organization. 2018. https://apps.who.int/iris/bitstream/handle/10665/275867/9789241565653-eng.pdf?ua=1. Accessed 9 Apr 2019.
Wassmer SC, Grau GER. Severe malaria: what’s new on the pathogenesis front? Int J Parasitol. 2017;47:145–52.
Recht J, Siqueira AM, Monteiro WM, Herrera SM, Herrera S, Lacerda MVG. Malaria in Brazil, Colombia, Peru and Venezuela: current challenges in malaria control and elimination. Malar J. 2017;16:273.
Rabinovich RN, Drakeley C, Djimde AA, Hall BF, Hay SI, Hemingway J, et al. malERA: an updated research agenda for malaria elimination and eradication. PLoS Med. 2017;14:e1002456.
Imwong M, Stepniewska K, Tripura R, Peto TJ, Lwin KM, Vihokhern B, et al. Numerical distributions of parasite densities during asymptomatic malaria. J Infect Dis. 2016;213:1322–9.
Tadesse FG, Van Den Hoogen L, Lanke K, Schildkraut J, Tetteh K, Aseffa A, et al. The shape of the iceberg: quantification of submicroscopic Plasmodium falciparum and Plasmodium vivax parasitaemia and gametocytaemia in five low endemic settings in Ethiopia. Malar J. 2017;16:99.
Lamptey H, Ofori MF, Kusi KA, Adu B, Yeboa EO, Baafour EK, et al. The prevalence of submicroscopic Plasmodium falciparum gametocyte carriage and multiplicity of infection in children, pregnant women and adults in a low malaria transmission area in Southern Ghana. Malar J. 2018;17:331.
Bejon P, Andrews L, Hunt-cooke A, Sanderson F, Gilbert SC, Hill AVS. Thick blood film examination for Plasmodium falciparum malaria has reduced sensitivity and underestimates parasite density. Malar J. 2006;5:104.
Haanshuus CG, Chandy S, Manoharan A, Vivek R. A high malaria prevalence identified by PCR among patients with acute undifferentiated fever in India. PLoS ONE. 2016;11:e0158816.
Cheng Q, Cunningham J, Gatton ML. Systematic review of sub-microscopic P. vivax infections: prevalence and determining factors. PLoS Negl Trop Dis. 2015;9:e3413.
Moreira CM, Shehada MA, Price RN, Drakeley CJ. A systematic review of sub-microscopic Plasmodium vivax infection. Malar J. 2015;14:360.
Almeida ACG, Kuehn A, Castro AJM, Vitor-Silva S, Figueiredo EFG, Brasil LW, et al. High proportions of asymptomatic and submicroscopic Plasmodium vivax infections in a peri-urban area of low transmission in the Brazilian Amazon. Parasit Vectors. 2018;11:194.
Alves FP, Gil LHS, Marrelli MT, Ribolla PEM, Camargo EP, Pereira Da Silva LH. Asymptomatic carriers of Plasmodium spp. as infection source for malaria vector mosquitoes in the Brazilian Amazon. J Med Entomol. 2005;42:777–9.
Schneider P, Bousema JT, Gouagna LC, Otieno S, van de Vegte-Bolmer M, Omar SASR. Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007;76:470–4.
Chen I, Clarke SE, Gosling R, Hamainza B, Killeen G, Magill A, et al. “Asymptomatic” malaria: a chronic and debilitating infection that should be treated. PLoS Med. 2016;13:e1001942.
Snounou G, Viriyakosol S, Zhu PZ, Jarra W, Pinheiro L, Rosario VE, et al. High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol. 1993;61:315–20.
Rougemont M, Van Saanen M, Sahli R, Hinrikson HP, Bille JJK. Detection of four Plasmodium species in blood from humans by 18S rRNA gene subunit-based and species-specific real-time PCR assays. J Clin Microbiol. 2004;42:5636–43.
Mangold KA, Manson RU, Koay ES, Stephens L, Regner M, Thomson RB Jr, Peterson LRKK. Real-time PCR for detection and identification of Plasmodium spp. J Clin Microbiol. 2005;43:2435–40.
Murphy SC, Prentice JL, Williamson K, Wallis CK, Fang FC, Fried M, et al. Real-time quantitative reverse transcription PCR for monitoring of blood-stage Plasmodium falciparum infections in malaria human challenge trials. Am J Trop Med Hyg. 2012;86:383–94.
Kamau E, Alemayehu S, Feghali KC, Saunders D, Ockenhouse CF. Multiplex qPCR for detection and absolute quantification of malaria. PLoS ONE. 2013;8:e71539.
Mercereau-Puijalon O, Barale J, Bischoff E. Three multigene families in Plasmodium parasites: facts and questions. Int J Parasitol. 2002;32:1323–44.
Costa DC, Madureira AP, Amaral LC, Sanchez BAM, Gomes LT, Fernandes Fontes CJ, et al. Submicroscopic malaria parasite carriage: how reproducible are polymerase chain reaction-based methods? Mem Inst Oswaldo Cruz. 2014;109:21–8.
Hofmann NE, Gruenberg M, Nate E, Ura A, Rodriguez-Rodriguez D, Salib M, et al. Assessment of ultra-sensitive malaria diagnosis versus standard molecular diagnostics for malaria elimination: an in-depth molecular community cross-sectional study. Lancet Infect Dis. 2018;18:1108–16.
Demas A, Oberstaller J, DeBarry J, Lucchi NW, Srinivasamoorthy G, Sumari D, et al. Applied genomics: data mining reveals species-specific malaria diagnostic targets more sensitive than 18S rRNA. J Clin Microbiol. 2011;49:2411–8.
Hofmann N, Mwingira F, Shekalaghe S, Robinson LJ, Mueller I, Felger I. Ultra-sensitive detection of Plasmodium falciparum by amplification of multi-copy subtelomeric targets. PLoS Med. 2015;12:e1001788.
Gruenberg M, Moniz CA, Hofmann NE, Wampfler R, Koepfli C, Mueller I, et al. Plasmodium vivax molecular diagnostics in community surveys: pitfalls and solutions. Malar J. 2018;17:55.
Kano FS, Sanchez BAM, Sousa TN, Tang ML, Saliba J, Oliveira FM, et al. Plasmodium vivax Duffy binding protein: baseline antibody responses and parasite polymorphisms in a well-consolidated settlement of the Amazon Region. Trop Med Int Health. 2012;17:989–1000.
Souza-silva A, Brito CFA, Adams JH, Kano FS, Carvalho LH. Duffy antigen receptor for chemokine (DARC) polymorphisms and its involvement in acquisition of inhibitory anti-duffy binding protein II (DBPII) immunity. PLoS ONE. 2014;9:e93782.
Kano FS, Souza-Silva FA, Torres LM, Lima BA, Sousa TN, Alves JR, Rocha RS, Fontes CJ, Sanchez BA, Adams JH, Brito CF, Pires DE, Ascher DB, Sell AMCL. The presence, persistence and functional properties of Plasmodium vivax duffy binding protein II antibodies are influenced by HLA class II allelic variants. PLoS Negl Trop Dis. 2016;10:e0005177.
Ladeia-Andrade S, Ferreira MU, de Carvalho ME, Curado I, Coura JR. Age-dependent acquisition of protective immunity to malaria in riverine populations of the Amazon Basin of Brazil. Am J Trop Med Hyg. 2009;80:452–9.
Pires CV, Alves JRS, Lima BAS, Paula RB, Sousa N, Soares IS, et al. Blood-stage Plasmodium vivax antibody dynamics in a low transmission setting: a nine year follow-up study in the Amazon region. PLoS ONE. 2018;13:e0207244.
Ministério da Saúde. Manual de Diagnóstico Laboratorial da Malária. Brasília: Secretaria de Vigilância em Saúde. 2009. http://bvsms.saude.gov.br/bvs/publicacoes/manual_diagnostico_laboratorial_malaria_2ed.pdf. Accessed 9 Apr 2019.
Ministério da Saúde. Guia prático de tratamento da malária no Brasil. Brasília: Secretaria de Vigilância em Saúde. 2010. http://bvsms.saude.gov.br/bvs/publicacoes/guia_pratico_malaria.pdf. Accessed 9 Apr 2019.
Olsson ML, Hansson C, Avent ND, Åkesson IE, Green CA, Daniels GL. A clinically applicable method for determining the three major alleles at the Duffy (FY) blood group locus using polymerase chain reaction with allele-specific primers. Transfusion. 1998;38:168–73.
Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193:673–5.
Lloyd YM, Esemu LF, Antallan J, Thomas B, Yunga ST, Obase B, et al. PCR-based detection of Plasmodium falciparum in saliva using mitochondrial cox3 and varATS primers. Trop Med Health. 2018;46:22.
Meerstein-Kessel L, Andolina C, Carrio E, Mahamar A, Sawa P, et al. A multiplex assay for the sensitive detection and quantification of male and female Plasmodium falciparum gametocytes. Malar J. 2018;17:441.
Andrews L, Andersen RF, Webster D, Dunachie S, Walther RM, Bejon P, et al. Quantitative real-time polymerase chain reaction for malaria diagnosis and its use in malaria vaccine clinical trials. Am J Trop Med Hyg. 2005;73:191–8.
Imwong M, Hanchana S, Malleret B, Rénia L, Day NPJ, Dondorp A, et al. High-throughput ultrasensitive molecular techniques for quantifying low-density malaria parasitemias. J Clin Microbiol. 2014;52:3303–9.
Lucchi NW, Oberstaller J, Kissinger JCUV. Malaria diagnostics and surveillance in the post-genomic era. Public Health Genomics. 2013;16:37–43.
Zimmerman PA, Howes RE. Malaria diagnosis for malaria elimination. Curr Opin Infect Dis. 2015;28:446–54.
Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.
Ehtesham R, Fazaeli A, Raeisi A, Keshavarz H, Heidari A. Detection of mixed-species infections of Plasmodium falciparum and Plasmodium vivax by nested PCR and rapid diagnostic tests in southeastern Iran. Am J Trop Med Hyg. 2015;93:181–5.
Singh US, Siwal N, Pande V, Das A. Can mixed parasite infections thwart targeted malaria elimination program in India? Biomed Res Int. 2017;2017:2847548.
Oliveira-ferreira J, Lacerda MVG, Brasil P, Ladislau JLB, Tauil PL, Daniel-ribeiro CT. Malaria in Brazil: an overview. Malar J. 2010;9:115.
Sampaio VS, Siqueira AM, Costa G, Paula M, Mourão G, Monteiro WM, et al. Malaria in the State of Amazonas: a typical Brazilian tropical disease influenced by waves of economic development. Rev Soc Bras Med Trop. 2015;48(Suppl 1):4–11.
Adams JH, Mueller I. The biology of Plasmodium vivax. Cold Spring Harb Perspect Med. 2017;7:9.
Jaiteh F, Masunaga Y, Okebe J, Alessandro UD, Balen J, Bradley J, et al. Community perspectives on treating asymptomatic infections for malaria elimination in The Gambia. Malar J. 2019;18:39.
Imwong M, Nguyen TN, Tripura R, Peto TJ, Lee SJ, Lwin KM, et al. The epidemiology of subclinical malaria infections in South-East Asia: findings from cross-sectional surveys in Thailand–Myanmar border areas, Cambodia, and Vietnam. Malar J. 2015;14:381.
Ferreira MU, Castro MC. Challenges for malaria elimination in Brazil. Malar J. 2016;15:284.
Nguyen T, Seidlein LV, Nguyen T, Truong P, Hung SD, Pham H, et al. The persistence and oscillations of submicroscopic Plasmodium falciparum and Plasmodium vivax infections over time in Vietnam: an open cohort study. Lancet Infect Dis. 2018;18:565–72.
Niang M, Diop F, Niang O, Sadio BD, Sow A, Faye O, et al. Unexpected high circulation of Plasmodium vivax in asymptomatic children from Kédougou, southeastern Senegal. Malar J. 2017;16:497.
Berzosa P, De Lucio A, Barja MR, Herrador Z, González V, García L, et al. Comparison of three diagnostic methods (microscopy, RDT, and PCR) for the detection of malaria parasites in representative samples from Equatorial Guinea. Malar J. 2018;17:333.
WHO. Nucleic acid amplification-based diagnostics. Geneva: World Health Organization. 2018. https://www.who.int/malaria/areas/diagnosis/nucleic-acid-amplification-tests/en/. Accessed 9 Apr 2019.
Biesheuvel C, Irwig L, Bossuyt P. Observed differences in diagnostic test accuracy between patient subgroups: is it real or due to reference standard misclassification? Clin Chem. 2007;53:1725–9.
Whiting PF, Rutjes AWS, Westwood ME, Mallett S, Deeks JJ, Reitsma JB, et al. QUADAS-2: a revised tool for the quality assessment of diagnostic accuracy studies. Ann Intern Med. 2011;155:529–36.
Siqueira AM, Mesones-lapouble O, Marchesini P, Sampaio VDS, Brasil P, Tauil PL, et al. Plasmodium vivax landscape in Brazil: scenario and challenges. Am J Trop Med Hyg. 2016;95(Suppl 6):87–96.
Zimmerman PA, Mehlotra RK, Kasehagen LJKJ. Why do we need to know more about mixed Plasmodium species infections in humans? Trends Parasitol. 2004;20:440–7.
Camargo M, León SCS, Del Río-ospina L, Páez AC, González Z, González E, et al. Micro-epidemiology of mixed-species malaria infections in a rural population living in the Colombian Amazon region. Sci Rep. 2018;8:5543.
LHC, CFAB, and TNS led the conception and study design. LCA contributed to the study design, performed PCR assays, and data analysis. DRR and LFFG performed PCR assays based on the amplification of 18S rRNA gene. JEL, CJFF, DBP and FSK participated in the field study. LCA, TNS and LHC wrote the manuscript. All authors read and approved the final manuscript.
The authors would like to thank the MR4/Bei Resources/NIAID/NIH/ATCC for providing the genomic DNA of the P. brasilianum Peruvian III strain (MRA-349). We also thank the Programme for Technological Development in Tools for Health-PDTIS-FIOCRUZ for use of the DNA sequencing (RPT01E) and Real-Time PCR (RPT09D) facilities. This work was supported by the Brazilian National Research Council (CNPq), the Coordination for the Improvement of Higher Education Personnel (CAPES), the Research Foundation of Minas Gerais State (FAPEMIG) and Program of Excellence in Research/Diagnosis (Proep/IRR/FIOCRUZ). Fellowships from CNPq (LHC, CFAB, TNS, CJFF, DRR) and CAPES (LCA, LFFG) are also acknowledged. The authors declare that the study was conducted in the absence of any financial relationships that could be construed as potential conflicts of interest.
The authors declare that they have no competing interests.
Availability of data and materials
The authors confirm that all data reported in the manuscript are publicly available.
Consent for publication
Ethics approval and consent to participate
Ethical and methodological aspects of this study were approved by the Ethical Committee of Research on Human Beings from the René Rachou Institute/Oswaldo Cruz Foundation (Protocols Nos. 24/2008, and 1.821.955/2016), according to the Brazilian National Council of Health (Resolutions 196/96 and 466/12). All participants were informed about the objectives and procedures of the study, with voluntary participation through written informed consent.
The study was supported by the National Council for Scientific and Technological Development (CNPq), the Research Foundation of Minas Gerais (FAPEMIG), and the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brasil (CAPES; Finance code 001)—Finance Code 001. Scholarships from the Coordination for the Improvement of Higher Education Personnel/CAPES (LCA, LFFG) and CNPq (DRR, LHC, CFAB, TNS, CJFF) are also acknowledged. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Analytical evaluation of NR-qPCR assay.
Amplification curves of NR-qPCR performed on field samples infected with P. vivax (n = 3; red), P. falciparum (n = 3; blue), P. malariae (n = 3; green), and P. brasilianum DNA (n = 1; orange). The amplification plots are shown for (A) Pvr47 and (B) Pfr364 assays.
Five-fold dilution of P. vivax-Pvr47 and P. falciparum-Pfr364 plasmids amplified by NR-qPCR.
Conditions of the Nested-PCR, R-qPCR, and NR-cPCR assays. The primers/probes used for P. vivax and P. falciparum targets were the original described. The products of nested-PCR and NR-cPCR were visualized by 2% agarose gel stained with ethidium bromide.
Limit of detection (LOD) for (A) Pvr47 and (B) Pfr364 targets amplified by NR-qPCR. Probit regression analysis was used on logarithmic scale using nine points of NR-qPCR standard curves (20,000 to 0.05 copies/μL). The calculated regression curves (blue lines) indicate the probability (y-axis) of obtaining positive results at any template concentration, and dashed brown lines shows 95% confidence intervals. Dashed black lines correspond to the lower DNA concentration in which 95% of positive samples were detected (0.66 copies/µL of Pvr47 and 3.27 copies/µL of Pfr364).
About this article
Cite this article
Amaral, L.C., Robortella, D.R., Guimarães, L.F.F. et al. Ribosomal and non-ribosomal PCR targets for the detection of low-density and mixed malaria infections. Malar J 18, 154 (2019). https://doi.org/10.1186/s12936-019-2781-3
- Molecular diagnosis
- Mixed-malaria infections