# Uptake of Plasmodium falciparum Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding

Arthur M. Talman, Dinkorma Ouologuem, Katie Love, Virginia Howick, Charles Mulamba, Aboubecrin Haidara, Niawanlou Dara, Daman Sylla, Adama Sacko, Mamadou Coulibaly, et al.

► To cite this version:

Arthur M. Talman, Dinkorma Ouologuem, Katie Love, Virginia Howick, Charles Mulamba, et al. Uptake of Plasmodium falciparum Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding. Frontiers in Microbiology, 2020, 11, pp.246. 10.3389/fmicb.2020.00246. hal-04973617

HAL Id: hal-04973617

https://hal.science/hal-04973617v1

Submitted on 3 Mar 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Check for
updates

# Uptake of *Plasmodium falciparum* Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding

Arthur M. Talman<fnref n="1" />,², Dinkorma T. D. Ouologuem<fnref n="3" />, Katie Love<fnref n="1" />, Virginia M. Howick<fnref n="1" />, Charles Mulamba<fnref n="1" />, Aboubecrin Haidara<fnref n="3" />, Niawanlou Dara<fnref n="3" />, Daman Sylla<fnref n="3" />, Adama Sacko<fnref n="3" />, Mamadou M. Coulibaly<fnref n="3" />, Francois Dao<fnref n="3" />, Cheick P. O. Sangare<fnref n="3" />, Abdoulaye Djimde<fnref n="3" /> and Mara K. N. Lawniczak<fnref n="1" />*

OPEN ACCESS

**Edited by:**
Brandon Keith Wilder,
Oregon Health & Science University,
United States

**Reviewed by:**
Robert Sinden,
Imperial College London,
United Kingdom
Chris Drakeley,
University of London, United Kingdom

*Correspondence:
Mara K. N. Lawniczak
mara@sanger.ac.uk

**Specialty section:**
This article was submitted to
Infectious Diseases,
a section of the journal
Frontiers in Microbiology

Received: 03 September 2019

Accepted: 03 February 2020

Published: 03 March 2020

**Citation:**
Talman AM, Ouologuem DTD, Love K, Howick VM, Mulamba C, Haidara A, Dara N, Sylla D, Sacko A, Coulibaly MM, Dao F, Sangare CPO, Djimde A and Lawniczak MKN (2020) *Uptake of Plasmodium falciparum Gametocytes During Mosquito Bloodmeal by Direct and Membrane Feeding*. Front. Microbiol. 11:246. doi: 10.3389/fmicb.2020.00246

*Plasmodium falciparum* remains one of the leading causes of child mortality, and nearly half of the world's population is at risk of contracting malaria. While pathogenesis results from replication of asexual forms in human red blood cells, it is the sexually differentiated forms, gametocytes, which are responsible for the spread of the disease. For transmission to succeed, both mature male and female gametocytes must be taken up by a female *Anopheles* mosquito during its blood meal for subsequent differentiation into gametes and mating inside the mosquito gut. Observed circulating numbers of gametocytes in the human host are often surprisingly low. A pre-fertilization behavior, such as skin sequestration, has been hypothesized to explain the efficiency of human-to-mosquito transmission but has not been sufficiently tested due to a lack of appropriate tools. In this study, we describe the optimization of a qPCR tool that enables the relative quantification of gametocytes within very small input samples. Such a tool allows for the quantification of gametocytes in different compartments of the host and the vector that could potentially unravel mechanisms that enable highly efficient malaria transmission. We demonstrate the use of our gametocyte quantification method in mosquito blood meals from both direct skin feeding on *Plasmodium* gametocyte carriers and standard membrane feeding assay. Relative gametocyte abundance was not different between mosquitoes fed through a membrane or directly on the skin suggesting that there is no systematic enrichment of gametocytes picked up in the skin.

**Keywords:** malaria, transmission, gametocyte, mosquito feeding, *Plasmodium falciparum*

## INTRODUCTION

Transmission of *Plasmodium falciparum* from humans to mosquitoes depends on the sexual phase of the parasite's life cycle. Both male and female gametocytes have to be picked up in the mosquito blood meal in order to mate and colonize the mosquito. Relatively low circulating gametocyte densities are typically observed (Bousema and Drakeley, 2011), and yet transmission remains efficient even at gametocyte densities for which random

<footnote n="1">¹ Wellcome Sanger Institute, Hinxton, United Kingdom, ² MIVEGEC, IRD, CNRS, University of Montpellier, Montpellier, France,</footnote>

<footnote n="3">³ Malaria Research and Training Centre, University of Science, Techniques and Technologies of Bamako, Bamako, Mali</footnote>
sampling of the host cannot explain successful infection of the mosquito (Lawniczak and Eckhoff, 2016). Sub-patently infected carriers that harbor parasites below the limit of detection by microscopy can readily infect mosquitoes (Schneider et al., 2007; Ouédraogo et al., 2009). Mosquitoes fed directly on the skin have also been found to be infected more readily than by a direct membrane feeding assay (Bonnet et al., 2000; Diallo et al., 2008). These findings have led to the hypothesis of gametocyte biology and behavior that enhance their transmission but have yet to be described (Pichon et al., 2000; Bousema and Drakeley, 2011). A pre-fertilization behavior, whereby malaria gametocytes associate in the circulating blood and/or adhere to sub-dermal capillaries could enhance their probability of being ingested in sufficient quantities. Intriguingly, two early studies revealed that gametocytes were on average three times more concentrated in skin biopsies than in the venous circulation (Chardome and Janssen, 1952; Van Den Berghe et al., 1954), although these experiments lacked appropriate controls. Whilst associative behaviors have been hypothesized, they have never been explicitly tested. Comparing gametocyte uptake in skin-fed and membrane-fed mosquitoes could determine if the patterns differ between these blood sources. Such an undertaking, however, requires a sensitive detection tool able to accurately measure gametocyte densities in individual mosquito blood meals. The most commonly used gametocyte molecular detection tools, albeit of good sensitivity, only detect female gametocytes in human blood sample (Schneider et al., 2015). In this study, we describe a molecular assay capable of measuring relative gametocyte densities directly in mosquito bloodmeals and investigate the uptake of gametocytes during natural and artificial blood feeding of symptomatic malaria patients.

# MATERIALS AND METHODS

## Gametocyte Specific Transcript Selection for Ultrasensitive Detection and Multiplex Rt qPCR

In order to find a suitable gametocyte biomarker, 107 transcripts displaying a 25-fold enrichment in gametocytes compared to rings in a published stage-specific RNAseq study were selected (López-Barragán et al., 2011; **Supplementary Table S1**). qPCR primers to each of these transcripts were designed using the IDT PrimerQuest® Tool. *P. falciparum* isolate NF54 was maintained in O+ blood in RPMI 1640 culture medium (GIBCO) supplemented with 25 mM HEPES (SIGMA), 10 mM D-Glucose (SIGMA), 50 mg/L hypoxanthine (SIGMA), and 10% human serum in a gas mix containing 5% O₂, 5% CO₂, and 90% N₂. Human O+ erythrocytes were obtained from NHS Blood and Transplant, Cambridge, United Kingdom. None of the blood products used contained identifying information from donors. *Plasmodium* culture using human serum and erythrocytes from donors has been approved by the NHS Cambridgeshire 4 Research Ethics Committee (REC reference 15/EE/0253) and the Wellcome Sanger Institute Human Materials and Data Management Committee. Pure ring stage parasites (10⁹) were

produced by double Percoll-sorbitol synchronization (Moll et al., 2008), followed by negative selection on a MACs LS column (Miltenyi Biotec). Pure stage V gametocyte samples (10⁸) were produced and purified on a MACs LS column (Ribaut et al., 2008). Parasite counts were established by Giemsa staining and hematocrit was measured with a hemocytometer. RNA was extracted with TRIzol according to the manufacturer's recommendations. RNAs were treated with DNA-free DNAse Turbo kit (Ambion) and reverse transcribed with the high capacity reverse transcriptase kit (Thermo Fisher), supplemented with oligo-dts (Thermo Fisher) at a final concentration of 2.5 µM. SYBR-green qPCR (Roche) was conducted on a Lightcycler 480 (Roche) for each of these transcripts on pure gametocyte and ring duplicate cDNAs. The DNA-free DNAse Turbo kit removes gDNA, a reverse transcriptase-less control was run with each extraction batch and was verified to be negative before inclusion of samples in the dataset. Enrichment of transcripts in gametocytes was estimated using the delta-delta *Ct* method for each 107 transcript and the housekeeping gene glucose-6-phosphate dehydrogenase-6-phosphogluconolactonase (PF3D7_1453800) as a control. A probe-based multiplex assay was designed to allow relative quantification and comparison of gametocytes within mosquito bloodmeals. A triplex qPCR assay was devised: it detected the gametocyte biomarker (*MDV-1*, PF3D7_1216500, cy5), an asexual transcript (*Mahrp2*, PF3D7_1353200, HEX) and a human transcript (*hsGAPDH*, FAM) (**Supplementary Table S2**). The *Mahrp2* transcript was selected because it was the most enriched transcript in rings compared to gametocytes by RNAseq; if one excludes *Hrp2* or *Hrp3*, two genes that have been shown to be deleted in some populations of parasites and are therefore not suitable (Watson et al., 2017). An alternative assay using the human transcript *UBA-1* was included to validate *hsGAPDH* as a suitable loading control for leukocytes. Multiplex primer efficiencies were calculated on qPCRs of gDNA serial dilutions for all primer/probe combinations (**Supplementary Table S3**). The qPCR conditions were as follows, the samples were first incubated for 10 min at 95°C, then 45 cycles were performed (95°C for 10 s, 60°C for 10 s, and 68°C for 20 s), followed by a melting curve step to ensure single product amplification. Single amplification was observed for all samples. The melting temperature were chosen as that recommended by the assay manufacturer without alteration. The *Ct* calculated by the relative quantification protocol of the Light cycle 480 software were used. Fold changes were calculated by the delta-delta *Ct* method corrected for primer efficiency (Pfaffl, 2004; p. 96 equation 3.5).

## Mock Blood Meals With a Gametocyte Serial Dilution

To test the triplex qPCR assay with differing concentrations of gametocytes, *in vitro* cultured asexual and sexual parasites were artificially fed to mosquitoes. A serial dilution of gametocytes (10–100,000 per blood meal) was spiked with 10,000 rings per µl, placed in whole blood and spun at 800 g for 5 min, heat-inactivated serum was used to resuspend the pellet to 50% hematocrit and aliquoted into a heated plastic membrane
feeder covered with stretched parafilm (Bemis NA). Eight female *Anopheles coluzzii* (N’Gousso strain) were allowed to feed on each dilution for 12 min and were immediately sacrificed in 70% ethanol. Fed mosquitoes were tapped briefly on absorbent paper to remove excess ethanol and immediately immersed in 50 µl of TriZol in RNAse-free 1.5 ml tubes and stored at -80°C. Upon thawing, 450 µl of fresh TriZol was added to each 50 µl sample, and mosquitoes were homogenized with a clean pestle and RNA extraction and cDNA generation were conducted as described above.

## Study Site and Ethical Approval

To test the assay in natural conditions, we recruited *P. falciparum* positive patients in Faladie, Mali. Faladie is situated in the Koulikoro region of Mali and is characterized by a seasonal hyperendemic transmission of mostly *P. falciparum* malaria. Patients, aged 6–14 years, with symptomatic non-severe malaria were recruited from November to December 2016 (Table 1). Patients were screened by thick smear microscopy during which the asexual parasite (rings) and gametocyte counts (stage V) were recorded for 1000 leukocytes. A standard concentration of 8000 leukocytes per µl was used to calculate parasitemia and gametocytaemias. The protocol was approved by an IRB from the Faculty Of Medicine, Pharmacy and Odontostomatology de Pharmacy; Université des Sciences, Techniques et Technologies of Bamako (IRB approval letter No. 2016/133/CE/FMPOS).

## Mosquito Feeding Assay

To compare the number of gametocytes acquired in the mosquito bloodmeal from feeding on the skin or peripheral bloodstream, a peripheral blood sample was first obtained from *P. falciparum* gametocyte carriers in a vacutainer with EDTA and immediately placed into two plastic membrane feeders covered with parafilm (600 µl each). 50 starved *A. coluzzii* females were allowed to feed on each feeder (n = 100 total). Membrane feedings were allowed to proceed for 12 min. Concomitantly 2 pots of 25 mosquitoes were allowed to direct skin-feed on each volunteer (back of left calf, right calf) until repletion (8 min approximatively). Following the feeds, all mosquitoes were immediately sacrificed with 70% ethanol. Mosquitoes were each transferred to 50 µl of TriZol and kept in liquid nitrogen or at -80°C until extraction.

## Mosquito Sample Processing for Parasite Quantification

For RNA extraction, upon thawing, 450 µl of fresh TriZol was added to each 50 µl sample, and mosquitoes were homogenized with a clean pestle. RNA extraction and cDNA generation were conducted as above. Two assays per sample were run: the triplex mentioned above (*MDV-1*, *hsGAPDH*, and *Mahrp2*) and a duplex with the alternative human transcript *hsUBA-1* (FAM) and a control mosquito transcript (ribosomal protein S7, HEX). Samples that yielded Ct > 38 for *MDV-1* or *hsGAPDH* were discarded. Fold changes were calculated using the delta-delta Ct method corrected for primer efficiency (Pfaffl, 2004; Supplementary Table S4); the test gene was always *MDV-1*

<table><caption>TABLE 1 | Feed data for six patients. Parasitemia and gametocytemia were established by microscopy.</caption><thead><tr><th rowspan="2">Patient number</th><th rowspan="2">Parasitemia-<br>day 0<br>(Microscopy)<br>(µl)</th><th rowspan="2">Gametocytemia<br>day 0<br>(Microscopy)<br>(µl)</th><th rowspan="2">Membrane<br>feeder 1<br>(positive<br>blood<br>meals/total)</th><th rowspan="2">Membrane<br>feeder 2<br>(positive<br>blood<br>meals/total)</th><th rowspan="2">Skin feed left<br>(positive<br>blood<br>meals/total)</th><th rowspan="2">Skin feed<br>right (positive<br>blood<br>meals/total)</th><th colspan="2">Variance (F-test)</th><th colspan="2">t-test</th></tr><tr><th>Variance in<br>skinfed than<br>membrane<br>fed</th><th>p-value</th><th>Mean fold<br>change in<br>skinfed vs<br>membrane<br>fed</th><th>p-value</th></tr></thead><tbody><tr><td>EGF006</td><td>14,680</td><td>520</td><td>16/17</td><td>15/17</td><td>16/17</td><td>12/13</td><td>Greater</td><td>0.03091</td><td>NS</td><td>0.1986</td></tr><tr><td>EGF007</td><td>37,640</td><td>48</td><td>17/18</td><td>6/8</td><td>19/20</td><td>14/15</td><td>Greater</td><td>0.006914</td><td>NS</td><td>0.05617</td></tr><tr><td>EGF009</td><td>4800</td><td>144</td><td>15/15</td><td>11/12</td><td>17/19</td><td>10/10</td><td>NS</td><td>0.1377</td><td>NS</td><td>0.617</td></tr><tr><td>EGF014</td><td>1440</td><td>360</td><td>16/17</td><td>9/9</td><td>15/15</td><td>19/20</td><td>Greater</td><td>0.01363</td><td>NS</td><td>0.2308</td></tr><tr><td>EGF017</td><td>512</td><td>296</td><td>21/21</td><td>20/21</td><td>23/23</td><td>23/23</td><td>Less</td><td>0.01764</td><td>NS</td><td>0.3374</td></tr><tr><td>EGF020</td><td>17,520</td><td>96</td><td>14/21</td><td>12/17</td><td>9/21</td><td>11/24</td><td>Less</td><td>2.545E-06</td><td>less</td><td>3.06E-05</td></tr><tr><td>Overall</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>Less</td><td>8.138e-06</td><td>NS</td><td>0.05251</td></tr></tbody><tfoot><tr><td colspan="12">Number of fed mosquitoes with a positive MDV-1 Ct (&lt;38) is given together with the total number of mosquitoes for each feed.</td></tr></tfoot></table>
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | The GQA, a multiplex assay capable of sensitively detecting gametocyte uptake by mosquitoes. (A) RNAseq (López-Barragán et al., 2011) and qPCR fold change of 107 transcripts in gametocytes compared to ring stage parasites. The color scale indicates sex-bias in expression based on (Lasonder et al., 2016) with high female:male ratios in red and balanced ratios in blue. (B) Relative abundance of MDV-1 in blood meals of mosquitoes fed on a serial dilution of gametocytes compared to the maximum gametocyte input (100,000/bloodmeal). Error bars represent the standard deviation. (C) Relationship between observed gametocyte numbers by microscopy and delta Ct of test and control gene in 10 µl blood samples from symptomatic malaria patients colored by asexual parasitemia (26 samples from 6 patients).

whilst the control genes used were either *hsUBA-1*, *Mahrp2* or *hsGAPDH*.

## Sampling and Statistical Evaluation

Statistical analyses were conducted in R and using package "car" for the ANOVA (Fox and Weisberg, 2018). F-tests were run with the var.test() function.

## RESULTS

### An Assay to Measure Gametocyte Uptake in the Bloodmeal

Because a mosquito bloodmeal is composed of only a few microliters of blood, we reasoned that maximum sensitivity in detecting bloodmeal gametocyte density will be achieved by detecting the most abundant transcripts found in gametocytes. We selected 107 transcripts that were at least 25-fold enriched in gametocytes over the ring form from a previous bulk RNAseq study (López-Barragán et al., 2011). We designed

individual RT-qPCR for the candidate transcripts in order to confirm their relative abundance. Previous bulk RNAseq work and our specific qPCR results on 107 transcripts were in good agreement (Figure 1A). We also overlaid a sex-specific transcriptome analysis (Lasonder et al., 2016), which revealed that the most abundant transcripts tend to be contributed by female gametocytes (Figure 1A). The most abundant transcript by both qPCR and RNAseq that did not show a sex bias was *MDV-1* (sex-bias = 1.01).

We tested the sensitivity of this potential gametocyte biomarker by qRT-PCR with a serial dilution of *in vitro* cultured gametocytes mixed with purified ring stage parasites in whole blood (spanning the range of 10–100 k per bloodmeal with a blood meal considered to be 2 µl). *A. coluzzii* females were fed on this serial dilution and immediately sacrificed, and whole mosquitoes were used to generate cDNA. *MDV-1* transcript abundance was normalized with either *Mahrp2*, which is a ring stage transcript, or with *hsGAPDH*, which is transcript present in human leukocyte. The gametocyte numbers estimated to be present in each dilution was significantly correlated to the fold
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Comparison of gametocyte density in natural and artificial mosquito blood meals. Gametocyte detected in mosquitoes fed on six patients either by direct skin feeding or membrane feeding. The left panel shows the data pooled across patients and the right displays results for each individual. Hinges correspond to the first and third quartiles, whiskers extend 1.5 interquartile range on either side of the hinges. Jitter points are overlaid, those beyond the whiskers are outlying points.</figcaption></figure>

changes normalized with *Mahrp2* (*r* = 0.79, *p*-values < 0.001 by Pearson's correlation) and *hsGAPDH* (*r* = 0.76, *p*-values < 0.001 by Pearson's correlation), indicating that the gametocytemia can be relatively quantified by these assays within 8 individual mosquito bloodmeals per condition over a wide dynamic range (10–100,000 gametocytes per blood meal) (**Figure 1B**). We found only 3 of 8 mosquitoes fed on a blood meal with 10 gametocytes had a signal below the limit of detection threshold (*Ct* < 38), suggesting the limit of reliable detection for the assay is between 10 and 100 gametocytes per bloodmeal. We named this assay gametocyte quantification assay (GQA).

We next tested the GQA on mock blood meals composed of 10 µl of venous blood from 6 symptomatic gametocyte carriers recruited in Faladje and one unfed mosquito (**Table 1**). We observed that peripheral gametocytemia as measured by microscopy is highly correlated with the delta Ct from the GQA (*r* = −0.82, *p* = 3.567e-07 by Pearson's correlation) and was independent of peripheral asexual parasitemia (**Figure 1C**). No single nucleotide polymorphism with a minor allele frequency superior to 0.002 was found in the MDV-1 region targeted by the GQA > 3000 genomes (MalariaGEN Plasmodium falciparum Community Project, 2016), indicating that genetic diversity in

natural populations is unlikely to affect the performance of the GQA. Altogether these results indicate that the GQA can be used to measure relative gametocytaemias in minute quantities of blood in both cultured parasites and for wild populations of parasites, it should be noted that the assay will be most useful for densities of gametocytes detectable by microscopy and not sensitive enough for subpatent gametocyte carriage.

## Uptake of Gametocytes During Natural and Artificial Bloodmeals

We next used the GQA to measure the relative uptake of gametocytes during mosquito blood feeding either directly from the skin or through an artificial membrane to test the potential role of the skin in gametocyte transmission biology. For the aforementioned six patients, we performed the GQA measurement after direct or membrane feeding. Control samples used in the quantification for each volunteer were mock blood meal samples (see above). We used the linear model: *fold change* ~ *feed* + *patient* + *feed* × *patient* (where *feed* is the feeding mode, skin or membrane). Relative gametocyte abundance was not different between mosquitoes fed through a membrane or directly on the skin when normalized to human *GAPDH* (type III ANOVA, *p* = 0.39331) (**Figure 2** and **Supplementary Table S4**), an alternative human transcript *UBA-1* (type III ANOVA, *p* = 0.18121), or the asexual transcript *Mahrp2* (type III ANOVA, *p* = 0.6999) (**Table 2**). This is indicative that there was no difference in density of gametocytes during natural or artificial blood meals. We did, however, note a patient-dependent effect on feeding mode but with inconsistent directionality (type III ANOVA, *feed* × *patient*: *p* = 0.01744). We, therefore, tested differences in the fold change of the gametocyte marker for each patient and found no difference in skin-fed vs membrane-fed except in the case of one patient (EGF020, **Table 1**). We also note that the variance observed between skin-fed and membrane-fed in each patient

<table><caption>**TABLE 2** | ANOVA full results.</caption><thead><tr><th>Reference gene</th><th>Model component</th><th>F-value</th><th>Pr(&gt;F)</th></tr></thead><tbody><tr><td>GAPDH</td><td>Feed</td><td>0.7305</td><td>0.39331</td></tr><tr><td>GAPDH</td><td>Patient</td><td>6.6312</td><td>6.461e-06</td></tr><tr><td>GAPDH</td><td>Feed × patient</td><td>2.7878</td><td>0.01744</td></tr><tr><td>UBA-1</td><td>Feed</td><td>1.7948</td><td>0.18121</td></tr><tr><td>UBA-1</td><td>Patient</td><td>5.5182</td><td>6.612e-05</td></tr><tr><td>UBA-1</td><td>Feed × patient</td><td>2.6143</td><td>0.02445</td></tr><tr><td>Mahrp2</td><td>Feed</td><td>0.1488</td><td>0.6999</td></tr><tr><td>Mahrp2</td><td>Patient</td><td>1.7555</td><td>0.1214</td></tr><tr><td>Mahrp2</td><td>Feed × patient</td><td>1.4370</td><td>0.2102</td></tr></tbody></table>
were different in 5 of the 6 patients but again with inconsistent
directionality (Table 1).

DISCUSSION

Malaria transmission from human to mosquito is a very efficient
process that is still not completely understood (Bousema and
Drakeley, 2011). In this study, we have established a new tool
to study *P. falciparum* transmission in natural infections. The
GQA detects a new biomarker of *P. falciparum* gametocytes
(*MDV-1*) which is both sensitive and detects both males and
female gametocytes. The GQA is particularly apt at comparing
blood samples originating from the same individual because
gametocyte densities can be normalized to either asexual parasite
or leukocyte transcripts, which should be good loading controls
within an infection. The GQA will, therefore, be extremely
useful to assess the spatial and temporal heterogeneity in
gametocyte densities within an infection. Moreover, this assay
requires very little initial input, therefore it can be used for
finger prick blood samples, skin punctures or mosquito blood
meals, as we have shown in this study. It should be noted
that the GQA is used at relatively high gametocytaemias
(>100 gametocyte/µl) compared to lower gametocytaemias that
may still yield successful mosquito infections.

Associative behaviors of gametocyte pre-fertilization have
been postulated as a mechanism for enhanced transmission
(Pichon et al., 2000; Lawniczak and Eckhoff, 2016; Nixon,
2016). In this study we have not seen a different density
of gametocytes in mosquito bloodmeals taken through skin
feeding versus membrane feeding on volunteers with patent
gametocytemia. Therefore, in our gametocyte concentration
window of observation (48-520 gametocytes/µl) there may
not be a differential uptake of gametocytes during skin
feeding. Another recent study comparing capillary and venous
gametocyte densities also found no difference (Sandeu et al.,
2017). Gametocyte clustering in the skin has previously been
hypothesized to potentially enhance transmission (Pichon et al.,
2000; Bousema and Drakeley, 2011; Lawniczak and Eckhoff,
2016). The effect of feeder type on the variance observed between
mosquitoes was also patient dependent. Therefore, if clustering
mechanisms do occur, they may not be specifically associated
with the skin or may be only observed post-fertilization.
There does not appear to be an association between higher
variance in skin feeding vs membrane feeding and parasitemia
or gametocytemia, noting that our study only examined six
patients. Our results warrant further examination of these
questions over a larger sample size both in terms of patients
and fed mosquitoes and also on volunteers with a wider range
of gametocytemia, including those without light microscopy
detectable gametocytemia. Indeed, pre-fertilization behaviors
might only be apparent at lower gametocyte concentration
when the likelihood of infection is low. Additionally, associative
behaviors between the sexes that ensure the presence of a male
and a female in a blood meal but don't alter local gametocyte
densities, for instance syzygy as observed in a species of
*Leucocytozoon* might also be occurring (Barraclough et al., 2008).

Overall the skin has been implicated as an organ that can
enhance vector-borne pathogens, including malaria transmission
to (Chardome and Janssen, 1952; Van Den Berghe et al., 1954)
and from the mosquito (Gueirard et al., 2010) as well as in
trypanosomes (Capewell et al., 2016). Conducting a full-scale
investigation of the role of this organ in pathogen transmission
is an important endeavor and will be facilitated by sensitive tools
such as the GQA and may allow a deeper understanding of the
infectious reservoir of malaria and how to eliminate it.

DATA AVAILABILITY STATEMENT

All datasets generated for this study are included in the
article/Supplementary Material.

ETHICS STATEMENT

The studies involving human participants were reviewed
and approved by the Faculty of Medicine, Pharmacy and
Odontostomatology de Pharmacy; Université des Sciences,
Techniques et Technologies of Bamako (IRB approval letter
no. 2016/133/CE/FMPOS). Written informed consent to
participate in this study was provided by the participants' legal
guardian/next of kin.

AUTHOR CONTRIBUTIONS

AT, DO, AD, and ML designed the study. AT, DO, CM, AH,
ND, DS, AS, MC, FD, and CS conducted the field work.
AT and KL performed the molecular biology. VH provided
statistical expertise. AT, VH, and ML wrote the manuscript
with contributions from other authors. All authors read and
approved the manuscript.

FUNDING

The Wellcome Sanger Institute is funded by the Wellcome Trust (grant 206194/Z/17/Z). ML was supported by an MRC Career Development Award (G1100339).

ACKNOWLEDGMENTS

We thank all the participants in the study.

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online
at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.
00246/full#supplementary-material

TABLE S1 | List of genes tested as gametocyte biomarkers.

TABLE S2 | Primer sequences in the multiplex assay.

TABLE S3 | Primer efficiencies.

TABLE S4 | Full qPCR data of skin and membrane feeds.
REFERENCES

Barraclough, R. K., Duval, L., Talman, A. M., Ariey, F., and Robert, V. (2008). Attraction between sexes: male-female gametocyte behaviour within a *Leucocytozoon toddi* (Haemosporida). *Parasitol. Res.* 102, 1321–1327. doi: 10.1007/s00436-008-0913-8

Bonnet, S., Gouagna, C., Safeukui, I., Meunier, J.-Y., and Boudin, C. (2000). Comparison of artificial membrane feeding with direct skin feeding to estimate infectiousness of *Plasmodium falciparum* gametocyte carriers to mosquitoes. *Trans. R. Soc. Trop. Med. Hyg.* 94, 103–106. doi: 10.1016/s0035-9203(00)90456-5

Bousema, T., and Drakeley, C. (2011). Epidemiology and infectivity of *Plasmodium falciparum* and *Plasmodium vivax* gametocytes in relation to malaria control and elimination. *Clin. Microbiol. Rev.* 5:17716.

Capewell, P., Cren-Travaillé, C., Marchesi, F., Johnston, P., Clucas, C., Benson, R. A., et al. (2016). The skin is a significant but overlooked anatomical reservoir for vector-borne African trypanosomes. *Elife* 5:e17716. doi: 10.7554/eLife.17716

Chardome, M., and Janssen, P. J. (1952). Enquete sur l'incidence malarienne par la methode dermique dans la region du Lubilash(Congo Belge). *Ann. Soc. Belg. Med. Trop.* 32, 209–211.

Diallo, M., Touré, A. M., Traoré, S. F., Niaré, O., Kassambara, L., Konaré, A., et al. (2008). Evaluation and optimization of membrane feeding compared to direct feeding as an assay for infectivity. *Malar. J.* 7:248.

Fox, J., and Weisberg, S. (2018). *An R Companion to Applied Regression*. Thousand Oaks, CA: SAGE Publications.

Gueirard, P., Tavares, J., Thiberge, S., Bernex, F., Ishino, T., Milon, G., et al. (2010). Development of the malaria parasite in the skin of the mammalian host. *Proc. Natl. Acad. Sci. U.S.A.* 107, 18640–18645. doi: 10.1073/pnas.1009346107

Lasonder, E., Rijpma, S. R., van Schaijk, B. C. L., Hoeijmakers, W. A. M., Kensche, P. R., Gresnigt, M. S., et al. (2016). Integrated transcriptomic and proteomic analyses of *P. falciparum* gametocytes: molecular insight into sex-specific processes and translational repression. *Nucleic Acids Res.* 44, 6087–6101. doi: 10.1093/nar/gkw536

Lawniczak, M. K. N., and Eckhoff, P. A. (2016). A computational lens for sexual-stage transmission, reproduction, fitness and kinetics in *Plasmodium falciparum*. *Malar. J.* 15, 487.

López-Barragán, M. J., Lemieux, J., Quiñones, M., Williamson, K. C., Molina-Cruz, A., Cui, K., et al. (2011). Directional gene expression and antisense transcripts in sexual and asexual stages of *Plasmodium falciparum*. *BMC Genom.* 12:587. doi: 10.1186/1471-2164-12-587

MalariaGEN Plasmodium falciparum Community Project, (2016). Genomic epidemiology of artemisinin resistant malaria. *eLife* 5:8714. doi: 10.7554/eLife.08714

Moll, K., Ljungström, I., Perlmann, H., Scherf, A., Wahlgren, M., and Manassas, V. (2008). *Methods IN Malaria Research*. Available at: https://www.researchgate.net/profile/Abhinav_Sinha/publication/312495167_Tight_synchronisation_protocol_for_in_vitro_cultures_of_Plasmodium_falciparum/links/587f233808ae9275d4ebaab5.pdf (accessed October 01, 2016).

Nixon, C. P. (2016). *Plasmodium falciparum* gametocyte transit through the cutaneous microvasculature: a new target for malaria transmission blocking vaccines? *Hum. Vaccin. Immunother.* 12, 3189–3195. doi: 10.1080/21645515.2016.1183076

Ouédraogo, A. L., Bousema, T., Schneider, P., de Vlas, S. J., Ilboudo-Sanogo, E., Cuzin-Ouattara, N., et al. (2009). Substantial contribution of submicroscopical *Plasmodium falciparum* gametocyte carriage to the infectious reservoir in an area of seasonal transmission. *PLoS One* 4:e8410. doi: 10.1371/journal.pone.0008410

Pfaffl, M. W. (2004). "Quantification strategies in real-time PCR," in *A-Z of Quantitative PCR*, ed. S. A. Bustin, (La Jolla, CA: International University Line), 87–112.

Pichon, G., Awono-Ambene, H. P., and Robert, V. (2000). High heterogeneity in the number of *Plasmodium falciparum* gametocytes in the bloodmeal of mosquitoes fed on the same host. *Parasitology* 121(Pt 2), 115–120. doi: 10.1017/S0031182099006277

Ribaut, C., Berry, A., Chevalley, S., Reybier, K., Morlais, I., Parzy, D., et al. (2008). Concentration and purification by magnetic separation of the erythrocytic stages of all human *Plasmodium* species. *Malar. J.* 7:45. doi: 10.1186/1475-2875-7-45

Sandeu, M. M., Bayibéki, A. N., Tchioffo, M. T., Abate, L., Gimonneau, G., Awono-Ambéné, P. H., et al. (2017). Do the venous blood samples replicate malaria parasite densities found in capillary blood? A field study performed in naturally-infected asymptomatic children in Cameroon. *Malar. J.* 16:345.

Schneider, P., Bousema, J. T., Gouagna, L. C., Otieno, S., van de Vegte-Bolmer, M., Omar, S. A., et al. (2007). Submicroscopic *Plasmodium falciparum* gametocyte densities frequently result in mosquito infection. *Am. J. Trop. Med. Hyg.* 76, 470–474. doi: 10.4269/ajtmh.2007.76.470

Schneider, P., Reece, S. E., van Schaijk, B. C. L., Bousema, T., Lanke, K. H. W., Meaden, C. S. J., et al. (2015). Quantification of female and male *Plasmodium falciparum* gametocytes by reverse transcriptase quantitative PCR. *Mol. Biochem. Parasitol.* 199, 29–33. doi: 10.1016/j.molbiopara.2015.03.006

Van Den Berghe, L., Chardome, M., and Peel, E. (1954). Superiority of preparations from skin scarification over preparations of peripheral blood for the diagnosis of malaria. *An. Inst. Med. Trop.* 9, 553–562.

Watson, O. J., Slater, H. C., Verity, R., Parr, J. B., Mwandagalirwa, M. K., Tshefu, A., et al. (2017). Modelling the drivers of the spread of *Plasmodium falciparum* hrp2 gene deletions in sub-Saharan Africa. *eLife* 6:25008. doi: 10.7554/eLife.25008

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2020 Talman, Ouologuem, Love, Howick, Mulamba, Haidara, Dara, Sylla, Sacko, Coulibaly, Dao, Sangare, Djimde and Lawniczak. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.