Comparison of molecular quantification of Plasmodium falciparum gametocytes by Pfs25 qRT-PCR and QT-NASBA in relation to mosquito infectivity
Malaria Journal volume 15, Article number: 539 (2016)
Quantifying gametocyte densities in natural malaria infections is important to estimate malaria transmission potential. Two molecular methods (Pfs25 mRNA quantitative reverse transcriptase PCR (qRT-PCR) and Pfs25 mRNA quantitative nucleic acid sequence based amplification (QT-NASBA)) are commonly used to determine gametocyte densities in clinical and epidemiological studies and allow gametocyte detection at densities below the microscopic threshold for detection. Here, reproducibility of these measurements and the association between estimated gametocyte densities and mosquito infection rates were compared.
To quantify intra- and inter-assay variation of QT-NASBA and qRT-PCR, a series of experiments was performed using culture-derived mature Plasmodium falciparum gametocytes from three different parasite isolates (NF54, NF135, NF166). Pfs25 mRNA levels were also determined in samples from clinical trials in Mali and Burkina Faso using both methods. Agreement between the two methods and association with mosquito infection rates in membrane feeding assays were assessed.
Intra- and inter-assay variability was larger in QT-NASBA compared to qRT-PCR, particularly at low gametocyte densities (< 1 gametocyte per μL). Logistic models, including log-transformed gametocytaemia estimated by QT-NASBA, explained variability in mosquito feeding experiment results as well as log-transformed gametocytaemia by qRT-PCR (marginal R2 0.28 and 0.22, respectively). Densities determined by both methods strongly correlated with mosquito infection rates [Spearman’s rank correlation coefficient, 0.59 for qRT-PCR and 0.64 for QT-NASBA (P < 0.001 for both)]. Gametocyte densities estimated by qRT-PCR were higher than levels estimated by QT-NASBA or light microscopy at high densities (>100 gametocyte per μL). Samples collected in one of the two transmission studies had extremely low gametocyte densities by both molecular methods, which is suggestive of RNA degradation due to an unknown number of freeze–thaw cycles and illustrates the reliance of molecular gametocyte diagnostics on a reliable cold-chain.
The experiments indicate that both qRT-PCR and QT-NASBA are of value for quantifying mature gametocytes in samples collected in field studies. For both assays, estimated gametocyte densities correlated well with mosquito infection rates. QT-NASBA is less reproducible than qRT-PCR, particularly for low gametocyte densities.
The transmission of malaria from humans to mosquitoes requires the presence of mature sexual stage parasites (gametocytes) in the peripheral blood. Upon ingestion by Anopheles mosquitoes and following sporogonic development, these gametocytes can render mosquitoes infectious to humans. Quantifying the density of these gametocytes in the peripheral blood is a fundamental part of estimating the infectiousness of individuals with malaria infections. For over a century, light microscopy was the only method available to quantify malaria parasites, including gametocytes. The development of sensitive molecular methods has uncovered some of the limitations of microscopy. In endemic areas, a substantial proportion of infected individuals carry parasites at levels below the microscopic threshold of detection [1–4], and many of these individuals also have low, sub-microscopic, densities of gametocytes [5, 6]. Together with early observations that individuals with no gametocytes detected by microscopy could be infectious to Anopheles mosquitoes , these studies suggest that the use of more sensitive methods is necessary to characterize the infectious reservoir of malaria [8, 9].
The most widely used target for molecular gametocyte detection and quantification is Pfs25 mRNA , which is gametocyte-specific and highly conserved among different parasite isolates . Although expression of Pfs25 mRNA is upregulated in female gametocytes , the female-biased gametocyte sex ratio and the high abundance of Pfs25 mRNA relative to the male specific marker Pfs230p make it an attractive target for sensitive gametocyte detection and quantification .
Various molecular techniques have been used for gametocyte detection and quantification. For Pfs25 mRNA there are two commonly used protocols, one based on quantitative reverse transcriptase PCR (qRT-PCR) and another based on quantitative nucleic acid sequence based amplification (QT-NASBA) [6, 12]. Although previous studies on other malaria-specific targets suggest that the accuracy, precision and operational attractiveness may differ between assays , QT-NASBA and qRT-PCR, as currently routinely used, have never been directly compared for Plasmodium falciparum gametocyte quantification. Additionally, differences in Pfs25 transcript levels between malaria parasite strains have been hypothesized but never directly examined .
Here, Pfs25 mRNA QT-NASBA and qRT-PCR were compared for intra- and inter-assay variation using gametocytes from three different isolates [6, 12]. The association between gametocyte densities estimated by these two molecular methods in natural malaria infections and infectiousness quantified by membrane feeding assays was also determined.
In this study, both parasite culture-derived samples and samples from naturally infected individuals were used. For the estimation of intra-assay variation for both QT-NASBA and qRT-PCR, culture-derived mature P. falciparum gametocytes at known densities were assayed in triplicates on the same 96-well plate with the same reaction mixtures. For the inter-assay variation analysis, gametocyte dilution series were tested on separate plates and with separately prepared reaction mixtures. Samples collected from naturally infected individuals living in malaria-endemic regions were used to assess the correlation between gametocyte levels measured by microscopy, qRT-PCR and QT-NASBA and human infectiousness to mosquitoes, determined by mosquito feeding experiments.
Plasmodium falciparum gametocytes from three laboratory parasite strains (NF54, NF135, NF166) were cultured in shaker flasks, as previously described [15–17]. Briefly, asexual blood stage parasites were reseeded at 5% haematocrit (erythrocytes provided by healthy Dutch donors with blood type A) and 0.5% parasitaemia in culture medium containing RPMI 1640 with HEPES (5.94 g/L), hypoxanthine (0.05 g/L), 10% (v/v) pooled human sera (also blood type A) obtained from malaria-naive individuals, and 0.2% (w/v) sodium bicarbonate. Culture medium was automatically changed twice a day, and on day 4, 50 mM of N-acetylglucosamine was added in order to prevent the next generation of merozoites from infecting new erythrocytes. Gametocytes were allowed to mature until day 14.
Gametocyte dilution series
On day 14, mature gametocytes were harvested and Percoll purified according to the protocol described in Kariuki et al. ; 63%, instead of 65%, Percoll gradient was used. Cultures and material used for purification were continuously kept at 37 °C to avoid gametocyte activation. Following Percoll purification, gametocytes were counted by a single reader using a haemocytometer and diluted in whole blood at densities of 103, 102, 101, 100, 10−1, and 10−2 or 5 × 103, 5 × 102, 5 × 101, 5 × 100, 5 × 10−1, and 5 × 10−2 gametocytes per µL. Fifty µL aliquots of each dilution were further diluted in 250 µL of RNAprotect Cell Reagent (Qiagen) when analysing standard curve triplicates and samples collected in Burkina Faso (see below) or in 450 µL of L6 buffer when analysing samples collected in Mali, and stored at −80 °C until further use .
Samples from naturally infected individuals in clinical trials
Blood samples from patent gametocyte carriers were collected during clinical trials in Mali and Burkina Faso [20, 21]. In the study undertaken in Mali, male participants with at least 32 gametocytes per μL of blood based on gametocyte enumeration against 500 white blood cells (WBC), assuming 8000 WBCs/µL, were recruited. The trial in Burkina Faso had as an inclusion criterion the presence of at least one microscopically detectable gametocyte at screening slides, where 100 microscopic fields were examined. In both studies, blood samples were taken at multiple time-points, before and after treatment initiation, to assess gametocyte clearance. Samples [venous blood in EDTA (Mali) or heparin (Burkina Faso) vacutainer] taken before initiation of anti-malarials were used in the current study. In Burkina Faso, 50 µL aliquots of blood were transferred into 250 µL of RNAprotect Cell Reagent. In Mali, 100 µL of blood were transferred into 900 µL of L6 buffer. Both sets of samples were kept at −80 °C until shipment, although the samples collected in Burkina Faso may have been exposed to freeze–thaw due to power-cut of unknown duration during political unrest in November 2014. For both studies, samples were shipped in dry ice at controlled temperature (−70 to −80 °C) to Radboud University Medical Centre (Nijmegen, The Netherlands), where the laboratory work was performed.
For all samples, total nucleic acids (NA) were extracted using a MagNAPure LC automatic extractor (Total Nucleic Acid Isolation Kit–High Performance, Roche Applied Science) and eluted into 50 µL of elution buffer; while all of the 300 µL of samples collected in Burkina Faso were extracted, only half of the volume (500 µL of the 1 mL) of samples collected in Mali were used for total NA extraction.
Mosquito membrane feeding experiments
All individuals whose samples were included in this analysis participated in mosquito feeding experiments: a small quantity of heparinized venous blood (400–500 μL) was offered to female Anopheles gambiae mosquitoes (180 mosquitoes in the Mali study and 60 in the Burkina Faso study) through a membrane [20–22]. Fully fed mosquitoes at the end of each experiment were selected, kept for one week on glucose at 27–29 °C to allow for parasite development, and dissected. In both studies, light microscopy was used to detect oocysts on the midgut wall, and their presence was confirmed by a second microscopist. The proportion of mosquitoes with at least one oocyst was used as outcome measure.
Pfs25 mRNA qRT-PCR
The use of reverse transcriptase PCR (RT-PCR) for the detection of Pfs25 mRNA was first described by Babiker et al. ; recently it was further developed into a qRT-PCR . Here, a protocol described by Schneider et al. was used  with primers designed by Wampfler et al. . Briefly, total NA extraction was followed by RQ1 DNaseI digest (Promega) and cDNA synthesis (High Capacity cDNA Reverse Transcription Kit, Applied Biosystems). The DNase treatment and the RT-step were done once for each sample; the qPCR step was done in triplicate for samples included in the intra-assay variation analysis. Success of the DNase treatment was tested by performing qPCR, with the same primers, on two to five randomly selected samples of each plate (in total 35 samples) before proceeding to the RT-step. No residual genomic DNA was detected in these assays, which indicates successful DNase treatment. PCR conditions were according to manufacturer’s instructions, described as standard cycling programme of GoTaq® qPCR Master Mix (Promega). qRT-PCR cycle of threshold (Ct) values (the cycle at which the fluorescence signal crosses a pre-determined threshold during the exponential phase of the amplification) were obtained using CFX96™ Real-Time PCR Detection System (BIO-RAD). For samples from naturally infected individuals, Ct values were converted into gametocyte density using plate-specific standard curve. Samples with estimated densities below 0.02 gametocytes per µL or 1 gametocyte per 50 µL of sample were considered negative. At the highest density of the gametocyte dilution series, cDNA corresponding to 363.6 gametocytes was estimated to have been included in qRT-PCR reactions (Additional file 1: Figure S1).
Pfs25 mRNA QT-NASBA
QT-NASBA was performed as in Schneider et al. . KCl concentrations in experiments performed for this study varied from 50 to 60 mM, and differed from the concentration used in the original protocol (80 mM) . For samples from the Burkina Faso study, the extracted total NA used for qRT-PCR was also used for QT-NASBA (before DNase treatment and cDNA production). For samples collected in the Mali trial, different aliquots of the same sample were extracted (total NA) and used for the different molecular assays; the QT-NASBA aliquot did not undergo DNAse treatment and cDNA production. QT-NASBA results are expressed as time to positivity (TTP), as in this method amplification happens at a constant temperature, not in cycles. Standardized manual thresholds to remove background noise are set at the begin value of 3 min and the end value of 15 min. As with qRT-PCR, gametocyte densities were assigned based on plate-specific gametocyte dilution series, which was diluted in whole blood before extraction of total NA. For samples from infected individuals, estimated gametocyte densities below 0.02 gametocytes per µL were considered to be negative. At the highest density of the gametocyte dilution curve, total NA corresponding to 2500 gametocytes was estimated to have been included in QT-NASBA reactions (Additional file 1: Figure S1).
Analyses were conducted using STATA version 12 (Stata Corporation, College Station, TX, USA) and R (R Foundation for Statistical Computing, Vienna, Austria).
Intra-assay variation was assessed by estimating coefficients of variation (CV, expressed in %) (i.e., standard deviation divided by mean, and multiplied by 100) for each level of gametocytaemia tested in triplicate, in culture-specific gametocyte dilution series. Gametocyte dilution series were also used to quantify inter-assay variation: CV was estimated for each gametocyte density for all series tested in different reaction plates. To analyse the agreement between qRT-PCR- and QT-NASBA-defined gametocyte levels in samples from infected individuals in the Mali study, a Bland–Altman plot was constructed, whereby differences between densities assigned by these two methods are presented by the average assigned level. Generalized linear mixed models, with study participant as random effect, were fit to assess the association between mosquito infection status in feeding assays and gametocyte densities (fixed effect) assigned by these two molecular methods. For each quantification method, Akaike information criterion was used to select gametocyte density scale (linear or logarithmic) in the models. Marginal R2, variance explained by fixed effects, was estimated using the ‘r.squaredGLMM’ function in MuMIn package in R .
Plasmodium falciparum gametocytes of the NF54, NF135 and NF166 strains were used to prepare dilution series with 103, 102, 101,100, 10−1, and 10−2 gametocytes per μL. Additionally, for the NF54 strain, samples originating from three different culture flasks were tested. For each gametocyte density of individual dilution series, Pfs25 mRNA expression was analysed in triplicate, using both qRT-PCR and QT-NASBA.
In Fig. 1 (panels a, b), TTP and Ct values for QT-NASBA and qRT-PCR, respectively, are presented. For both assays, Pfs25 mRNA was not detected at the lowest gametocyte densities (0.01 gametocytes per µL) of the NF135 and NF166 dilution series. Although results from dilution series appeared comparable between all three P. falciparum strains, the failure to detect the lowest gametocyte concentration for NF135 and NF166 might indicate lower Pfs25 transcript levels compared to the NF54 strain that only affects detection and quantification at extremely low densities. Additionally, in 5/15 0.10-gametocytes-per-µL samples, Pfs25 mRNA was not detected by qRT-PCR. In samples where amplified Pfs25 transcript was detected CV were inversely related to gametocyte levels and much lower for qRT-PCR compared to QT-NASBA at all gametocyte densities.
Seven (including the following densities 5 × 103, 5 × 102, 5 × 101, 5 × 100, 5 × 10−1, and 5 × 10−2 gametocytes per μL) and ten (including the densities 103, 102, 101,100, 10−1, and 10−2 gametocytes per μL) gametocyte dilution series, each tested on a different plate, were used to assess QT-NASBA and qRT-PCR inter-assay variation, respectively. Only NF54 strain parasite cultures were used in these experiments. Samples with the lowest density in each series were often negative for Pfs25 mRNA. Despite the higher densities used for the QT-NASBA versus qRT-PCR dilution series, a similar pattern to the intra-assay variation experiments was observed (Fig. 1, panels g, h): CV were considerably higher for QT-NASBA compared to qRT-PCR and decreased with increasing gametocyte densities.
Gametocyte densities in naturally infected individuals
Blood samples collected prior to treatment from gametocytaemic individuals participating in mosquito feeding experiments in Mali were analysed by microscopy, qRT-PCR and QT-NASBA. For the two molecular assays, gametocyte densities were estimated based on dilution series of culture-derived P. falciparum gametocytes of the NF54 strain tested in the same plate, as routinely done [5, 12]. As a consequence of enrolment criteria, all individuals (N = 73) had microscopically detectable gametocytes at the time of the feeding assay (median (IQR) 80 (64–112) gametocytes per μL), and all their samples were assigned gametocyte levels above the threshold of positivity for both molecular assays (minima gametocyte levels were 0.04 and 6.1 gametocytes/µL for qRT-PCR and QT-NASBA, respectively). Median (IQR) gametocytaemias estimated by QT-NASBA and qRT-PCR were similar: 94.8 (41.6–180.7) and 69.8 (29.8–263.0) gametocytes per μL, respectively. The agreement between microscopy and QT-NASBA, between microscopy and qRT-PCR and between the two molecular methods was assessed by Bland–Altman plots (Fig. 2). At high gametocyte concentrations (>100 gametocytes per μL), gametocyte densities estimated by qRT-PCR were higher than densities estimated by microscopy or QT-NASBA.
In the same study in Mali, 73 membrane feeding assays were performed. The median (IQR) number of mosquitoes fed and dissected was 176 (174–177) and 142 (130–152), respectively. Overall 58/73 (79.4%) individuals infected at least one mosquito, and the median (IQR) mosquito infection rate (proportion of infected mosquitoes, i.e., proportion of mosquitoes with ≥ one oocyst) was 12.4% (3.4–32.6). Gametocyte levels were lower in non-infectious versus infectious individuals (Additional file 2: Figure S2), and correlated with mosquito infection rates (Spearman’s rank correlation coefficient 0.59, 0.62 and 0.64 (P < 0.001 for all assays) for qRT-PCR, microscopy and QT-NASBA, respectively). Statistical models including log-transformed gametocyte levels estimated by QT-NASBA or microscopy explained a slightly higher proportion (marginal R2 = 0.28 and 0.26, respectively) of infectiousness variability compared to models including log-transformed densities by qRT-PCR (marginal R2 = 0.22) (Fig. 3).
Lower-than-expected gametocyte densities in molecular assays
Venous samples collected from 77/79 children participating in mosquito feeding experiments in Burkina Faso were also analysed. Of these 77 children, 27 did not carry microscopically detectable gametocyte densities on the day of their membrane feeding assay, despite having gametocytes detected by microscopy during the screening phase of the study (one to three days earlier). Median (IQR) gametocyte density determined by microscopy was 12 (0–24) gametocytes per μL. Twenty of 77 did not have Pfs25 mRNA detected by qRT-PCR, and four had assigned gametocyte levels between 0 and 0.02 gametocytes per μL, and were considered gametocyte negative. However, ten of these 24 children with no gametocytes detected by qRT-PCR had gametocytes detected by microscopy. Gametocyte densities estimated by QT-NASBA were also very low for most samples: 55 children, including 31 with microscopically detectable gametocytes, had assigned levels below 0.02 gametocytes per μL (and were thus classified as gametocyte negative) and for one child of 76 Pfs25 mRNA was not detected. All 23 samples with no Pfs25 mRNA amplification or with densities below 0.02 gametocytes per μL by qRT-PCR and that were also tested by QT-NASBA had assigned densities below 0.02 gametocytes per μL by this method. Considering that all these children had microscopically detectable gametocytes one to three days before sample collection, it is likely that RNA material was of insufficient quality. It is highly plausible that freeze–thaws with unknown duration or frequency have affected RNA integrity of this set of samples, all kept in the same box throughout the storage period and shipment. The infectivity of study participants in the clinical trial in Burkina Faso was lower than the infectivity of individuals enrolled in the Malian study: 29/77 (37.6%) children infected at least one mosquito in feeding experiments. Although differences in the prevalence of host factors that influence infectivity, such as transmission-blocking immunity  or haemoglobinopathies , could have contributed to this observation, this disparity in infectiousness prevalence is probably due to real differences in gametocyte levels as suggested by microscopy-based quantification and differences in the number of dissected mosquitoes (median number of mosquitoes dissected per assay, 142 and 45 in Mali and Burkina Faso, respectively). The facts that five of 29 and 17 of 29 individuals who infected mosquitoes did not have Pfs25 mRNA detected by qRT-PCR or had very low assigned levels (<0.02 gametocytes per μL) by QT-NASBA, respectively, and that the gametocyte density-mosquito infection rate plot is shifted to the left compared to the corresponding graph for the Malian study (Additional file 3: Figure S3) support the freeze–thaw explanation.
In this study, Pfs25 mRNA QT-NASBA and qRT-PCR were evaluated using culture-derived gametocyte samples and samples from naturally infected individuals. In in vitro experiments, both intra- and inter-assay variations were lower for qRT-PCR compared to QT-NASBA. When analysing samples from naturally infected individuals, gametocyte densities estimated by qRT-PCR were higher than densities estimated by QT-NASBA and microscopy in individuals with high gametocyte levels. Gametocyte densities estimated by the two molecular assays, as well as by microscopy, correlated well with mosquito infection rates. Overall, the findings indicate that both assays are useful for quantifying sexual stage parasite densities in samples from naturally infected individuals, but qRT-PCR showed better reproducibility.
Determining the variability of quantitative assays is important for the interpretation of results. The experiments here suggest that qRT-PCR has lower intra- and inter-assay variation compared to QT-NASBA. The higher precision of qRT-PCR is evident for the entire range of gametocyte densities used in the dilution series, despite involving extra reaction steps, such as DNase treatment and cDNA production, which are known to be a source of variation . These two steps were performed once for each dilution series density, for the three different NF54 cultures and the NF135 and the NF166 cultures. There was considerably less between-strain and -culture variation for similar densities in qRT-PCR versus QT-NASBA. In QT-NASBA reactions, from the time of primer depletion, amplification will mostly depend on the initial number of transcripts, resulting in a linear increase of amplicons. This exact moment when amplification starts may differ between samples depending on when and how the enzyme is added into the reaction wells, which adds a level of uncertainty to this methodology . The range of TTP values in dilution series is also sensitive to the KCl concentration in the mix, which requires regular optimization with new targets, batches of primers or molecular beacons .
Despite qRT-PCR being a more precise method for gametocyte quantification, QT-NASBA might detect a larger number of very low gametocyte densities (<0.1 gametocytes per μL). QT-NASBA was originally developed for highly sensitive qualitative detection. The apparent higher sensitivity of QT-NASBA compared to qRT-PCR may be explained by methodological differences: for QT-NASBA it is not necessary to dilute samples, in order to complete DNase treatment and cDNA production (Additional file 1: Figure S1) and the equivalent of a larger volume of the original sample is therefore added into the final reaction in QT-NASBA compared to qRT-PCR. This additional dilution in the qRT-PCR protocol could be avoided by using column-based RNA extraction with on-column DNase digestion . Also, while qRT-PCR, in the protocol used here, carries on for a set amount of 40 cycles, the QT-NASBA reaction is continuous at a stable temperature for the duration of 90 min. This might allow for more sensitive amplification of low numbers of transcripts, but could also lead to the detection of single stranded DNA, in cases in which genomic DNA has been degraded and exists in the sample in a single stranded form. A stringent cut-off density below which samples are considered negative must therefore be applied in QT-NASBA, as with all other molecular assays.
Whether qRT-PCR overestimates true gametocyte levels or QT-NASBA and microscopy underestimate densities is not clear. When qRT-PCR is compared with the other methods, high mean densities correspond to the qRT-PCR measurement being larger than the other measurements. There are two possible reasons for observing this pattern: (1) qRT-PCR overestimates high gametocyte concentrations; or, (2) qRT-PCR measurements have greater variability . Since the in vitro experiments suggest that qRT-PCR has lower variability than QT-NASBA, it is likely that qRT-PCR has an upward bias at high gametocyte concentrations (>100 gametocytes per μL). The overestimation is unlikely to be caused by residual DNA in the RNA samples used in the qRT-PCR analysis, as the success of the DNase treatment was tested for two to five samples per plate, with no amplification detected.
The value of molecular methods in detecting epidemiologically relevant malaria infections is becoming increasingly well established. A recent meta-analysis comparing PCR, microscopy and rapid diagnostic test (RDT) sensitivities showed that on average RDTs detect less than half of PCR-positive P. falciparum infections, and microscopy detects slightly less infections than RDTs . The density of these sub-microscopic infections and their concurrent gametocyte levels are of relevance in estimating the contribution to malaria transmission. For P. falciparum parasite density (including asexual parasites and gametocytes), RNA-based QT-NASBA was previously compared to DNA-based qPCR for the quantification of P. falciparum 18S RNA and DNA, respectively . Another recent study, this time on Leishmania parasites, compared 18S rRNA-based quantification by qRT-PCR and QT-NASBA to 18S rRNA gene DNA-based detection with qPCR . In the first study, QT-NASBA was found to be more convenient and equally applicable for quantification purposes as real-time PCR, with strong correlation with quantification by microscopy and similar inter-assay variation to both methods . In the second study, intra- and inter-assay CV were deemed equal for all three molecular methods and qRT-PCR was preferred over the other two methods for convenience reasons . One explanation for these conflicting conclusions is the method used for detection of the amplified target in QT-NASBA: in the study comparing Leishmania quantification, the electrochemiluminescence (ECL) read-out added to the general workload of the assay. Additionally, co-amplification with a known amount of quantitative (Q)-RNA, as was originally done for the quantification of the different gametocyte stages as well as the Leishmania parasites, added an extra factor to the analysis of results [12, 31]. The direct comparison to a standard curve with microscopically determined density of parasites, in addition to the use of molecular beacons, simplified the quantification process of QT-NASBA results . As mentioned earlier in the context of assay sensitivity, also convenience-wise, an important advantage of QT-NASBA is that it can be performed directly on extracted NA.
Several factors, such as the type of anticoagulant used in sample collection and the storage conditions, might influence the outcome of molecular assays on samples [9, 32]. The analysis of samples from Burkina Faso provides indirect evidence for this: the extremely low assigned gametocyte levels, including in samples collected from children with microscopically detectable gametocytes, suggest that these samples might have gone through more than one freeze-thaw cycle. In these situations, sample quality may be compromised: mRNA may be degraded and, although abundant transcripts, such as Pfs25 mRNA, might still be detectable, quantification becomes less reliable . Quantification of constitutively expressed human RNA targets and comparison of transcript levels among samples is one method that could have confirmed whether RNA degradation occurred and is recommended for future studies where there is uncertainty about sample integrity.
Ultimately, the goal in quantifying gametocytes is to indirectly estimate human infectivity, as mosquito feeding assays are logistically complex and only a handful of research institutes currently have the infrastructure to perform these experiments in sub-Saharan Africa. In this study, it was observed that QT-NASBA- and microscopy-defined densities explain a slightly higher proportion of the variation in infectivity compared to gametocyte levels assigned by qRT-PCR but this difference should be interpreted with caution. The statistical method used here to evaluate this relationship was recently employed to assess the association between viraemia and mosquito infection risk in feeding assays involving dengue-infected individuals . While the sigmoidal curves estimated by these models seem to fit mosquito infection data well for QT-NASBA and microscopy, it is possible that qRT-PCR-defined high densities were overestimated (Fig. 2b) and that could possibly explain the poorer fit of the model for that assay. These comparisons are based on a statistical model. More complex models , allowing for more flexible sigmoidal curves, would be required to formally assess the shape of the association between gametocyte density and mosquito infection rates. The quantification of other parasite (e.g., sex ratio) or host factors (e.g., transmission-blocking immune responses and haematological factors) that might influence infectivity is likely to improve the predictive value of these models . Of note, the enrolment criterion of multiple gametocytes being observed by microscopy (two or more gametocytes per 500 WBC) makes it difficult to directly extrapolate these results to sub-microscopic infections. By definition, microscopy has no quantitative value in sub-patent infections and molecular methods will need to be used to characterize the gametocytaemia-infectivity curve at low gametocyte densities.
Despite being the most commonly used amplification target for estimating mature P. falciparum gametocytes counts, the use of Pfs25 transcripts for quantifying gametocyte densities has drawbacks. Pfs25 mRNA levels are much lower in male compared to female gametocytes . Although it is generally assumed, due to the female bias in sex ratios (4–5:1) in natural infections, that female gametocyte levels are a good surrogate for all-sex gametocytes counts, quantification of the less abundant male gametocyte-specific Pfs230p mRNA or other male targets would provide more accurate estimates of the total number of mature gametocytes as well as sex ratios , a factor known to influence infectiousness . In addition to these targets, transcribed only in mature gametocytes, transcripts of other genes have been used for the study of sexual stage malaria: an alternative marker of female gametocytes is pfg377, the expression of which begins in the sequestered stage III ; a traditional marker for commitment to gametocytogenesis is Pfs16 [12, 39], which is present in all gametocyte stages. Another early marker of commitment is PfGEXP5, the expression of which starts in ring-stage parasites committed to gametocytogenesis, is a potentially useful tool to be used in studies aiming to identify factors influencing commitment to sexual stage . For all of these targets, transcript detection is possible both by the more convenient QT-NASBA and by the more laborious but more reproducible qRT-PCR. Design of intron-spanning primers may close the gap between the attractiveness of both methods and retain the advantage of the more reproducible qRT-PCR .
Estimating the contribution of sub-microscopic infections to malaria transmission is a priority in the malaria elimination era. Molecular methods that can quantify gametocytes levels in microscopically undetectable infections come with a promise to improve understanding of malaria epidemiology . qRT-PCR of Pfs25 mRNA is more reproducible compared to QT-NASBA. Gametocyte densities estimated by both methods, and by microscopy, correlate well with infectiousness in untreated individuals. Although this analysis of field samples suggest that these assays are suitable to quantify patent gametocyte levels, their use to understand the detectability and infectivity of sub-microscopic gametocyte densities is even more important. Future studies thus need to be designed to include individuals with sub-patent gametocyte levels to confirm the sensitivity, precision and relative merits of these assays. Precise quantification of gametocytes at very low levels is of relevance to determine if a threshold density is associated with infectivity provided molecular targets are informative of gametocyte viability and not merely of density . Importantly, neither of the two molecular methods is routinely used in areas where field studies are undertaken and both assays are highly dependent on sample RNA quality.
Okell LC, Ghani AC, Lyons E, Drakeley CJ. Submicroscopic infection in Plasmodium falciparum-endemic populations: a systematic review and meta-analysis. J Infect Dis. 2009;200:1509–17.
Okell LC, Bousema T, Griffin JT, Ouedraogo AL, Ghani AC, Drakeley CJ. Factors determining the occurrence of submicroscopic malaria infections and their relevance for control. Nat Commun. 2012;3:1237.
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.
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.
Ouedraogo AL, Goncalves BP, Gneme A, Wenger EA, Guelbeogo MW, Ouedraogo A, et al. Dynamics of the human infectious reservoir for malaria determined by mosquito feeding assays and ultrasensitive malaria diagnosis in Burkina Faso. J Infect Dis. 2016;213:90–9.
Wampfler R, Mwingira F, Javati S, Robinson L, Betuela I, Siba P, et al. Strategies for detection of Plasmodium species gametocytes. PLoS ONE. 2013;8:e76316.
Muirhead-Thomson RC. Low gametocyte thresholds of infection of Anopheles with Plasmodium falciparum; a significant factor in malaria epidemiology. Br Med J. 1954;1:68–70.
Wu L, van den Hoogen LL, Slater H, Walker PG, Ghani AC, Drakeley CJ, et al. Comparison of diagnostics for the detection of asymptomatic Plasmodium falciparum infections to inform control and elimination strategies. Nature. 2015;528:S86–93.
Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.
Lasonder E, Rijpma SR, van Schaijk BC, Hoeijmakers WA, Kensche PR, Gresnigt MS, et al. Integrated transcriptomic and proteomic analyses of P. falciparum gametocytes: molecular insight into sex-specific processes and translational repression. Nucleic Acids Res. 2016;44(13):6087–101.
Schneider P, Reece SE, van Schaijk BC, Bousema T, Lanke KH, Meaden CS, et al. Quantification of female and male Plasmodium falciparum gametocytes by reverse transcriptase quantitative PCR. Mol Biochem Parasitol. 2015;199:29–33.
Schneider P, Schoone G, Schallig H, Verhage D, Telgt D, Eling W, et al. Quantification of Plasmodium falciparum gametocytes in differential stages of development by quantitative nucleic acid sequence-based amplification. Mol Biochem Parasitol. 2004;137:35–41.
Schneider P, Wolters L, Schoone G, Schallig H, Sillekens P, Hermsen R, et al. Real-time nucleic acid sequence-based amplification is more convenient than real-time PCR for quantification of Plasmodium falciparum. J Clin Microbiol. 2005;43:402–5.
Schneider P. Submicroscopic Plasmodium falciparum gametocytaemia and the contribution to malaria transmission. PhD thesis, Nijmegen: Radboud University Department of Medical Microbiology; 2006.
Ponnudurai T, Meuwissen JH, Leeuwenberg AD, Verhave JP, Lensen AH. The production of mature gametocytes of Plasmodium falciparum in continuous cultures of different isolates infective to mosquitoes. Trans R Soc Trop Med Hyg. 1982;76:242–50.
Teirlinck AC, Roestenberg M, van de Vegte-Bolmer M, Scholzen A, Heinrichs MJ, Siebelink-Stoter R, et al. NF135.C10: a new Plasmodium falciparum clone for controlled human malaria infections. J Infect Dis. 2013;207:656–60.
Eldering M, Morlais I, van Gemert GJ, van de Vegte-Bolmer M, Graumans W, Siebelink-Stoter R, et al. Variation in susceptibility of African Plasmodium falciparum malaria parasites to TEP1 mediated killing in Anopheles gambiae mosquitoes. Sci Rep. 2016;6:20440.
Kariuki MM, Kiaira JK, Mulaa FK, Mwangi JK, Wasunna MK, Martin SK. Plasmodium falciparum: purification of the various gametocyte developmental stages from in vitro-cultivated parasites. Am J Trop Med Hyg. 1998;59:505–8.
Boom R, Sol CJ, Salimans MM, Jansen CL, Wertheim-van Dillen PM, van der Noordaa J. Rapid and simple method for purification of nucleic acids. J Clin Microbiol. 1990;28:495–503.
Dicko A, Brown JM, Diawara H, Baber I, Mahamar A, Soumare HM, et al. Primaquine to reduce transmission of Plasmodium falciparum malaria in Mali: a single-blind, dose-ranging, adaptive randomised phase 2 trial. Lancet Infect Dis. 2016;16(6):674–84.
Goncalves BP, Tiono AB, Ouedraogo A, Guelbeogo WM, Bradley J, Nebie I, et al. Single low dose primaquine to reduce gametocyte carriage and Plasmodium falciparum transmission after artemether-lumefantrine in children with asymptomatic infection: a randomised, double-blind, placebo-controlled trial. BMC Med. 2016;14:40.
Ouédraogo AL, Guelbéogo WM, Cohuet A, Morlais I, King JG, Gonçalves BP, Bastiaens GJ, Vaanhold M, Sattabongkot J, Wu Y, Coulibaly M. A protocol for membrane feeding assays to determine the infectiousness of P. falciparum naturally infected individuals to Anopheles gambiae. Malar World J. 2013;4(16):1–4.
Babiker HA, Abdel-Wahab A, Ahmed S, Suleiman S, Ranford-Cartwright L, Carter R, et al. Detection of low level Plasmodium falciparum gametocytes using reverse transcriptase polymerase chain reaction. Mol Biochem Parasitol. 1999;99:143–8.
Nakagawa S, Schielzeth H. A general and simple method for obtaining R2 from generalized linear mixed-effects models. Methods Ecol Evol. 2013;4:133–42.
Stone WJ, Dantzler KW, Nilsson SK, Drakeley CJ, Marti M, Bousema T, et al. Naturally acquired immunity to sexual stage P. falciparum parasites. Parasitology. 2016;143:187–98.
Gouagna LC, Bancone G, Yao F, Yameogo B, Dabire KR, Costantini C, et al. Genetic variation in human HBB is associated with Plasmodium falciparum transmission. Nat Genet. 2010;42:328–31.
Bustin S, Dhillon HS, Kirvell S, Greenwood C, Parker M, Shipley GL, et al. Variability of the reverse transcription step: practical implications. Clin Chem. 2015;61:202–12.
Akkermans ADL, Elsas JDv, Bruijn FJd. Molecular microbial ecology manual. Dordrecht: Boston Kluwer Academic Publishers; 1995.
Mahony JB, Song X, Chong S, Faught M, Salonga T, Kapala J. Evaluation of the NucliSens Basic Kit for detection of Chlamydia trachomatis and Neisseria gonorrhoeae in genital tract specimens using nucleic acid sequence-based amplification of 16S rRNA. J Clin Microbiol. 2001;39:1429–35.
Bartlett JW, Frost C. Reliability, repeatability and reproducibility: analysis of measurement errors in continuous variables. Ultrasound Obstet Gynecol. 2008;31:466–75.
van der Meide W, Guerra J, Schoone G, Farenhorst M, Coelho L, Faber W, et al. Comparison between quantitative nucleic acid sequence-based amplification, real-time reverse transcriptase PCR, and real-time PCR for quantification of Leishmania parasites. J Clin Microbiol. 2008;46:73–8.
Pritsch M, Wieser A, Soederstroem V, Poluda D, Eshetu T, Hoelscher M, et al. Stability of gametocyte-specific Pfs25-mRNA in dried blood spots on filter paper subjected to different storage conditions. Malar J. 2012;11:138.
Waltmann A, Karl S, Chiu C, Mueller I. Limited Degradation of the Plasmodium falciparum Gametocyte Marker pfs25 mRNA exposed to tropical temperatures: considerations for malaria transmission field studies. Am J Trop Med Hyg. 2016;94:886–9.
Duong V, Lambrechts L, Paul RE, Ly S, Lay RS, Long KC, et al. Asymptomatic humans transmit dengue virus to mosquitoes. Proc Natl Acad Sci USA. 2015;112:14688–93.
Bousema T, Churcher TS, Morlais I, Dinglasan RR. Can field-based mosquito feeding assays be used for evaluating transmission-blocking interventions? Trends Parasitol. 2013;29:53–9.
Bousema T, Drakeley C. Epidemiology and infectivity of Plasmodium falciparum and Plasmodium vivax gametocytes in relation to malaria control and elimination. Clin Microbiol Rev. 2011;24:377–410.
Robert V, Read AF, Essong J, Tchuinkam T, Mulder B, Verhave JP, et al. Effect of gametocyte sex ratio on infectivity of Plasmodium falciparum to Anopheles gambiae. Trans R Soc Trop Med Hyg. 1996;90:621–4.
Wampfler R, Timinao L, Beck HP, Soulama I, Tiono AB, Siba P, et al. Novel genotyping tools for investigating transmission dynamics of Plasmodium falciparum. J Infect Dis. 2014;210:1188–97.
Moelans II, Meis JF, Kocken C, Konings RN, Schoenmakers JG. A novel protein antigen of the malaria parasite Plasmodium falciparum, located on the surface of gametes and sporozoites. Mol Biochem Parasitol. 1991;45:193–204.
Tiburcio M, Dixon MW, Looker O, Younis SY, Tilley L, Alano P. Specific expression and export of the Plasmodium falciparum Gametocyte EXported Protein-5 marks the gametocyte ring stage. Malar J. 2015;14:334.
Joice R, Narasimhan V, Montgomery J, Sidhu AB, Oh K, Meyer E, et al. Inferring developmental stage composition from gene expression in human malaria. PLoS Comput Biol. 2013;9:e1003392.
Babiker HA, Schneider P, Reece SE. Gametocytes: insights gained during a decade of molecular monitoring. Trends Parasitol. 2008;24:525–30.
White NJ, Ashley EA, Recht J, Delves MJ, Ruecker A, Smithuis FM, et al. Assessment of therapeutic responses to gametocytocidal drugs in Plasmodium falciparum malaria. Malar J. 2014;13:483.
HEP and BPG contributed equally to the completion of this work, taking part in interpretation and analysis of data as well as writing. AD, IC, JoelleB, TB, and RG contributed to the study design for the trial in Mali. AD and SFT oversaw the data collection in Mali. HD, AM, HMS, and IB contributed to the conduct of the study and the data collection in Mali. BPG, IN, ABT, SBS, TB, and CD took part in coordinating the trial in Burkina Faso. KL completed the laboratory work, taking part in the development of methodologies. JohnB contributed to the analysis and interpretation of data. RB provided the facilities for the laboratory work and took part in useful discussions on data interpretation. IF provided critical assessment and input on the content of this work during completion process. TB contributed to and supervised interpretation and analysis of data, as well as writing. All authors read and approved the final manuscript.
The authors acknowledge Karin Koolen, Marga van de Vegte-Bolmer, Wouter Graumans, Rianne Siebelink-Stoter and Sanna Rijpma for providing in vitro produced P. falciparum gametocytes, and Tom Churcher for his comments on a previous version of this manuscript.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
The study in Mali was approved by the Ethics Committee of the Malaria Research and Training Centre Faculty of Medicine, Pharmacy and Dentistry of the University of Science, Techniques and Technologies of Bamako, and the Committee on Human Research at the University of California, San Francisco (UCSF) and was registered with ClinicalTrials.gov, Number NCT01743820 . The study in Burkina Faso was approved by the London School of Hygiene and Tropical Medicine ethics committee, and the Comité d’Ethique pour la Recherche en Santé (Ministère de la Santé du Burkina Faso) and was registered with ClinicalTrials.gov, Number NCT01935882 .
This work was supported by the Bill and Melinda Gates Foundation (AFIRM OPP1034789); TB is supported by a VIDI fellowship from The Netherlands Organization for Scientific Research (NWO; Project Number 016.158.306).
Helmi Pett and Bronner P. Gonçalves contributed equally to this work
Additional file 1: Figure S1. Estimated gametocyte nucleic acid quantities in each reaction step based on highest dilution of gametocyte standard curve (103 gametocytes/μL). For dilution series used to assign gametocyte levels to field samples, the RNA stabilising buffer was L6 Buffer (450 μL) for the Mali study and RNAprotect (250 μL) for the Burkina Faso study.
Additional file 2: Figure S2. Gametocyte levels by infectiousness. Individuals who infected less than the median proportion of infected mosquitoes in positive assays (12.4%) were considered to have low infectiousness. Three samples with qRT-PCR levels lower than 1 gametocyte/μL had their values set to 1 gametocyte/μL so they could be represented in this graph.
Additional file 3: Figure S3. Gametocyte levels and mosquito infection rates in the Burkina Faso (A) and Mali (B) studies. Individuals who infected at least one mosquito are presented in this graph. In (A), data from the Burkina Faso transmission study is presented; light blue circles with donut hole represent samples with assigned gametocyte levels by QT-NASBA below 0.02 gametocytes/μL. Red circles with donut hole represent 5 samples for which the Pfs25 amplicon was not detected by qRT-PCR; in this graph, these qRT-PCR negative samples were assigned the value of 0.02 gametocytes/μL. One sample had levels quantified by qRT-PCR below 0.02 gametocytes/μL. In (B), similar data from the Mali study is presented using the same x-axis range for comparison; black circles represent gametocyte levels estimated by microscopy.
About this article
Cite this article
Pett, H., Gonçalves, B.P., Dicko, A. et al. Comparison of molecular quantification of Plasmodium falciparum gametocytes by Pfs25 qRT-PCR and QT-NASBA in relation to mosquito infectivity. Malar J 15, 539 (2016). https://doi.org/10.1186/s12936-016-1584-z