- Open Access
Detection of Plasmodium falciparum male and female gametocytes and determination of parasite sex ratio in human endemic populations by novel, cheap and robust RTqPCR assays
Malaria Journal volume 16, Article number: 468 (2017)
The presence of Plasmodium falciparum gametocytes in peripheral blood is essential for human to mosquito parasite transmission. The detection of submicroscopic infections with gametocytes and the estimation of the gametocyte sex ratio are crucial to assess the human host potential ability to infect mosquitoes and transmit malaria parasites.
Aim and objectives
The aim of this work was to develop sensitive and cheap Real Time qPCR assays for large-scale epidemiological surveys, based on detection and amplification of gametocyte sex specific transcripts selected from the literature: the female-specific pfs25 and pf glycerol kinase (pfGK) and the male-specific pfs230p and pf13 transcripts.
RTqPCR assays were used to test the gametocyte- and sex-specific expression of the target genes using asexual stages of the gametocyteless parasite clone F12 and FACS purified male and female gametocytes of the PfDynGFP/P47mCherry line. Assays were performed on 50 blood samples collected during an epidemiological survey in the Soumousso village, Burkina Faso, West-Africa, and amplification of the human housekeeping gene 18S rRNA was employed to normalize RNA sample variability.
SYBR Green assays were developed that showed higher sensitivity compared to Taqman assays at a reduced cost. RTqPCR results confirmed that expression of pfs25 and pfs230p are female and male-specific, respectively, and introduced two novel markers, the female-specific pfGK and the male-specific pf13. A formula was derived to calculate the ratio of male to female gametocytes based on the ratio of male to female transcript copy number. Use of these assays in the field samples showed, as expected, a higher sensitivity of RTqPCR compared to microscopy. Importantly, similar values of gametocyte sex-ratio were obtained in the field samples based on the four different target combinations.
Novel, sensitive, cheap and robust molecular assays were developed for the detection and quantification of female and male P. falciparum gametocytes. In particular, the RTqPCR assays based on the female-specific pfs25 and the newly described male gametocyte-specific pf13 transcripts, including normalization by the human 18S, reliably assess presence and abundance of female and male gametocytes and enable to determine their sex-ratio in human subjects in endemic areas.
Despite a remarkable decrease in mortality and morbidity over the last 10 years thanks to the successful application of control measures, malaria still imposes a major public health burden: in 2015 malaria caused 212 million clinical cases and 429,000 deaths worldwide, 90% of which occurred in sub-Saharan Africa, where Plasmodium falciparum is the most prevalent species .
Transmission of malaria parasites from the human host to the mosquito vector requires the presence of mature gametocytes in the peripheral blood. Mature (stage V) gametocytes are morphologically and physiologically distinguishable as males or females and circulate in peripheral blood for an average of 4.6–6.5 days [2,3,4]. Parasite transmission to the vector is a very efficient process as it has been shown that mosquitoes can be infected even when feeding on subjects carrying sub-microscopic gametocyte densities, corresponding to < 4 gametocytes/μl [5, 6]. As the infection to the mosquito requires the mating of two gametes of opposite sex, the gametocyte sex ratio, i.e. the proportion of gametocytes that are male, is a critical parameter to evaluate/predict infectiousness. In malaria parasites, gametocyte sex ratio is female biased, can vary during the course of individual infections and is affected by parasite density and transmission level [7, 8].
Highly sensitive molecular tools relying on mRNA amplification, such as Reverse Transcriptase quantitative PCR (RTqPCR) and Quantitative Nucleic Acid Sequence Based Amplification (qNASBA), have been developed to detect gametocytes in peripheral blood targeting sexual stage-specific transcripts [9, 10]. The higher sensitivity of these assays compared to microscopy has been consistently described in the last years [11,12,13,14].
However, existing assays pose important cost constraints for use in the field as they are based on quite expensive qPCR Taqman and qNASBA technologies. Furthermore, they do not include an internal control gene, making them susceptible to variation in RNA quantity and quality among samples. Finally, the significant variability between reports on the absolute number of gametocyte-specific transcript copies corresponding to one gametocyte prevents to consistently translate the assay readouts into actual numbers of gametocytes.
The aim of this work was to develop sensitive, cheap and robust RTqPCR assays based on amplification of sex specific transcripts for use in large-scale epidemiological surveys. Based on previously published assays  and on a recent transcriptomic analysis of male and female gametocytes , four target genes were selected: the female-specific pfs25 and pfGK, and the male-specific pfs230p and pf13 transcripts. SYBR Green assays were developed for the detection of the above transcripts that showed higher sensitivity and reduced costs compared to existing Taqman assays, and introduced amplification of a reference gene (the human 18S rRNA) to normalize transcript copy number in field samples. Finally, a protocol to derive gametocyte sex ratio from the ratio of the male to female transcript copy number was described and tested on a panel of blood samples collected during an epidemiological survey in the Soumousso village, Burkina Faso, West-Africa.
The P. falciparum lines 3D7A , DynGFP/P47mCherry  and F12  were cultured in human 0+ erythrocytes at 5% haematocrit under 5% CO2, 2% O2, 93% N2 . Cultures were grown in RPMI 1640 medium (Gibco) supplemented with 25 mM Hepes, 50 µg/ml hypoxanthine, 0.25 mM NaHCO3, 50 µg/ml gentamicin sulfate, 10% pooled heat-inactivated O+ human serum. The PfDynGFP/P47mCherry line was cultured under selection of 5 µg/ml blasticidin.
Gametocyte culture and purification
Parasites from the 3D7A and the DynGFP/P47mCherry lines were induced to produce gametocytes by parasite overgrowth. After a 5-day treatment with 50 mM N-acetyl-glucosamine (NAG) to eliminate residual asexual parasites, stage III/IV gametocytes were partially purified from uninfected erythrocytes on MACS Separation Columns CS (Miltenyi Biotec). Gametocytes were then put back in culture and allowed to mature to stage V over the following 7 days. Aliquots of 5 × 105 mature gametocytes were centrifuged and pellets were stored at − 80 °C for further analyses.
Isolation of trophozoites and ring stages from the F12 line
Asynchronous parasites from the P. falciparum F12 line were grown to high parasitaemia (8%). Trophozoites were purified from uninfected erythrocytes and from ring stage parasites over a 60% Percoll cushion , counted in a haemocytometer chamber and aliquots were centrifuged and resuspended in 50 µl of medium and then in 450 µl of RNAlater (Applied Biosystem). Other aliquots of purified trophozoites were further grown at a low haematocrit (2%) for additional 26 h to obtain synchronous, newly invaded ring stage parasites. Aliquots of ring stage parasites were centrifuged and resuspended in 50 µl of medium and then in 450 µl of RNAlater.
Culture and purification of male and female gametocytes of the PfDynGFP/P47mCherry line
Parasites from the PfDynGFP/P47mCherry line were induced to produce gametocytes as above and stage III/IV gametocytes were partially purified from uninfected erythrocytes on MACS Separation Columns CS (Miltenyi Biotec). Gametocytes were resuspended (10 × 106/ml) in 1× PBS and male and female gametocytes were separated by sorting at room temperature using a BD FACSAria flow cytometer (Becton–Dickinson). First, gametocytes were separated from uninfected erythrocytes using forward and sideward scatter and then the GFP+ (male) gametocytes were separated from the GFP− (female) gametocytes. Aliquots of the sorted samples were analysed by UV fluorescence microscopy and Giemsa stained to determine the purity of the gametocyte populations. The GFP+ and GFP− gametocytes were put back in separate cultures and were allowed to mature to stage V over the following 7 days. Mature gametocytes were counted in a haemocytometer chamber and aliquots of male and female gametocytes were centrifuged and resuspended in 50 µl medium and then in 450 µl of RNAlater.
Epidemiological survey and study sample
An epidemiological survey was conducted during the peak of the 2013 malaria high transmission season (July–September) in the Soumousso village (11°00′46″N, 4°02′45″W), Burkina Faso, West Africa. Whole blood samples were collected from 500 subjects aged from 2 to 20 years of both sexes in order to investigate the P. falciparum infection reservoir in the population and how it varies with host factors (manuscript in preparation). Blood was collected by finger-prick in a 1.5 ml tube containing EDTA and 50 µl were immediately added to 450 µl of RNALater solution and stored at − 20 °C. Venous blood from Italian healthy donors (N = 3) were collected and a pool was used as negative control.
A representative sample subset from 50 subjects, belonging to different age groups and of both sexes, was selected from the original collection and used for the present work using cluster random sampling according to microscope positivity for P. falciparum gametocytes. Random sampling was used to select 14 (1/3) samples from subjects microscopy-negative for P. falciparum gametocytes and 36 (2/3) samples from subjects microscopy positive for P. falciparum gametocytes.
Thick and thin blood smears were prepared according to the WHO guidelines for the microscopic diagnosis of malaria , for the detection and quantification of P. falciparum gametocytes density by light microscopy (LM). Gametocyte density was estimated per microlitre of blood, by reading 100 microscopic fields (~ 20 leukocytes per field at × 1000 corresponding to ~ 0.25 μl of blood) of the thick blood smear. Microscopic readings were performed in duplicate by two experienced microscopists and the mean value of the reading was used in the data analysis.
RNA extraction and cDNA synthesis
500 μl blood/RNA later (50 μl blood, 450 μl RNALater) were centrifuged and resuspended in 250 μl of DEPC-treated water. RNA extraction was then performed using TRIzol (Life Technologies) reagent following the manufacturer’s protocol and the RNA pellet was finally resuspended in 20 μl DEPC-treated water. Quality and concentration of isolated RNA was evaluated using the microplate spectrophotometer Take3 (Synergy HT, BIOTEK). RNA samples were treated with TURBO™ DNase (Ambion) to remove contaminant genomic DNA and cDNA synthesis was performed using the High Capacity cDNA Reverse Transcription kit (Life Technologies), according to the manufacturer’s instruction. The efficacy of DNAse treatment was tested for each sample by RT minus PCR.
SYBR Green assay primers design
Primers (pfs25, pfGK, pf13, pf230p, Table 1) for SYBR Green assays were designed using the Primer Express software (Applied Biosystems) following recommended guidelines for qPCR primer design. Primers do not span introns, and gDNA contamination was avoided as described above. In order to prevent non-P. falciparum-specific amplification, primers were checked for homology by performing BLASTn analysis against genomes of Plasmodium species infecting humans other than P. falciparum (PlasmoDB, Orthology and Synteny tool), as well as against the human genome (NCBI). To experimentally check the specificity of each primer pair a melting curve analysis was added at the end of the PCR cycles. Each primer pair produced a single amplicon and showed a clean amplification as revealed by single sharp peaks (Additional file 1: Figure S1).
Standard curves of target constructs
Genomic DNA was extracted from a culture of P. falciparum using the QIAamp DNA Mini Kit (Qiagen) following manufacturer’s instruction. DNA amplicons of each P. falciparum gametocyte target (pfs25, pfGK, pf13, pf230p) and of the human reference gene (18S rRNA) were produced by end point PCR, using specific primers tailed with T7 promoter sequences to enable in vitro RNA transcription (Additional file 1: Table S2). Total RNA was therefore transcribed using the T7 RNA Polymerase Kit (Promega). After DNAse treatment, the concentration of RNA was determined by spectrophotometer analysis and copy numbers were estimated using ENDMEMO tool (http://www.endmemo.com/bio/dnacopynum.php). Finally, cDNA was prepared as described above and tenfold dilutions were used to obtain standard curves using three replicates per dilution. The amplification efficiency (E) of the qPCR assays is estimated on the basis of the equation E = (10 − 1/slope − 1) × 100. R2 is the coefficient of correlation obtained for the standard curve and should be > 0.99. The Limit of Quantification (LOQ) is the lowest concentration of construct (copies/µl) that can be quantified by a given qPCR assay.
The qPCR reactions were performed using a StepOnePlus Real-Time PCR system (Applied Biosystems). Each SYBR Green qPCR reaction (final volume of 20 μl) contained 10 μl of 2× Master mix SYBR Green (Applied Biosystems), 225 nM of Forward and Reverse primers and 2 μl of cDNA. The thermal cycle conditions were as follows: 10 min at 95 °C, followed by 40 cycles at 95 °C for 15 s, 60 °C for 1 min. A melting curve analysis was added at the end of each reaction, with fluorescence detection every 0.2 °C. pfs25 Taqman protocol was performed following the protocol described by Schneider et al. . All samples were analyzed in duplicate. In each plate, a negative control (uninfected sample), a negative amplification control and a standard curve for each target were included. A Ct value of < 40 was used as threshold of positivity. Absolute quantification protocol was applied to evaluate P. falciparum target genes copy numbers, using the standard curves of the target constructs previously described. Relative quantification was used to assess male and female gametocyte numbers, using standard curves generated from serial dilutions of male and female gametocytes isolated from cultures as described above. P. falciparum target genes copy numbers extrapolated for field samples were further normalised using the endogenous control gene by expressing them as a ratio to the to the concentration of the h18S transcript (ng/μl) measured in the same qPCR experiment.
Statistical analyses and graphs
Statistical analysis and graphics were performed using GraphPad Prism 6.0 (GraphPad Software, Inc. La Jolla, CA). Nonlinear standard curve interpolations were applied to fit data from RTqPCR experiments: Cq values from marker’s absolute quantifications were correlated to Cq values from purified female or male gametocytes relative quantifications. Nonlinear regression analysis was also employed to fit data from field samples in a straight line in log–log graphs (ratio of Male/Female gametocytes on Y-axis, ratio of copies of Male Marker/Female Marker on X-axis). Slope, intercept and R2 parameters were calculated for each equation. Resulting equation formulas were used to interpolate ratios of Male/Female gametocytes from RTqPCR values (copies of Male/Female Markers) in human blood samples using all four primers combinations.
Absolute (counts) and relative (percentages) frequencies of gametocyte-positive samples in the field collection were calculated and showed in bar plots using Microsoft Office Excel 2007.
Linear correlation between transcript copy number/μl and gametocytes/μl in field samples was assessed by correlation coefficient (r2) and Spearman test (rho, p value) using STATA v.10 (StataCorp. 2007. Stata Statistical Software: Release 10. College Station, TX: StataCorp LP.).
Results and discussion
Selection of target and reference genes
The P. falciparum pfs25 (PF3D7_1031000) and pf230p (PF3D7_0209000) genes were selected respectively as markers of female and male gametocytes in RTqPCR assays according to Schneider et al. . The P. falciparum pfGK gene (PF3D7_1351600) was initially selected to detect and quantify both male and female gametocytes, based on the evidence that this gene is both gametocyte-specific and expressed at very high levels [22, 23]. Since the work by Lasonder et al. showed that pfGK is a female-specific transcript (Additional file 1: Table S1), as it turned out to be the case for virtually all highly expressed gametocyte-specific transcripts, pfGK was used as a further female-specific assay. As the analysis of Lasonder et al.  showed that pf230p mRNA is not a top ranking male-specific transcript, a new male gametocyte-specific marker was developed using gene pf13 (PF3D7_1311100), as in that analysis this mRNA ranked first as male gametocyte-specificity, with the lowest female to male expression ratio (Additional file 1: Table S1).
To overcome the inherent sample variability in RNA quantity and quality expected in field samples, a human reference gene was selected to normalize gametocyte- and sex-specific target quantification. To choose a reference gene evenly expressed across the samples, qPCR was analysed for four human genes: 18S rRNA , B2M (β-2-microglobulin) , ACTB (β-actin) , HPRT (hypoxanthine phosphoribosyl-transferase) . The four genes were tested in 20 of the 50 cDNA samples collected in the epidemiological survey in the Soumousso village, from males and females of different age groups. According to the uniformity of expression 18S rRNA has been chosen as reference gene.
Analytical validation of the assays
Specific parameters (melting temperatures, standard curve slope and efficiency) relative to PCR amplification for each primer pair were evaluated and resulted satisfactory for downstream application (Additional file 1: Figure S1 and S2, Table S3). The LOQ, i.e. the lowest concentration of construct (copies/reaction) that can be quantified with acceptable precision and accuracy, is equal to 100 copies/μl for pfs25-Taqman, 10 copies/μl for pfs25-SYBR Green, 100 copies/μl for pfGK, 73.6 copies/μl for pf13, and 100 copies/μl for pf230p (Additional file 1: Table S2).
Stage and sex specific expression of target genes
Stage-specific expression of the selected target genes was assessed performing RTqPCR assays on asexual ring stages and trophozoites of the F12 P. falciparum strain and from gametocytes of the PfDynGFP/P47mCherry strain. The F12 line was used because it is unable to produce morphologically recognisable gametocytes [18, 28] as it is mutated in the AP2-G transcription factor, a master regulator of gametocyte differentiation ; these parasites were used as source of asexual stages to minimize the possibility of amplifying sexual stage transcripts from the low level of gametocyte contamination virtually unavoidable in asexual stage samples prepared from wild type parasite lines.
The pfs25 and pfGK female markers showed a high level of expression in the gametocyte stage, a very low level of expression in trophozoites and an undetectable level of expression in ring stages (Fig. 1). The pf13 and pf230p male markers showed comparable expression levels in gametocytes and in trophozoites, but no detectable level of expression in ring stage parasites (Fig. 1). These results indicate that the expression of female gametocyte markers pfs25 and pfGK is highly specific for the sexual stages. In contrast, the male gametocyte markers pf13 and pfs230p could be amplified also from asexual trophozoites. This indicates that the latter assays are not suited to quantify male gametocytes in in vitro cultures of asynchronous sexual and asexual blood stages. However, importantly, these assays can be confidently used with a high specificity for the sexual stages in field blood samples, as ring stages are the only asexual parasites found in circulation in peripheral blood.
The sex-specific expression of the target genes was assessed by RTqPCR on samples of purified female and male gametocytes, which were obtained by FACS sorting of GFP+ fluorescent male gametocytes and GFP− female gametocytes of the P. falciparum line PfDynGFP/P47mCherry . Sorted sub-populations of female and male gametocytes were counted and used as templates to evaluate the expression intensity of each target. Two independent sorting experiments were considered. Differences in copy number of each target gene between the two sorting experiments were observed (Additional file 1: Figure S3), likely due to variability accumulated between samples in the experimental steps from cell sorting to cDNA production. However, when transcript copy numbers were used to calculate for each target the fold difference of expression (FD) in female gametocytes compared to males, the female/male FD values were consistent across biological and technical replicates (Fig. 2). These results therefore confirmed the expected sex-specificity of all target genes both at high and low gametocyte density (250 gam/reaction; 25 gam/reaction; 2.5 gam/reaction): female markers pfs25 and pfGK are strongly enriched in the female gametocytes samples (logFDpfs25 = 1.97 and logFDpfGK = 0.86), whereas male markers pf13 and pfs230p are strongly enriched in the male gametocytes samples (logFDpf13 = − 0.33 and logFDpfs230p = − 1.23).
Estimation of gametocyte sex ratio
Determination of gametocyte sex ratio (male/female) is a useful tool to evaluate Plasmodium infectivity. Here, a method is proposed to extrapolate gametocyte sex ratio in a blood sample based on the number of male and female transcript copies quantified by RTqPCR. Results of assays performed on two independent purifications of female and male gametocytes of the P. falciparum line PfDynGFP/P47mCherry were used to build a linear correlation between the ratio of male to female transcripts copy number and the ratio of male to female gametocytes. Figure 3 shows the correlation curves for each targets pair. The slope and intercept of a curve can be used to estimate the gametocyte sex ratio in field samples (see below) using the formula.
Using this formula researchers can obtain gametocyte sex-ratios by inputing their own RTqPCR data.
Performance of RTqPCR assays on blood samples collected in the field
A pilot test of the performance of the proposed protocol, including RTqPCR of the four gametocyte targets and the human reference gene, was conducted on 50 RNA samples extracted from whole peripheral blood collected during an epidemiological survey in Soumousso, Burkina Faso. The results obtained by RTqPCR assays were compared with those of microscopy, which identified gametocyte presence in 34 samples out of 50. The pfs25 assay was performed with both Taqman and SYBR Green technologies to compare performance of the two methods.
As expected, all gametocyte assays identified a higher number of gametocyte-positive samples compared to microscopic examination (Fig. 4). Among female gametocyte markers, the pfs25-SYBR Green assay identified the highest number of positive samples (N = 44, 88%), 15 out of 16 microscopy positive samples (94%), and 29 out of 34 negative samples (85%). Among male gametocyte markers, the pf13 assay identified the highest number of positive samples (N = 43), 14 out of 16 microscopy positive samples (88%), and 29 out of 34 negative samples (85%). In the majority of the above samples (n = 37, 74%) gametocytes were detected by both assays. Interestingly, in one case a microscopy-positive sample was negative for both pfs25 (SYBR Green and Taqman) and pfGK assays but positive for both the pfs230p and pf13 assays, raising the intriguing possibility that this infection contained in that moment a vast excess of male gametocytes. Similarly, in one of the two instances where the pf13 assay was negative among the microscopy positive samples, this sample was also negative for pfs230p while it was positive for both pfs25 and pfGK assays, suggesting that this infection contained only female gametocytes.
In pfs25 assays, SYBR Green resulted in a higher number of positive samples (N = 44), compared to Taqman (N = 34, of which 32 or 94% also positive by SYBR Green). The lower sensitivity of pfs25-Taqman with respect to pfs25-SYBR Green in this experiment could be at least partly explained by a poorer performance of the Taqman assay when RNA quantity/quality is low (i.e. h18S quantity less than 1 ng/μl, Additional file 1: Figure S4).
The impact of RNA quality/quantity and normalization by human 18S rRNA was further assessed by measuring the correlation between the number of gametocytes/μl as assessed by microscopy and the number of gametocyte transcript copies/μl assessed by RTqPCR, in absence and presence of normalization. For pfs25, linear correlation coefficient increased from r2 = 0.031 prior to normalization (Spearman rho = − 0.278; p = 0.315) to r2 = 0.200 after normalization (Spearman rho = 0.522; p = 0.046); similarly, for pf13, linear correlation coefficient increased from r2 = 0.028 prior to normalization (Spearman rho = 0.155; p = 0.598) to r2 = 0.356 after normalization (Spearman rho = 0.568; p = 0.034) (Fig. 5).
Finally, the formulas described above were applied to derive gametocyte sex-ratio in each field sample from RTqPCR data, using the four combination of male and female specific transcripts. Results obtained using different target combinations were very similar (Fig. 6).
Analysis of costs
In order for researchers to make a cost-informed decision as to the protocols of choice for their studies, cost per sample for each major step of the work were compiled: (i) collection of blood from survey participant in the field (including EDTA tubes, RNAlater, tips and tubes) amounted to 0.60€ per sample; (ii) obtaining cDNA from blood (RNA extraction, DNAse treatment, cDNA synthesis) amounted to 3.10€ per sample; (iii) RTqPCR of P. falciparum gametocyte targets and human reference target cost 0.55€ per assay per sample using SYBR Green technology, and 1.35€ per sample per assay using Taqman technology. The total cost per sample to obtain both pfs25 and pf13 data normalised for h18S using SYBR Green assays is 5.35€, with a reduction of 0.80€ per sample compared to using the pfs25-Taqman assay.
The aim of this work was to develop molecular, RNA-based assays to detect female and male P. falciparum gametocytes at sub-microscopic densities and to estimate gametocyte sex ratio in human blood samples. Based on the recent literature, pfs25 and pfGK were selected as female gametocyte markers and pf13 and pfs230p as male gametocyte markers. Sensitive (LOQ ≤ 100 copies/μl) RTqPCR assays were developed for all markers using SYBR Green technology, which results in a substantial reduction of costs compared to previously published protocols based on Taqman or qNASBA. The gametocyte stage specificity of expression was confirmed for all markers, as well as the female-sex specificity of expression for pfs25 and pfGK and the male-sex specificity of expression for pf13 and pfs230p. A robust correlation was observed between the ratio of male to female transcript copy number and the ratio of male to female gametocytes, and derived formulas to calculate gametocyte sex ratio in a blood sample based on RTqPCR data. The gametocyte assays, along with a human reference gene assay for data normalization, were tested in blood samples collected during a malaria epidemiological survey in the village of Soumousso, Burkina Faso. All assays were able to detect a higher number of gametocyte positive samples compared to microscopy, as expected, with those for pfs25 and pf13 as the most sensitive assays. The pfs25 SYBR Green assay was more sensitive and less affected by RNA availability compared to the pfs25 Taqman assay. Normalization with the human reference gene enabled to establish a positive correlation of transcript copy number with microscopy gametocyte counts. Gametocyte sex-ratio based on different male and female marker pairs showed consistent results in field samples. These assays can therefore be used in epidemiological studies and in clinical trials on transmission blocking interventions to detect P. falciparum sub-microscopic densities of male and female gametocytes and to estimate gametocyte sex-ratio in human infected blood. These are crucial parameters to assess the infectivity of individuals to mosquitoes and to investigate factors affecting malaria transmission at population level.
WHO. World Malaria Report 2016. Geneva: World Health Organization; 2016. http://www.who.int/malaria/publications/world-malaria-report-2016/report/en/.
Eichner M, Diebner HH, Molineaux L, Collins WE, Jeffery GM, Dietz K. Genesis, sequestration and survival of Plasmodium falciparum gametocytes: parameter estimates from fitting a model to malaria therapy data. Trans R Soc Trop Med Hyg. 2001;95:497–501.
Bousema T, Okell L, Shekalaghe S, Griffin JT, Omar S, Sawa P, et al. Revisiting the circulation time of Plasmodium falciparum gametocytes: molecular detection methods to estimate the duration of gametocyte carriage and the effect of gametocytocidal drugs. Malar J. 2010;9:136.
Bousema T, Okell L, Felger I, Drakeley C. Asymptomatic malaria infections: detectability, transmissibility and public health relevance. Nat Rev Microbiol. 2014;12:833–40.
Schneider P, Bousema JT, Gouagna LC, Otieno S, van de Vegte-Bolmer M, Omar SA, et al. Submicroscopic Plasmodium falciparum gametocyte densities frequently result in mosquito infection. Am J Trop Med Hyg. 2007;76:470–4.
Churcher TS, Bousema T, Walker M, Drakeley C, Schneider P, Ouédraogo AL, et al. Predicting mosquito infection from Plasmodium falciparum gametocyte density and estimating the reservoir of infection. Elife. 2013;2:e00626.
Reece SE, Duncan AB, West SA, Read AF. Host cell preference and variable transmission strategies in malaria parasites. Proc R Soc B Biol Sci. 2005;272:511–7.
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.
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.
Menegon M, Severini C, Sannella A, Paglia MG, Sangaré D, Abdel-Wahab A, et al. Genotyping of Plasmodium falciparum gametocytes by reverse transcriptase polymerase chain reaction. Mol Biochem Parasitol. 2000;111:153–61.
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.
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.
Wampfler R, Timinao L, Beck H-P, 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.
Pett H, Gonçalves BP, Dicko A, Nébié I, Tiono AB, Lanke K, et al. Comparison of molecular quantification of Plasmodium falciparum gametocytes by Pfs25 qRT-PCR and QT-NASBA in relation to mosquito infectivity. Malar J. 2016;15:539.
Schneider P, Reece SE, van Schaijk BCL, Bousema T, Lanke KHW, Meaden CSJ, et al. Quantification of female and male Plasmodium falciparum gametocytes by reverse transcriptase quantitative PCR. Mol Biochem Parasitol. 2015;199:29–33.
Lasonder E, Rijpma SR, van Schaijk BL, Hoeijmakers WM, 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:6087–101.
Walliker D, Quakyi IA, Wellems TE, McCutchan TF, Szarfman A, London WT, et al. Genetic analysis of the human malaria parasite Plasmodium falciparum. Science. 1987;236:1661–6.
Alano P, Roca L, Smith D, Read D, Carter R, Day K. Plasmodium falciparum: parasites defective in early stages of gametocytogenesis. Exp Parasitol. 1995;81:227–35.
Trager W, Jensen JB. Human malaria parasites in continuous culture. Science. 1976;193:673–5.
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.
WHO. Bench Aids for malaria microscopy. Geneva: World Health Organization; 2009.
Le Roch KG, Zhou Y, Blair PL, Grainger M, Moch JK, Haynes JD, et al. Discovery of gene function by expression profiling of the malaria parasite life cycle. Science. 2003;301:1503–8.
Young JA, Fivelman QL, Blair PL, de la Vega P, Le Roch KG, Zhou Y, et al. The Plasmodium falciparum sexual development transcriptome: a microarray analysis using ontology-based pattern identification. Mol Biochem Parasitol. 2005;143:67–79.
Schmittgen TD, Zakrajsek BA. Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods. 2000;46:69–81.
Shen Y, Li Y, Ye F, Wang F, Lu W, Xie X. Identification of suitable reference genes for measurement of gene expression in human cervical tissues. Anal Biochem. 2010;405:224–9.
Liu L-L, Zhao H, Ma T-F, Ge F, Chen C-S, Zhang Y-P. Identification of valid reference genes for the normalization of RT-qPCR expression studies in human breast cancer cell lines treated with and without transient transfection. PLoS ONE. 2015;10:e0117058.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002;3:research0034.
Tiburcio M, Dixon MWA, 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.
Kafsack BFC, Rovira-Graells N, Clark TG, Bancells C, Crowley VM, Campino SG, et al. A transcriptional switch underlies commitment to sexual development in malaria parasites. Nature. 2014;507:248–52.
FS, PaA, GS and VM performed experimental work. FY performed the epidemiological survey. JO, DM and VM coordinated the epidemiological survey. FS, PaA, PiA and VM designed the experimental work. FS, PaA, FL and VM analyzed the data. FS, PaA, GS, FL, PiA and VM wrote the paper. DM, PiA and VM contributed equipment and reagents. All authors read and approved the final manuscript.
The authors wish to thank P. Schneider, School of Biological Sciences, University of Edinburgh, for useful assistance in setting up the pfs25 Taqman protocol, and in providing the pfs25 construct to set up the standard curve, S. Rijpma, T. Bousema and R. Sawerwein, Radboud University Medical Centre, Nijmegen, The Netherlands for kindly providing the P. falciparum DynGFP/P47mCherry line and M. Sanchez and V. Tirelli, Istituto Superiore di Sanità, Rome, for precious technical assistance in the FACS purification of male and female gametocytes. Human 0+ erythrocytes were kindly provided by G. Girelli (Sapienza University of Rome, Italy).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Consent for publication
The manuscript does not include details, images, or videos relating to an individual study participant.
Ethics approval and consent to participate
The study (Reference Number 009-2011/CE-CM) received approval by the Ethical Committee of Centre Muraz, Ministry of Health, Burkina Faso. Study subjects or their guardians gave written informed consent to participation.
This work was supported by the European Virtual Institute for Malaria Research (EVIMalaR, European Community’s Seventh Framework Programme; Grant Agreement 242095 to D.M., Pi.A., V.M.) and by Grant OPP1040394 of the Bill & Melinda Gates Foundation (to Pi.A.). The funding bodies had no role in the design of the study, in the collection, analysis, and interpretation of data, and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1. Plasmodium falciparum genes showing the 10 lowest female to male mRNA ratio (male specific) and the 10 highest female to male mRNA ratio (female specific), from Lasonder et al. . Table S2. Primer sequences for amplification of target constructs. Table S3. Limit of Quantification and quality parameters of RTqPCR assays. Table S4. Parameters of the formula to extrapolate gametocyte sex ratio from RTqPCR data. Figure S1. Melting curves of RTqPCR assays. Figure S2. Standard Curves of RTqPCR assays. Figure S3. Expression of target genes in male and female P. falciparum gametocytes. Figure S4. Sensitivity of pfs25 Taqman and pfs25 SYBR Green assays according to RNA quality/quantity.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Santolamazza, F., Avellino, P., Siciliano, G. et al. Detection of Plasmodium falciparum male and female gametocytes and determination of parasite sex ratio in human endemic populations by novel, cheap and robust RTqPCR assays. Malar J 16, 468 (2017). https://doi.org/10.1186/s12936-017-2118-z
- Malaria transmission
- Plasmodium falciparum
- Real Time qPCR