# The Rare, the Best: Spread of Antimalarial-Resistant Plasmodium falciparum Parasites by Anopheles Mosquito Vectors

Antoine Berry, Sandie Menard, Sandrine ENsango, Luc Abate, Didier

Concordet, Majoline Tchioffo Tsapi, Xavier Iriart, Parfait H Awono-Ambéné,

Benjamin Roche, Isabelle Morlais

► To cite this version:

Antoine Berry, Sandie Menard, Sandrine E Nsango, Luc Abate, Didier Concordet, et al.. The Rare, the Best: Spread of Antimalarial-Resistant Plasmodium falciparum Parasites by Anopheles Mosquito Vectors. Microbiology Spectrum, 2021, 9 (2), 10.1128/Spectrum.00852-21. hal-03510741

HAL Id: hal-03510741

https://hal.inrae.fr/hal-03510741

Submitted on 4 Jan 2022

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>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

# The Rare, the Best: Spread of Antimalarial-Resistant *Plasmodium falciparum* Parasites by *Anopheles* Mosquito Vectors

Antoine Berry,ª,ᵇ Sandie Menard,ª Sandrine E. Nsango,ᶜ,>d Luc Abate,ᵉ Didier Concordet,f Majoline Tchioffo Tsapi,ᵉ Xavier Iriart,ª,ᵇ
Parfait H. Awono-Ambéné,g Benjamin Roche,ᵉ,ᵇ,ᵢ Isabelle Morlaisᵉ

**ABSTRACT** The emergence of resistance to antimalarials has prompted the steady switch to novel therapies for decades. Withdrawal of antimalarials, such as chloroquine in sub-Saharan Africa in the late 1990s, led to rapid declines in the prevalence of resistance markers after a few years, raising the possibility of reintroducing them for malaria treatment. Here, we provide evidence that the mosquito vector plays a crucial role in maintaining parasite genetic diversity. We followed the transmission dynamics of *Plasmodium falciparum* parasites through its vector in natural infections from gametocytes contained in the blood of asymptomatic volunteers until sporozoites subsequently developed in the mosquito salivary glands. We did not find any selection of the mutant or wild-type *pfcrt* 76 allele during development in the *Anopheles* mosquito vector. However, microsatellite genotyping indicated that minority genotypes were favored during transmission through the mosquito. The analysis of changes in the proportions of mutant and wild-type *pfcrt* 76 alleles showed that, regardless of the genotype, the less-represented allele in the gametocyte population was more abundant in mosquito salivary glands, indicating a selective advantage of the minority allele in the vector. Selection of minority genotypes in the vector would explain the persistence of drug-resistant alleles in the absence of drug pressure in areas with high malaria endemicity and high genetic diversity. Our results may have important epidemiological implications, as they predict the rapid re-emergence and spread of resistant genotypes if antimalarials that had previously selected resistant parasites are reintroduced for malaria prevention or treatment.

**IMPORTANCE** Drug selection pressure in malaria patients is the cause of the emergence of resistant parasites. Resistance imposes a fitness cost for parasites in untreated infections, so withdrawal of the drug leads to the return of susceptible parasites. Little is known about the role of the malaria vector in this phenomenon. In an experimental study conducted in Cameroon, an area of high malaria transmission, we showed that the vector did not favor the parasites based on sensitivity or resistance criteria, but it did favor the selection of minority clones. This finding shows that the vector increases the diversity of plasmodial populations and could

**Citation** Berry A, Menard S, Nsango SE, Abate L, Concordet D, Tchioffo Tsapi M, Iriart X, Awono-Ambéné PH, Roche B, Morlais I. 2021. The rare, the best: spread of antimalarial-resistant *Plasmodium falciparum* parasites by *Anopheles* mosquito vectors. *Microbiol Spectr* 9:e00852-21. https://doi.org/10.1128/Spectrum.00852-21.

**Editor** Laura A. Kirkman, Weill Cornell Medicine

**Copyright** © 2021 Berry et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

Address correspondence to Antoine Berry, berry.a@chu-тулouse.fr, or Isabelle Morlais, isabelle.morlais@ird.fr.

**Received** 12 July 2021

**Accepted** 31 August 2021

**Published** 20 October 2021

<footnote>ªInstitut Toulousain des Maladies Infectieuses et Inflammatoires (Infinity), Université Toulouse, CNRS UMR5051, INSERM UMR1291, UPS, Toulouse, France</footnote>

<footnote>ᵇCentre Hospitalier Universitaire de Toulouse, Hôpital Purpan, Service de Parasitologie-Mycologie, Toulouse, France</footnote>

<footnote>ᶜDepartment of Biological Sciences, Faculté de Médecine et des Sciences Pharmaceutiques, Université de Douala, Douala, Cameroon</footnote>

<footnote>ᵈMalaria Research Unit, Centre Pasteur du Cameroun, Yaoundé, Cameroon</footnote>

<footnote>ᵉMIVEGEC, IRD, CNRS, Université de Montpellier, Montpellier, France</footnote>

<footnote>fINTHERES, Université de Toulouse, INRA, ENVT, Toulouse, France</footnote>

<footnote>ᵉLaboratoire d'Entomologie Médicale, Organisation de Coordination pour la Lutte contre les Endémies en Afrique Centrale, Yaoundé, Cameroon</footnote>

<footnote>ᵇInternational Center for Mathematical and Computational Modeling of Complex Systems (UMI IRD/UPMC UMMISCO), Bondy, France</footnote>

<footnote n="1">¹Departamento de Etología, Fauna Silvestre y Animales de Laboratorio, Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México (UNAM), Ciudad de México, México</footnote>

<footnote>Antoine Berry and Sandie Menard contributed equally to this article as first authors. Author order was determined both alphabetically and in order of increasing seniority.</footnote>
play an important role in *falciparum* malaria epidemiology by maintaining resistant clones despite the absence of therapeutic pressure.

**KEYWORDS** malaria, *Anopheles*, transmission, *Plasmodium falciparum*, resistance, antimalarials, genetic diversity, mosquito

Antimalarial drug resistance represents one of the main obstacles for the control of malaria (1, 2). Indeed, *Plasmodium falciparum* has an outstanding ability to develop resistance to numerous drug regimens, including the most efficient ones, like artemisinins. Drug resistance to *P. falciparum* is associated with specific mutations or gene amplification that confer a fitness advantage to mutant parasites when exposed to drugs. In contrast, drug resistance imposes a fitness cost in untreated infections, and withdrawal of the drug leads to the return of sensitive parasites (3–6). Chloroquine (CQ) resistance (CQR) is related to point mutations in the *P. falciparum* chloroquine resistance transporter (*pfcrt*) gene, and the *pfcrt* K76T mutation (encoding a change of K to T at position 76) is the best-known CQR marker (7).

Malaria infections are complex under natural conditions and generally harbor multiple genotypes of parasites (8, 9). Several lines of evidence suggest that genetic diversity within infections correlates with transmission intensity (10, 11), but whether the genetic composition influences the transmission dynamics remains elusive. Transmission of malaria parasites depends on the successful development of the sexual stages, the gametocytes, within the mosquito vectors. Parasite population structure and dynamics in the insect vector are complex and influenced by numerous factors from human or vector hosts (12–17). We have previously shown that mixed-genotype infections in humans lead to lower parasite burdens in mosquitoes, supporting the existence of competitive interactions between coinfecting genotypes (17, 18). However, the importance of mosquitoes in the transmission of drug-resistant parasites in naturally acquired multiple-genotype infections is still poorly understood.

In this study, we investigated the genetic diversity of *P. falciparum* parasites in both gametocytes isolated from the blood of asymptomatic donors and salivary glands dissected from mosquitoes experimentally fed with blood from the same gametocyte donors. Samples were genotyped at the *pfcrt* 76 codon and at 7 microsatellite loci. We then assessed the role of mosquitoes in the spread of *Plasmodium* resistance to antimalarial drugs within the population. Our results shed light on the key role of the mosquito vector in shaping the antimalarial-resistant parasite population.

## RESULTS

**Characteristics of the infections.** *P. falciparum* gametocyte-containing blood from 67 naturally infected carriers was used to isolate gametocyte samples and to perform experimental feedings. The median gametocyte density of blood donors was 53 parasites/µl (minimum = 11; maximum = 2,304) (Table S1 in the supplemental material). A total of 723 mosquito salivary gland samples were *P. falciparum* positive.

**No selection of the mutant or wild-type *pfcrt* 76 alleles by the *Anopheles* mosquito vector.** We excluded from the analysis all blood donors carrying a single genotype, either the mutant or the wild-type *pfcrt* 76 allele, in the gametocyte samples and all the paired sporozoite samples. Thus, 57 gametocyte samples and 643 sporozoite (salivary gland) samples, with a mean of 11 (minimum = 1; maximum = 24) mosquitoes per blood donor, were analyzed for *pfcrt* 76 alleles (Table S1). No selection of the mutant (or wild-type) allele was observed, as the median of the proportion of mutants (67.9%; interquartile range [IQR], 33.4% to 100.0%) in gametocytes was not different from that of paired sporozoite samples (mutant allele: 60.2%; IQR, 41.3% to 79.2%; Wilcoxon signed-rank test, *P* = 0.268) (Fig. 1). By symmetry, the same observation was made for the wild-type allele (Fig. S3).

**The *pfcrt* 76 allele had an advantage for transmission from humans to the vector when present at very low frequencies.** For each parasite donor, the fraction of the *pfcrt* 76 allele obtained for the sporozoite samples was plotted against the one
<figure><img src="image_4.png" /><figcaption>**FIG 1** Box plot of mutant *pfcrt* 76 alleles in gametocyte and sporozoite samples. The dots represent the observed proportions of mutant *pfcrt* 76 alleles in the gametocyte samples and the observed average proportions of mutant *pfcrt* 76 alleles in the sporozoite samples. Box plots represent the median, 25th, and 75th percentiles and the minimum and maximum.</figcaption></figure>

obtained for the corresponding gametocyte samples (Fig. 2 and Fig. S4). The *pfcrt* 76 proportion in the sporozoite samples is expected to be similar to the one in the gametocyte sample in the absence of selection through the mosquito vector (blue line in Fig. 2 and Fig. S4). We tested this hypothesis using a logistic regression where the proportion of the mutant *pfcrt* 76 allele in the gametocyte sample was compared to that of the paired sporozoite samples. The probability of observing a higher proportion of the *pfcrt* 76 alleles in sporozoite samples than in gametocyte samples (red line in Fig. 3) was significantly different (P = 5.29e-05) from the expected 50% probability (blue line in Fig. 3) for the mutant *pfcrt* 76 allele. As the proportion of mutant *pfcrt* 76

<figure><img src="image_5.png" /><figcaption>**FIG 2** Dispersion of the proportions of mutant *pfcrt* 76 alleles in gametocyte samples and in paired sporozoite samples. The dots represent the observed proportion of mutant *pfcrt* 76 alleles in each sporozoite sample against the observed proportion of mutant *pfcrt* 76 alleles in the paired gametocyte sample. The blue line represents the expected and theoretical distribution of dots if the vector has no influence in the transmission of *pfcrt* 76 alleles. The dot size is weighted by the number of samples.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**FIG 3** Logistic regression of the probability of a higher proportion of mutant *pfcrt* 76 alleles in sporozoite samples as a function of the proportion of mutant *pfcrt* 76 alleles in gametocyte samples. The proportion of mutant *pfcrt* 76 alleles in sporozoite samples was compared to that in the paired gametocyte samples. The probability of having a higher proportion of mutant *pfcrt* 76 alleles in sporozoite samples was calculated for each gametocyte sample using a logistic regression; the donor effect was taken into account. This probability (expressed as a percentage) as a function of the proportion of mutant *pfcrt* 76 alleles in gametocyte samples was plotted (red line). The blue line represents the expected probability if the variation of the proportion of mutant *pfcrt* 76 alleles in sporozoite samples is randomly distributed (50% of a higher proportion and 50% of a lower proportion). By symmetry, the result was identical for the wild-type allele (see Fig. S3).</figcaption></figure>

alleles decreased in gametocyte samples, the probability of having a higher proportion of this mutant allele in sporozoite samples increased, while as the proportion of mutant *pfcrt* 76 alleles increased in gametocyte samples, the probability of having a lower proportion of this mutant allele in sporozoite samples increased. By symmetry and logically, the same observation was made for the wild-type allele (Fig. S5). These results suggest that the minority alleles, whether mutant or wild type, had a selective advantage within a parasite population in the transmission from humans to the vector.

Next, we restricted our analysis to experimental feedings where gametocyte samples had a proportion of mutant *pfcrt* alleles below 20%, defining these parasites as a minority variant. We computed the difference of proportions between the gametocyte and sporozoite stages, where a positive value meant an increase in the proportion of parasites with the mutant *pfcrt* allele at the sporozoite stage. Through a simple chi-square test, we showed that the proportion of parasites with the minority allele at the sporozoite stage was significantly higher than zero (mean = 40%, P < 1e-16), suggesting that minority variants had a competitive advantage (Fig. 4). A similar pattern was observed for the wild-type minority allele (mean = 21%, P = 3.5e-9) (Fig. S6).

**Contribution of parasitological factors influencing the advantage of the minority variant.** To quantify the contribution of each factor influencing the advantage of this minority variant, we designed a generalized linear model through a forward approach (increasing sequentially the complexity of the model to find the most-parsimonious one) based on an improvement of the Akaike information criterion (AIC). This model aims at explaining the difference between the proportions of parasites with the mutant *pfcrt* allele at the gametocyte and sporozoite stages according to gametocytemia, oocyst count, and allelic diversity. It turned out that the most-parsimonious model involved a positive relationship with gametocytemia (coefficient = 0.046, P = 4.49e-6) and negative relationships with the oocyst count (coefficient = -0.331, P = 0.016) and genetic diversity at the gametocyte stage (coefficient = -6.692, P = 5.68e-6) (Table 1).

**Advantage of the minority variant: a general phenomenon?** DNAs from *P. falciparum*-positive salivary glands and from gametocyte samples were successfully
<figure><img src="image_7.png" /><figcaption>**FIG 4** Distribution of the increase in the frequency of parasites with mutant *pfcrt* 76 alleles between the gametocyte and sporozoite stages for experimental infections where the proportion of parasites harboring the considered allele did not exceed 20% at the gametocyte stage. The comparison of the proportion of mutant *pfcrt* 76 alleles was restricted to experimental feedings where gametocyte samples had a proportion of mutant *pfcrt* alleles below 20%, defining these parasites as a minority variant. The difference of proportions between the gametocyte and sporozoite stages was computed, where a positive value meant an increase in the proportion of parasites with the mutant *pfcrt* allele at the sporozoite stage. The dashed line represents the average of this increase.</figcaption></figure>

genotyped at 7 microsatellite loci for 49 blood donors and 560 salivary glands (Table S1). The multiplicity of infection (MOI) in the gametocyte samples was not different from that of the paired sporozoite bulks: the median MOI was 3 (mean = 3.47) for gametocyte samples versus 4 (mean = 3.96) for the sporozoite samples (Wilcoxon signed-rank test, *P* = 0.120). At each microsatellite locus, the number of alleles varied from 4 to 16 (mean = 11.4) and was not different between gametocyte and sporozoite samples (chi-square test, *P* = 0.998). Furthermore, the distribution of alleles at each microsatellite locus was not different (Table 2). All these data indicate that the allelic diversity between the gametocyte and sporozoite samples was not different at the meta-population level. In contrast, when considering the blood donor level, we observed

<table><caption>**TABLE 1** Results of the most-parsimonious generalized linear model explaining the increase in the frequency of *pfcrt* 76 mutant alleles between the gametocyte and sporozoite stagesª</caption><thead><tr><td>Parameter</td><td>Coefficient</td><td>P value</td></tr></thead><tbody><tr><td>Intercept</td><td>9.579</td><td>1.41e-4</td></tr><tr><td>Gametocytemia</td><td>0.046</td><td>4.49e-6</td></tr><tr><td>Oocyst count</td><td>-0.331</td><td>0.016</td></tr><tr><td>MOI</td><td>-6.692</td><td>5.68e-5</td></tr></tbody></table>

<footnote>ªThe results of the most-parsimonious generalized linear model explaining the increase in the frequency of *pfcrt* 76 mutant alleles between the gametocyte and sporozoite stages where the percentage of parasites harboring *pfcrt* 76 mutant alleles does not exceed 20% at the gametocyte stage. MOI, multiplicity of infection in the gametocyte population.</footnote>
<table><caption>TABLE 2 Allele distribution in the gametocyte or sporozoite samples for each microsatellite locus</caption><thead><tr><th rowspan="2">Microsatellite locus</th><th rowspan="2">Type of sample<sup>a</sup></th><th colspan="16">Number of observed alleles for the considered allele</th><th rowspan="2">χ<sup>2</sup> test</th></tr><tr><th>132</th><th>139</th><th>142</th><th>145</th><th>147</th><th>150</th><th>153</th><th>156</th><th>159</th><th>162</th><th>165</th><th>168</th><th>171</th><th>174</th><th>177</th><th>180</th><th>183</th><th>186</th></tr></thead><tbody><tr><td rowspan="2">P20</td><td><strong>Allele</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>5</td><td>1</td><td>4</td><td>9</td><td>0</td><td>9</td><td>25</td><td>17</td><td>19</td><td>15</td><td>14</td><td>4</td><td>4</td><td>5</td><td>8</td><td>0</td><td>1</td><td>1</td><td>0.707</td></tr><tr><td></td><td>Sporo</td><td>4</td><td>2</td><td>6</td><td>5</td><td>2</td><td>7</td><td>26</td><td>19</td><td>19</td><td>15</td><td>10</td><td>5</td><td>5</td><td>6</td><td>12</td><td>1</td><td>0</td><td>0</td><td></td></tr><tr><td rowspan="2">P24</td><td><strong>Allele</strong></td><td>157</td><td>160</td><td>163</td><td>166</td><td>169</td><td>172</td><td>175</td><td>177</td><td>180</td><td>183</td><td>186</td><td>189</td><td>192</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>0</td><td>15</td><td>25</td><td>25</td><td>18</td><td>9</td><td>8</td><td>5</td><td>8</td><td>8</td><td>3</td><td>4</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>2</td><td>11</td><td>21</td><td>25</td><td>24</td><td>18</td><td>13</td><td>6</td><td>13</td><td>6</td><td>3</td><td>4</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">P21</td><td><strong>Allele</strong></td><td>65</td><td>68</td><td>71</td><td>74</td><td>77</td><td>80</td><td>83</td><td>86</td><td>92</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>3</td><td>1</td><td>26</td><td>16</td><td>5</td><td>24</td><td>14</td><td>3</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>3</td><td>3</td><td>31</td><td>19</td><td>5</td><td>29</td><td>17</td><td>5</td><td>2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">P23</td><td><strong>Allele</strong></td><td>93</td><td>96</td><td>99</td><td>102</td><td>105</td><td>108</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>6</td><td>21</td><td>42</td><td>3</td><td>0</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>7</td><td>30</td><td>45</td><td>5</td><td>1</td><td>1</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">P26</td><td><strong>Allele</strong></td><td>159</td><td>162</td><td>165</td><td>168</td><td>171</td><td>174</td><td>177</td><td>180</td><td>183</td><td>186</td><td>189</td><td>192</td><td>198</td><td>201</td><td>213</td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>2</td><td>33</td><td>27</td><td>3</td><td>1</td><td>23</td><td>22</td><td>2</td><td>2</td><td>7</td><td>3</td><td>3</td><td>0</td><td>1</td><td>1</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>3</td><td>33</td><td>32</td><td>3</td><td>3</td><td>26</td><td>26</td><td>4</td><td>2</td><td>8</td><td>5</td><td>3</td><td>3</td><td>2</td><td>0</td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">P22</td><td><strong>Allele</strong></td><td>54</td><td>57</td><td>60</td><td>63</td><td>66</td><td>69</td><td>72</td><td>75</td><td>78</td><td>81</td><td>84</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>0</td><td>1</td><td>11</td><td>9</td><td>29</td><td>21</td><td>13</td><td>17</td><td>3</td><td>6</td><td>0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>1</td><td>1</td><td>14</td><td>7</td><td>29</td><td>21</td><td>14</td><td>21</td><td>0</td><td>8</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">P25</td><td><strong>Allele</strong></td><td>93</td><td>96</td><td>99</td><td>102</td><td>105</td><td>108</td><td>111</td><td>114</td><td>117</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Gameto</td><td>1</td><td>1</td><td>14</td><td>17</td><td>29</td><td>11</td><td>17</td><td>7</td><td>3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>Sporo</td><td>6</td><td>5</td><td>10</td><td>19</td><td>36</td><td>13</td><td>17</td><td>11</td><td>5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

changes in allelic composition in sporozoite samples compared to their paired gametocyte samples (Fig. 5). We detected additional alleles in sporozoite samples for 43 of 49 feedings (87.6%), whereas for 35 feedings (71.4%), alleles identified in the gametocyte samples were not recovered in sporozoite samples. Details for each locus can be found in Fig. 5.

DISCUSSION

In this study, we investigated the transmission dynamics of malaria parasites through the mosquito vector from genetically diverse natural infections. Our results reveal that (i) the mosquito vector contributes to the spread of CQ-resistant parasites by maintaining subpatent genotypes and (ii) the minority genotype, whether wild-type or mutated, in the gametocyte population is favored during transmission through the vector.

Gametocyte carriers were identified in an area of high malaria transmission where infections are mostly asymptomatic and the complexity of infection is great (17, 18). In such infections, a high proportion of individuals carry gametocytes and contribute to parasite transmission (14, 19, 20). In this study, up to 10 coinfecting genotypes were identified in the gametocyte samples, which reflects high parasite exposure and a high level of acquired immunity (21, 22). Microscopy and even standard PCR do not allow the accurate determination of the relative densities of coinfecting genotypes, and this often represents a limitation in studies on transmission dynamics. Here, we identified microsatellite alleles in sporozoites that were not detected in gametocyte samples. These newly detected alleles were obviously already present in gametocytes but remained below the detection threshold of the quantitative PCR (qPCR). Even if we cannot preclude misamplification of low-density genotypes, our work confirms that gametocytes present at undetectable densities are infectious for the mosquito vector (17, 19).

Our data indicate that development within the mosquito favored minority, low-density genotypes, and this higher transmissibility of minority genotypes could represent a trade-off

<footnote><sup>a</sup>Gameto, gametocyte samples; Sporo, sporozoite samples. The text in bold corresponds to all the alleles observed for the considered microsatellite locus.</footnote>
<figure><img src="image_8.png" /></figure>

<table><caption>**FIG 5** Frequencies of newly detected alleles and nondetected alleles in paired sporozoites for each feeding as a function of microsatellites (x axis). The allelic composition in sporozoite samples was compared to that in their paired gametocyte samples. Each bar represents the newly detected (frequency of >0%) or the nondetected (frequency of <0%) alleles in sporozoite samples for each paired gametocyte sample. The table shows, for each microsatellite, the frequency of newly detected or nondetected alleles in all sporozoite samples compared to their frequency in the paired gametocyte samples.</caption><thead><tr><th></th><th>P20</th><th>P24</th><th>P21</th><th>P23</th><th>P26</th><th>P22</th><th>P25</th><th>All microsatellites</th></tr></thead><tbody><tr><td>Newly detected alleles (%)</td><td>27/49 (55.1%)</td><td>31/49 (63.2%)</td><td>15/49 (30.6%)</td><td>17/49 (34.7%)</td><td>24/49 (49.0%)</td><td>21/49 (42.9%)</td><td>26/49 (53.0%)</td><td>43/49 (87.6%)</td></tr><tr><td>Non-detected alleles (%)</td><td>25/49 (51.0%)</td><td>18/49 (36.7%)</td><td>5/49 (10.2%)</td><td>3/49 (6.1%)</td><td>17/49 (34.7%)</td><td>16/49 (32.6%)</td><td>17/49 (34.7%)</td><td>35/49 (71.4%)</td></tr></tbody></table>

to maintain genetic diversity in the population (23, 24). A specific selection of genotypes by the mosquito cannot explain our results, as no differences in MOIs, allelic compositions, or allelic diversities were observed between gametocytes and sporozoites. Even if the microsatellite analysis is not able to formally demonstrate that low-density genotypes are favored by vectors, because the data produced are only qualitative (presence or absence of a given allele), these results strongly strengthen this hypothesis. As observed in the human host (25), multiple balances likely determine the outcome of malaria parasite infections in the vector. Indeed, we previously reported from natural *P. falciparum* infections that genetic complexity allows parasites to escape the mosquito's immune responses (18) and that gametocytes are capable of sensing the genetic content within the infection to adapt inbreeding levels (17).

In Cameroon, chloroquine was used as a first-line treatment of uncomplicated malaria until 2002, when CQ resistance was reported throughout the country and the drug was deemed inefficient (26). Artemisinin-based combination therapies (ACTs) were adopted starting in 2004 and have been used nationwide since 2007. Molecular monitoring of CQ resistance reported the reemergence of wild-type parasites a few years after chloroquine removal; 45% of blood isolates carried the K76 allele in 2009 and 75% in 2012 (27). Here, the wild-type genotype was identified in 69% of the gametocyte samples. Our study may indicate stability in the prevalence of the wild-type genotype, but it would be necessary to genotype the asexual blood stages to further investigate this, as the genetic diversity might differ between the sexual and asexual populations (19). In any case, the increase in wild-type prevalence in Cameroon occurred to a lesser extent than in malaria-endemic countries in East
Africa (3–6). Other ecological parameters determine the distribution of resistant and sensitive parasites, even if they are mostly linked to drug pressure. The prevalence of resistant genotypes decreases during the dry, low-transmission season, when drug pressure is limited (28, 29). CQ-resistant parasites are more frequent in rural areas, where self-medication is a common practice and chloroquine is still available from drug vendors (27, 30). Moreover, a higher prevalence of wild-type genotypes has been reported in clinical infections than in asymptomatic parasite carriers, which could be suggestive of higher virulence in wild-type parasites (27, 31). Alternatively, wild-type parasites are expected to have better fitness in more-immune individuals, such as asymptomatic infections in areas of high malaria transmission (32, 33). Nonetheless, the heterogeneity of the epidemiological context in different malaria settings renders predictions about the spread of drug selection difficult.

The role of mosquitoes in the transmission dynamics of drug-resistant genotypes has already been investigated (34, 35). In field studies among sympatric human and mosquito populations from Zambia, *P. falciparum* dihydrofolate reductase (*pfdhfr*) and *pfcrt*-resistant alleles were found at much lower frequencies in mosquitoes than in blood samples, and the authors suggested that the mosquito vector contributes to recovery of drug-sensitive parasites (34, 35). In our mosquito infection experiments, the frequencies of wild-type and mutant *pfcrt* alleles were similar over all gametocyte and sporozoite samples, which indicates that the cost of drug resistance observed in the human host does not occur in the mosquito vector. As we quantified the proportions of the *pfcrt* K76 allele in our samples, we further explored the dynamics of mutant and wild-type genotypes from gametocytes to sporozoites for each experiment and found that the least represented allele in the sexual stage population is more abundant in mosquito salivary glands, regardless of its identity, K76 or 76T. This result thus indicates a selective advantage of the minority allele in the mosquito vector. Negative frequency-dependent selection that maintains genetic diversity within populations has been largely described in the literature (36, 37). This pattern of selection would explain the persistence of drug-resistant genotypes in the absence of drug pressure in our study area, where malaria endemicity is high and genetic diversity great. Unfortunately, this may also have important epidemiological consequences, as it predicts the rapid reemergence and spread of resistant genotypes if antimalarial drugs that had previously selected resistant parasites are reintroduced for malaria prevention or treatment. Chloroquine has been massively used in Africa, including in Cameroon, since the coronavirus disease 2019 (COVID-19) pandemic (38), and thus, it will be important to follow the dynamics of the *pfcrt* K76T mutation and examine whether this novel change in drug pressure has induced an adaptive response by increasing the frequency of the *pfcrt* 76T allele.

The positive association we observed between gametocytemia and the increase in the minority *pfcrt* genotype strengthens the scenario of a competitive advantage of the minority allele in the mosquito vector. Conversely, the negative relationship between the MOI at the gametocyte stage and the increase in the minority *pfcrt* allele frequency suggests cooperative interactions, as more genotypes are at low densities in complex infections. Mixed infections lead to lower parasite burdens in mosquitoes (lower oocyst counts) but contribute to infecting more mosquitoes (higher infection prevalence), thereby increasing the chances of malaria parasite transmission to the next host (17, 18). Our results based on the *pfcrt* K76T mutation suggest that genetically diverse infections create a favorable environment for a genotype that would otherwise be less fit. We then provide evidence that heterogeneous interactions in mixed-genotype infections and the outcome of an infection depend on the interplay between competitive and cooperative behaviors.

Our study showed that the vector plays an important role in maintaining parasite genetic diversity and that the transmission of the minority *pfcrt* alleles by mosquitoes is possibly an evolutionary mechanism for the survival of *P. falciparum* genotypes with lower fitness in the human host. Further studies with a larger sample size are required to confirm our findings, and the genotyping of a broad set of polymorphisms using deep sequencing will be necessary to investigate whether the selective advantage of minority genotypes
occurs. Indeed, our results may have important implications for malaria control, particularly
in the persistence of resistant strains in the absence of drug pressure.

MATERIALS AND METHODS

**Study sites and origin of samples.** The study was carried out over 2 years, from 2007 to 2008, during high-transmission (rainy) seasons. The recruitment sites were in primary schools from the Mfou district (3°43′N, 11°38′E), 26 km from the center of Yaoundé, Cameroon. This study was reviewed and approved by the Cameroonian National Ethics Committee (protocol number 039/CNE/MP/06). Naturally infected *P. falciparum* gametocyte carriers were identified by microscopy of a thick blood smear obtained by finger prick and stained in 10% Giemsa stain. Gametocyte densities were expressed as the number of gametocytes observed per 1,000 leukocytes, assuming a standard concentration of 8,000 leukocytes per microliter. Venous blood samples were drawn and used for mosquito feedings and gametocyte isolation. Children with trophozoite densities of >500 parasites per microliter were then treated with an artemisinin-based combination therapy according to national recommendations.

**Gametocyte isolation.** Gametocytes were isolated from 1 ml of serum-free blood using a MACS column (Miltenyi Biotec, Germany) as previously described (39). DNA from purified gametocytes was extracted with DNAzol (Molecular Research Center, Inc., USA) and was subjected to whole-genome amplification (WGA) using the GenomiPhi V2 DNA amplification kit (GE Health Care, Sweden). All DNA samples were frozen at -20°C until gametocyte genotyping.

**Experimental feedings.** For each mosquito feeding experiment, 2 batches of about 60 females of the local Ngousso laboratory strain of *Anopheles coluzzii* were fed with 400 µl of gametocyte-containing blood using the serum replacement procedure (18). Blood-fed mosquitoes were kept under standard insectary conditions (26°C, 80% relative humidity) with a daily 6% sucrose solution until dissections. At day 8 postinfection (p.i.), mosquitoes from 1 batch were dissected. The midguts were removed and stained in a 0.4% Mercurochrome solution, and the developed oocysts were counted by light microscopy. The prevalence of infection was defined as the proportion of infected mosquitoes among the total number of dissected mosquitoes and the infection intensity as the number of oocysts per *P. falciparum*-positive mosquito. At day 14 p.i., salivary glands from mosquitoes in the second batch were dissected and transferred into 200 µl of DNAzol. DNA extractions were performed according to the manufacturer's protocol. A *P. falciparum*-specific PCR (PF1 5'-GGAATGTTATTGCTAACAC-3' and PF2 5'-AATGAAGAGCTGTGTATC-3') was carried out on salivary gland DNA to identify sporozoite-positive samples. DNA from positive salivary glands was frozen at -20°C until sporozoite genotyping.

**Genotyping of gametocyte and sporozoite samples.** Genotyping of the *pfcrt* 76 codon (chromosome 7) was performed on both gametocyte and sporozoite DNA with a real-time PCR assay using fluorescence resonance energy transfer (FRET) hybridization probes and a melting curve analysis as previously described (40). Each run included two control DNA samples of *P. falciparum*, i.e., the DNA from the CQ-susceptible F32/Tanzania strain corresponding to the *pfcrt* K76 wild-type allele and that of the CQ-resistant FCM29/Cameroon clone carrying the *pfcrt* 76T mutant allele.

DNA from gametocytes and DNA from positive salivary glands, after being pooled for each donor, were amplified to assess genetic polymorphism at 7 microsatellite loci located on 5 different chromosomes according to Anderson et al. (41) and Annan et al. (12). PCR products were resolved on an ABI Prism 3100 DNA genetic analyzer (Applied Biosystems, Foster City, CA, USA) using 500-LIZ as the internal size standard. Alleles were read under GeneMapper software (Applied Biosystems). The maximum number of alleles at the more polymorphic locus provided the minimum number of clones per isolate and determined the multiplicity of the infection.

**Quantification of mutant and wild-type alleles.** To assess the proportions of mutant and wild-type *pfcrt* 76 alleles in each sample, we used a hybridization probe-based FRET qPCR. We first produced standards with known proportions of mutant and wild-type *pfcrt* 76 alleles from laboratory strains (Fig. S1 and S2). All assays were performed in duplicate. The relative quantification of mutant/wild-type alleles of *pfcrt* 76 was determined using the plot of the negative first derivative of the fluorescence with respect to temperature (−dF/dT) versus temperature as previously described (42, 43). The relative amounts of mutant or wild-type alleles were determined by calculating the ratio of the melting peak heights of the high (mutant) or low (wild-type) melting temperature (*T*<sub>m</sub>) curves against the sum of the heights of both mutant and wild-type melting peaks, respectively, as follows: mutant % = [mutant *T*<sub>m</sub> peak height/(mutant *T*<sub>m</sub> peak height + wild-type *T*<sub>m</sub> peak height)] × 100, and wild-type % = [wild-type *T*<sub>m</sub> peak height/(mutant *T*<sub>m</sub> peak height + wild-type *T*<sub>m</sub> peak height)] × 100.

The proportions of *pfcrt* 76 alleles were assessed in field-derived gametocyte and sporozoite samples after validating the reproducibility of the quantification method (Fig. S1 and S2). Samples were tested only once due to the small amount of DNA and the number of assays to be performed.

**Statistical analysis.** Data that were not normally distributed were displayed as median values along with interquartile ranges and were compared using the Wilcoxon signed-rank test for two-group comparisons. Proportions were compared using the χ² test or Fisher's exact test, as appropriate. The relationship between the ratio of mutant to wild-type *pfcrt* 76 alleles and premixed standards was calculated by simple linear regression, and the strength of the relationship was assessed by correlation analysis using GraphPad Prism software (GraphPad, Inc., San Diego, CA, USA). Differences in comparisons were considered statistically significant if the *P* value was ≤0.05.

All the following statistical analyses were carried out using R (version 3.2) (44) and the lme4 package (45). To analyze the variations of proportions of *pfcrt* 76 alleles, we compared the proportion of wild-type *pfcrt* 76 alleles observed in the gametocyte samples to the proportion of wild-type *pfcrt* 76 alleles
in the paired sporozoite samples, respectively. The same comparison was performed for the mutant *pfcrt* 76 alleles. The variations of proportions of the wild-type or mutant *pfcrt* 76 allele were analyzed using a logistic regression that took into account the donor effect because data were not independent (1 game-
tocyte sample corresponded to *n* sporozoite samples). The generalized linear mixed model fit by maxi-
mum likelihood (Laplace approximation) and the binomial family used the following formula: variations
of proportions of *pfcrt* 76 alleles ~ 1 + percentage in gametocyte sample + (1 | gametocyte sample).

The dynamics of *pfcrt* 76 alleles was assessed by considering infections where an allele, either mu-
tant or wild type, was present below the threshold of 20% in the gametocyte sample, i.e., it was consid-
ered a “minority” variant. Differences in the proportions of the minority parasite among the two devel-
opmental stages were checked using the chi-square test. The increase in the proportion of the minority
variant was modeled through a generalized linear model in which the infection characteristics (gameto-
cytemia, MOI, and oocyst count) were included as explanatory variables. We then used the function glm
to identify the most-parsimonious model through a forward approach and test for significant interac-
tions between explanatory variables and the minority allele increase.

SUPPLEMENTAL MATERIAL

Supplemental material is available online only.

SUPPLEMENTAL FILE 1, PDF file, 0.9 MB.

ACKNOWLEDGMENTS

We are grateful to the volunteers from Mfou for their participation in this study, and
to the teams from Mfou Hospital, IRD, OCEAC, and University Hospital of Toulouse for
fieldwork, mosquito rearing, and *P. falciparum* genotyping.

This work was supported by the Institut de Recherche pour le Développement
under scholarships awarded, by the Institut de Médecine et d'Epidémiologie Appliquée
under grant no. MEA-AO-07, and by the Agence Nationale de la Recherche under grant
no. ANR-11-BSV7-009-01.

We declare no conflicts of interest.

REFERENCES

1. Petersen I, Eastman R, Lanzer M. 2011. Drug-resistant malaria: molecular mechanisms and implications for public health. FEBS Lett 585:1551-1562. https://doi.org/10.1016/j.febslet.2011.04.042.

2. White NJ. 2004. Antimalarial drug resistance. J Clin Invest 113:1084-1092.
https://doi.org/10.1172/JCI21682.

3. Frosch AE, Laufer MK, Mathanga DP, Takala-Harrison S, Skarbinski J, Claassen CW, Dzinjalamala FK, Plowe CV. 2014. Return of widespread chloroquine-sensitive *Plasmodium falciparum* to Malawi. J Infect Dis 210: 1110-1114. https://doi.org/10.1093/infdis/jiu216.

4. Kublin JG, Cortese JF, Njunju EM, Mukadam RA, Wirima JJ, Kazembe PN, Djimde AA, Kouriba B, Taylor TE, Plowe CV. 2003. Reemergence of chloroquine-sensitive *Plasmodium falciparum* malaria after cessation of chloroquine use in Malawi. J Infect Dis 187:1870-1875. https://doi.org/10.1086/375419.

5. Laufer MK, Thesing PC, Eddington ND, Masonga R, Dzinjalamala FK, Takala SL, Taylor TE, Plowe CV. 2006. Return of chloroquine antimalarial efficacy in Malawi. N Engl J Med 355:1959-1966. https://doi.org/10.1056/NEJMoa062032.

6. Mita T, Kaneko A, Lum JK, Bwijo B, Takechi M, Zungu IL, Tsukahara T, Tanabe K, Kobayakawa T, Bjorkman A. 2003. Recovery of chloroquine sensitivity and low prevalence of the *Plasmodium falciparum* chloroquine resistance transporter gene mutation K76T following the discontinuance of chloroquine use in Malawi. Am J Trop Med Hyg 68:413-415. https://doi.org/10.4269/ajtmh.2003.68.413.

7. Fidock DA, Nomura T, Talley AK, Cooper RA, Dzekunov SM, Ferdig MT, Ursos LM, Sidhu AB, Naude B, Deitsch KW, Su XZ, Wootton JC, Roepe PD, Welllems TE. 2000. Mutations in the *P. falciparum* digestive vacuole transmembrane protein PfCRT and evidence for their role in chloroquine resistance. Mol Cell 6:861-871. https://doi.org/10.1016/S1097-2765(05)00077-8.

8. Read AF, Taylor LH. 2001. The ecology of genetically diverse infections. Science 292:1099-1102. https://doi.org/10.1126/science.1059410.

9. Zwetyenga J, Rogier C, Tall A, Fontenille D, Snounou G, Trape JF, Mercereau-Puijalon O. 1998. No influence of age on infection complexity and allelic distribution in *Plasmodium falciparum* infections in Ndiop, a Senegalese village with seasonal, mesoendemic malaria. Am J Trop Med Hyg 59:726-735. https://doi.org/10.4269/ajtmh.1998.59.726.

10. Bruce MC, Donnelly CA, Alpers MP, Galinski MR, Barnwell JW, Walliker D, Day KP. 2000. Cross-species interactions between malaria parasites in humans. Science 287:845-848. https://doi.org/10.1126/science.287.5454.845.

11. Nkhoma SC, Nair S, Al-Saai S, Ashley E, McGready R, Phyo AP, Nosten F, Anderson TJ. 2013. Population genetic correlates of declining transmission in a human pathogen. Mol Ecol 22:273-285. https://doi.org/10.1111/mec.12099.

12. Annan Z, Durand P, Ayala FJ, Arnathau C, Awono-Ambene P, Simard F, Razakandrainibe FG, Koella JC, Fontenille D, Renaud F. 2007. Population genetic structure of *Plasmodium falciparum* in the two main African vectors, *Anopheles gambiae* and *Anopheles funestus*. Proc Natl Acad Sci U S A 104:7987-7992. https://doi.org/10.1073/pnas.0702715104.

13. Boissiere A, Tchioffo MT, Bachar D, Abate L, Marie A, Nsango SE, Shahbazkia HR, Awono-Ambene PH, Levashina EA, Christen R, Morlais I. 2012. Midgut microbiota of the malaria mosquito vector *Anopheles gam*
*biae* and interactions with *Plasmodium falciparum* infection. PLoS Pathog 8:e1002742. https://doi.org/10.1371/journal.ppat.1002742.

14. Bousema T, Drakeley C. 2011. Epidemiology and infectivity of *Plasmodium falciparum* and *Plasmodium vivax* gametocytes in relation to malaria control and elimination. Clin Microbiol Rev 24:377-410. https://doi.org/10.1128/CMR.00051-10.

15. Churcher TS, Bousema T, Walker M, Drakeley C, Schneider P, Ouedraogo AL, Basanez MG. 2013. Predicting mosquito infection from *Plasmodium falciparum* gametocyte density and estimating the reservoir of infection. *Elife* 2:e00626. https://doi.org/10.7554/eLife.00626.

16. Eldering M, Morlais I, Van Gemert GJ, Van de Vegte-Bolmer M, Graumans W, Siebelink-Stoter R, Vos M, Abate L, Roeffen W, Bousema T, Levashina EA, Sauerwein RW. 2016. Variation in susceptibility of African *Plasmodium falciparum* malaria parasites to TEP1 mediated killing in *Anopheles gam*
*biae* mosquitoes. Sci Rep 6:20440. https://doi.org/10.1038/srep20440.

17. Morlais I, Nsango SE, Toussile W, Abate L, Annan Z, Tchioffo MT, Cohuet A, Awono-Ambene PH, Fontenille D, Rousset F, Berry A. 2015. *Plasmodium falciparum* mating patterns and mosquito infectivity of natural isolates of gametocytes. PLoS One 10:e0123777. https://doi.org/10.1371/journal.pone.0123777.
18. Nsango SE, Abate L, Thoma M, Pompon J, Fraiture M, Rademacher A, Berry A, Awono-Ambene PH, Levashina EA, Morlais I. 2012. Genetic clonal-ity of *Plasmodium falciparum* affects the outcome of infection in *Anopheles gambiae*. Int J Parasitol 42:589-595. https://doi.org/10.1016/j.ijpara.2012.03.008.

19. Grignard L, Goncalves BP, Early AM, Daniels RF, Tiono AB, Guelbeogo WM, Ouedraogo A, van Veen EM, Lanke K, Diarra A, Nebie I, Sirima SB, Targett GA, Volkman SK, Neafsey DE, Wirth DF, Bousema T, Drakeley C. 2018. Transmission of molecularly undetectable circulating parasite clones leads to high infection complexity in mosquitoes post feeding. Int J Parasitol 48:671-677. https://doi.org/10.1016/j.ijpara.2018.02.005.

20. Stone W, Goncalves BP, Bousema T, Drakeley C. 2015. Assessing the infectious reservoir of *falciparum* malaria: past and future. Trends Parasitol 31: 287-296. https://doi.org/10.1016/j.pt.2015.04.004.

21. Eldh M, Hammar U, Arnot D, Beck HP, Garcia A, Liljander A, Mercereau-Puijalon O, Migot-Nabias F, Mueller I, Ntoumi F, Ross A, Smith T, Sonden K, Vafa Homann M, Yman V, Felger I, Farnert A. 2020. Multiplicity of asymptomatic *Plasmodium falciparum* infections and risk of clinical malaria: a systematic review and pooled analysis of individual participant data. J Infect Dis 221:775-785. https://doi.org/10.1093/infdis/jiz510.

22. Galatas B, Bassat Q, Mayor A. 2016. Malaria parasites in the asymptomatic: looking for the hay in the haystack. Trends Parasitol 32:296-308. https://doi.org/10.1016/j.pt.2015.11.015.

23. Mackinnon MJ, Marsh K. 2010. The selection landscape of malaria parasites. Science 328:866-871. https://doi.org/10.1126/science.1185410.

24. Nkhoma SC, Trevino SG, Gorena KM, Nair S, Khoswe S, Jett C, Garcia R, Daniel B, Dia A, Terlouw DJ, Ward SA, Anderson TJC, Cheeseman IH. 2020. Co-transmission of related malaria parasite lineages shapes within-host parasite diversity. Cell Host Microbe 27:93-103.e4. https://doi.org/10.1016/j.chom.2019.12.001.

25. Deroost K, Pham TT, Opdenakker G, Van den Steen PE. 2016. The immunological balance between host and parasite in malaria. FEMS Microbiol Rev 40:208-257. https://doi.org/10.1093/femsre/fuv046.

26. Niba PTN, Nji AM, Evehe MS, Ali IM, Netongo PM, Ngwafor A, Moyeh MN, Ngum LN, Ndum OE, Acho FA, Mbu'u CM, Fosah DA, Atogho-Tiedeu B, Achonduh-Atijegbe O, Djokam-Dadjeu R, Chedjou JPK, Bigoga JD, Moukoko CEE, Ajua A, Achidi E, Tallah E, Leke RGF, Tourgordi A, Ringwald P, Alifrangis M, Mbacham WF. 2021. Drug resistance markers within an evolving efficacy of anti-malarial drugs in Cameroon: a systematic review and meta-analysis (1998-2020). Malar J 20:32. https://doi.org/10.1186/s12936-020-03543-8.

27. Menard S, Morlais I, Tahar R, Sayang C, Mayengue PI, Irart X, Benoit-Vical F, Lemen B, Magnaval JF, Awono-Ambene P, Basco LK, Berry A. 2012. Molecular monitoring of *Plasmodium falciparum* drug susceptibility at the time of the introduction of artemisinin-based combination therapy in Yaounde, Cameroon: implications for the future. Malar J 11:113. https://doi.org/10.1186/1475-2875-11-113.

28. Abdel-Muhsin AM, Mackinnon MJ, Ali E, Nassir el KA, Suleiman S, Ahmed S, Walliker D, Babiker HA. 2004. Evolution of drug-resistance genes in *Plasmodium falciparum* in an area of seasonal malaria transmission in Eastern Sudan. J Infect Dis 189:1239-1244. https://doi.org/10.1086/382509.

29. Ord R, Alexander N, Dunyo S, Hallett R, Jawara M, Targett G, Drakeley CJ, Sutherland CJ. 2007. Seasonal carriage of *pfcrt* and *pfmdr1* alleles in Gambian *Plasmodium falciparum* imply reduced fitness of chloroquine-resistant parasites. J Infect Dis 196:1613-1619. https://doi.org/10.1086/522154.

30. Mawili-Mboumba DP, Ndong Ngomo JM, Maboko F, Guiyedi V, Mourou Mbina JR, Kombila M, Bouyou Akotet MK. 2014. *pfcrt* 76T and *pfmdr1* 86Y allele frequency in *Plasmodium falciparum* isolates and use of self-medication in a rural area of Gabon. Trans R Soc Trop Med Hyg 108:729-734. https://doi.org/10.1093/trstmh/tru147.

31. Tukwasibwe S, Mugenyi L, Mbogo GW, Nankoberanyi S, Maiteki-Sebuguzi C, Joloba ML, Nsobya SL, Staedke SG, Rosenthal PJ. 2014. Differential prevalence of transporter polymorphisms in symptomatic and asymptomatic *falciparum* malaria infections in Uganda. J Infect Dis 210:154-157. https://doi.org/10.1093/infdis/jiu044.

32. Escalante AA, Smith DL, Kim Y. 2009. The dynamics of mutations associated with anti-malarial drug resistance in *Plasmodium falciparum*. Trends Parasitol 25:557-563. https://doi.org/10.1016/j.pt.2009.09.008.

33. Klein EY, Smith DL, Boni MF, Laxminarayan R. 2008. Clinically immune hosts as a refuge for drug-sensitive malaria parasites. Malar J 7:67. https://doi.org/10.1186/1475-2875-7-67.

34. Mharakurwa S, Kumwenda T, Mkulama MA, Musapa M, Chishimba S, Shiff CJ, Sullivan DJ, Thuma PE, Liu K, Agre P. 2011. Malaria antifolate resistance with contrasting *Plasmodium falciparum* dihydrofolate reductase (DHFR) polymorphisms in humans and *Anopheles* mosquitoes. Proc Natl Acad Sci U S A 108:18796-18801. https://doi.org/10.1073/pnas.1116162108.

35. Mharakurwa S, Sialumano M, Liu K, Scott A, Thuma P. 2013. Selection for chloroquine-sensitive *Plasmodium falciparum* by wild *Anopheles arabiensis* in Southern Zambia. Malar J 12:453. https://doi.org/10.1186/1475-2875-12-453.

36. Carius HJ, Little TJ, Ebert D. 2001. Genetic variation in a host-parasite association: potential for coevolution and frequency-dependent selection. Evolution 55:1136-1145. https://doi.org/10.1554/0014-3820(2001)055[1136:GVIΑΗΡ]2.0.CO;2.

37. Fitzpatrick MJ, Feder E, Rowe L, Sokolowski MB. 2007. Maintaining a behaviour polymorphism by frequency-dependent selection on a single gene. Nature 447:210-212. https://doi.org/10.1038/nature05764.

38. Belayneh A. 2020. Off-label use of chloroquine and hydroxychloroquine for COVID-19 treatment in Africa against WHO recommendation. Res Rep Trop Med 11:61-72. https://doi.org/10.2147/RRTM.S269936.

39. 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. Malar J 7:45. https://doi.org/10.1186/1475-2875-7-45.

40. Vessiere A, Berry A, Fabre R, Benoit-Vical F, Magnaval JF. 2004. Detection by real-time PCR of the *Pfcrt* T76 mutation, a molecular marker of chloroquine-resistant *Plasmodium falciparum* strains. Parasitol Res 93:5-7. https://doi.org/10.1007/s00436-004-1096-6.

41. Anderson TJ, Haubold B, Williams JT, Estrada-Franco JG, Richardson L, Molinedo R, Bockarie M, Mokili J, Mharakurwa S, French N, Whitworth J, Velez ID, Brockman AH, Nosten F, Ferreira MU, Day KP. 2000. Microsatellite markers reveal a spectrum of population structures in the malaria parasite *Plasmodium falciparum*. Mol Biol Evol 17:1467-1482. https://doi.org/10.1093/oxfordjournals.molbev.a026247.

42. Karadag A, Riminucci M, Bianco P, Cherman N, Kuznetsov SA, Nguyen N, Collins MT, Robey PG, Fisher LW. 2004. A novel technique based on a PNA hybridization probe and FRET principle for quantification of mutant genotype in fibrous dysplasia/McCune-Albright syndrome. Nucleic Acids Res 32:e63. https://doi.org/10.1093/nar/gnh059.

43. Nie H, Evans AA, London WT, Block TM, Ren XD. 2011. Quantification of complex precore mutations of hepatitis B virus by SimpleProbe real time PCR and dual melting analysis. J Clin Virol 51:234-240. https://doi.org/10.1016/j.jcv.2011.05.020.

44. R Core Team. 2020. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.

45. Bates D, Mächler M, Bolker B, Walker S. 2015. Fitting linear mixed-effects models using lme4. J Stat Soft 67:1-48. https://doi.org/10.18637/jss.v06i01.