# Measuring malaria morbidity in an area of seasonal transmission: Pyrogenic parasitemia thresholds based on a 20-year follow-up study

Marion Dollat, Cheikh Talla, Cheikh Sokhna, Fatoumata Diene Sarr,
Jean-Francois Trape, Vincent Richard

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

► To cite this version:

Marion Dollat, Cheikh Talla, Cheikh Sokhna, Fatoumata Diene Sarr, Jean-Francois Trape, et al. Measuring malaria morbidity in an area of seasonal transmission: Pyrogenic parasitemia thresholds based on a 20-year follow-up study. PLoS ONE, 2019, 14 (6), pp.e0217903. 10.1371/journal.pone.0217903. hal-02508073

HAL Id: hal-02508073

https://amu.hal.science/hal-02508073

Submitted on 17 Mar 2020

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_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Measuring malaria morbidity in an area of seasonal transmission: Pyrogenic parasitemia thresholds based on a 20-year follow-up study

Marion Dollat<sup>1,2,*</sup>, Cheikh Talla<sup>1*</sup>, Cheikh Sokhna<sup>3</sup>, Fatoumata Diene Sarr<sup>1</sup>, Jean-François Trape<sup>3</sup>, Vincent Richard<sup>1</sup>

1. Unité d'Epidémiologie des Maladies Infectieuses, Institut Pasteur de Dakar, Dakar, Sénégal, 2. Service de Maladies Infectieuses et Tropicales, Hôpital Avicenne, Assistance Publique-Hôpitaux de Paris (AP-HP), Paris, France, 3. Laboratoire de Paludologie, Institut de Recherche pour le Développement, Dakar, Sénégal

* These authors contributed equally to this work.

† dollatmarion@gmail.com

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

## Abstract

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

### OPEN ACCESS

**Citation:** Dollat M, Talla C, Sokhna C, Diene Sarr F, Trape J-F, Richard V (2019) Measuring malaria morbidity in an area of seasonal transmission: Pyrogenic parasitemia thresholds based on a 20-year follow-up study. PLoS ONE 14(6): e0217903. https://doi.org/10.1371/journal.pone.0217903

**Editor:** Thomas A. Smith, Swiss Tropical & Public Health Institute, SWITZERLAND

**Received:** July 5, 2018

**Accepted:** May 22, 2019

**Published:** June 27, 2019

**Copyright:** © 2019 Dollat et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the manuscript and its Supporting Information files.

**Funding:** The Dielmo and Ndiop Cohort was implemented by "Institut Pasteur de Dakar" and "Institut de recherche pour le développement" through their own budget funded by different sources (Institut Pasteur, Ministry of research...) since 1990. The combination of different funding was able to maintain the follow-up activities included in the annual budget of each institution.

## Introduction

Asymptomatic carriage of *P. falciparum* is frequent in areas endemic for malaria and individual diagnosis of clinical malaria attacks is still difficult. We investigated the impact of changes in malaria endemicity on the diagnostic criteria for malaria attacks in an area of seasonal malaria transmission.

## Methods

We analyzed the longitudinal data collected over 20 years from a daily survey of all inhabitants of Ndiop, a rural community in central Senegal, in a logistic regression model to investigate the relationship between the level of *Plasmodium falciparum* parasitemia and the risk of fever, with the aim of determining the best parasitemia thresholds for attributing to malaria a fever episode.

## Results

A total of 34,136 observations recorded from July 1993 to December 2013 from 850 individuals aged from 1 day to 87 years were included. *P. falciparum* asymptomatic carriage declined from 36% to 1% between 1993 and 2013. A total of 9,819 fever episodes were associated with a positive blood film for *P. falciparum*. Using age-dependent parasitemia thresholds for attributing to malaria a fever episode, we recorded 6,006 malaria attacks during the study period. Parasitemia thresholds seemed to be lower during the low-to-zero transmission season and tended to decrease with changes in control policies. The number of clinical malaria attacks was overestimated for all age groups throughout the study when all fever episodes associated with *P. falciparum* parasitemia were defined as malaria attacks.
The staff of the IRD and Pasteur Institutes of Dakar and Paris contributed to the design, healthcare, data collection, and data treatment during the 20 years of the project.

**Competing interests:** The authors have declared that no competing interests exist.

## Conclusion

Pyrogenic thresholds are particularly sensitive to changes in malaria epidemiology and are therefore an interesting tool to accurately assess the burden of malaria in the context of declining transmission.

## Introduction

In recent years, the scale-up of new treatments and effective prevention tools has led to major advances in the fight against malaria [1,2]. However, it is still difficult to precisely assess the population-scale impact of these various strategies on malaria morbidity, due to the lack of reliable surveillance data, the varying diagnosis criteria, and the limited epidemiological methods used to estimate the disease burden. Indeed, most individuals in areas endemic for malaria progressively acquire partial and labile immunity, which allows them to tolerate low to moderate levels of parasitemia without experiencing clinical symptoms [3]. Thus, the detection of parasites in the blood film from a febrile individual is not sufficient to distinguish a malaria attack from other causes of fever.

The measurement of parasite density has long been the cornerstone of the approaches to assess malaria morbidity in research and clinical trials in endemic areas [4,5]. Methods based on parasite density were developed in the early 1990s to estimate the fraction of fever cases attributable to malaria in a population [6–8]. Several studies provided evidence for an age-dependent threshold effect in the relationship between the level of parasitemia and the occurrence of fever at the individual level and showed that such pyrogenic parasitemia thresholds can be used to confirm or rule out the diagnosis of clinical malaria attack in a given area and population [7,9–11].

The level of endemicity has also been shown to critically influence the pyrogenic parasitemia thresholds [12,13]. In the current context of declining malaria in many parts of the world, changes in the acquisition of immunity and thus in the resulting levels of parasitemia associated with malaria attacks are expected. Previous studies in Dielmo, an area with intense and perennial malaria transmission in Senegal, have shown that the supervised introduction of combination therapy for first-line treatment of malaria attacks and the deployment of long-lasting insecticide-treated nets (LLINs) were associated with a dramatic decrease in parasite density levels in asymptomatic individuals and altered pyrogenic thresholds for *P. falciparum* malaria attacks in all age-groups [13]. Here, we analyze the longitudinal data collected uninterruptedly during a 20-year period in the neighboring community of Ndiop, Senegal, an area with seasonal malaria transmission typical of most Sahelo-Sudanian savannah areas of West Africa. Our objectives were to investigate morbidity evolution according to transmission decrease, and to determine the trend of pyrogenic parasitemia thresholds for diagnosing *P. falciparum* malaria attack on the road towards malaria elimination.

## Population and methods

### Study area

The study was performed in Ndiop, a village located in central Senegal approximately 290 km southeast of Dakar and 10 km from the Gambian border (15°95′N, 16°35′W). The area is characterized by an average annual rainfall of 750 millimeters, concentrated during the rainy season between June and October, followed by a dry season of seven to eight months.
A prospective longitudinal study was carried out from July 1993 to December 2013 to investi-
gate the relationship between malaria host, vector, and parasite. Both active and passive detection
of fever cases was performed and regular cross-sectional surveys of malaria prevalence conducted.

Procedures

The procedures of medical, parasitological, entomological, and epidemiological surveillance
were the same as those used in the village of Dielmo and have already been described [13,14].

Briefly, inclusion in the cohort was offered to all villagers at the beginning of the study, and
to any newcomer thereafter. During twenty years from July 1993 to December 2013, all house-
holds included in the study were visited daily except Sundays and nominative information on
the inhabitants, including the presence of fever or other symptoms (allegation of fever, asthe-
nia, headache, vomiting, diarrhea, abdominal pain, cough) were recorded. Blood tests were
performed for all suspected or confirmed cases of fever. Asymptomatic malaria carriage was
investigated by performing cross-sectional surveys for all included individuals at regular inter-
vals: every week during the first six months of the study, then every month unless rare excep-
tions, and finally three times a year from 2004 onwards: at the beginning and end of the dry
season (January or February, and May or June), and at the end of the rainy season (October or
November). Blood was taken by finger stick and 200 oil-immersion fields were examined. The
parasite/leukocyte ratio for each plasmodial species was measured. For readability and clinical
relevance reasons, parasitemia was expressed in trophozoites/µl: to overcome the lack of simul-
taneous measurement of leukocytemia, we adopted a mean standard leukocyte count of 8,000
per µl of blood for all age groups, following WHO recommendations and in accordance with
data from a preliminary investigation in the study area [15,16].

Four first-line drug regimens were successively used for malaria attack treatment during the study period: oral quinine (July 1993 to December 1994), chloroquine (January 1995 to October 2003), sulfadoxine-pyrimethamine plus amodiaquine (SP+AQ: November 2003 to May 2006), and artesunate plus amodiaquine (Artemisinin-based combination therapy or ACT: from June 2006). Until 2011, antimalarial drugs were systematically given to young children (less than 5 years old) with fever associated with a parasitemia $\ge$ 2400 trophozoites/µl, whereas the choice to treat or not with antimalarial drugs was left to the appreciation of the physician or nurse when parasitemia was lower. Only symptomatic treatment was generally given to adults and older children permanently living in the village (except for pregnant women). This attitude was justified by the willing to limit parasitic resistance to antimalarials at that time, and by the fact that a number of fever cases with positive parasitemia did not necessarily correspond to malaria, in a context where the close surveillance of the villagers allowed immediate reactivity in case of clinical worsening. The treatment policy was modified in 2011 to limit malaria transmission and since then, ACT has been systematically given to all febrile patients showing the presence of trophozoites in blood films, independent of age and parasite density.

LLINs were deployed in Ndiop in July 2008, at the time of their distribution throughout the
country. Malaria transmission averaged 43 infected bites per person per year before the intro-
duction of combination therapy, with high annual variations (maximum 169 in 2001), and
steeply decreased to very low levels after the deployment of LLINs (average of five for the
period from 2008 to 2013) (personal data).

Parasite prevalence

Parasite prevalence was measured during cross-sectional studies in all villagers enrolled in the project. Comparisons of the prevalence rate between seasons were analyzed using the Fisher exact test to investigate seasonal variations.
**Definition of case and control observations**

Febrile cases were identified by active and passive surveillance and compared with asymptomatic controls, identified during cross-sectional investigations. We included only individuals who had participated in at least four consecutive cross-sectional surveys in the period from 1993 to 1997 in the analysis, when these studies were performed monthly, or at least three surveys in one year during the subsequent years. We excluded observations performed during pregnancy and those showing the presence of two or more parasitic species. We established a baseline during the first six months of follow-up to have a reference period, keeping only observations performed before the first antimalarial intake for each individual. Definitions of case and control were the same as those used for the previous work about the Dielmo cohort [13]. Briefly, fever cases corresponded to observations for which rectal temperature was ≥ 38°C or axillary temperature ≥ 37.5°C, and control observations were those recorded during cross-sectional surveys with rectal/axillary temperature < 38°C/37.5°C and no episode of illness within 15 days before and seven days after the thick smear was performed.

**Pyrogenic threshold calculations**

The same previous method for pyrogenic threshold calculations applied to the Dielmo cohort was used [13]. Risk of fever was analyzed using logistic regression of age and *P. falciparum* parasitemia. Odds ratio estimated by the model could be used to measure the association between risk of fever and parasitemia variation. Parasitemia could affect the risk of fever as a continuous variable and also as a binary one: the existence of a threshold effect had been previously demonstrated [9,13]. Our data suggested a variation of the threshold level according to age: following the method applied to Dielmo cohort, we aimed to estimate this pyrogenic threshold in order to attribute a fever episode to malaria. To define precisely the shape and the model of this age-dependent threshold, we define it as a function of age and five parameters (a, b, c, d, e) whose values could be estimated by successive fits, as previously described [9,13]. The section below describes the steps of modeling leading us to estimate these parameters, and so the pyrogenic thresholds.

At first, bivariate analysis was performed to explore the association between risk of fever, parasitemia and age. Comparing several age groups with each other, the best fit based on deviance criterion was obtained using a series of *k* = 5 dummy variables for the following age groups: 0–23 months, 24–59 months, 5–9 years, 10–14 years, and ≥15 years. Previous studies had shown that parasitemia thresholds for attributing fever episodes to malaria decreased in relation to control policies and decrease in transmission [12,13]. The treatment period (i.e., baseline, quinine, chloroquine, SP+AQ, ACT, and ACT+LLINs periods) was therefore expected to affect the relationship between parasite density and fever, due to the impact on the reservoir and the changes in malaria transmission. This was further supported by the model: considering the probability of fever associated with parasitemia and treatment period in bivariate analysis, stratification on the treatment period revealed significantly different ORs. This effect modifier was also found in the logistic regression model. Consequently, we separately analyzed the six treatment periods. We increased the power of analysis by pooling the two shortest periods (SP+AQ and ACT), which had a similar profile, in a global period called “Bitherapy”. The January to June semester was defined as the low-to-zero transmission season and the July to December semester as the high transmission season based on entomological and climatic data (unpublished). Since Ndiop is a seasonal transmission area, the relationship between fever and parasitemia may differ from one season to another. Using the same method used for treatment period (ie stratification in bivariate analysis then integration of interaction between the effect of age and that of parasitemia in the logistic regression model), we found
that the season significantly impacted the relationship between parasitemia and fever, which led us to analyze our data for each season separately. At each modeling step an individual random effect was tested quantifying unmeasured or unmeasurable inter-individual variability. Akaike criterion (AIC) was used to compare models with and without random effect [17]. Analyses were performed using the lme4 package of R software version 3.4 [18,19].

As previously shown, the logit of the probability $\pi_{ij}$ that the individual i presents a fever episode during the observation j can be expressed as a linear function of age $z_{ik}$ (with k representing the five age-groups) and parasitemia $x_{ij}$ (model A):

$$ \operatorname{logit}(\pi_{ij}) = \beta_0 + \sum_{k=1}^{5} \beta_{1k} z_{ik} + \beta_2 f(x_{ij}) + \alpha_i \quad (\text{model A}) $$

In this model, $\beta_0$ was a constant, $\beta_1$ and $\beta_2$ the regression coefficients, and $\alpha_i$ the random-effects individual term. Looking for the best manner to describe the fever risk $f(x_{ij})$ as a continuous function of parasitemia, different functions were tested: linear, log and power of x (parasitemia). The best criterion of model selection (deviance) was obtained for the rth power function of parasitemia. The parsimonious exponent r was then tested and obtained based on the model deviance for each treatment period and within each treatment period for both transmission seasons for different values with a precision of 0.01 (S1 and S2 Figs).

The existence of a threshold effect, in addition to the previous continuous effect of parasitemia, was previously demonstrated [9,13]. It was introduced as a binary variable $s_{ij}$ in the model (model B):

$$ \operatorname{logit}(\pi_{ij}) = \beta_0 + \sum_{k=1}^{5} \beta_{1k} z_{ik} + \beta_2 (x_{ij})^r + \beta_3 s_{ij} + \alpha_i \quad (\text{model B}) $$

$s_{ij}$ took the value 0 when the jth parasitemia of the individual i was below the tested threshold and the value 1 when it was higher. We tested constant age-independent thresholds at different values and then age-dependent thresholds. LOESS (locally weighted smoothing) was used in regression analysis to describe the relationship between age and parasitemia (S1 and S2 Figs). The trend of the loess curve could be approximated by two different equations whose shape depends on five parameters (a, b, c, d, e) as previously described [9,13]. We identified a as age of maximum parasitemia, b as the highest parasitemia, c as the parasitemia at age 0, d as the level of parasitemia in the oldest adults, and e as the shape of the decrease. The first equation ($h_1$) concerned the youngest children before age a (age of maximal parasitemia), whereas the second ($h_2$) was applied to older children and adults after age a:

$$ h_1(z_i) = [z_i(2a) - z_i^2][(b-c)/a^2] + c $$

$$ h_2(z_i) = \{[a(2a) - a^2][(b-c)/a^2] + c - d\}\{\exp[-e(z_i - a)]\} + d $$

By varying these five parameters (a, b, c, d, e) on a defined set, we obtained many combinations enabling us to define a binary threshold variable. This variable was coded in 1 or 0 according to whether the parasite density was below or above the curve of $h_1$ and $h_2$, and introduced into the model B as a qualitative variable. All models were compared with each other using the AIC criterion. For each study period, the model with the smallest AIC was retained with the five corresponding estimated parameters [17].

## Definition of *P. falciparum* malaria attack

A *P. falciparum* clinical malaria attack was defined as any case with fever or fever-related symptoms for which parasitemia was higher than the threshold derived from the above model. Cases were counted separately if they occurred 15 days or more apart.
## Comparison of malaria attack definitions

We then tested four parasitemia definitions and compared the number of malaria attacks determined by each: (A) episode of illness (fever or fever-related symptoms) with parasitemia higher than the age-dependent threshold measured for the corresponding season and treatment period; (B) episode of illness with parasitemia higher than the age-dependent threshold measured for the corresponding treatment period, regardless of the season; (C) episode of illness with parasitemia higher than the constant threshold of 5,000 trophozoites/µl, frequently used in malaria endemic areas; and (D) episode of illness associated with the presence of malaria parasites, regardless of the level of parasitemia.

## Ethical considerations

The project was initially approved by the Ministry of Health of Senegal and the assembled village population. Approval was then renewed on a yearly basis with written informed consent from individuals enrolled in the project and the parents or guardians of the children enrolled. The National Ethics Committee of Senegal and *ad-hoc* committees of the Ministry of Health, The Pasteur Institutes (Dakar and Paris), and the Institut de Recherche pour le Développement (IRD) regularly performed audits.

## Results

### *P. falciparum* prevalence

We measured parasite prevalence in 41,334 blood films collected during 159 cross-sectional surveys from all present villagers, irrespective of clinical symptoms. Before the beginning of the project, a preliminary survey carried out in June 1993 (end of the dry season) showed a parasite prevalence of 17%. We observed that prevalence reached 39% on average during the quinine period (which includes two rainy seasons and one dry season) and remained high during the chloroquine period. It started to decrease after the beginning of the combination therapy period (average prevalence of 14% between November 2003 and July 2008), and fell further after the deployment of LLINs (1%). There were consistently wide variations between seasons: parasite prevalence was significantly lower during the low-to-zero malaria transmission season. This difference was found for each treatment period, except at the last years of the study upon the distribution of LLINs (Table 1).

## Pyrogenic thresholds

We excluded 176 individuals (453 observations) because of insufficient follow up, 1,246 observations performed during pregnancy, and 3,201 observations with a thick blood film positive

<table><caption>Table 1. *P. falciparum* prevalence rate (sexual forms and gametocytes) during cross-sectional surveys in Ndiop according to season and treatment period.</caption><thead><tr><th></th><th>Prevalence rate during high transmission season (%)</th><th>Prevalence rate during low-to-zero transmission season (%)</th><th>P value*</th></tr></thead><tbody><tr><td>Quinine<br>(07/1993-12/1994)</td><td>41.5</td><td>26.5</td><td>&lt;0.001</td></tr><tr><td>Chloroquine<br>(01/1995-10/2003)</td><td>23.5</td><td>16</td><td>&lt;0.001</td></tr><tr><td>Bitherapy<br>(11/2003-07/2008)</td><td>19.9</td><td>10.2</td><td>&lt;0.001</td></tr><tr><td>ACT+LLINs<br>(08/2008-12/2013)</td><td>1.6</td><td>1.3</td><td>0.41</td></tr></tbody></table>

<footnote>* *p* values were calculated using the Fisher exact test</footnote>

<footnote>https://doi.org/10.1371/journal.pone.0217903.t001</footnote>
for two or more parasitic species when determining the parasitemia thresholds used to define
a *P. falciparum* attack. A total of 34,136 observations recorded from July 1993 to December
2013 from 850 individuals aged from 1 day to 87 years were included: 22,827 observations
matched the definition for a case of fever and 11,309 for that of a control observation.

P. falciparum parasitemia during fever episodes

The proportion of fever cases harboring *P. falciparum* trophozoites and the mean level of para-
sitemia during fever episodes decreased markedly during the 20 years of the study (Fig 1). Par-
asite prevalence was lower in children below three years of age than in older children and
adults, regardless of the treatment period. The highest levels of parasitemia were observed in
three-year-old children, except at the end of the study, in which children between seven and
nine years old had the highest levels. During the 1993 high malaria transmission season, 71%
of thick blood films from febrile patients were positive for *P. falciparum*. This proportion
remained high during the quinine and chloroquine periods (84% and 77%, respectively), then
decreased during the bitherapy period (62%) and more so after the introduction of LLINs
(26%). During the low-to-zero malaria transmission season, the prevalence of fever associated
with positive parasitemia was lower: 32%, 44%, and 28% during the quinine, chloroquine, and
bitherapy periods respectively, and 4% in years following introduction of LLINs.

Asymptomatic *P. falciparum* infections

Parasite prevalence and the mean level of parasitemia during asymptomatic infections are shown for each period, season, and age group in Fig 2. At the beginning of the study, asymptomatic carriage concerned 36% of the Ndiop villagers, with a maximal prevalence of 56% among young adults (15–19 years). The proportion in the whole population gradually decreased afterwards during the quinine (27%), chloroquine (17%), and bitherapy (118%) periods. Asymptomatic carriage was very low (1%) during the most recent period (Fig 3). We observed the highest asymptomatic parasitemia levels in children between 4 and 10 years old, the highest mean being 1,200 trophozoites/µl during the chloroquine period in the rainy season for the four to six-year age group. Asymptomatic carriage was lower during the low-to-zero transmission season during the first 15 years of the study: 22% of villagers were asymptomatic carriers in the low-to-zero transmission season vs 32% in the high transmission season during the quinine period, 16% vs 19% during the chloroquine period, and 8% vs 20% during the bitherapy period. Asymptomatic carriage further declined after the introduction of LLINs and was low in both seasons: 1% during the low-to-zero transmission season vs 1% during the high transmission season.

Attributing fever episodes to *P. falciparum* malaria

Our data showed the existence of two distinct groups based on parasite density, which sug-
gested a possible existence of a parasitemia threshold between fever and no fever cases (Fig 4).
For each treatment period and within each treatment period for both transmission seasons, we
tested this threshold effect by introducing the binary variable s<sub>ij</sub> in the model: we found a sig-
nificant association for each study period. Parameters based on the lowest deviance that define
the shape and level of pyrogenic thresholds during the baseline and each season for the four
treatment periods by age are given in Table 2.

Modeling failed to identify a specific threshold for the low-to-zero transmission season of the
ACT+LLINs period, due to the small number of observations with positive parasitemia (38 obser-
vations) such that the relationship between parasitemia and risk of fever was no longer significant.
There was thus a single threshold for the ACT+LLINs period, independent of the season.
<figure><img src="image_5.png" /><figcaption>Fig 1. Age distribution of parasite prevalence rate, class of parasite density, and mean *P. falciparum* parasitemia observed during all causes of fever episodes for each study period. *Geometric mean incalculable.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217903.g001

The Fig 5 presents the age-dependent thresholds across the periods. The level of thresholds tended to decrease from the first periods to the ACT+LLINs period. The highest thresholds were observed among children between 5 and 8 years old, except in the last years of the study (ACT+LLINs period) during which the highest threshold shifted to older children. The thresholds seemed to be higher during the high transmission season than the low-to-zero transmission season especially among children and young adults.
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Age distribution of parasite prevalence rate, class of parasite density, and mean asymptomatic *P. falciparum* parasitemia in control observations for each study period. Data for the ACT+LLINs period are not displayed because levels were too low.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217903.g002

We compared the four definitions of clinical malaria attack. Considering definition A as the reference (episode of illness with parasitemia higher than the age-dependent threshold measured for the corresponding season and treatment period), definition D (episode of illness

<figure><img src="image_7.png" /><figcaption>**Fig 3.** Changes in asymptomatic carriage, Ndiop, July 1993—December 2013.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217903.g003
<figure><img src="image_8.png" /><figcaption>Fig 4. Parasite density by age, season, and period of treatment, with pyrogenic threshold determined by the model.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217903.g004

with positive parasitemia, whatever its level) considerably overestimated the number of clinical
malaria attacks for all age groups and all periods (Table 3). During the course of the study, a
total of 9,819 fever episodes were associated with a positive blood film for *P. falciparum*, within
which 6,006 could be attributed to malaria using the age-dependent parasitemia thresholds.
Overestimation was greatest during the quinine and chloroquine periods during the low-to-
zero transmission season (264% and 138%, respectively), and lowest during the ACT+LLINs
period, regardless of the season (25% and 21%). When a constant threshold of 5,000 trophozo-
ites/µl was used (definition C), the number of clinical malaria attacks for the quinine and chlo-
roquine periods was underestimated during the low-to-zero transmission season (-9% and
-6%, respectively), but overestimated during the high transmission season (+19% and +5%,
respectively), with disparities depending on age. This threshold led to an underestimation of
the number of attacks for almost all age groups for the bitherapy period, regardless of the sea-
son. Finally, the use of a season-independent threshold led to a substantial overestimation of
malaria attacks for the quinine and chloroquine periods (7% and 10%, respectively).
<table><caption>Table 2. Estimates of the parameters defining the age-dependent pyrogenic threshold according to season for each study period.</caption><thead><tr><th></th><th>Baseline period</th><th>Quinine period, Low-to-zero transmission season</th><th>Quinine period, high transmission season</th><th>Chloroquine period, Low-to-zero transmission season</th><th>Chloroquine period, high transmission season</th><th>Bitherapy period, Low-to-zero transmission season</th><th>Bitherapy period, high transmission season</th><th>ACT +LLINs period</th></tr></thead><tbody><tr><td>Number of observations</td><td>3007</td><td>1524</td><td>2202</td><td>8097</td><td>9608</td><td>2753</td><td>2132</td><td>4813</td></tr><tr><td>Number of individuals</td><td>327</td><td>312</td><td>333</td><td>580</td><td>590</td><td>453</td><td>418</td><td>472</td></tr><tr><td>Total number of case observations/case observations with positive parasitemia†</td><td>349/248 (71.1%)</td><td>80/26 (32.5%)</td><td>944/792 (83.9%)</td><td>1489/659 (44.3%)</td><td>5555/4295 (77.3%)</td><td>355/100 (28.2%)</td><td>1528/953 (62.4%)</td><td>1009/183 (18.1%)</td></tr><tr><td>Total number of control observations/case observations with positive parasitemia†</td><td>2658/957 (36.0%)</td><td>1444/322 (22.3%)</td><td>1258/405 (32.2%)</td><td>6608/1056 (16.0%)</td><td>4053/754 (18.6%)</td><td>2398/203 (8.5%)</td><td>604/122 (20.2%)</td><td>3804/41 (1.1%)</td></tr><tr><td>Exponent for the function of parasitemia (r)</td><td>0.26</td><td>0.35</td><td>0.42</td><td>0.53</td><td>0.28</td><td>0.41</td><td>0.42</td><td>0.26</td></tr><tr><td>Age of maximum parasitemia (a) ††</td><td>7</td><td>8</td><td>8</td><td>5</td><td>6</td><td>8</td><td>6</td><td>10</td></tr><tr><td>Maximum parasitemia threshold (b) †††</td><td>12 000</td><td>8500</td><td>16500</td><td>6000</td><td>10500</td><td>3000</td><td>8500</td><td>1500</td></tr><tr><td>Parasitemia threshold at year 0 (c) †††</td><td>7000</td><td>5000</td><td>11000</td><td>4500</td><td>3500</td><td>2500</td><td>200</td><td>100</td></tr><tr><td>Lowest parasitemia threshold in adults (d) †††</td><td>2500</td><td>1000</td><td>6500</td><td>1000</td><td>1000</td><td>200</td><td>100</td><td>100</td></tr><tr><td>Shape of the decrease (e)</td><td>0.09</td><td>0.09</td><td>0.05</td><td>0.09</td><td>0.07</td><td>0.06</td><td>0.09</td><td>0.07</td></tr><tr><td>Threshold effect OR (95% CI)</td><td>3.36** (1.57–7.18)</td><td>23.19** (2.63–261.59)</td><td>6.05* (1.61–39.59)</td><td>2.66*** (1.56–4.56)</td><td>2.27*** (1.58–3.29)</td><td>5.01* (1.44–17.69)</td><td>2.55* (1.10–6.20)</td><td>27.49** (3.46–606.8)</td></tr><tr><td>Continuous effect of parasitemia OR (95% CI)</td><td>1.34*** (1.27–1.43)</td><td>1.05 (0.98–1.12)</td><td>1.08*** (1.07–1.09)</td><td>1.02*** (1.02–1.02)</td><td>1.31*** (1.28–1.33)</td><td>1.05*** (1.02–1.08)</td><td>1.06*** (1.04–1.08)</td><td>1.45*** (1.24–1.70)</td></tr><tr><td>0–23 months OR (95% CI)</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>24–59 months OR (95% CI)</td><td>0.69 (0.41–1.14)</td><td>0.73 (0.38–1.40)</td><td>0.71 (0.46–1.12)</td><td>0.49*** (0.40–0.59)</td><td>0.59*** (0.48–0.72)</td><td>0.56** (0.39–0.79)</td><td>0.72 (0.45–1.16)</td><td>0.76* (0.60–0.95)</td></tr><tr><td>5–9 years OR (95% CI)</td><td>0.47** (0.28–0.80)</td><td>0.22*** (0.10–0.47)</td><td>0.52** (0.34–0.81)</td><td>0.23*** (0.18–0.28)</td><td>0.33*** (0.27–0.40)</td><td>0.19*** (0.13–0.28)</td><td>0.37*** (0.24–0.55)</td><td>0.44*** (0.35–0.56)</td></tr><tr><td>10–14 years OR (95% CI)</td><td>0.15*** (0.08–0.28)</td><td>0.13*** (0.04–0.36)</td><td>0.32*** (0.19–0.53)</td><td>0.15*** (0.12–0.19)</td><td>0.20*** (0.16–0.25)</td><td>0.12*** (0.07–0.19)</td><td>0.21*** (0.13–0.32)</td><td>0.34*** (0.26–0.45)</td></tr><tr><td>≥15 years OR (95% CI)</td><td>0.18*** (0.12–0.27)</td><td>0.12*** (0.06–0.24)</td><td>0.44*** (0.30–0.66)</td><td>0.10*** (0.09–0.12)</td><td>0.18*** (0.15–0.21)</td><td>0.11*** (0.07–0.15)</td><td>0.14*** (0.10–0.21)</td><td>0.20*** (0.16–0.25)</td></tr></tbody></table>

Abbreviations: ACT, artemisinin-based combination therapy; LLINs, long-lasting-insecticide-treated nets; CI, confidence interval; OR, odds ratio.

† *P. falciparum*, asexual forms

†† age in years

†††† parasitemia in trophozoites/µL.

* p<0.05

** p<0.01

*** p<0.001
<figure><img src="image_9.png" /><figcaption>Fig 5. Random-effect logistic regression model-derived threshold levels of *P. falciparum* parasitemia for attributing fever episodes to *P. falciparum* malaria by age, season, and period of treatment.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0217903.g005

Discussion

In this cohort of villagers living in an initially moderate and seasonal malaria transmission area, we observed that malaria epidemiology and the pyrogenic parasitemia thresholds guiding the definition of malaria attacks tend to change over time. These results corroborate those observed in Dielmo, a neighboring locality with initial intense and perennial malaria transmission and where the same treatment policies and control measures were implemented in parallel [13,19].

Parasite prevalence progressively decreased with the implementation of new malaria control strategies. It decreased in the early years of the study, possibly due to the project itself, which resulted in more administered treatments. Thereafter, it decreased when combination therapy replaced chloroquine as first-line treatment. Finally, the most important decline followed the deployment of LLINs: parasite prevalence fell to 1% by the end of the study, and gametocyte carriage to 0.2%. These changes mimic those observed in Dielmo, where mean *P. falciparum* prevalence was 69% during the first year of the study in 1990 (vs 26% during the first full year of study in Ndiop), and fell to 0.3% in 2012 [19].

Asymptomatic carriage was more frequent during the first two periods than during subsequent periods, but parasite densities were lower: high prevalence of low density carriers in all age groups could reflect robust immunity acquired before the beginning of the survey. During the study, the peak prevalence of fever with presence of trophozoites on blood film (regardless of the level of parasitemia) shifted gradually towards older children and young adults. In correlation with the dramatic decline observed in asymptomatic carriage, this result highlights the
<table><caption>Table 3. Number of *P. falciparum* malaria attacks by age, season, and treatment period, according to four definitions of malaria attacks.</caption><thead><tr><th rowspan="2">Period</th><th rowspan="2">Malaria definitions</th><th colspan="8">Age (in years)</th><th rowspan="2">Total (children)</th><th rowspan="2">Total (adults)</th><th rowspan="2">Total</th></tr><tr><th>[0-1]</th><th>[2-3]</th><th>[4-6]</th><th>[7-9]</th><th>[10-14]</th><th>[15-19]</th><th>[20-29]</th><th>[≥ 30]</th><th>[0-14]</th><th>[≥15]</th></tr></thead><tbody><tr><td rowspan="4">Baseline</td><td>A</td><td>28</td><td>33</td><td>43</td><td>30</td><td>29</td><td>21</td><td>28</td><td>41</td><td>163</td><td>90</td><td>253</td></tr><tr><td>B</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>C</td><td>29 (+3.6%)</td><td>33 (+0.0%)</td><td>45 (+4.7%)</td><td>33 (+10.0%)</td><td>31 (+6.9%)</td><td>22 (+4.8%)</td><td>27 (-3.6%)</td><td>36 (-12.2%)</td><td>171 (+4.9%)</td><td>85 (-5.6%)</td><td>256 (+1.2%)</td></tr><tr><td>D</td><td>47 (+67.9%)</td><td>45 (+36.4%)</td><td>60 (+39.5%)</td><td>45 (+50%)</td><td>50 (+72.4%)</td><td>53 (+152.4%)</td><td>68 (+142.9%)</td><td>95 (+131.7%)</td><td>247 (+51.5%)</td><td>216 (+140%)</td><td>463 (+83.0%)</td></tr><tr><td rowspan="4">Quinine period (low-to-zero transmission season)</td><td>A</td><td>1</td><td>1</td><td>2</td><td>2</td><td>0</td><td>1</td><td>3</td><td>1</td><td>6</td><td>5</td><td>11</td></tr><tr><td>B</td><td>1 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (-0.0%)</td><td>0 (-100%)</td><td>0 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (-33.3%)</td><td>1 (-0.0%)</td><td>4 (-33.3%)</td><td>4 (-20%)</td><td>8 (-27.3%)</td></tr><tr><td>C</td><td>1 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (-0.0%)</td><td>2 (-0.0%)</td><td>0 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (-33.3%)</td><td>1 (-0.0%)</td><td>6 (-0.0%)</td><td>4 (-20%)</td><td>10 (-9.1%)</td></tr><tr><td>D</td><td>6 (+500%)</td><td>5 (+400%)</td><td>3 (+50%)</td><td>3 (+50%)</td><td>5 (+∞%)</td><td>4 (+300%)</td><td>11 (+266.7%)</td><td>3 (+200%)</td><td>22 (+266.7%)</td><td>18 (+260%)</td><td>40 (+263.6%)</td></tr><tr><td rowspan="4">Quinine period (high transmission season)</td><td>A</td><td>59</td><td>102</td><td>154</td><td>104</td><td>79</td><td>20</td><td>16</td><td>17</td><td>498</td><td>53</td><td>551</td></tr><tr><td>B</td><td>59 (-0.0%)</td><td>102 (+0.0%)</td><td>154 (+0.0%)</td><td>104 (+0.0%)</td><td>82 (+3.8%)</td><td>25 (+25%)</td><td>27 (+68.7%)</td><td>38 (+123.5%)</td><td>501 (+0.6%)</td><td>90 (+69.8%)</td><td>591 (+7.3%)</td></tr><tr><td>C</td><td>70 (+18.6%)</td><td>114 (+11.8%)</td><td>171 (+11.0%)</td><td>126 (+21.2%)</td><td>96 (+21.5%)</td><td>27 (+35%)</td><td>27 (+68.7%)</td><td>25 (+47.1%)</td><td>577 (+15.9%)</td><td>79 (+49.1%)</td><td>656 (+19.1%)</td></tr><tr><td>D</td><td>93 (+57.6%)</td><td>131 (+28.4%)</td><td>198 (+28.6%)</td><td>147 (+41.3%)</td><td>129 (+63.3%)</td><td>60 (+200%)</td><td>57 (+256.2%)</td><td>80 (+252.9%)</td><td>698 (+40.2%)</td><td>197 (+271.7%)</td><td>895 (+62.4%)</td></tr><tr><td rowspan="4">Chloroquine period (low-to-zero transmission season)</td><td>A</td><td>26</td><td>73</td><td>131</td><td>77</td><td>92</td><td>31</td><td>13</td><td>14</td><td>399</td><td>58</td><td>457</td></tr><tr><td>B</td><td>26 (-0.0%)</td><td>73 (-0.0%)</td><td>131 (-0.0%)</td><td>76 (-1.3%)</td><td>91 (-1.1%)</td><td>30 (-3.2%)</td><td>13 (-0.0%)</td><td>14 (-0.0%)</td><td>397 (-0.5%)</td><td>57 (-1.7%)</td><td>454 (-0.7%)</td></tr><tr><td>C</td><td>26 (-0.0%)</td><td>76 (+4.1%)</td><td>133 (+1.5%)</td><td>75 (-2.6%)</td><td>86 (-6.5%)</td><td>20 (-35.5%)</td><td>9 (-30.8%)</td><td>6 (-57.1%)</td><td>396 (-0.8%)</td><td>35 (-39.7%)</td><td>431 (-5.7%)</td></tr><tr><td>D</td><td>113 (+334.6%)</td><td>163 (+123.3%)</td><td>250 (+90.8%)</td><td>166 (+115.6%)</td><td>205 (+122.8%)</td><td>83 (+167.7%)</td><td>43 (+246.2%)</td><td>65 (+364.3%)</td><td>897 (+124.8%)</td><td>191 (+229.3%)</td><td>1088 (+138.1%)</td></tr><tr><td rowspan="4">Chloroquine period (high transmission season)</td><td>A</td><td>273</td><td>509</td><td>777</td><td>576</td><td>664</td><td>298</td><td>166</td><td>216</td><td>2799</td><td>680</td><td>3479</td></tr><tr><td>B</td><td>281 (+2.9%)</td><td>538 (+5.7%)</td><td>845 (+8.8%)</td><td>660 (+14.6%)</td><td>763 (+14.9%)</td><td>342 (+14.8%)</td><td>204 (+22.9%)</td><td>238 (+10.2%)</td><td>3087 (+10.3%)</td><td>784 (+15.3%)</td><td>3871 (+11.3%)</td></tr><tr><td>C</td><td>281 (+2.9%)</td><td>549 (+7.9%)</td><td>856 (+10.2%)</td><td>659 (+14.4%)</td><td>726 (+9.3%)</td><td>301 (+1.0%)</td><td>153 (-7.8%)</td><td>124 (-3.7%)</td><td>3071 (+9.7%)</td><td>578 (-15.0%)</td><td>3649 (+4.9%)</td></tr><tr><td>D</td><td>403 (+47.6%)</td><td>664 (+30.5%)</td><td>1049 (+35.0%)</td><td>861 (+49.5%)</td><td>1063 (+60.1%)</td><td>550 (+84.6%)</td><td>393 (+136.7%)</td><td>516 (+138.9%)</td><td>4040 (+44.3%)</td><td>1459 (+114.6%)</td><td>5499 (+58.1%)</td></tr><tr><td rowspan="4">Bitherapy period (low-to-zero transmission season)</td><td>A</td><td>10</td><td>15</td><td>25</td><td>15</td><td>22</td><td>6</td><td>2</td><td>5</td><td>87</td><td>13</td><td>100</td></tr><tr><td>B</td><td>10 (-10.0%)</td><td>13 (-13.3%)</td><td>22 (-12.0%)</td><td>13 (-13.3%)</td><td>17 (-22.7%)</td><td>4 (-33.3%)</td><td>1 (-50.0%)</td><td>5 (-0.0%)</td><td>75 (-13.8%)</td><td>10 (-23.1%)</td><td>85 (-15.0%)</td></tr><tr><td>C</td><td>10 (-10.0%)</td><td>13 (-13.3%)</td><td>23 (-8.0%)</td><td>13 (-13.3%)</td><td>17 (-22.7%)</td><td>3 (-50.0%)</td><td>1 (-50.0%)</td><td>1 (-80.0%)</td><td>75 (-13.8%)</td><td>5 (-61.5%)</td><td>80 (-20.0%)</td></tr><tr><td>D</td><td>16 (+60.0%)</td><td>26 (+73.3%)</td><td>33 (+32.0%)</td><td>23 (+53.3%)</td><td>29 (+31.8%)</td><td>15 (+150%)</td><td>9 (+350%)</td><td>15 (+200%)</td><td>127 (+46.0%)</td><td>39 (+200%)</td><td>166 (+66.6%)</td></tr><tr><td rowspan="2">Bitherapy period(high transmission season)</td><td>A</td><td>69</td><td>122</td><td>218</td><td>176</td><td>168</td><td>79</td><td>49</td><td>79</td><td>753</td><td>207</td><td>960</td></tr><tr><td>B</td><td>68 (-1.4%)</td><td>121 (-0.8%)</td><td>218 (-0.0%)</td><td>175 (-0.6%)</td><td>161 (-4.2%)</td><td>76 (-3.8%)</td><td>48 (-2.0%)</td><td>74 (-6.3%)</td><td>743 (-1.3%)</td><td>198 (-4.3%)</td><td>941 (-2.0%)</td></tr></tbody></table>

(Continued)
<table><caption>**Table 3.** (Continued)</caption><thead><tr><th rowspan="2">Period</th><th rowspan="2">Malaria definitions</th><th colspan="8">Age (in years)</th><th rowspan="2">Total (children)</th><th rowspan="2">Total (adults)</th><th rowspan="2">Total</th></tr><tr><th>63 (-8.7%)</th><th>125 (+2.4%)</th><th>240 (+10.1%)</th><th>187 (+6.2%)</th><th>166 (-1.2%)</th><th>70 (-11.4%)</th><th>35 (-28.6%)</th><th>30 (-62.0%)</th></tr></thead><tbody><tr><td></td><td>C</td><td>63 (-8.7%)</td><td>125 (+2.4%)</td><td>240 (+10.1%)</td><td>187 (+6.2%)</td><td>166 (-1.2%)</td><td>70 (-11.4%)</td><td>35 (-28.6%)</td><td>30 (-62.0%)</td><td>781 (+3.7%)</td><td>135 (-34.8%)</td><td>916 (-4.6%)</td></tr><tr><td></td><td>D</td><td>100 (+44.9%)</td><td>157 (+28.7%)</td><td>290 (+33.0%)</td><td>267 (+51.7%)</td><td>261 (+55.4%)</td><td>138 (+74.7%)</td><td>87 (+77.6%)</td><td>131 (+65.8%)</td><td>1075 (+42.8%)</td><td>356 (+72.0%)</td><td>1431 (+49.1%)</td></tr><tr><td>ACT+LLINs period</td><td>A</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>(low-to-zero transmission season)</td><td>B</td><td>0</td><td>1</td><td>1</td><td>1</td><td>5</td><td>1</td><td>2</td><td>1</td><td>8</td><td>4</td><td>12</td></tr><tr><td></td><td>C</td><td>0 (-0.0%)</td><td>0 (-100%)</td><td>1 (-0.0%)</td><td>1 (-0.0%)</td><td>5 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (-0.0%)</td><td>1 (-0.0%)</td><td>7 (-12.5%)</td><td>4 (-0.0%)</td><td>11 (-8.3%)</td></tr><tr><td></td><td>D</td><td>0 (-0.0%)</td><td>1 (-0.0%)</td><td>2 (+100.0%)</td><td>1 (-0.0%)</td><td>6 (+16.7%)</td><td>2 (+100.0%)</td><td>2 (-0.0%)</td><td>1 (-0.0%)</td><td>10 (+25.0%)</td><td>5 (+25.0%)</td><td>15 (+25.0%)</td></tr><tr><td>ACT+LLINs period</td><td>A</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>(high transmission season)</td><td>B</td><td>4</td><td>12</td><td>18</td><td>21</td><td>38</td><td>33</td><td>22</td><td>35</td><td>93</td><td>90</td><td>183</td></tr><tr><td></td><td>C</td><td>4 (-0.0%)</td><td>9 (-25%)</td><td>15 (-16.7%)</td><td>20 (-4.8%)</td><td>34 (-10.5%)</td><td>20 (-39.4%)</td><td>15 (-31.8%)</td><td>15 (-57.1%)</td><td>82 (-11.8%)</td><td>50 (-44.4%)</td><td>132 (-27.9%)</td></tr><tr><td></td><td>D</td><td>13 (+225%)</td><td>14 (+16.7%)</td><td>23 (+27.8%)</td><td>24 (+14.3%)</td><td>43 (+13.2%)</td><td>40 (+21.2%)</td><td>30 (+36.4%)</td><td>35 (+0.0%)</td><td>117 (+25.8%)</td><td>105 (+16.7%)</td><td>222 (+21.3%)</td></tr></tbody></table>

A = episode of illness with parasitemia higher than the age-dependent threshold measured for the corresponding season and treatment period

B = episode of illness with parasitemia higher than the age-dependent threshold measured for the corresponding treatment period, regardless of the season

C = episode of illness with parasitemia higher than the constant threshold of 5,000 trophozoites/µl, frequently used in malaria endemic areas

D = all episodes of illness associated with the presence of malaria parasites, regardless of the level of parasitemia

https://doi.org/10.1371/journal.pone.0217903.t003

fact that reduced exposure consecutive to malaria control intensification is reflected in waning immunity in adults. It illustrates that was observed in serological studies [20].

The pyrogenic parasitemia thresholds seemed to vary substantially during the study period and to be highest during the high transmission season in all periods except the last years, although the method used here do not allow us to demonstrate a significant difference between periods and between seasons. Pyrogenic threshold is an old concept in malariology [21]. It involves a discontinuous relationship between fever and parasitemia, of which the existence was confirmed under holoendemic perennial malaria conditions in the 1990s with the development of efficient statistical methods [9]. Our results show the existence of an age-dependent threshold effect of parasitemia on the risk of fever in a mesoendemic malaria setting with seasonal transmission, but at lower levels for all age group than in holoendemic area as in Dielmo [13] (S3 Fig). Thresholds were higher in children aged 5 to 10 years than in younger children for all periods, unlike in Dielmo, where the highest levels were observed in children aged 1 to 2 years [13]. This difference could be explained by a slower acquisition of immunity, due to lower cumulated exposure to infected Anopheles bites, as the level of the parasite threshold reflects the balance between the parasite and the individuals' immune response.

Malaria control strategies seem to modify pyrogenic parasitemia threshold levels, as previously observed in Dielmo [13]. The only exception concerned the quinine period, where the pyrogenic threshold determined for the high transmission season was higher than that determined for the baseline. Several hypotheses can be put forward. At first, the short duration of the quinine period (18 months) make its analysis particularly sensitive to interannual variations. Secondly, the baseline period has been built keeping only the first treated fever episode
for each individual: subsequent episodes, with higher parasitaemia, were recorded in the quinine period. Thus, the threshold levels at baseline may be artificially low.

Pyrogenic threshold levels varied depending on the season: they tended to be lower during the low-to-zero transmission season. Similar results have been reported in areas of similar endemicity in several studies with various approaches to determine the parasitemia thresholds used for case definition of malaria attack [10,11,22,23]. The decrease in threshold levels during the dry season may be linked to a rapid decrease in anti-parasitic immunity following the interruption of transmission. Indeed, some studies have shown a rapid decrease in antibody levels against *P. falciparum* during the dry season in areas of seasonal malaria transmission [24–26].

Our study shows the persistence of a pyrogenic threshold during the ACT+LLINs period, although malaria became hypoendemic and asymptomatic carriage quite rare. This may reflect, in part, the long-term persistence of some immunity as observed in immigrants exposed to malaria after years or even decades without contact with the parasite [27]. The role of memory B cells and long-term persistence of IgG antibody-secreting cells has been proposed [27,28]. It also appears that there is a persistent antigenic stimulus that is undetectable by microscopy, requiring more sensitive tools, such as PCR, to be measured. However, a recent study in Ndiop indicated quite unfrequent sub-microscopic asymptomatic carriage in 2013 [20]. Thus, the role of such sub-microscopic infections in maintaining immunity is yet to be demonstrated [29].

Our results confirm that parasite density is a key determinant of malaria morbidity in endemic areas with seasonal transmission. Defining malaria attacks as all fever cases associated with any level of parasitemia overestimated the malaria burden for all age groups in an endemic area. This was observed during all study periods and remained true during the last years of the study, although to a lesser degree, when malaria became hypoendemic.

This study has several limitations. First, it is well known that there are rapid changes in parasite density during the same fever episode [30]. Following the methodology previously applied to the Dielmo data, we considered in the analysis only the highest level of parasitaemia when several thick blood smears were performed during the same episode to reduce this measurement bias, but we are aware that this assumption might tend to push the pyrogenic threshold up. However, since this criterion of highest parasitaemia was applied throughout the study, it does not prevent us from studying the evolution of the thresholds according to the treatment and prevention strategy. Second, the existence of genetic factors in the immune response to the parasite, as well health of individual, micro-heterogeneity in transmission, and recent exposure, is likely to result in variability of the threshold between individuals. Nevertheless, for each period, the five parameters *a*, *b*, *c*, *d*, and *e* estimated by the model remained unchanged whether the random-effect was included or excluded. It must also be acknowledged that the choice to group data by periods did not take into account dynamic changes within treatment periods, such as variations observed in annual parasite prevalence, wear and tear of the nets, or changing in clinicians appreciation to treat or not with antimalarial drugs a fever case with low parasitemia. However, implementation of new strategies at the village level seemed to be the most relevant breaking point for study epidemiological changes. The thresholds determined for Ndiop may not be extrapolated to other settings with similar endemicity: malaria transmission is spatially heterogeneous and a recent study showed the complexity of such spatial distribution in eight neighboring villages of Ndiop [31]. Finally, asymptomatic carriage may have been underestimated particularly in the last years of the study, first because optical microscopy used here is less sensitive than molecular biology tools and second because it was estimated three times a year during the 10 last years of the study rather than monthly at the beginning.
**Conclusion**

The pyrogenic parasitemia threshold model is applicable to cohort studies in the context of seasonal malaria transmission, as previously documented for intense and perennial transmission settings. Pyrogenic thresholds are not fixed and particularly sensitive to the evolution of the epidemiological profile, and are therefore an interesting tool to accurately assess the burden of malaria in the context of declining transmission.

**Supporting information**

S1 Fig. A: Parasite density by age (loess curves in solid red line), B: zoom of A on low values.
(TIF)

S2 Fig. Fitting of regression model by r. The parameter r is the exponent of the power function of parasite density used for modeling the relationship between parasitemia and fever risk as a continuous function.
(TIF)

S3 Fig. Comparison of parasitemia thresholds for attributing fever episodes to *P. falciparum* malaria between the villages of Ndiop and Dielmo, by period of treatment. The thresholds presented for Ndiop are those of the high transmission season.
(TIF)

**Acknowledgments**

We thank the villagers of Ndiop for their involvement in the project. We thank all the staff of the IRD and Pasteur Institutes of Dakar and Paris who contributed to the design, healthcare, data collection, and data treatment during the 20 years of the project, especially the field-workers and Joseph Faye who managed the data. We are grateful to Christophe Rogier and Odile Mercereau-Puijalon for useful comments on the manuscript.

**Author Contributions**

**Conceptualization:** Cheikh Sokhna, Jean-François Trape, Vincent Richard.

**Formal analysis:** Marion Dollat, Cheikh Talla, Vincent Richard.

**Investigation:** Cheikh Sokhna, Fatoumata Diene Sarr, Jean-François Trape.

**Methodology:** Marion Dollat, Cheikh Talla, Vincent Richard.

**Project administration:** Fatoumata Diene Sarr, Vincent Richard.

**Resources:** Cheikh Sokhna, Jean-François Trape.

**Software:** Marion Dollat, Cheikh Talla.

**Supervision:** Jean-François Trape, Vincent Richard.

**Writing – original draft:** Marion Dollat, Cheikh Talla, Jean-François Trape, Vincent Richard.

**References**

1. Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015 Oct 8; 526(7572):207–11. https://doi.org/10.1038/nature15535 PMID: 26375008

2. WHO | World malaria report 2017 [Internet]. WHO. [cited 2018 Jan 12]. Available from: http://www.who.int/malaria/publications/world-malaria-report-2017/report/en/
3. Laishram DD, Sutton PL, Nanda N, Sharma VL, Sobti RC, Carlton JM, et al. The complexities of malaria disease manifestations with a focus on asymptomatic malaria. Malar J. 2012 Jan 31; 11:29. https://doi.org/10.1186/1475-2875-11-29 PMID: 22289302

4. Miller MJ. Observations on the natural history of malaria in the semi-resistant West African. Trans R Soc Trop Med Hyg. 1958 Mar; 52(2):152–68. https://doi.org/10.1016/0035-9203(58)90036-1 PMID: 13543904

5. Trape JF, Peelman P, Morault-Peelman B. Criteria for diagnosing clinical malaria among a semi-immune population exposed to intense and perennial transmission. Trans R Soc Trop Med Hyg. 1985; 79(4):435–42. https://doi.org/10.1016/0035-9203(85)90054-9 PMID: 3909550

6. Richard A, Lallemant M, Trape JF, Carnevale P, Mouchet J. [Malaria in the forest region of Mayombe, People's Republic of the Congo. III. The role of malaria in general morbidity]. Ann Soc Belg Med Trop. 1988 Dec; 68(4):317–29. PMID: 3069060

7. Smith T, Schellenberg JA, Hayes R. Attributable fraction estimates and case definitions for malaria in endemic areas. Stat Med. 1994 Nov 30; 13(22):2345–58. PMID: 7855468

8. Schellenberg JR, Smith T, Alonso PL, Hayes RJ. What is clinical malaria? Finding case definitions for field research in highly endemic areas. Parasitol Today Pers Ed. 1994 Nov; 10(11):439–42.

9. Rogier C, Commenges D, Trape JF. Evidence for an age-dependent pyrogenic threshold of Plasmodium falciparum parasitemia in highly endemic populations. Am J Trop Med Hyg. 1996 Jun; 54(6):613–9. https://doi.org/10.4269/ajtmh.1996.54.613 PMID: 8686780

10. McGuinness D, Koram K, Bennett S, Wagner G, Nkrumah F, Riley E. Clinical case definitions for malaria: clinical malaria associated with very low parasite densities in African infants. Trans R Soc Trop Med Hyg. 1998 Oct; 92(5):527–31. https://doi.org/10.1016/s0035-9203(98)90902-6 PMID: 9861370

11. Dicko A, Mantel C, Kouriba B, Sagara I, Thera MA, Doumbia S, et al. Season, fever prevalence and pyrogenic threshold for malaria disease definition in an endemic area of Mali. Trop Med Int Health TM IH. 2005 Jun; 10(6):550–6. https://doi.org/10.1111/j.1365-3156.2005.01418.x PMID: 15941418

12. Mwangi TW, Ross A, Snow RW, Marsh K. Case Definitions of Clinical Malaria under Different Transmission Conditions in Kilifi District, Kenya. J Infect Dis. 2005 Jun 1; 191(11):1932–9. https://doi.org/10.1086/430006 PMID: 15871128

13. Roucher C, Rogier C, Dieye-Ba F, Sokhna C, Tall A, Trape J-F. Changing malaria epidemiology and diagnostic criteria for Plasmodium falciparum clinical malaria. PloS One. 2012; 7(9):e46188. https://doi.org/10.1371/journal.pone.0046188 PMID: 23029433

14. Trape J-F, Tall A, Sokhna C, Ly AB, Diagne N, Ndiath O, et al. The rise and fall of malaria in a West African rural community, Dielmo, Senegal, from 1990 to 2012: a 22 year longitudinal study. Lancet Infect Dis. 2014 Jun; 14(6):476–88. https://doi.org/10.1016/S1473-3099(14)07012-1 PMID: 24813159

15. Trape JF, Rogier C, Konate L, Diagne N, Bouganali H, Canque B, et al. The Dielmo project: a longitudinal study of natural malaria infection and the mechanisms of protective immunity in a community living in a holoendemic area of Senegal. Am J Trop Med Hyg. 1994 Aug; 51(2):123–37. https://doi.org/10.4269/ajtmh.1994.51.123 PMID: 8074247

16. World Health Organization, editor. Basic malaria microscopy. 2nd ed. Geneva: WHO; 2010. 2 p.

17. Akaike H (1974) A new look at the statistical model identification. IEEE Trans Automat Contr AC-19: 716–723.

18. Bates D, Mächler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw [Internet]. 2015 [cited 2018 Jan 28]; 67(1). Available from: http://www.jstatsoft.org/v67/i01/

19. R: The R Project for Statistical Computing [Internet]. [cited 2018 Jan 28]. Available from: https://www.r-project.org/

20. Perraut R, Varela M-L, Loucoubar C, Niass O, Sidibé A, Tall A, et al. Serological signatures of declining exposure following intensification of integrated malaria control in two rural Senegalese communities. PloS One. 2017; 12(6):e0179146. https://doi.org/10.1371/journal.pone.0179146 PMID: 28609450

21. Bruce-Chwatt LJ. Malaria in African infants and children in Southern Nigeria. Ann Trop Med Parasitol. 1952 Sep; 46(2):173–200. https://doi.org/10.1080/00034983.1952.11685522 PMID: 12986705

22. Ouédraogo A, Tiono AB, Diarra A, Sanon S, Yaro JB, Ouedraogo E, et al. Malaria morbidity in high and seasonal malaria transmission area of Burkina Faso. PloS One. 2013; 8(1):e50036. https://doi.org/10.1371/journal.pone.0050036 PMID: 23320064

23. Bouvier P, Rougemont A, Breslow N, Doumbo O, Delley V, Dicko A, et al. Seasonality and malaria in a west African village: does high parasite density predict fever incidence? Am J Epidemiol. 1997 May 1; 145(9):850–7. https://doi.org/10.1093/oxfordjournals.aje.a009179 PMID: 9143216

24. Akpogheneta OJ, Duah NO, Tetteh KKA, Dunyo S, Lanar DE, Pinder M, et al. Duration of naturally acquired antibody responses to blood-stage Plasmodium falciparum is age dependent and antigen specific. Infect Immun. 2008 Apr; 76(4):1748–55. https://doi.org/10.1128/IAI.01333-07 PMID: 18212081
25. Cavanagh DR, Elhassan IM, Roper C, Robinson VJ, Giha H, Holder AA, et al. A longitudinal study of type-specific antibody responses to Plasmodium falciparum merozoite surface protein-1 in an area of unstable malaria in Sudan. J Immunol Baltim Md 1950. 1998 Jul 1; 161(1):347-59.

26. Perraut R, Mercereau-Puijalon O, Diouf B, Tall A, Guillotte M, Le Scanf C, et al. Seasonal fluctuation of antibody levels to Plasmodium falciparum parasitized red blood cell-associated antigens in two Senegal- lese villages with different transmission conditions. Am J Trop Med Hyg. 2000 Jun; 62(6):746-51.
https://doi.org/10.4269/ajtmh.2000.62.746 PMID: 11304067

27. Moncunill G, Mayor A, Jiménez A, Nhabomba A, Casas-Vila N, Puyol L, et al. High antibody responses against Plasmodium falciparum in immigrants after extended periods of interrupted exposure to malaria. PloS One. 2013; 8(8):e73624. https://doi.org/10.1371/journal.pone.0073624 PMID: 23967347

28. White MT, Griffin JT, Akpogheneta O, Conway DJ, Koram KA, Riley EM, et al. Dynamics of the antibody response to Plasmodium falciparum infection in African children. J Infect Dis. 2014 Oct 1; 210(7):1115-22. https://doi.org/10.1093/infdis/jiu219 PMID: 24719471

29. Fowkes FJI, Boeuf P, Beeson JG. Immunity to malaria in an era of declining malaria transmission. Parasitology. 2016 Feb; 143(2):139-53. https://doi.org/10.1017/S0031182015001249 PMID: 26741253

30. Delley V, Bouvier P, Breslow N, Doumbo O, Sagara I, Diakite M, et al. What does a single determination of malaria parasite density mean? A longitudinal survey in Mali. Trop Med Int Health TM IH. 2000 Jun; 5 (6):404-12. PMID: 10929139

31. Espié E, Diene Sarr F, Diop F, Faye J, Richard V, Tall A, et al. Spatio-Temporal Variations in Malaria Incidence in Children Less than 10 Years Old, Health District of Sokone, Senegal, 2010-2013. PloS One. 2015; 10(9):e0137737. https://doi.org/10.1371/journal.pone.0137737 PMID: 26381623