RESEARCH

Open Access

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

# Exploring the heterogeneity of human exposure to malaria vectors in an urban setting, Bouaké, Côte d'Ivoire, using an immuno-epidemiological biomarker

Dipomin F. Traoré¹,²,³*, André B. Sagna¹,², Akré M. Adja¹,⁴, Dounin D. Zoh¹,⁴, Kouassi A. Adou¹,⁵, Kouassi N. Lingué¹, Issa Coulibaly¹, N'Cho Bertin Tchiekoi¹, Serge B. Assi¹, Anne Poinsignon¹,², Mamadou Dagnogo³ and Franck Remoue¹,²

**Abstract**

**Background:** In some African cities, urban malaria is a threat to the health and welfare of city dwellers. To improve the control of the disease, it is critical to identify neighbourhoods where the risk of malaria transmission is the highest. This study aims to evaluate the heterogeneity of malaria transmission risk in one city (Bouaké) in a West African country (Côte d'Ivoire) that presents several levels of urbanization.

**Methods:** Two cross-sectional studies were conducted in three neighbourhoods (Dar-es-Salam, Kennedy and N'gattakro) in Bouaké during both the rainy and dry seasons. Data on insecticide-treated net (ITN) use and blood samples were collected from children aged between 6 months and 15 years to determine the parasite density and the prevalence of *Plasmodium falciparum* and the level of IgG against the *Anopheles* gSG6-P1 salivary peptide, used as the biomarker of *Anopheles* bite exposure.

**Results:** The specific IgG levels to the gSG6-P1 salivary peptide in the rainy season were significantly higher compared to the dry season in all neighbourhoods studied (all *p* < 0.001). Interestingly, these specific IgG levels did not differ between neighbourhoods during the rainy season, whereas significant differences in IgG level were observed in the dry season (*p* = 0.034). ITN use could be a major factor of variation in the specific IgG level. Nevertheless, no difference in specific IgG levels to the gSG6-P1 salivary peptide was observed between children who declared “always” versus “never” sleeping under an ITN in each neighbourhood. In addition, the prevalence of *P. falciparum* in the whole population and immune responders was significantly different between neighbourhoods in each season (*p* < 0.0001).

**Conclusion:** This study highlights the high risk of malaria exposure in African urban settings and the high heterogeneity of child exposure to the *Anopheles* vector between neighbourhoods in the same city. The *Anopheles* gSG6-P1 salivary peptide could be a suitable biomarker to accurately and quantitatively assess the risk of malaria transmission in urban areas.

**Keywords:** Urban malaria, salivary biomarker of exposure, *Anopheles*, Transmission risk

<footnote>*Correspondence: francois-dipomin.traore@ird.fr;
traodipominfrancois@yahoo.fr
¹Institut Pierre Richet (IPR), Institut National de la Santé Publique (INSP),
Bouaké, Côte d'Ivoire
Full list of author information is available at the end of the article</footnote>
## Background

The fast urbanization of large African cities, the anarchic occupation of urban space and various socioeconomic conditions have major implications in the epidemiology of urban malaria. The presence of shallow waters, rice cultivation and gardening in an urban environment could lead to variation in malaria transmission [1]. In Bouaké, the second largest city of Côte d'Ivoire, several shallows have been transformed for rice farming and vegetable farming. This factor increases mosquito proliferation, in particular *Anopheles gambiae*, the major malaria vector in Africa, which is adapted to this urban environment and ensures a continuous transmission of malaria in several of the city's neighbourhoods [2]. The presence of *An. gambiae* mosquitoes depends on local conditions, which explain the considerable variability of malaria distribution. Malaria transmission may vary from one region, neighbourhood and household to another, reflecting the concept of transmission hot spots [1, 3–5]. Even though malaria transmission in urban settings is generally considered low compared to rural areas, city dwellers could be considered at high risk of severe malaria because of their low acquired immunity specific to malaria, highlighting the particular health problem of urban malaria [6, 7].

The evaluation of malaria transmission is currently based on entomological methods (human-landing catch) and on parasitological assessments in human populations. However, these methods are labour-intensive and difficult to sustain on a large scale, especially when transmission and exposure levels are low (dry season, high altitude, urban settings or after vector control) [8, 9]. The entomological methods commonly used to assess human exposure to mosquito bites do not provide a measure of the individual exposure in a given area. In addition, these methods inevitably increase the hazard of the participants' exposure to mosquito-borne infections and, therefore, cannot be used in children [10, 11].

To improve the evaluation of malaria transmission/ exposure according to the World Health Organization (WHO) recommendations, much effort is being made to develop new indicators and methods at the individual level. Over the past few decades, several studies have shown that the measurement in human populations of antibody (Ab) responses to saliva molecules of vector insect was an adequate method to assess the human exposure level to vector bites and the risk of vector-borne disease [12, 13]. Specifically, the gSG6-P1 peptide (*An. gambiae* Salivary Gland Protein-6 peptide 1) of *Anopheles* saliva has been identified as a pertinent biomarker of *Anopheles* bites [14]. This salivary peptide is specific to the *Anopheles* genus, antigenic, easy to synthesize and highly conserved between *Anopheles* mosquitoes [14]. In

particular, the human IgG response to the gSG6-P1 peptide was especially relevant as a biomarker in a context of low exposure to *Anopheles* bites, for example in urban settings and during the dry season [15, 16]. In one study, carried out in northern Senegal in 2013, this salivary biomarker was used to observe a considerable heterogeneity of human exposure to *Anopheles* between neighbouring villages in a low-transmission setting [16]. In addition, this biomarker was recently used to identify hotspots of malaria transmission in severe areas in Thailand [17]. It could therefore be applied to malaria surveillance and control (i) by assessing the level of heterogeneity of human exposure to *Anopheles* bites [18] and (ii) by evaluating the efficacy (Phase 3 study) and the effectiveness (operational level, post-implementation) of vector control strategies [19, 20]. Indeed, in low urban transmission areas in Dakar, Senegal, it was shown that human IgG responses to the gSG6-P1 peptide could be, at both the population and individual levels, a credible new alternative tool to assess the heterogeneity of exposure levels to *Anopheles* bites and malaria risk [15]. Moreover, by proving the usefulness of this biomarker for assessing the effectiveness of anti-malaria vector control in populations, it was shown that this biomarker could be used as a potential alternative to the standard entomological methods, especially in low-endemic areas and urban settings [20]. Previous results in the same areas, using this immunological biomarker, indicated that human exposure to *Anopheles* bites remained similar in both urban and rural areas, whatever the season [21]. Surprisingly, urban populations could therefore be as highly exposed to *Anopheles* bites as populations living in rural areas.

The major aim of the present study was, therefore, to explore the heterogeneity of human exposure to *Anopheles* bites in different neighbourhoods in Bouaké, Côte d'Ivoire, using the *An. gambiae* salivary biomarker (gSG6-P1). Additionally, in this context of urban exposure to the *Anopheles* vector, the potential impact of the declared use of insecticide-treated nets (ITNs) was evaluated on human-vector contact.

## Methods

### Study area

In Côte d'Ivoire, malaria transmission is stable throughout the country and presents peaks during the rainy season [22]. This study was conducted in three neighbourhoods of Bouaké: Kennedy (KEN), N'gattakro (NGA) and Dar-es-Salam (DAR). The population of Bouaké is estimated at 536,719 inhabitants: 5000 inhabitants in KEN, 14,684 in NGA and 54,992 in DAR [23]. The Bouaké region is located in a climatic transition area that shows two seasons: the dry season and the rainy season. The annual average rainfall varies between
1000 and 1200 mm [24]. The urban area is crossed by many small watercourses located 500–800 m from each other. Apart from the city centre, all the Bouaké neighbourhoods studied are crossed by humid shallows

(Fig. 1). For some years, most of these shallows have been used for gardening and as paddy fields. The neighbourhoods of NGA and DAR are both known for their paddy field areas and vegetable gardens. In contrast, in

<figure><img src="image_2.png" /><figcaption>**Fig. 1** Characterization of Bouaké neighbourhood environment. This figure shows maps of different bio-ecological aspects of Bouaké neighbourhoods (a N'gattakro; b Kennedy; c Dar-es-salam), which could be potential breeding sites of *Anopheles*. People have transformed shallows to grow rice and vegetable crops and use well water for their field activities</figcaption></figure>
KEN, both types of urban farming are found in the centre of the neighbourhood.

## Population studied

Blood samples and a sociological questionnaire concerning epidemiological data and ITN use were collected during two periods: (1) in the rainy season (August 2014) and (2) in the dry season (April 2015). Households and children (6 months to 15 years) were randomly selected by the project's sociological team. At each survey, children with an axillary temperature of 37.5 °C or higher or a recent history of fever (within 24 h of the survey) with a positive malaria rapid diagnostic test (RDT) were treated with anti-malarial drugs, as recommended by National Malaria Control Programme (NMCP). A total of 801 thick blood smear samples were collected from the children studied for parasitological measurements. Parasite density (parasitaemia) was calculated as the number of *P. falciparum* parasites per microliter of blood (PD = Number of trophozoites counted × 8000/Number of leukocytes counted); the geometric mean of parasitaemia was also calculated. In addition, capillary blood samples were collected in Microvettes® microtubes (Sarstedt®, Marnay, France) and centrifuged to obtain serums, which were stored at -20 °C for immunological study.

## Salivary peptide gSG6-P1

The gSG6-P1 peptide was designed using bioinformatics to maximize its *Anopheles* specificity and its immunogenicity, as previously described [14]. It was synthesized and purified (> 95%) by Genepep SA (Saint Jean de Védas, France). The peptide was shipped in lyophilized form and then resuspended in 0.22-µm ultra-filtered water and stored at -20 °C for later use.

## Evaluation of human IgG antibody level

ELISAs (Enzyme-Linked Immunosorbent Assay) were carried out on individual sera to measure the IgG level to the gSG6-P1 peptide as previously described [16, 21]. Briefly, Maxisorp plates (Nunc, Roskilde, Denmark) were coated with gSG6-P1 (20 µg/ml) in PBS (phosphate buffered saline). After washing (demineralized water + Tween 0.1%), each serum was incubated in duplicate at 4 °C overnight at a 1/320 dilution (in PBS with 1% Tween 20). A biotinylated mouse anti-human IgG (BD Pharmingen, San Diego, CA, USA) was incubated at a 1/4000 dilution in PBS with 1% Tween (1.5 h at 37 °C) and peroxidase conjugated ExtrAvidin (Sigam, St. Louis, MO, USA) was then added (1/20,000; 1 h at 37 °C). Colorimetric analysis was carried out using ABTS (2.2-azino-bis (3 ethylbenzthiazoline 6-sulfonic acid) diammonium; Sigma) in 50 mM citrate buffer (Sigma, pH=4, containing 0.003% H₂O₂) and the absorbance (OD) was measured at 405 nm.

Individual results were expressed as: $ΔOD = ODx - ODn$, where $ODx$ represents the mean of the individual optical density (OD) value in both wells with gSG6-P1 antigen and $ODn$ the individual OD value in a blank well containing no gSG6-P1 antigen.

The positivity threshold (PT) of the IgG level to anti-gSG6-P1 was calculated using the following formula: $PT = mean(ΔODneg) + 3SD$. The ΔODneg mean of non-*Anopheles*-exposed individuals from Bordeaux (southern France) was zero. Consequently, $PT=0$ and an exposed individual was then classified as an immune responder if the ΔOD was greater than zero.

## Statistical analysis

Data analysis was carried out using Graph Pad Prism® (Graph Pad Software, San Diego, CA, USA). Values in each group did not assume a Gaussian distribution. The nonparametric Mann-Whitney U test was then used for the comparison of IgG levels in children between age groups and also between individuals who always slept under an ITN and those who never slept under an ITN. The nonparametric Kruskal-Wallis test was used to compare the three neighbourhoods. The Dunn post-test was used for multiple comparisons between neighbourhoods. All differences were considered as significant at *p* < 0.05. Excel software was used to calculate the mean geometrical parasitaemia of *P. falciparum*. The Chi square test was used to compare *P. falciparum* prevalence between neighbourhoods in both seasons.

## Ethics statement

This study followed the ethics principles recommended by the Edinburgh revision of the Helsinki Declaration. The present study was approved by the Ethics Committee of the Côte d'Ivoire Ministry of Health (June 2014; No. 41/MSLS/CNER-dkn). Written informed consent of all the parents or guardians of children who participated in the study was obtained before inclusion.

## Results

### IgG level against gSG6-P1 salivary peptide according to neighbourhoods

To explore if different bio-ecological environments between neighbourhoods would influence individual exposure to *Anopheles* bites (Fig. 1), the specific IgG response in children was compared between the three neighbourhoods studied during the rainy and dry seasons, classically known for the periods of high and low exposure to *Anopheles* vectors, respectively. The first analysis showed that the level of specific IgG varied according to seasons and was higher in the rainy season compared to the dry season in all neighbourhoods (all *p* < 0.0001; Mann-Whitney test, data not shown).
The second analysis compared the levels in specific IgG responses of children between neighbourhoods in the rainy season (Fig. 2a) and the dry season (Fig. 2b). In the rainy season, high specific IgG responses were observed in all the three neighbourhoods studied and no statistically significant difference was observed between them (p = 0.691; Kruskal–Wallis test). In contrast, specific IgG levels in the dry season differed between the three neighbourhoods (p = 0.034; Kruskal–Wallis test) despite the very low level of IgG medians. Specific IgG levels were

significantly higher in NGA compared to KEN and DAR (p = 0.019 and p = 0.012, respectively). No significant difference of specific IgG level was observed between KEN and DAR. When taking into account the whole population studied in the dry season (responders and non-responders to the gSG6-P1 peptide), the median levels of specific IgG responses, even if very close to zero, remained different between neighbourhoods. Nevertheless, as it appeared to be difficult to differentiate the specific IgG medians at the population level between the

<figure><img src="image_3.png" /><figcaption>**Fig. 2** Comparison of IgG level with gSG6-P1 salivary peptide between neighbourhoods. Dot plots show the individual specific IgG level (ΔOD value) to gSG6-P1 in the children of the neighbourhoods studied between the seasons: in the rainy season (n = 281; a) and the dry season (n = 161; b) for the whole "immune responders and not responders to the peptide gSG6-P1" population and for the only immune responders to gSG6-P1 (n = 59; c). Bars indicate the median value in each group. The number in parentheses on the figure above the dot plots indicates values above ΔOD = 2.5. The statistical results using the nonparametric Kruskal–Wallis test are indicated between the groups studied</figcaption></figure>
different neighbourhoods, a second analysis was then set up considering only the IgG responder individuals to gSG6-P1 (i.e., IgG immune responders; Fig. 2c). A high median IgG level was observed in NGA, showing high individual values for some children, whereas very low IgG responses were detected in DAR and KEN. The specific IgG level remained significantly higher in NGA compared to KEN and DAR (p=0.021 and P=0.014, respectively) in the subpopulation of IgG immune responders.

## IgG levels against gSG6-P1 peptide in children according to age group

To evaluate the potential impact of children's age on the specific IgG level between neighbourhoods, three age groups were first defined in both the rainy and dry seasons: 0–5, 6–10 and 11–15 years old. No significant difference in specific IgG level according to the three neighbourhoods was observed in the two seasons. However, there were few children aged from 11 to 15 years (n=3; 14; 16 for DAR, KEN, NGA respectively, in the dry season) and, for this reason, they were grouped into one age group (6–15 years). Two age groups were then defined for immunological analysis: ≤5 and 6–15 years old (Additional file 1). The age-dependent IgG levels were then compared between the neighbourhoods studied in the dry and rainy seasons. The anti-gSG6-P1 IgG levels in children under 5 years old were higher in NGA than in DAR and KEN. This difference was only significant in the dry season (p=0.0145). No significant difference was observed in children over 5 years old (6–15 years) between neighbourhoods in the two seasons. These age-dependent IgG levels were also compared within neighbourhoods in the dry and rainy seasons. No significant difference was observed (Additional file 1).

## Parasitological data

The prevalence of *P. falciparum* and the geometric mean of parasitaemia were calculated in the whole population and in individuals who were positive in specific IgGs to the gSG6-P1 peptide (i.e., immune responders) during both seasons (Table 1). In the whole population and in immune responders, the prevalence differed between neighbourhoods for each season (p<0.0001, Chi square test) (Table 1). The prevalence was higher in KEN (90.58%; 90.36% in the whole population and immune responders, respectively) and NGA (92.71%; 91.40%) compared to DAR (50.56%; 56.70%) in the rainy season. Similar results were obtained in the dry season (Table 1). Moreover, the prevalence of *P. falciparum* did not vary between the two seasons in each neighbourhood. Similar trends between neighbourhoods and seasons were observed for the geometric mean of *P. falciparum* parasitaemia (Table 1).

## IgG level against gSG6-P1 peptide according to the use of insecticide-treated nets

The use of ITNs by children in the urban malaria context (low exposure to *Anopheles* bites but high risk of malaria transmission) could have a major impact on the level of human-vector contact. Indeed, the results showed that, when the different neighbourhoods were considered, the proportion of individuals who had declared that they slept under an ITN significantly varied between the neighbourhoods studied in the rainy season (Chi square test: X-squared=24, 98; df=2; p<0.0001) but not in the dry season (Table 2). Therefore, the use of ITNs could be one of the major factors of variation of specific IgG levels between the different neighbourhoods. Specific IgG responses were compared between children depending on whether

<caption>**Table 1 Characterization of prevalence and parasite to neighbourhood during the rainy season and dry season**</caption>

density of *Plasmodium falciparum* average according

<table><thead><tr><th rowspan="2">Neighbourhood</th><th colspan="3">Rainy season</th><th rowspan="2">p-value</th><th colspan="3">Dry season</th><th rowspan="2">p-value</th></tr><tr><th>Dar-es-salam</th><th>Kennedy</th><th>N'gattakro</th><th>Dar-es-salam</th><th>Kennedy</th><th>N'gattakro</th></tr></thead><tbody><tr><td>Total population</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Prevalence</td><td>56/99 (56,56%)</td><td>77/85 (90,58%)</td><td>89/96 (92,71%)</td><td>***</td><td>24/47 (51%)</td><td>53/56 (94,64%)</td><td>52/59 (88,1%)</td><td>***</td></tr><tr><td>MPD (geometric mean)</td><td>16</td><td>541</td><td>483</td><td></td><td>3</td><td>76</td><td>66</td><td></td></tr><tr><td>Responders to gSG6-P1 peptide</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Prevalence</td><td>55/97 (56,70%)</td><td>75/83 (90,36%)</td><td>85/93 (91,40%)</td><td>***</td><td>7/14 (50%)</td><td>17/18 (94,44%)</td><td>26/27 (96,30%)</td><td>***</td></tr><tr><td>MPD (geometric mean)</td><td>17</td><td>549</td><td>489</td><td></td><td>5</td><td>69</td><td>65</td><td></td></tr></tbody></table>

MPD mean parasite density

P-value: χ² test, *** very significant
<table><caption>**Table 2** Comparison of the "declared or not" use of insecticide-treated nets (ITN) between neighbourhoods according to seasons</caption><caption>Insecticide-treated Net (ITN)</caption><thead><tr><th colspan="5">Rainy season</th><th colspan="4">Dry season</th></tr><tr><th>Net use</th><th>Dar-es-salam</th><th>Kennedy</th><th>N'gattakro</th><th>p-value</th><th>Dar-es-salam</th><th>Kennedy</th><th>N'gattakro</th><th>p-value</th></tr></thead><tbody><tr><td>Yes</td><td>8 (8.88%)</td><td>32 (42.1%)</td><td>21 (23.59%)</td><td>***</td><td>26 (55.31%)</td><td>27 (48.21%)</td><td>38 (64.40%)</td><td>NS</td></tr><tr><td>No</td><td>82 (91.12%)</td><td>44 (57.9%)</td><td>68 (76.41%)</td><td></td><td>21 (44.69%)</td><td>29 (51.79%)</td><td>21 (35.6%)</td><td></td></tr></tbody></table>

NS not significant

P-value: χ² test, *** p < 0.0001

or not they slept (the night before sampling) under an ITN. This comparison was done in the whole population and in the two age groups (≤ 5 and > 5 years old). Surprisingly, the results indicated that the anti-gSG6-P1 IgG levels were similar in ITN users compared to non-ITN users in the rainy (p = 0.337 Mann-Whitney) and dry seasons (p = 0.094 Mann-Whitney) (Additional file 2). No significant results were observed in the two age groups during the two seasons between children depending on whether or not they declared they slept under an ITN (children ≤ 5: p = 0.376 in the rainy season; p = 0.435 in the dry season, and children > 5: p = 0.569 in the rainy season; p = 0.93 in the dry season). When comparing the age groups ≤ 5 versus > 5 years old between the two seasons, the results showed that there was no difference in specific IgG levels, regardless of the season (p = 0.407 in the rainy and p = 0.288 the dry season).

## Discussion

To improve malaria surveillance and control in an urban context, it is essential to evaluate the risk of transmission between neighbourhoods in the same urban area. The results of the present study indicated that human exposure to *Anopheles* bites varied between the three neighbourhoods studied depending on the season. In the rainy season, the level of specific IgG was higher than in the dry season, but it did not vary according to the neighbourhoods. In contrast, the specific IgG level varied significantly according to the neighbourhoods in the dry season. It was higher in NGA compared to KEN and DAR. The semi-urban environment of Bouaké is particular in its urban agricultural practices (rice farming and gardening), the presence of shallows, different habitat types and different population densities between neighbourhoods. Household drinking water storage practices and discharge of sewage into the streets could also create potential breeding sites of *Anopheles*, with an intensity

that can vary according to the neighbourhood. This could explain the difference observed in human exposure to malaria vectors between neighbourhoods and highlights the necessity of assessing the differences of exposure to *Anopheles* bites in this context.

In the rainy season, several new temporary breeding sites of *Anopheles* (puddles, ponds, etc.) arise with rainwater. These add to the number of permanent breeding sites that already exist in neighbourhoods. Therefore, *Anopheles* densities increase and individuals receive more bites, which would explain the higher levels of specific IgG response observed in populations in the rainy season compared to the dry season. A recent study has shown that the populations of the city of Bouaké are highly exposed to *Anopheles* bites, similar to populations of rural areas, probably linked to the semi-urban context of Bouaké city [21]. The present study confirms these results and shows high specific IgG responses in Bouaké's population during the rainy season, regardless of the neighbourhood studied. In a recently urbanized area in Kenya, a study showed a strong preference of *An. gambiae* for man-made, aquatic sites such as permanent habitats in the rainy season [25]. Individuals living in the KEN neighbourhood appear to be slightly more exposed to *Anopheles* bites compared to the DAR and NGA neighbourhoods, even if this difference is not significant. This could be explained by the fact that in KEN, rice farming and vegetable gardening are more prevalent in the centre of the neighbourhood in contrast to DAR and NGA where they are located at the periphery.

Interestingly, in the dry season, even though the median of the specific IgG level sharply decreased compared to the rainy season, it was observed that the specific IgG levels were only significantly different between neighbourhoods in the dry season. Many factors could explain this heterogeneity of human exposure in an urban area. In the dry season, individuals in NGA, and especially children under 5 years old, were more exposed to *Anopheles* bites and were at higher risk of malaria transmission than those of DAR and KEN. This could be explained by
the fact that most children in this age group usually still
sleep with their mother. Indeed, it has been shown in
Kenya that the greater exposure to malaria vector bites
occurs indoors in the early evening when LLINs (long-
lasting insecticide nets) are not used [26]. Thus, children
staying indoors with their parents a long time before
sleeping under the nets could be exposed to *Anopheles*
bites. This early biting habit of some malaria vectors
could then reduce the protective efficacy of LLINs. In
addition, the decrease of mosquito bites received by indi-
viduals and the heat during the dry season could encour-
age some people to stay outdoors longer before going
to bed. All these phenomena in the urban environment
could explain the differences in the risk of malaria trans-
mission between neighbourhoods during this hot period.
However, the present results show that the frequency of
ITN use (sleeping under ITNs) did not differ between
neighbourhoods in the dry season, suggesting that ITN
use is similar between neighbourhoods in this particu-
lar season. This difference in exposure to *Anopheles* bites
between the individuals in each neighbourhood could
be explained by an inadequate ITN use by the popula-
tions living in NGA or use of damaged ITNs (presence
of holes). In addition, it could be due to the fact that indi-
viduals go to sleep late under ITNs because of the heat or
for other sociological reasons in this neighbourhood. In
areas where individuals typically stay outside later in the
evenings, the protective effect of using an ITN or protec-
tion inside houses is nullified if they are exposed during
the evening hours outside the home [4].

Interestingly, a trend according to the density of human
population was observed. DAR's specific IgG response
appeared the lowest in both the rainy and dry seasons
(highly evident in responders to the gSG6-P1 peptide)
and was significantly lower compared to KEN and NGA.
This observation could be explained by a kind of dilution
effect of *Anopheles* bites when population size increases.
For example, individuals living in KEN could be exposed
more to *Anopheles* bites because of the low human den-
sity [4], in contrast to DAR presenting a high population
density. Several factors related to rapid and uncontrolled
population and/or household growth can have major
implications for the disease transmission patterns in sub-
Saharan African cities [4, 27], such as Bouaké.

The high level of specific IgG response observed
between neighbourhoods in the rainy season seemed to
be in accordance with the parasitological results. High
*Plasmodium* prevalence in the populations of the three
neighbourhoods was observed in the rainy season and
remained stable during the dry season. The decrease of
the specific IgG level and parasite density in the dry sea-
son compared to the rainy season was related to a con-
siderable decrease in the level of exposure to *Anopheles*

bites. This could be explained by the disappearance of the
temporary breeding sites and the decrease of *Anopheles*
density [21, 28]. The very high prevalence of *P. falcipa-*
*rum* in the neighbourhood populations studied during
the two seasons can be explained by the presence of per-
manent sites (shallows, rice crops, vegetable gardens)
during the dry season, which continue to serve as breed-
ing sites for *Anopheles*. Khaemba et al. [25] showed that
dams and swamps remained the preferred sites of *An.
gambiae* during the dry season.

The present study is a first step to explore the neigh-
bourhood-dependent heterogeneity of human exposure
to the *Anopheles* vector in an urban area. It cannot be
conclude which pertinent socioepidemiological and/or
environmental factors explain the differences of human
exposure to *Anopheles* bites between the neighbour-
hoods, specifically observed in the dry season. Multi-
ple factors may be implicated and future studies will be
necessary to precisely indicate the environmental and/
or sociological factors involved, some of which may be
specific to urban contexts, responsible for this hetero-
geneity of *Anopheles* exposure between neighbourhoods
in the same urban area. In addition, it could be interest-
ing to conduct similar studies in different African cities
presenting various population sizes (megalopolis versus
medium cities, for example) and different degrees and
histories of urbanization.

Conclusion

The use of the *An. gambiae* salivary biomarker (gSG6-P1)
showed that the human exposure to malaria vectors was
very low in the dry season compared to the rainy sea-
son, although malaria prevalence remained stable during
both seasons. In the rainy season, the level of exposure
to malaria vectors was high, but did not vary according
to neighbourhoods. Interestingly, it varied according
to neighbourhoods in the dry season. The salivary bio-
marker of human exposure to *Anopheles* bites could be
a relevant tool to orient malaria control strategies aim-
ing to target malaria control in urban neighbourhoods at
high risk of malaria.

**Additional files**

**Additional file 1.** IgG level with gSG6-P1 salivary peptide according to age. Fig A and B: IgG levels to gSG6-P1 peptide of children under 5 years according districts (Fig. A rainy season, Fig. B dry season). Fig. C and D: IgG levels to gSG6-P1 peptide of children over 5 years according districts (Fig. C rainy season, Fig. D dry season).

**Additional file 2.** IgG level with gSG6-P1 salivary peptide according to the use of insecticide-treated nets (ITN) in the whole population (A) and in age groups (B and C). Fig. A: IgG level to gSG6-P1 peptide in the whole population according to the ITN use in the rainy (Fig. A1) and dry (Fig. 1B)
seasons. Fig. B: IgG level to gSG6-P1 peptide in "under 5 years" age group according to the ITN use in the rainy (Fig. A1) and dry (Fig. 1B) seasons. Fig. C: IgG level to gSG6-P1 peptide in "over 5 years" age group according to the ITN use in the rainy (Fig. A1) and dry (Fig. 1B) seasons.

**Abbreviations**

gSG6-P1: *An. gambiae* Salivary Gland Protein-6 peptide 1; IgG: immunoglobulin G; ELISAs: Enzyme-Linked Immunosorbent Assay; ΔOD: delta optical density; ITN: insecticide-treated net; NGA: N'gattakro; KEN: Kennedy; DAR: Dar-es-salam.

**Authors' contributions**

FR, AMA, AP, SBA and MD conceived the study. DFT, ABS, TBN, DDZ, LNK, IC collected data in field surveys; DFT, ABS carried out the immunological assessments and analyzed data. AAK, DFT designed the different geographical maps of the neighbourhoods. DFT drafted the first manuscript. FR participated in preparing and writing the manuscript and revised the manuscript. ABS, AP, AMA, MD provided guidance on revision of the manuscript. FR and DFT are guarantors of the paper. All authors read and approved the final manuscript.

**Author details**

¹ Institut Pierre Richet (IPR), Institut National de la Santé Publique (INSP), Bouaké, Côte d'Ivoire. ² MIVEGEC, Univ. Montpellier, CNRS, IRD, Montpellier, France. ³ Unité de Formation et de Recherche en Sciences de la nature (UFR SN), Université Nangui Abrogoua, Abidjan, Côte d'Ivoire. ⁴ Unité de Formation et de Recherche en Biosciences (UFR Biosciences), Université Félix Houphouët Boigny, Abidjan, Côte d'Ivoire. ⁵ Institut de Géographie Tropicale (IGT), Université Félix Houphouët Boigny, Abidjan, Côte d'Ivoire.

**Acknowledgements**

We gratefully acknowledge the populations of Bouaké's health neighbourhoods (Kennedy, Dar-es-Salaam and N'Gattakro) and Pétéssou and Allokokro villages, especially householders and guardians of children for their kind support and collaboration. We thank Dr. Linda Northrup, English Solutions (France) for her editing of the manuscript.

**Competing interests**

The authors declare that they have no competing interests.

**Availability of data and materials**

The data sets generated and analyzed during the current study are not publicly available due to confidentiality concerns (data collected from patients) but are available from the corresponding author on reasonable request.

**Consent for publication**

Not applicable.

**Ethics approval and consent to participate**

This study followed the ethical principles recommended by the Edinburgh revision of the Declaration of Helsinki. The present study was approved by the Ethics Committee of the Côte d'Ivoire Ministry of Health (June 2014; No. 41/MSLS/CNER-dkn). Written informed consent of all parents or guardians of children who participated in the study was obtained before inclusion.

**Funding**

This research was integrated into the "Jeune Equipe Associée à l'IRD" EVAPALCI multidisciplinary study funded by the IRD (Institut de Recherche pour le Développement; Département des Programmes de Recherche et de la Formation au Sud; convention of February 18, 2013). D. F. Traoré and D. D. Zoh were supported by PhD fellowships provided by the IRD (Département des Programmes de Formation au Sud) and A. B. Sagna by a fellowship provided by the "Méditerranée Infection" Foundation (Marseille, France).

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**References**

1. Pages F, Orlandipradines E, Corbel V. Vecteurs du paludisme: biologie, diversité, contrôle et protection individuelle. Med Mal Infect. 2007;37:153–61.

2. Dossou-Yovo J, Doannio JMC, Diarrassouba S, Chauvancy G. Impact d'aménagements de rizières sur la transmission du paludisme dans la ville de Bouaké, Côte d'Ivoire. Bull Soc Path Exot. 1998;91:327–33.

3. Mouchet J, Carnevale P, Coosemans M, Julvez J, Manguin S, Richard-Lenoble D, et al. Biodiversité du paludisme dans le monde [Internet]. Paris: J. Libbey; 2004. http://www/documentation.ird.fr/hor/fdi:010035112. Accessed 24 Nov 2018.

4. Robert V, Macintyre K, Keating J, Trape J-F, Duchemin J-B, Warren M, et al. Malaria transmission in urban sub-Saharan Africa. Am J Trop Med Hyg. 2003;68:169–76.

5. Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malar J. 2013;12:133.

6. Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW. Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005;3:81–90.

7. Donnelly MJ, McCall P, Lengeler C, Bates I, D'Alessandro U, Barnish G, et al. Malaria and urbanization in sub-Saharan Africa. Malar J. 2005;4:12.

8. Lengeler C. Insecticide-treated bed nets and curtains for preventing malaria. Cochrane Database Syst Rev. 2004;2:CD000363.

9. Maxwell CA, Msuya E, Sudi M, Njunwa KJ, Carneiro IA, Curtis CF. Effect of community-wide use of insecticide-treated nets for 3–4 years on malarial morbidity in Tanzania. Trop Med Int Health. 2002;7:1003–8.

10. Mboera LEG. Sampling techniques for adult Afrotropical malaria vectors and their reliability in the estimation of entomological inoculation rate. Tanzan Health Res Bull. 2005;7:117–24.

11. Badu K, Siangla J, Larbi J, Lawson BW, Afrane Y, Ong'echa J, et al. Variation in exposure to *Anopheles gambiae* salivary gland peptide (gSG6-P1) across different malaria transmission settings in the western Kenya highlands. Malar J. 2012;11:318.

12. Remoue F, Cornelie S, Ngom A, Boulanger D, Simondon F. Immune responses to arthropod bites during vector-borne diseases. In: Garraud O, editor. Update tropical Immunology. Kerala: Transworld research Network; 2005. p. 377–400.

13. Schwartz BS, Ribeiro JMC, Goldstein MD. Anti-tick antibodies: an epidemiologic tool in Lyme disease research. Am J Epidemiol. 1990;132:58–66.

14. Poinsignon A, Cornelie S, Mestres-Simon M, Lanfrancotti A, Rossignol M, Boulanger D, et al. Novel peptide marker corresponding to salivary protein gSG6 potentially identifies exposure to *Anopheles* bites. PLoS ONE. 2008;3:e2472.

15. Drame PM, Machault V, Diallo A, Cornélie S, Poinsignon A, Lalou R, et al. IgG responses to the gSG6-P1 salivary peptide for evaluating human exposure to *Anopheles* bites in urban areas of Dakar region, Sénégal. Malar J. 2012;11:72.

16. Sagna AB, Sarr JB, Gaayeb L, Drame PM, Ndiath MO, Senghor S, et al. gSG6-P1 salivary biomarker discriminates micro-geographical heterogeneity of human exposure to *Anopheles* bites in low and seasonal malaria areas. Parasit Vectors. 2013;6:68.

17. Ya-umphan P, Cerqueira D, Parker DM, Cottrell G, Poinsignon A, Remoue F, et al. *Anopheles* salivary biomarker to assess malaria transmission risk along the Thailand-Myanmar Border. J Infect Dis. 2016. https://doi.org/10.1093/infdis/jiw543.

18. Drame PM, Poinsignon A, Marie A, Noukpo H, Doucoure S, Cornelie S, et al. New salivary biomarkers of human exposure to malaria vector bites. In: Manguin S (ed.) *Anopheles* mosquitoes—new insights into malaria vectors. InTech; 2013. http://www.intechopen.com/books/anopheles-mosquitoes-new-insights-into-malaria-vectors/new-salivary-biomarkers-of-human-exposure-to-malaria-vector-bites. Accessed 24 Nov 2018.

19. Drame PM, Poinsignon A, Besnard P, Cornelie S, Le Mire J, Toto J-C, et al. Human antibody responses to the *Anopheles* salivary gSG6-P1 peptide: a novel tool for evaluating the efficacy of ITNs in malaria vector control. PLoS ONE. 2010;5:e15596.
20. Drame PM, Diallo A, Poinsignon A, Boussari O, Dos Santos S, Machault V, et al. Evaluation of the effectiveness of malaria vector control measures in urban settings of Dakar by a specific *Anopheles* salivary biomarker. PLoS ONE. 2013;8:e66354.

21. Traoré DF, Sagna AB, Adja AM, Zoh DD, Lingué KN, Coulibaly I, et al. Evaluation of malaria urban risk using an immuno-epidemiological biomarker of human exposure to *Anopheles* bites. Am J Trop Med Hyg. 2018;98:1353-9.

22. Hien S, Angbo K, Dasse S, Yeboah O, N'guessan K, Kouacou A, et al. L'haptoglobine chez les enfants atteints de paludisme grave a *Plasmo-dium falciparum* de 0 à 15 ans : relations avec l'âge, la parasitémie et le taux d'hémoglobine. J Appl Biosci. 2014;76:6425.

23. RGPH. Recensement Général de la Population et de l'Habitat. 2014 [Internet]. 2014. http://www.ins.ci/n/documents/RGPH2014_expo_dg.pdf. Accessed 03 July 2018

24. Dossou-Yovo J, Doannio JMC, Diarrassouba S. Préférences trophiques des vecteurs du paludisme dans la ville de Bouaké et dans les villages environnants de Côte d'Ivoire. Bull Soc Path Exot. 1998;91:257-8.

25. Khaemba BM, Mutani A, Bett MK. Studies of anopheline mosquitoes transmitting malaria in a newly developed highland urban area: a case study of Moi University and its environs. East Afr Med J. 1994;71:159-64.

26. Cooke MK, Kahindi SC, Oriango RM, Owaga C, Ayoma E, Mabuka D, et al. 'A bite before bed': exposure to malaria vectors outside the times of net use in the highlands of western Kenya. Malar J. 2015;14:259.

27. Vanek MJ, Shoo B, Mtasiwa D, Kiama M, Lindsay SW, Fillinger U, et al. Community-based surveillance of malaria vector larval habitats: a baseline study in urban Dar es Salaam, Tanzania. BMC Public Health. 2006;6:154.

28. Rossi P, Belli A, Sabatinelli G. Etude sur la dispersion d'Anopheles gambiae s.l. dans une zone urbaine à Ouagadougou, Burkina Faso. Parassitologia. 1986;28:33-9.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

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

BMC