Real-time PCR assays for detection and quantification of early P. falciparum gametocyte stages
==============================================================================================

* Amal A.H. Gadalla
* Giulia Siciliano
* Ryan Farid
* Pietro Alano
* Lisa Ranford-Cartwright
* James S McCarthy
* Joanne Thompson
* Hamza Babiker

## Abstract

**Introduction** The use of reverse transcription, quantitative qRT-PCR assays for detection and quantification of late gametocyte stages has revealed the high transmission capacity of the human malaria parasite, *Plasmodium falciparum*. A full understanding how the parasite adjusts its transmission in response to varying in-host environmental conditions during natural infections requires simultaneous quantification of early and late gametocytes. Here, we describe qRT-PCR assays that are specific for detection and quantification of early-stage gametocytes of *P. falciparum*.

**Methods** The assays are based on expression of known early gametocyte genes (*pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5*). The specificity of the qRT-PCR assays was tested using purified stage II and stage V gametocytes. These validated assays were used with qRT-PCR assays targeting late stage (*pfs25)* and all-stage (*pfs16*) gametocyte-specific transcripts to quantify gametocytes in natural *P. falciparum* infections and in a controlled human clinical infection study.

**Results** The relative expression of *pfpeg4, pfg27* and *pfge3*, but not of *pfge1* and *pfgexp5*, was significantly higher in purified stage II compared to stage V gametocytes, indicating early gametocyte specificity. In natural infections, 71.2% of individuals had both early and late gametocyte transcripts (*pfpeg4*/*pfg27* plus *pfs25*), 12.6% harboured only early gametocytes transcripts (*pfpeg4*/*pfg27)*, and 15.2% had only late gametocytes transcripts (*pfs25*). In natural infections, the limit of detection was equivalent to 190 and 390 gametocytes/mL blood for *pfpeg4* and *pfg27*, respectively. In infected volunteers, transcripts of *pfpeg4* and *pfg27* were detected shortly after the onset of blood stage infection, demonstrating the specificity of the assays.

**Conclusion** The *pfpeg4* and *pfg27* qRT-PCR assays can be used specifically to quantify circulating immature gametocytes. Quantification of early gametocytes will improve understanding of epidemiological processes that modulate *P. falciparum* transmission and enhance the evaluation of transmission blocking interventions.

Keywords
*   Malaria
*   *Plasmodium falciparum*
*   early gametocytes qPCR
*   *pfpeg4*
*   *pfg27*

## 1. Introduction

The malaria parasite *Plasmodium falciparum* continues to be a major cause of human disease despite mounting control measures. Central to this success is its ability to undergo sexual development and efficient transmission between human hosts through its mosquito vector. In *P. falciparum*, sexual development occurs when asexual blood stage parasites exit the proliferative cycle and develop into schizonts containing merozoites committed to sexual differentiation (Bruce et al., 1990). These merozoites invade erythrocytes and differentiate into the sexual forms, gametocytes, over a period of approximately 10 days, passing through five morphologically distinct stages (I to V) (Carter and Miller, 1979). A second direct path to commitment to gametocytogenesis may also occur shortly after invasion of the erythrocyte by an uncommitted merozoite (Bancells et al., 2019). While very early and late-stage gametocytes (stages I and V) are present in the peripheral circulation, immature gametocytes (stage II-IV) withdraw from peripheral circulation and sequester in tissue microvasculature (Smalley et al., 1981, Joice et al., 2014, Aguilar et al., 2014, Obaldia et al., 2018).

Stage I gametocytes are morphologically indistinguishable from pigmented asexual trophozoites, and thus cannot be identified in peripheral blood by microscopic examination. Transcriptome analyses has led to the identification of transcripts upregulated before morphological differentiation into stage II gametocytes. These include *pfmdv-1 (pfpeg3), pfgdv1, pfge1*-*3* (Silvestrini et al., 2005, Eksi et al., 2005, Eksi et al., 2012) and *pfpeg4* (Joice et al., 2013), which encodes the membrane protein ETRAMP 10.3 (Silvestrini et al., 2005, MacKellar et al., 2010, Buchholz et al., 2011). Further analyses of the transcriptomes of *P. falciparum* isolates have demonstrated that transcripts of *pfge1, pfge2* and *pfg27* cluster away from those of late stage-specific genes *pfs25, pfs28* and *pfs47*, suggesting that they are expressed exclusively in early stage gametocytes (Eksi et al., 2012). In addition, Tiburcio et al., (2015) described *pfgexp5* as the earliest gametocyte-specific gene, expressed 14 hours post-invasion. This is consistent with the detection of *pfgexp5* in blood samples taken at early time points following experimental infection of volunteers with *P. falciparum* (Farid et al., 2017).

The development of RT-PCR assays based on late sexual stage-specific mRNA transcripts has enabled epidemiological surveys and revealed the transmission potential of *P. falciparum* that drives its success in the face of continued control efforts (Babiker et al., 1999, Nassir et al., 2005, Schneider et al., 2006, Menegon et al., 2000). These surveys have identified a range of epidemiological correlates that modulate transmission, including drug treatment, anaemia, mixed species infection and multiplicity of genotypes (Bousema et al., 2004, Babiker and Schneider, 2008, Mawili-Mboumba et al., 2013, Ouédraogo et al., 2009). These epidemiological findings are consistent with the predictions of the evolutionary hypotheses that the interplay between asexual replication and transmission is altered by a changing in-host environment (Reece et al., 2008, Carter et al., 2013). However, to disentangle the effect of these factors on transmission success, there is a need for assays that can detect and quantify very early gametocytes close to the point of commitment.

Here, we describe sensitive qRT-PCR assays for quantification of early stage gametocytes, that complement existing late-stage assays, to enable the analysis of epidemiological factors that drive *P. falciparum* transmission success, and robust assessment of control strategies targeting sexual stages. The assays were validated using *in vitro* purified early and late *P. falciparum* gametocytes and their robustness was tested using samples from natural infections and from experimentally-infected volunteers.

## 2. Methods

### 2.1. Asexual and purified stage II and V gametocyte samples

*P. falciparum* clone 3D7A was cultivated as described by Trager and Jensen (1976) with minor modifications. Briefly, parasites were maintained in human type O-positive RBCs at 5% haematocrit (Hct) in RPMI 1640 medium supplemented with 25 mM HEPES (Sigma), 50 µg/ml hypoxanthine and with the addition of 10% (v/v) naturally-clotted heat-inactivated 0+ human serum (Interstate Blood Bank, Inc.). The cultures were maintained at 37°C in a standard gas mixture consisting of 3% O2, 5% CO2 and 92% N2. To obtain synchronous gametocytes, an asexual culture (0.5-1% starting parasitaemia) was grown without further dilution as described (D’Alessandro et al., 2013). Induced gametocyte cultures were supplemented with 50 mM N-acetylglucosamine (NAG; Sigma-Aldrich) to clear residual asexual parasites and obtain a virtually pure gametocyte culture (Gupta et al., 1985). To obtain the stage II gametocyte sample, two days after NAG treatment gametocytes were inspected on Giemsa-stained smears to assess maturation, and purified by centrifugation through a discontinuous Percoll gradient, as described in Silvestrini et al. (2010). Gametocytes were examined in a counting chamber for quantification and control of contamination. Stage IV and V gametocytes were not detectable. A small fraction of residual unhealthy/dead trophozoites and uninfected RBCs was detectable. To obtain the stage V gametocyte sample, gametocytes were induced as described before (D’Alessandro et al., 2013), and 12 days after NAG treatment gametocytes were purified by MACS Separation Columns CS (Miltenyi Biotec) (Ribaut et al., 2008) and examined as above. Immature sexual stages and asexual stages could not be detected. In both cases the enriched parasite preparations were pelleted and frozen in liquid nitrogen. Samples containing mixed-stage asexual parasites, purified stage II gametocytes, and purified stage V gametocytes were used to examine the stage specificity of the candidate early gametocyte genes and to establish the mathematical relationship between transcript numbers of the candidate genes and gametocyte counts (more details are given in section 2.6).

### 2.2. Natural *P. falciparum* infections

A total of 250 samples were examined from an existing collection, stored at −80°C as packed RBCs. The samples were obtained from patients with uncomplicated *P. falciparum* malaria from Asar village, eastern Sudan, with ethical clearance from the Ministry of Health, Sudan (Ali et al., 2006). The samples were employed in the current study to investigate whether the expression of early gametocyte markers is: 1) detectable in natural infections; 2) detectable specifically in a subset of samples that were *pfs16* positive (early and late gametocyte markers) and *pfs25* negative (late gametocyte marker); 3) associated with total parasitaemia.

### 2.3. Experimentally infected human volunteers

Samples from experimentally-infected human volunteers, who had participated in a previously published clinical trial (Farid et al., 2017, Pasay et al., 2016), were utilized to study the stage-specificity of the gametocyte markers *in vivo*. Ethical approval of the study was obtained from QIMR Berghofer Human Research Ethics Committee. The volunteer cohort has been previously described. In brief, it comprised malaria-naïve, healthy males and non-pregnant females, aged 18-50 years. Infections were induced using approximately 1800 viable *P. falciparum*–infected human erythrocytes of clone 3D7. Volunteers were then treated with 480 mg of piperaquine (which affects only asexual stages) on Day 8 (D08) post-infection, when parasitaemia reached a predetermined threshold of >1,000 parasites/mL. In the current study, samples were examined from two of these volunteers who showed a relatively high post-treatment parasitaemia, and high levels of *pfs25* expression at or beyond D11 post-infection (Farid et al., 2017). Ten time points were selected for each volunteer, between D07 to D09 post-infection. This represents a suitable window for interrogation for early gametocytes, as parasitaemia started to increase at D07, with limited possibility for the presence of late gametocytes in the initial inoculum as discussed (Farid et al., 2017).

### 2.4. Total parasite density quantification

Parasite density was quantified as parasites/mL of blood using absolute quantification of the *18S rRNA* gene by qPCR assay (Hermsen et al., 2001, Gadalla et al., 2016). *18S rRNA* copy numbers were converted to parasite numbers using a calibration curve from *P. falciparum* clone 3D7 parasite DNA with a range of 0.14 to 138938 parasites/μL of DNA (Nwakanma et al., 2009). The *18S rRNA* qPCR amplification efficiency was 97.7% (se 0.01%). Quantification was carried out in duplicate with 0.35 standard deviation between replicates. Quantification of parasite density in the volunteers’ samples has been described elsewhere (Farid et al., 2017).

### 2.5. DNA, RNA extraction, purification and cDNA preparation

For samples from naturally-infected patients, DNA and RNA were extracted from 100 µL and 50 µL of packed RBCs, using the QIAamp DNA mini kit (Qiagen) and SV Total RNA Isolation System (Promega, UK) respectively. For RNA purification, all samples were treated with a unified concentration of 1 unit of RQ1 RNase-Free DNase (Promega, UK) per 8 µL of RNA sample to remove any genomic DNA (gDNA) carryover. Then, purified RNA samples were checked for gDNA by*18S rRNA* qPCR assays. Pure RNA samples were converted to cDNA using the High-Capacity cDNA Reverse Transcription Kit with random primers (Thermo Fisher, UK). For the volunteer cohort, RNA was extracted from 800 µL of whole blood and cDNA was prepared as described (Farid et al., 2017).

### 2.6. qRT-PCR assays for early and late gametocytes

#### 2.6.1. Oligonucleotide design

Full gene names and accession numbers are provided in supplementary Table S1. qRT-PCR primers and TaqMan dual-labelled probes were designed within exon sequences, avoiding polymorphic regions to ensure reproducibility of the experiment for field isolates (Table S2). BLAST alignments (Ensembl.org, 2014, PlasmoDB.org, 2014) against sequences available in PlasmoDB.org (2014) showed that the primers and probe sets had 100% identity with *P. falciparum* target genes, and no high identity alignment with other human *Plasmodium* species or other *P. falciparum* genes. This was further confirmed by the absence of amplification when primers and probes were used in PCR assays using *Plasmodium vivax, Plasmodium malariae* and *Plasmodium ovale* gDNA. Probe design for *pfgexp5* was described previously by Farid et al. (2017).

#### 2.6.2. Optimization

A Taqman assay was designed and optimised for each of the putative early gametocyte-specific genes (*pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5*), for one late gametocyte-specific gene (*pfs25*), for one pan-gametocyte marker (*pfs16*), and for one reference gene (*pf40S*). Primer concentrations for the optimised qRT-PCR reactions were 0.3 µM (*pfpeg4, pfg27, pfge1* and *pfs25*), 0.6 µM (*pfgexp5*), and 1.2 µM (*pfge3*). Probe concentrations used were 0.1 µM in the *pfpeg4, pfge1, pfge3, pfgexp5* and *pfs25* qRT-PCR assays and 0.2 µM in the *pfg27* and *pf40S* assays. The final reaction contained 1X TaqMan Universal PCR Master Mix, No AmpErase UNG, ABI. The temperature profile was; 2 min at 50°C, 10 min at 95°C, and then 45 cycles of 15 sec at 95°C and 1 min at 60°C. Quantification was carried out in duplicate with 0.35 standard deviation between replicates.

#### 2.6.3. Generation of qRT-PCR standard curves

To provide enough material for qRT-PCR optimization, establishment of sensitivity and for sample quantification assays, standard curves were based on cDNA that was transcribed *in vitro* from DNA templates, hereafter named ivcDNA. DNA templates for each gene were prepared by PCR amplification of the target regions from *P. falciparum* clone 3D7 incorporating primers containing SP6 or T7 RNA polymerase promoter sequences (Table S2). Amplicons were purified from excess reagent and primers using Wizard SV Gel and PCR Clean-Up System (Promega) as described by the manufacturer, and the concentration was determined using nanodrop (NanoDrop Spectrophotometer, ND1000, Thermo Fisher Scientific). The number of molecules of DNA amplicon was calculated using the equation (ng DNA × 6.022 × 1023) / length (bp) of the DNA template × 109 x 650. Purified DNA amplicon was then transcribed to RNA using Riboprobe Combination System-SP6/T7 RNA Polymerase (Promega) as described by the manufacturer. gDNA was removed from RNA using 1 unit of DNase treatment (RQ1 DNase, Promega) per 8 µL of RNA and qPCR was carried out to confirm the purity of the RNA prior to cDNA preparation. ivcDNA was prepared by RNA reverse transcription using High-Capacity cDNA Reverse Transcription Kit, Applied Biosystems (ABI), USA as described by the manufacturer. The concentration of the ivcDNA was calculated based on the initial concentration of DNA amplicon and incorporating the dilution factors occurred during the processes above (DNA transcription, RNA purification and RNA reverse transcription). ivcDNA was then used to generate standard curves to assess qRT-PCR efficiency and sensitivity and to quantify the transcript number of the early gametocyte candidate genes in clinical samples.

#### 2.6.4. Quantification of gametocyte gene expression and gametocyte numbers

The limit of quantification (LoQ) for the qRT-PCR assays was defined as the lowest early or late gametocyte concentration (serially diluted cDNA from the purified stage II or stage V, respectively) detectable in all experiments that fell within the log-linear relationship of the qRT-PCR assay (Boyer et al., 2013). There was insufficient purified early gametocyte cDNA available to run a quantitative standard curve against all clinical samples, thus an ivcDNA standard curve was used. Detection of ivcDNA concentrations below the LoQ was possible in >50% of the standard curves. Although points below the LoQ were outside the log-linear relationship of the qRT-PCR, and could not be reliably quantified, they were considered to represent reliable detection because they fell within a clear amplification curve, in contrast to the cDNA-free sample (negative control) (Klymus et al., 2020).

The transcript number of the stage-specific gametocyte genes was quantified in clinical samples based on the log-linear mathematical relationship between qRT-PCR cycle threshold (CT) and the log10 of the concentration of the ivcDNA standard curve (supplementary Figure S1). Then, the log-linear mathematical relationship between the transcript numbers and the early gametocyte numbers in the purified stage II gametocytes culture was established (supplementary Figure S2). The number of gametocytes/mL of blood in natural infections was estimated based on the average number of transcripts obtained from known numbers of purified stage II or stage V gametocytes, assuming 100% efficiency of the processes of reverse transcription, RNA purification and RNA extraction.

### 2.7. Statistical analysis

Intra- and inter-assay coefficients of variation of the standard curves were calculated as the standard deviation of the CT values divided by their mean CT, and then multiplied by 100. The levels of the reference gene expression between different gametocyte stages were compared using Student’s T test. A Wilcoxon test was used to compare the relative expression of early and late gametocyte markers among stage II, stage V and asexual samples. Logistic regression was used to investigate the probability of detecting early gametocytes at variable levels of parasitaemia. Spearman’s correlation was used to assess the correlation between *pfpeg4* and *pfg27* expression in field samples. It was also used to assess the correlation between early and late gametocyte gene expression in field samples. R version 3.2.3 (2015-12-10) was used for statistical analysis.

## 3. Results

### 3.1. *In silico* analysis of gametocyte markers and the reference gene

To identify putative early-stage gametocyte markers, we searched published *P. falciparum* stage-specific gametocyte transcription analyses and transcriptome data (Bruce et al., 1994, Lobo et al., 1994, Furuya et al., 2005, Silvestrini et al., 2005, Farid et al., 2017, Eksi et al., 2012, PlasmoDB.org, 2014). Candidate genes were prioritised by comparing the transcript levels of fragments per kilobase of exon model per million mapped reads (FPKM) as described (Lopez-Barragan et al., 2011). The genes *pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5* were selected for evaluation in Taqman qRT-PCR assays. The gene *pfs16* was used as a pan-gametocyte stage marker, since Pfs16 protein expression starts 24-30h post-invasion of a sexually-committed merozoite and continues throughout gametocyte development (Bruce et al., 1994). *P. falciparum* transcription data available in PlasmoDB were also searched to identify a candidate reference gene, and *pfs40S* was selected as it constitutively expressed, at the same level, in asexual parasites and gametocytes and could be used to estimate the relative expression of early gametocyte candidate genes.

### 3.2. Efficiency, reproducibility and sensitivity of qRT-PCR assay

The efficiency, sensitivity and reproducibility of the early-stage gametocyte qRT-PCR assays was assessed using ivcDNA of *pfpeg4, pfg27, pfge1, pfge3* and *pfs16*, and a late-stage gametocyte culture was used to test the *pfs25* qRT-PCR assay. In 4-6 independent experiments, the average qRT-PCR efficiency for the genes ranged between 85.6% and 92.2% and the maximum intra- and inter-assay coefficient of variation (CV) was 0.9% - 6.4% and 1.9% - 6.6% (Table S3 and Figure S1). In all cases, the low CV (<10%) demonstrates the high reproducibility of the assays. The LoQ of the *pfpeg4, pfg27, pfge1, pfge3, pfs16* and *pfs25* qPCR assays was 0.14, 0.28, 0.2, 0.14, 0.68 and 0.15 gametocyte/µL of cDNA, which is equivalent to 0.56, 1.12, 0.8, 0.56, 2.72 and 0.60 gametocytes per qRT-PCR assay, respectively.

### 3.3. Validation of early markers in purified stage II and V gametocytes

To establish that the expression level of the reference *pf40S* is similar across early and late gametocytes, qRT-PCR was performed. The reference gene, *pf40S*, showed similar expression levels (t-test T= 1.25, P= 0.28) in samples containing equal numbers of purified stage II (mean CT= 38.73, SD= 0.67, in 4 technical replicates of 1 sample) and stage V gametocytes (mean CT = 37.60, SD= 1.43, in 4 technical replicates of 1 sample).

The stage specificity of the early gametocyte genes was assessed by quantifying their expression in comparison to the reference gene *pf40S* in three samples of *P. falciparum* culture: (i) mixed asexual stages (ii) purified stage II gametocytes and (iii) purified stage V gametocytes. The relative expression of *pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5* was examined in purified stage II, stage V gametocytes (4-6 technical replicates each), and in mixed asexual stages (3-4 technical replicates). Significantly higher expression of *pfpeg4* (∼24-fold) and *pfg27* (∼43-fold) was observed in stage II compared to stage V gametocytes (Wilcoxon signed-rank test, *pfpeg4*; P=7.9×10−3; *pfg27*, P=0.01), with negligible expression in asexual stages (Figure 1). *Pfge1* and *pfge3* were expressed at higher levels in stage II compared to stage V, but their relative expression in stage II was considerably lower than *pfpeg4* (125-fold and 187-fold lower, respectively) and *pfg27* (50-fold and 75-fold lower, respectively). In contrast, *pfgexp5* was highly expressed in stage II gametocytes compared to asexual stages, but at higher levels in stage V gametocytes (Wilcoxon signed-rank test, P= 0.02). As expected, a significantly higher relative expression of *pfs25* (∼11-fold) was observed in stage V compared to stage II gametocytes (Wilcoxon signed-rank test, P = 7.9×10−3, Figure 1). Therefore, *pfpeg4* and *pfg27* were selected as early gametocyte markers for further validation in clinical samples.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/04/2021.03.28.21254192/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2021/04/04/2021.03.28.21254192/F1)

Figure 1: Validation of early gametocyte markers.
Figure shows fold change (y axis) in expression of early (*pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5*) and late (*pfs25*) markers relative to the reference gene (*pf40S*) expression in parasites obtained from *in vitro* culture at different parasite developmental stages. Points represents actual data points, boxplots represent median, first and third quartiles. Error bars represent the minimum and the maximum fold change value.

### 3.4. Detection of early and late gametocytes in natural *P. falciparum* infections

Two hundred and fifty samples from patients with natural *P. falciparum* infections (Ali et al., 2006) were used to investigate the specificity of the early gametocyte assays, across a range of parasitaemias, and in the presence of late gametocytes. Of the 250 samples, 79.2% (n=198) were positive for *pfs16*, indicating the presence of gametocytes (supplementary Table S4). Of those 198 samples, 25 (12.6%) were positive for *pfpeg4* and/or *pfg27* but not for *pfs25*, indicating the presence of early gametocytes only, and 30 (15.2%) were positive for *pfs25* and not *pfpeg4* and *pfg27*, indicating the presence of only late gametocytes (Figure 2A). The majority of the samples contained mixed early and late gametocyte stages (n = 141, 71.2%). However, for 2 (1%) samples out of the 198 that were *pfs16-*positive no transcripts for either early or late gametocyte genes were detected.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/04/2021.03.28.21254192/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2021/04/04/2021.03.28.21254192/F2)

Figure 2: Early gametocytes in natural infections.
(**A**) Number of samples with different gametocyte stages in field samples (n= 198). Early gametocyte detection is based on both *pfpeg4* and *pfg27* expression. (**B**) Correlation between log10 early gametocyte densities as deteted by *pfpeg4* and *pfg27*. The vertical and horizontal dotted lines represent the limit of quantification (LoQ) for *pfpeg4* (equivalent to 190 gametocytes/ml) and *pfg27* (equivalent to 390 gametocytes/ml), respectively. The larger panel shows correlation in all samples and the smaller panel shows correlations in samples where densities were above the LoQ (represents the right upper square on the larger panel). (**C**) Correlation between log10 early (*pfpeg4* or *pfg27*) and late stage gametocytes/ ml of blood among samples with mixed stages. The figure shows that early stage gametocytes densities are independent, particularly for *pfpeg4*, to the density of late stage gametocytes. The dotted line in both panels represents the limit of quantification. Correlation line, Spearman correlation coefficient (rho) and P value were calculated for points above the limit of quantification. Points with lower density (below LoQ line) are robustly detectable however with higher uncertainty around the estimated quantities.

The *pfpeg4* and *pfg27* qRT-PCR assays detected early gametocytes in patients with a wide range of total parasitaemia, 79 to 2.6×108 parasites/mL of blood as determined by 18s rRNA qPCR. The presence of early stage gametocytes was associated with parasitaemia: for each 10-fold increase in parasite density the probability of detecting early gametocytes increased by a factor of 2.22 (95% CI= 1.63 – 3.01) for *pfpeg4* and 2.12 (95% CI= 1.67 – 2.70) for *pfg27*. Of note, >80% of samples containing early gametocytes were from patients with parasite density of >104 parasites/mL of blood Figure S3.

The *pfpeg4* and *pfg27* qRT-PCR assays detected variable densities of early gametocytes, with a range of 1 to 5.6×104 and 18.4 to 4.3×105 early gametocytes/mL of blood, respectively. However, the densities of 104 and 19 samples were below LoQ for *pfpeg4* and *pfg27*, respectively (Section 2.6.4). The LoQ of the assays was estimated as 0.14 and 0.28 early gametocyte/µL of cDNA for *pfpeg4* and *pfg27*, equivalent to 190 and 390 early gametocytes/mL of blood, respectively, assuming 100% efficiency of the RNA extraction, purification and cDNA processes. A moderate correlation was seen between the early gametocyte densities quantified by *pfpeg4* and *pfg27* qRT-PCR assays (Spearman’s correlation coefficient (rho) =0.63, P < 0.01, Figure 2B). However, where early gametocyte densities were above the LoQ, this correlation was lower (rho =0.40, P < 0.008, Figure 2B).

Among samples containing a mix of early and late gametocytes and where early gametocyte densities were above the LoQ, early gametocyte densities did not correlate (*pfpeg4*: rho = 0.06, P= 0.70) or weakly correlated (*pfg27*: rho= −0.21, P= 0.03) with densities of late stage gametocytes and detected by *pfs25* (Figure 2C).

### 3.5. Detection of early gametocytes in experimentally infected volunteers

The specificity of the *pfg27* and *pfpeg4* assays was further tested in two volunteers experimentally infected with *P. falciparum* (Farid et al., 2017, Pasay et al., 2016). The parasite density ranged between 1.3×103 – 1.6×105 (median 4.4×104) parasites/mL blood in volunteer S035 and 3.3×105 – 8.4×102 (median 5.4×104) parasites/mL blood in volunteer S042 (Figure 3) over days 7 to 10 post-infection. *pfpeg4* and *pfg27* transcripts were detectable at an early stage of the infection, day 8 and day 9 post infection (Figure 3). Early gametocyte density (as identified by the *pfpeg4* and *pfg27* assays) fluctuated at low levels, often below the LoQ and with overlapping 95% CI throughout the follow up time points (*pfpeg4* range: 106 – 850; *pfg27* range: 113 – 371 early gametocytes /mL of blood). However, late gametocytes (80 to 247 gametocytes /mL) were detected below the LoQ (using the *pfs25* assay) in 4 and 5 out of the 10 time points tested in volunteer S035 and volunteer S042, respectively. The detection of late gametocytes was associated with wide 95% CI often crossing the zero due to detection occurring in only one of the sample duplicates (Figure 3).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/04/04/2021.03.28.21254192/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2021/04/04/2021.03.28.21254192/F3)

Figure 3: Detection of early gametocytes *in vivo* in human experimental infections.
Left Y axis represents number of gametocytes/ml of blood, quantified by *pfpeg4* or *pfg27* RT-qPCR. On the right Y axis plotted total parasite density as quantified by *18S rRNA* qPCR. Horizontal grey dotted line represents limit of quantification of early gametocytes and reads to the left Y axis. Data shown is from 2 volunteers (S035 and S042) during follow up time-points between D07 and D09 post-infection. The x-axis represents day post-infection and digits after the decimal points are fractions of that day. Volunteers were treated with piperaquine at D08 (grey shaded column). Data is provided in supplementary file (Table S5).

## Discussion

The present study describes sensitive qRT-PCR assay for detection and quantification of early-stage gametocytes of *P. falciparum*. We initially examined the relative expression of five early gametocyte markers (*pfpeg4, pfg27, pfge1, pfge3* and *pfgexp5*) in purified early (stage II) and late (stage V) gametocyte and demonstrated that *pfpeg4* and *pfg27* are expressed predominantly in the early stages. In contrast, *pfgexp5* is expressed at high levels in both stage II and stage V gametocytes, and therefore is unable to distinguish between early- and late-stage gametocytes. Transcript quantification of gametocyte genes was estimated relative to that of the reference gene *pf40S*, which is expressed at the same level in stage II and V gametocytes of *P. falciparum*. Thus, variations in the amount of the starting transcripts were controlled for. The sensitivity of the *pfpeg4* and *pfg27* qRT-PCR assays were equivalent to detection thresholds of 190 and 390 gametocytes per mL of blood, respectively.

The specificity of the *pfpeg4 and pfg27* qRT-PCR assays was evident by analysis of 198 field *P. falciparum* isolates from eastern Sudan (Ali et al., 2006). The *pfpeg4* and/or *pfg27* transcripts were detected in 25 (12.6%) isolates in the absence of the *pfs25* transcripts, suggesting the presence of a cohort of circulating early gametocytes. Conversely, *pfs25* transcripts were detected in 30 (15.2%) isolates in the absence of *pfpeg4* and *pfg27* transcripts, indicating the presence of a circulating population comprising only late gametocytes (Figure 2A). Evidence of circulating early stage *P. falciparum* gametocytes in natural infections is scarce (Aguilar et al., 2014, Pelle et al., 2015), but has been documented in cases of exceedingly high parasite burdens, when asexual stages that usually sequester are also visible in circulation, as well as in splenectomized infected individuals (Bachmann et al., 2009).

Recent transcriptome analysis identified a number of early gametocyte-specific transcripts from sexually-committed ring stages (Poran et al., 2017), which were also found in peripheral blood samples in natural infections (Usui et al., 2019). The authors reported a high proportion of infections with circulating sexually-committed ring stage (76%), but the ratio of gametocytes to sexually committed rings (*in vitro*) varied dramatically, ranging from 78% to absent. This pattern is consistent with the pattern of detection of early gametocytes using *pfpeg4* and/or *pfg27* assays in our field setting in eastern Sudan. In this setting, 25 (12,6%) of 198 gametocyte-producing *P. falciparum* isolates were positive for *pfpeg4* and/or *pfg27*, but not for *pfs25*, indicating the presence of early gametocytes alone. Nonetheless, the majority of the samples (n = 141, 71.2%) contained transcripts of early (*pfpeg4* and/or *pfg27*) and late gametocyte stages (*pfs25*). These data suggest that, in addition to sexually committed ring stages, slightly older stage I gametocytes may also be found in the peripheral circulation. Further longitudinal analysis of asymptomatic infections, in the absence of re-infection in a highly seasonal endemic setting, such as that in eastern Sudan (Abdel-Wahab et al., 2002) will elucidate epidemiological factors that drive gametocytogenesis. A previous longitudinal survey of asymptomatic infection in this region demonstrated an increase in *pfs25* transcripts, indicative of mature gametocytes, following resurgence of mosquitoes at the start of the rainy season. However, no corresponding increase in the density or prevalence of total parasites or gametocytes was seen (Gadalla et al., 2016). Analysis of early markers would verify the hypothesis that *P. falciparum* may respond to environmental cues, such as mosquito biting, to modulate its transmission strategy (Reece et al., 2010).

The robustness of the *pfpeg4* and *pfg27* assays is demonstrated by the ability to detect transcripts at a wide range of parasitaemia in field samples, from very high (2.6×108 parasites/mL blood), always associated with clinical presentation, to sub-microscopic levels (79 parasites/mL blood), often present in asymptomatic infections (Rogier et al., 1996, Roucher et al., 2012, Gadalla et al., 2016). The detectability of *pfpeg4* and/or *pfg27* transcripts is strongly linked to the level of total parasitaemia, in line with the performance of other qRT-PCR assays for gametocyte quantification; *pfgexp5* (Farid et al., 2017), *pfs230p* (Schneider et al., 2015) and the *pfs25* (Tadesse et al., 2017).

The specificity of the *pfpeg4* and *pfg27* qRT-PCR assays was further tested by analysis of samples from volunteers experimentally infected with *P. falciparum*, taken at D07 - D09 post-infection, when only early but not late gametocyte stages are present. The volume of the initial inoculum was 1,800 parasites; equivalent to approximately 1/10−14 parasite/ mL, precluding the detection of gametocytes at early time points. Transcripts of *pfpeg4* and *pfg27* were detected at significant levels (above the LoQ) at several time points between D7 and D9 (Figure 3), demonstrating the reliability of the assay, even at a low density of circulating early gametocytes at this early stage of infection. In contrast, *pfs25* transcripts were detected, but below the reliable limit of quantification (LoQ), at only two time points in each volunteer, with a wide 95% CI of replica including zero (Figure 3).

In summary, both *pfpeg4* and *pfg27* qRT-PCR assays are specific and sensitive, and can quantify early *P. falciparum* gametocytes, as low as 190 and 390 gametocytes per mL of whole blood, respectively. The late-stage specificity (*pfs25* assay), currently used for field surveys, quantifies late gametocytes (Babiker et al., 1999). These early stage-specific qRT-PCR assays will facilitate the study of epidemiological factors that influence transmission, such as the impact of drugs and the multiplicity of infection, as well as climatic variables that govern fluctuations of the mosquito vector. A better understanding of the gametocyte reservoir in natural infections is essential for novel approach for malaria elimination strategies and for the assessment of transmission blocking strategies.

## Supporting information

Supplementary file [[supplements/254192_file03.pdf]](pending:yes)

## Data Availability

Data is available as supplementary file.

## Acknowledgments

We thank Department of Biochemistry staff, College of Medicine and Health Sciences, Sultan Qaboos University for support. Amal Gadalla was supported by SQU scholarship and International Atomic Energy Agency (IAEA).

*   Received March 28, 2021.
*   Revision received March 28, 2021.
*   Accepted April 4, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  Abdel-Wahab, A., Abdel-Muhsin, A. M., Ali, E., Suleiman, S., Ahmed, S., Walliker, D. & Babiker, H. A. 2002. Dynamics of gametocytes among Plasmodium falciparum clones in natural infections in an area of highly seasonal transmission. J Infect Dis, 185, 1838–42.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/340638&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12085337&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000175940800025&link_type=ISI) 

2.  Aguilar, R., Magallon-Tejada, A., Achtman, A. H., Moraleda, C., Joice, R., Cistero, P., Li Wai Suen, C. S., Nhabomba, A., Macete, E., Mueller, I., Marti, M., Alonso, P. L., Menendez, C., Schofield, L. & Mayor, A. 2014. Molecular evidence for the localization of Plasmodium falciparum immature gametocytes in bone marrow. Blood, 123, 959–66.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJsb29kam91cm5hbCI7czo1OiJyZXNpZCI7czo5OiIxMjMvNy85NTkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNC8wNC8yMDIxLjAzLjI4LjIxMjU0MTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

3.  Ali, E., Mackinnon, M. J., Abdel-Muhsin, A. A., Ahmed, S., Walliker, D. & Babiker, H. A. 2006. Increased density but not prevalence of gametocytes following drug treatment of Plasmodium falciparum. Trans R Soc Trop Med Hyg, 100, 176–83.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.trstmh.2005.04.021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16257026&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

4.  Babiker, H. A., Abdel-Wahab, A., Ahmed, S., Suleiman, S., Ranford-Cartwright, L., Carter, R. & Walliker, D. 1999. Detection of low level Plasmodium falciparum gametocytes using reverse transcriptase polymerase chain reaction. Mol Biochem Parasitol, 99, 143–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0166-6851(98)00175-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10215031&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079510300013&link_type=ISI) 

5.  Babiker, H. A. & Schneider, P. 2008. Application of molecular methods for monitoring transmission stages of malaria parasites. Biomed Mater, 3, 034007.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/1748-6041/3/3/034007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18708712&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

6.  Bachmann, A., Esser, C., Petter, M., Predehl, S., Von Kalckreuth, V., Schmiedel, S., Bruchhaus, I. & Tannich, E. 2009. Absence of Erythrocyte Sequestration and Lack of Multicopy Gene Family Expression in Plasmodium falciparum from a Splenectomized Malaria Patient. PLOS ONE, 4, e7459.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0007459&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19826486&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

7.  Bancells, C., Llorà-Batlle, O., Poran, A., NöTzel, C., Rovira-Graells, N., Elemento, O., Kafsack, B. F. C. & Cortés, A. 2019. Revisiting the initial steps of sexual development in the malaria parasite Plasmodium falciparum. Nature Microbiology, 4, 144–154.
    
    

8.  Bousema, J., Gouagna, L., Drakeley, C., Meutstege, A., Okech, B., Akim, I., Beier, J., Githure, J. & Sauerwein, R. 2004. Plasmodium falciparum gametocyte carriage in asymptomatic children in western Kenya. Malaria Journal, 3, 18.
    
    

9.  Boyer, T. C., Hanson, T. & Singer, R. S. 2013. Estimation of Low Quantity Genes: A Hierarchical Model for Analyzing Censored Quantitative Real-Time PCR Data. PLOS ONE, 8, e64900.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0064900&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23741414&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

10. Bruce, M. C., Alano, P., Duthie, S. & Carter, R. 1990. Commitment of the malaria parasite Plasmodium falciparum to sexual and asexual development. Parasitology, 100, 191–200.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0031182000061199&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2189114&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DB88900004&link_type=ISI) 

11. Bruce, M. C., Carter, R. N., Nakamura, K.-I., Aikawa, M. & Carter, R. 1994. Cellular location and temporal expression of the Plasmodium falciparum sexual stage antigen Pfs16. Molecular and Biochemical Parasitology, 65, 11–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0166-6851(94)90111-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7935618&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994NW86200002&link_type=ISI) 

12. Buchholz, K., Burke, T. A., Williamson, K. C., Wiegand, R. C., Wirth, D. F. & Marti, M. 2011. A High-Throughput Screen Targeting Malaria Transmission Stages Opens New Avenues for Drug Development. The Journal of Infectious Diseases, 203, 1445–1453.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jir037&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21502082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

13. Carter, L. M., Kafsack, B. F., LlináS, M., Mideo, N., Pollitt, L. C. & Reece, S. E. 2013. Stress and sex in malaria parasites: Why does commitment vary? Evol Med Public Health, 2013, 135–47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/emph/eot011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24481194&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

14. Carter, R. & Miller, L. H. 1979. Evidence for environmental modulation of gametocytogenesis in Plasmodium falciparum in continuous culture. Bull World Health Organ, 57 Suppl 1, 37–52.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=397008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

15. D’Alessandro, S., Silvestrini, F., Dechering, K., Corbett, Y., Parapini, S., Timmerman, M., Galastri, L., Basilico, N., Sauerwein, R., Alano, P. & Taramelli, D. 2013. A Plasmodium falciparum screening assay for anti-gametocyte drugs based on parasite lactate dehydrogenase detection. J Antimicrob Chemother, 68, 2048–58.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jac/dkt165&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23645588&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000323424100015&link_type=ISI) 

16. Eksi, S., Haile, Y., Furuya, T., Ma, L., Su, X. & Williamson, K. C. 2005. Identification of a subtelomeric gene family expressed during the asexual–sexual stage transition in Plasmodium falciparum. Molecular and Biochemical Parasitology, 143, 90–99.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molbiopara.2005.05.010&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15996767&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231335800011&link_type=ISI) 

17. Eksi, S., Morahan, B. J., Haile, Y., Furuya, T., Jiang, H., Ali, O., Xu, H., Kiattibutr, K., Suri, A., Czesny, B., Adeyemo, A., Myers, T. G., Sattabongkot, J., Su, X.-Z. & Williamson, K. C. 2012. Plasmodium falciparum Gametocyte Development 1 (Pfgdv1) and Gametocytogenesis Early Gene Identification and Commitment to Sexual Development. PLoS Pathog, 8, e1002964.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1002964&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23093935&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

18. ENSEMBL.ORG. 2014. EnsEMBL release 75-February 2014 [Online]. Available: [http://www.ensembl.org/Help/Faq?id=392](http://www.ensembl.org/Help/Faq?id=392) [Accessed].
    
    

19. Farid, R., Dixon, M. W., Tilley, L. & Mccarthy, J. S. 2017. Initiation of gametocytogenesis at very low parasite density in Plasmodium falciparum infection. The Journal of Infectious Diseases, 215, 1167–1174.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jix035&link_type=DOI) 

20. Furuya, T., Mu, J., Hayton, K., Liu, A., Duan, J., Nkrumah, L., Joy, D. A., Fidock, D. A., Fujioka, H., Vaidya, A. B., Wellems, T. E. & Su, X.-Z. 2005. Disruption of a Plasmodium falciparum gene linked to male sexual development causes early arrest in gametocytogenesis. Proceedings of the National Academy of Sciences of the United States of America, 102, 16813–16818.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTAyLzQ2LzE2ODEzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMDQvMjAyMS4wMy4yOC4yMTI1NDE5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

21. Gadalla, A. A. H., Schneider, P., Churcher, T. S., Nassir, E., Abdel-Muhsin, A.-M. A., Ranford-Cartwright, L. C., Reece, S. E. & Babiker, H. A. 2016. Associations between Season and Gametocyte Dynamics in Chronic Plasmodium falciparum Infections. PloS one, 11, e0166699–e0166699.
    
    

22. Gupta, S. K., Schulman, S. & Vanderberg, J. P. 1985. Stage-Dependent Toxicity of N-Acetyl-Glucosamine to Plasmodium falciparum1. The Journal of Protozoology, 32, 91–95.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1550-7408.1985.tb03020.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3886901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

23. Hermsen, C. C., Telgt, D. S. C., Linders, E. H. P., Van De Locht, L. A. T. F., Eling, W. M. C., Mensink, E. J. B. M. & Sauerwein, R. W. 2001. Detection of Plasmodium falciparum malaria parasites in vivo by real-time quantitative PCR. Molecular and Biochemical Parasitology, 118, 247–251.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0166-6851(01)00379-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11738714&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173251600013&link_type=ISI) 

24. Joice, R., Narasimhan, V., Montgomery, J., Sidhu, A. B., Oh, K., Meyer, E., Pierre-Louis, W., Seydel, K., Milner, D., Williamson, K., Wiegand, R., Ndiaye, D., Daily, J., Wirth, D., Taylor, T., Huttenhower, C. & Marti, M. 2013. Inferring developmental stage composition from gene expression in human malaria. PLoS Comput Biol, 9, e1003392.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1003392&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24348235&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

25. Joice, R., Nilsson, S. K., Montgomery, J., Dankwa, S., Egan, E., Morahan, B., Seydel, K. B., Bertuccini, L., Alano, P., Williamson, K. C., Duraisingh, M. T., Taylor, T. E., Milner, D. A. & Marti, M. 2014. Plasmodium falciparum transmission stages accumulate in the human bone marrow. Science translational medicine, 6, 244re5–244re5.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6InNjaXRyYW5zbWVkIjtzOjU6InJlc2lkIjtzOjEyOiI2LzI0NC8yNDRyZTUiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNC8wNC8yMDIxLjAzLjI4LjIxMjU0MTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

26. Klymus, K. E., Merkes, C. M., Allison, M. J., Goldberg, C. S., Helbing, C. C., Hunter, M. E., Jackson, C. A., Lance, R. F., Mangan, A. M., Monroe, E. M., Piaggio, A. J., Stokdyk, J. P., Wilson, C. C. & Richter, C. A. 2020. Reporting the limits of detection and quantification for environmental DNA assays. Environmental DNA, 2, 271–282.
    
    

27. Lobo, C.-A., Ruud, N. H. Konings & Kumar, N. 1994. Expression of early gametocyte-stage antigens Pfg27 and Pfs16 in synchronized gametocytes and non-gametocyte producing clones of Plasmodium falciparum. Molecular and Biochemical Parasitology, 68, 151–154.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0166-6851(94)00155-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7891740&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

28. Lopez-Barragan, M. J., Lemieux, J., Quinones, M., Williamson, K. C., Molina-Cruz, A., Cui, K., Barillas-Mury, C., Zhao, K. & Su, X. Z. 2011. Directional gene expression and antisense transcripts in sexual and asexual stages of Plasmodium falciparum. BMC Genomics, 12, 587.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2164-12-587&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22129310&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

29. Mackellar, D. C., O’Neill, M. T., Aly, A. S. I., Sacci, J. B., Cowman, A. F. & Kappe, S. H. I. 2010. Plasmodium falciparum PF10_0164 (ETRAMP10.3) Is an Essential Parasitophorous Vacuole and Exported Protein in Blood Stages. Eukaryotic Cell, 9, 784–794.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZXVrY2VsbCI7czo1OiJyZXNpZCI7czo3OiI5LzUvNzg0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDQvMDQvMjAyMS4wMy4yOC4yMTI1NDE5Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

30. Mawili-Mboumba, D., Nikiema, R., Bouyou-Akotet, M., Bahamontes-Rosa, N., Traore, A. & Kombila, M. 2013. Sub-microscopic gametocyte carriage in febrile children living in different areas of Gabon. Malaria Journal, 12, 375.
    
    

31. Menegon, M., Severini, C., Sannella, A., Paglia, M. G., Sangaré, D., Abdel-Wahab, A., Abdel-Muhsin, A.-M. A., Babiker, H., Walliker, D. & Alano, P. 2000. Genotyping of Plasmodium falciparum gametocytes by reverse transcriptase polymerase chain reaction. Molecular and Biochemical Parasitology, 111, 153–161.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0166-6851(00)00314-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11087925&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000165702100014&link_type=ISI) 

32. Nassir, E., Abdel-Muhsin, A. M., Suliaman, S., Kenyon, F., Kheir, A., Geha, H., Ferguson, H. M., Walliker, D. & Babiker, H. A. 2005. Impact of genetic complexity on longevity and gametocytogenesis of Plasmodium falciparum during the dry and transmission-free season of eastern Sudan. Int J Parasitol, 35, 49–55.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpara.2004.10.014&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15619515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000226570700006&link_type=ISI) 

33. Nwakanma, D. C., Gomez-Escobar, N., Walther, M., Crozier, S., Dubovsky, F., Malkin, E., Locke, E. &Conway, D. J. 2009. Quantitative detection of Plasmodium falciparum DNA in saliva, blood, and urine. J Infect Dis, 199, 1567–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/598856&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19432544&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266000600003&link_type=ISI) 

34. Obaldia, N., 3rd, Meibalan, E., Sa, J. M., Ma, S., Clark, M. A., Mejia, P., Moraes Barros, R. R., Otero, W., Ferreira, M. U., Mitchell, J. R., Milner, D. A., Huttenhower, C., Wirth, D. F., Duraisingh, M. T., Wellems, T. E. & Marti, M. 2018. Bone Marrow Is a Major Parasite Reservoir in Plasmodium vivax Infection. mBio, 9.
    
    

35. Oué Draogo, A. L., Bousema, T., Schneider, P., De Vlas, S. J., Ilboudo-Sanogo, E., Cuzin-Ouattara, N., Né Bié, I., Roeffen, W., Verhave, J. P., Luty, A. J. & Sauerwein, R. 2009. Substantial contribution of submicroscopical Plasmodium falciparum gametocyte carriage to the infectious reservoir in an area of seasonal transmission. PLoS One, 4.
    
    

36. Pasay, C. J., Rockett, R., Sekuloski, S., Griffin, P., Marquart, L., Peatey, C., Wang, C. Y. T., O’Rourke, P., Elliott, S., Baker, M., Möhrle, J. J. & Mccarthy, J. S. 2016. Piperaquine Monotherapy of Drug-Susceptible Plasmodium falciparum Infection Results in Rapid Clearance of Parasitemia but Is Followed by the Appearance of Gametocytemia. The Journal of Infectious Diseases, 214, 105–113.
    
    

37. Pelle, K. G., Oh, K., Buchholz, K., Narasimhan, V., Joice, R., Milner, D. A., Brancucci, N. M. B., Ma, S., Voss, T. S., Ketman, K., Seydel, K. B., Taylor, T. E., Barteneva, N. S., Huttenhower, C. & Marti, M. 2015. Transcriptional profiling defines dynamics of parasite tissue sequestration during malaria infection. Genome Medicine, 7, 19.
    
    

38. PLASMODB.ORG. 2014. PlasmoDB: a functional genomic database for malaria parasites. Nucleic Acids Res. 2008 Oct 31. Aurrecoechea C, et al. [Online]. Available: [http://plasmodb.org/plasmo/home.do](http://plasmodb.org/plasmo/home.do) [Accessed].
    
    

39. Poran, A., Nötzel, C., Aly, O., Mencia-Trinchant, N., Harris, C. T., Guzman, M. L., Hassane, D. C., Elemento, O. & Kafsack, B. F. C. 2017. Single-cell RNA sequencing reveals a signature of sexual commitment in malaria parasites. Nature, 551, 95–99.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature24280&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29094698&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

40. Reece, S. E., Ali, E., Schneider, P. & Babiker, H. A. 2010. Stress, drugs and the evolution of reproductive restraint in malaria parasites. Proc Biol Sci, 277, 3123–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rspb.2010.0564&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20484242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000282645700009&link_type=ISI) 

41. Reece, S. E., Drew, D. R. & Gardner, A. 2008. Sex ratio adjustment and kin discrimination in malaria parasites. Nature, 453, 609–14.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature06954&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18509435&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000256185200034&link_type=ISI) 

42. Ribaut, C., Berry, A., Chevalley, S., Reybier, K., Morlais, I., Parzy, D., Nepveu, F., Benoit-Vical, F. & Valentin, A. 2008. Concentration and purification by magnetic separation of the erythrocytic stages of all human Plasmodium species. Malaria Journal, 7, 45.
    
    

43. Rogier, C., Commenges, D. & Trape, J. F. 1996. Evidence for an age-dependent pyrogenic threshold of Plasmodium falciparum parasitemia in highly endemic populations. Am J Trop Med Hyg, 54, 613–9.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoidHJvcG1lZCI7czo1OiJyZXNpZCI7czo4OiI1NC82LzYxMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA0LzA0LzIwMjEuMDMuMjguMjEyNTQxOTIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

44. Roucher, C., Rogier, C., Dieye-Ba, F., Sokhna, C., Tall, A. & Trape, J.-F. 2012. Changing Malaria Epidemiology and Diagnostic Criteria for Plasmodium falciparum Clinical Malaria. PLOS ONE, 7, e46188.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0046188&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23029433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

45. Schneider, P., Bousema, T., Omar, S., Gouagna, L., Sawa, P., Schallig, H. & Sauerwein, R. 2006. (Sub)microscopic Plasmodium falciparum gametocytaemia in Kenyan children after treatment with sulphadoxine-pyrimethamine monotherapy or in combination with artesunate. International Journal for Parasitology, 36, 403–408.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijpara.2006.01.002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16500657&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237861800003&link_type=ISI) 

46. Schneider, P., Reece, S. E., Van Schaijk, B. C., Bousema, T., Lanke, K. H., Meaden, C. S., Gadalla, A., Ranford-Cartwright, L. C. & Babiker, H. A. 2015. Quantification of female and male Plasmodium falciparum gametocytes by reverse transcriptase quantitative PCR. Mol Biochem Parasitol, 199, 29–33.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molbiopara.2015.03.006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25827756&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

47. Silvestrini, F., Bozdech, Z., Lanfrancotti, A., Giulio, E. D., Bultrini, E., Picci, L., Derisi, J. L., Pizzi, E. & Alano, P. 2005. Genome-wide identification of genes upregulated at the onset of gametocytogenesis in Plasmodium falciparum. Molecular and Biochemical Parasitology, 143, 100–110.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.molbiopara.2005.04.015&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16026866&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231335800012&link_type=ISI) 

48. Silvestrini, F., Lasonder, E., Olivieri, A., Camarda, G., Van Schaijk, B., Sanchez, M., Younis Younis, S., Sauerwein, R. & Alano, P. 2010. Protein Export Marks the Early Phase of Gametocytogenesis of the Human Malaria Parasite Plasmodium falciparum. Molecular & Cellular Proteomics, 9, 1437–1448.
    
    

49. Smalley, M. E., Abdalla, S. & Brown, J. 1981. The distribution of Plasmodium falciparum in the peripheral blood and bone marrow of Gambian children. Trans R Soc Trop Med Hyg, 75, 103–5.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0035-9203(81)90019-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7022784&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F04%2F04%2F2021.03.28.21254192.atom) 

50. Tadesse, F. G., Van Den Hoogen, L., Lanke, K., Schildkraut, J., Tetteh, K., Aseffa, A., Mamo, H., Sauerwein, R., Felger, I., Drakeley, C., Gadissa, E. & Bousema, T. 2017. The shape of the iceberg: quantification of submicroscopic Plasmodium falciparum and Plasmodium vivax parasitaemia and gametocytaemia in five low endemic settings in Ethiopia. Malaria journal, 16, 99–99.
    
    

51. Tibúrcio, M., Dixon, M. W. A., Looker, O., Younis, S. Y., Tilley, L. & Alano, P. 2015. Specific expression and export of the Plasmodium falciparum Gametocyte EXported Protein-5 marks the gametocyte ring stage. Malaria Journal, 14, 334.
    
    

52. Trager, W. & Jensen, J. B. 1976. Human malaria parasites in continuous culture. Science, 193, 673–5.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIxOTMvNDI1NC82NzMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wNC8wNC8yMDIxLjAzLjI4LjIxMjU0MTkyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

53. Usui, M., Prajapati, S. K., Ayanful-Torgby, R., Acquah, F. K., Cudjoe, E., Kakaney, C., Amponsah, J. A., Obboh, E. K., Reddy, D. K., Barbeau, M. C., Simons, L. M., Czesny, B., Raiciulescu, S., Olsen, C., Abuaku, B. K., Amoah, L. E. & Williamson, K. C. 2019. Plasmodium falciparum sexual differentiation in malaria patients is associated with host factors and GDV1-dependent genes. Nat Commun, 10, 2140.