Conditions of malaria transmission in Dakar from 2007 to 2010

Background Previous studies in Dakar have highlighted the spatial and temporal heterogeneity of Anopheles gambiae s.l. biting rates. In order to improve the knowledge of the determinants of malaria transmission in this city, the present study reports the results of an extensive entomological survey that was conducted in 45 areas in Dakar from 2007 to 2010. Methods Water collections were monitored for the presence of anopheline larvae. Adult mosquitoes were sampled by human landing collection. Plasmodium falciparum circumsporozoïte (CSP) protein indexes were measured by ELISA (enzyme-linked immunosorbent assay), and the entomological inoculation rates were calculated. Results The presence of anopheline larvae were recorded in 1,015 out of 2,683 observations made from 325 water collections. A water pH of equal to or above 8.0, a water temperature that was equal to or above 30°C, the absence of larvivorous fishes, the wet season, the presence of surface vegetation, the persistence of water and location in a slightly urbanised area were significantly associated with the presence of anopheline larvae and/or with a higher density of anopheline larvae. Most of the larval habitats were observed in public areas, i.e., freely accessible. A total of 496,310 adult mosquitoes were caught during 3096 person-nights, and 44967 of these specimens were identified as An.gambiae s.l. The mean An. gambiae s.l. human-biting rate ranged from 0.1 to 248.9 bites per person per night during the rainy season. Anopheles arabiensis (93.14%), Anopheles melas (6.83%) and An. gambiae s.s. M form (0.03%) were the three members of the An. gambiae complex. Fifty-two An. arabiensis and two An. melas specimens were CSP-positive, and the annual CSP index was 0.64% in 2007, 0.09% in 2008-2009 and 0.12% in 2009-2010. In the studied areas, the average EIR ranged from 0 to 17.6 infected bites per person during the entire transmission season. Conclusion The spatial and temporal heterogeneity of An. gambiae s.l. larval density, adult human-biting rate (HBR) and malaria transmission in Dakar has been confirmed, and the environmental factors associated with this heterogeneity have been identified. These results pave the way for the creation of malaria risk maps and for a focused anti-vectorial control strategy.


Background
Urban malaria is considered to be an emerging problem in Africa because the populations of most large African cities have grown exponentially over the last 30 years [1]. Furthermore, it has been estimated that by 2030, 54% of the African population are expected to live in cities [2,3]. Many studies have reported evidence of malaria transmission in urban areas, and although levels are usually lower than in peri-urban and rural areas [2,4], urban malaria is considered to be an emerging health problem of major importance in Africa [1]. In African cities, transmission is spatially heterogeneous and occurs in areas where conditions are favourable for malaria vectors [5][6][7][8]. Indeed, in urban settings, malaria risk heterogeneity is recorded over small distances due to the degrees of diversity of types of urbanization, density of the human population, quality of water and waste management, vector control measures, household factors, access to health care and patterns of human migration, which could import parasites from rural areas [4,9,10]. The occurrence of malaria in African towns has also been linked to agricultural practices [11][12][13][14], to the distance from breeding sites [15][16][17][18][19] and to the vegetation cover [19].
To optimise the use of human and financial resources of the programmes for malaria control and because of the heterogeneity of malaria transmission in towns, it is appropriate to guide vector control interventions in only the areas of transmission. When drawn at appropriate scales, entomological maps of the spatial and temporal distribution of malaria vector larvae and adults could provide valuable information for targeted malaria control and the selective allocation of resources. Drawing such maps could theoretically be based on exhaustive entomological studies that cover the entire city over several years, but identification of the environmental factors that are associated with malaria transmission and the remotely sensed measurement of those factors could lead the way to an operational malaria risk mapping. Indeed, over the past several decades, remote sensing (RS) and geographic information systems (GIS) have become tools for evaluating the environmental, meteorological and climatic factors that influence the geographical and temporal risk distribution of malaria [20][21][22][23][24][25]. From 2005 to 2007 in Dakar, entomological data have been collected from several areas [6,8], and a portion of these data have allowed for a preliminary step towards malaria risk mapping, which is based on remotely sensed data, by setting up two risk maps for the years 1996 and 2007 [26]. For more accurate and precise malaria risk mapping, the objective of the present study was to investigate the environmental conditions associated with malaria transmission in-depth in 45 studied areas of Dakar and its suburbs from the years 2007 to 2010.

Study site
Dakar (14°40'20" North, 17°25'22" West), the capital city of Senegal, is part of the Cap-Vert peninsula, which is located at the western-most point of Africa. The estimated population of Dakar was 1,030,594 inhabitants in 2005, which amounted to approximately 20% of the country's population, and the population density was 12,233 inhabitants per km 2 . The altitude of this area peaks at 104 m above sea level (Mamelles). The present study was conducted in 45 different areas of downtown Dakar and in Pikine, Thiaroye and Guediawaye, three of Dakar's satellite cities. The studied zones were chosen to cover as many diverse environments as possible in terms of type of urbanization, road network, vegetation and socio-economic level. Each zone was delimited on the ground to cover an area of approximately 200 × 200 m that depended on the technical and logistical limitations that were presented by the landscape. The studied zones are presented and listed in Figure 1, and each point is the centre of the area.

Climate and study period
The Cap-Vert peninsula has a mild Sahelian climate. The hot and wet season occurs from June to November, the average temperature ranges from 24 to 30°C and the average rainfall amount is about 400 mm. The cool and dry season occurs from December to May, and the average temperature ranges from 19°C to 25°C. In 2007, 2008 and 2009, the annual rainfall was 178 mm, 510 mm and 565 mm, respectively. Field studies were undertaken in September and October of 2007 and from July of 2008 to June of 2010; therefore, three rainy seasons and two dry seasons were covered. In summary, the number of zones studied per year were as follows: 10 zones in September-October 2007 (see [8]), 30 zones in July 2008-June 2009 and 30 zones in July 2009-June 2010. Each zone was followed for one or two years, and the periods of study are presented in Figure 1. For those periods, daily precipitation values that were recorded at the international airport were available (data from the national weather agency). Two of the 45 zones (Bima and BA 160) were also studied in 2005 and 2006, but these studies followed a slightly different sampling methodology (see [6]).

Mosquito larval sampling Characterization of water collections
Each of the studied areas were visited every 10 days during the three rainy seasons and every month during the two dry seasons, in September-October 2007 and from June 2008 to April 2010 without any interruption. The physical, biological and chemical characteristics of all of the open water collections found in the studied areas were recorded. For visual classifications, consistency was maintained by having the same team perform the field work for the full duration of the project (two persons out of four rotated to perform this work). Subsequently, all of the collected information was transferred to a Microsoft Excel spreadsheet.
The type of habitat was categorised as ditches or puddles, swamp areas, marshes, ponds or lakes, "céanes" (non cemented wells used in market gardens), cemented wells or basins, waterproof containers, and canals. A water collection was considered to be temporary when it was found to be dry at least once during the followup; otherwise, it was recorded as permanent. The presence of larvivorous fishes, such as Guppies, Gambusia or Tilapia, was assessed visually. In 2007, the temperature of the water was measured using a mercury-in-glass thermometer that was immerged into the water for 60 seconds. From 2008 to 2010, a waterproof pH meter (PHSCAN 30, Hanna) provided the temperature, pH and conductivity of the water. Turbidity was estimated by using a graduated, transparent bottle that had black letters written on the bottom. The bottle was filled with water from the collection, and turbidity was evaluated by measuring the level that the water reached before the letters were no longer visible. Graduations ranged from 0 cm, at the top of the bottle, to 26 cm, at the bottom of the bottle; therefore, a higher value indicated greater turbidity. The proportion of the water surface that was covered by vegetation, e.g., water lettuce, water lentils, grass or algae, was estimated visually, and the proportion of the shadow on the water at midday, which was dependent on the surroundings, was estimated visually. The salinity of the water was proportional to the conductivity.

Field larval sampling
All of the water collections in the 45 studied areas were examined for larvae. Larvae and pupae were sampled using a standard dipping method [27], which entailed performing a minimum of 20 dips for small collections and up to 100 dips for large water bodies. When anopheline specimens were found, the larval density was calculated as the number of larvae (all instars) and pupae, which were further emerged and identified at the laboratory, per dip and recorded for each water collection. The presence of Culicinae larvae was also recorded.

Larval mosquito molecular analysis
Random samples of the Anopheles specimen, which were found in the studied water bodies, were stored by date and breeding site in numbered vials that were filled with 70% alcohol. In the studied areas where adult An. melas had been caught, the species of the larvae specimens were identified by PCR (polymerase chain reaction). The PCR analyses were done on pools of specimens, and each pool corresponded to a sample that was collected at one breeding site on one date between September and October of 2007 and between July and September of 2008.

Larval sampling inside private properties during the 2009 rainy season
A transversal field study was undertaken in September of 2009 to investigate larval collections inside private properties. In each of the 45 studied zones, a random geographical sampling was performed using a GIS to choose the places that were to be visited. The sampling scheme included 30 private properties in each zone, which were numbered from 1 to 30, and the places were visited in the order of their numbers. When a private property was closed or if access was denied, the following property was visited, until a maximum of approximately 10 private properties were accessed. At every acceded property, the gardens, courtyards, balconies and flat roofs were searched for the presence of water bodies. A description of these water collections was recorded, and the presence of larvae was researched visually for small collections or by using a standard dipping method. These data were analysed separately from those that were collected outside private properties.

Adult mosquitoes
Field adult mosquito sampling Adult mosquito sampling was carried out by human landing catch, once every two weeks during September-October of 2007, during both extended 2008 and 2009 wet seasons (July to December), and once every month during both 2009 and 2010 dry seasons (January to June). For every studied zone, one catching point was located indoors and two were located outdoors except for Point E and Nord Foire. In the latter two areas, all three catching points were located outdoors, respectively from October of 2009 (and on the following 10 nights) and from September of 2009 (and on the following 10 nights), because it was not possible to find indoors catching points. Indoor captures were conducted by leaving a window or door slightly ajar. The three catching points were located around the centre of each studied area. Within each area, the distance between each of the three catching points was approximately 30 meters. Collectors gave prior informed consent and received yellow fever immunizations and anti-malarial chemoprophylaxis for the duration of the study and one month afterward. Two collectors were contracted to collect mosquitoes at each catching point from 7:00 p.m. to 7:00 a.m. or from 8:00 p.m. to 7:00 a.m., depending on the hours of dawn and dusk in Dakar, and each collector rested every two hours. Collectors were rotated among the catching points on different collection nights to minimise sampling bias. The mosquitoes were recorded by catching point, date and hour of capture, and they were then sorted by genera. The anopheline mosquitoes were identified morphologically following the keys of Gillies and Coetzee [28], and Culicinae were identified morphologically following the keys of Edwards [29]. All anopheline mosquitoes were individually stored in numbered vials with desiccant and preserved at -20°C until processing.

Molecular identification of Anopheles gambiae s.l
Depending on the number of Anopheles caught per studied area, all specimens, or a random sample of approximately 100 specimens, belonging to the Anopheles gambiae complex were selected from each study site for species identification by PCR [30]. All Plasmodium falciparum circumsporozoite protein (CSP)-positive anopheline mosquitoes were also tested.

Biting patterns and sporozoïtes rates
The human biting rate (HBR) was expressed as the average number of female An. gambiae s.l. bites per person per night for both outdoor and indoor catching points. The entomological inoculation rates (EIR) were calculated to be the products of the HBR and CSP indexes and provided the number of infected Anopheles bites per person per night. The EIR for the rainy seasons (henceforth referred to as Rainy EIR) was derived using the Rainy HBR and the CSP indexes. The Rainy EIR was calculated following two methods: by using the average CSP indexes for all areas together or by using the CSP indexes that were actually evaluated in each area and that took into account the number of infected mosquitoes that were caught in each area. The Rainy EIR was multiplied by 60 to provide the total number of infected Anopheles bites that could be received during the peak of transmission (September and October). For the studied area that was followed during more than one year, the CSP indexes were compared.

Statistical analysis
The statistical analyses for the larval study aimed to identify the following: the determinants of the presence or absence of Anopheles larvae in the water collections, andthe factors associated with the density of Anopheles larvae at the breeding sites. Continuous independent variables were dichotomised at the mean or median, and the statistical unit was defined as the measurement of variables for a given water collection at one date. In longitudinal studies, some correlation may exist between observations made for the same water collection; therefore, to take this possible interdependence of observations into account, random effect models were used. The presence or absence of larvae was analysed using a logistic regression model, and the larval density was analysed using a negative binomial regression model with the number of larvae per 100 dips as the dependant variable. The variables associated with the presence or the density of the larvae with a p-value < 0.25 in univariate analysis were retained for multivariate analysis, and a backward, stepwise selection procedure was applied in the final model to keep variables with a pvalue < 0.05 in the analysis. Proportions were compared using a chi-squared test or Fisher exact test, as appropriate, and all analyses were performed using STATA 9.0 software (Stata-Corp LP). Some of the 2007 results have been presented in [8], but all data from 2007 through 2010 were taken into account in the analysis to increase statistical power.

Correlation between larvae and adults densities
The adult density fraction was estimated for each of the studied areas by calculating the total number of Anopheles that were caught during one decade, and these numbers were divided by the total number of Anopheles that were caught during the entire studied year. The larval density fraction at each studied area was estimated as follows: the product of the larval density was multiplied by the water surface for each breeding site for a decade, the sum of these products for all of the water collections was calculated for the entire studied area, the resulting number was then divided by the yearly value, and correlations between the larval density fraction and adult density fraction were graphically examined.

Mosquito larval sampling Characterization of water collections
A total of 325 water collections were examined. Depending on the persistence of the water bodies and on their characteristics of draining, emptying and drying out, each water body was investigated between one to 37 times. The types of open water collections included the following: 203 ditches or puddles; 24 swamp areas, marshes, ponds or lakes; 70 "céanes", cemented wells or basins; 11 large artificial holes (type underpinning); nine small, waterproof containers (tyres or buckets); and eight canals. Among them, 261 (80%) were temporary and 64 (20%) were permanent. The description of the quantitative physical, biological and chemical parameters that were recorded for the water collections are presented in Additional file 1.

Field larval sampling
A total of 228 of the 325 (70%) water bodies were found to be breeding sites for Anopheles at least once during the follow-up period, and 57 were breeding sites at every date when water was present. Of the temporary collections, 191 (73%) were habitats for anopheline larvae at least once during the follow-up, while for the permanent collections, 37 (58%) were Anopheles breeding sites at least once. Several of the studied areas did not harbour any breeding site during the full duration of their follow-up (Grand Yoff, HLM, Nord Foire and Point E). In breeding habitats, the density of Anopheles larvae and pupae ranged from 0.01 per dip (in large permanent collections in Golf and Hann IRD), to 42 per dip (in a concrete basin of an abandoned market-garden in Yarakh in 2008).

Larval mosquito molecular analysis
In 2007, PCR amplification undertaken on Anopheles larvae highlighted the presence of An. arabiensis only (no An. melas) among An. gambiae s.l. larvae [8]. During the 2008 rainy season, a total of 164 pools of larvae were tested for species identification by PCR amplification, and 147 pools were found to be positive for An. Arabiensis only, 15 pools were positive for both An. arabiensis and An. melas, and two pools were positive for An. melas only. All of the pools that contained An. melas were collected from three areas (Cafetériat, Pikine or Zone A), and those samples were taken from a large marshy area (Niaye), artificial basins or puddles. Due to laboratory constraints, no pools from Bima, Cambérène, Golf, or Potou where tested. In the breeding sites containing An. melas, the conductivity of the water was found to be 5.03 (range: 0.3 to > 20; 95% CI: 1.43 -8.62), while this value was found to be 1.75 (range: 0.06 to 12.03; 95% CI: 1.40 -2.09) in the breeding sites that did not contain An. melas. This difference was significant (Student t-test p < 0.0001).

Longitudinal survey of breeding sites
A total of 2,903 observations were recorded in the 325 sampled water collections in October-September of 2007 and between July of 2008 and April of 2010. The presence or absence of Anopheles larvae was known for 2,683 of those observations (1015 observations, 38%, were positive for larvae), and the Anopheles larval density was known for 1008 of the observations. The results of the univariate and multivariate analyses of the presence of Anopheles larvae and of the Anopheles larval density are presented in Additional file 2 and Table 1. A water pH equal to or above 8.0, the absence of larvivorous fishes, the wet season, the presence of surface vegetation, the persistence of water collection and a location in a lowly urbanised area were significantly associated with the presence of anopheline larvae in open water collections. Among the water collections that contained anopheline larvae, a water temperature equal to or above 30°C, a water pH equal to or above 8.0, the absence of larvivorous fishes and the wet season were significantly associated with a higher density of anopheline larvae.

Larval sampling inside private properties during the 2009 rainy season
Between five to 12 private properties or buildings that were not usually accessible for monitoring water collections were visited during September of 2009 in each of the 45 studied areas, totalling 355 private properties. In 153 of the private properties, no water collection or water container was found, but in the remaining 202 private properties, water was found in small containers (e.g., buckets, bowls, saucepans, and cans), in pools or basins, or in puddles that were located in gardens or on flat roofs or balconies. Culicidae larvae were found in 80 (23%) of the private properties, which were located in 27 of the study areas, and Anopheles larvae were found in 11 (3%) of the 355 private properties, which were located in 8 study areas ( (Table 3). Cumulatively for the 45 studied areas, the peak biting time was between 1:00 a.m. and 5:00 a.m. (Figure 2), and 77.4% of the bites occurred after midnight.
Molecular identification of An. gambiae s.l Of the 3,775 specimens of An. gambiae s.l. that were successfully tested by PCR, 3,516 were identified as An. arabiensis (93.14%), 258 were identified as An. melas Biting patterns and sporozoïte rates  Figure 4 provides the mean Rainy HBR for every studied area, Figure 5 shows the spatial distribution of these Rainy HBR, and Figure 6 provides the mean Dry HBR for each zone. The highest Anopheles Dry HBR value reached more than 15 bites per person per night (Golf area). A total of 44,628 An. gambiae s.l. were analysed by ELISA for the detection of the P. falciparum antigen, and 54 were positive (52 An. arabiensis and two An. melas). Table 4 Figure 7 shows a very close parallel between larval and adult density fractions for the three studied years. In fact, the correlation coefficient between these fractions was 0.87.

Anopheles larval habitats in Dakar
Several factors were found to be associated with the occurrence and abundance of anopheline larvae that  were collected in Dakar from 2007 to 2010. Using multivariate analyses, the probability of the presence of anopheline larvae and the larval density were positively associated with numerous factors that could be classified into the following four categories: climatic factors (rainy season), environmental factors (lower level of urbanization), morphological factors (temporary collections) and biological factors (increased surface vegetation, higher water temperature, higher pH, and the absence of larvivorous fishes). These results were consistent with the known preference of An. gambiae s.l. for breeding in temporary pools [28,32]. Higher water temperature was not related to the presence of anopheline larvae but was linked to the density of larvae in a collection. The water temperature is perhaps not a key factor in the laying of eggs by parous females; however, it impacts the development of An. gambiae immature stages. Low-floating vegetation was previously found to be a determinant of the presence of anopheline larvae [12,33]; however, in the present study, the level of floating vegetation was positively associated with the presence of larvae. This difference could be explained by the difficulties of measuring the surface of vegetation that covers a water collection. Also, the surface vegetation can be a proxy for the presence of food, which can favour the presence of anopheline larvae. The pH was significantly associated with the occurrence or the abundance of larvae, and this factor could also be considered to be an indicator of the presence of food for larvae. The existence of a muddy bottom in and around the water collection was favourable for the presence of larvae, and this observation was consistent with the possibility for An. gambiae to lay their eggs on soil around puddle larval habitats [34,35]. However, a muddy bottom could also be a proxy for an increased persistence of water, which is different from asphalt or sand. As was previously described, the presence of predator fishes was associated with a lower probability of larvae or lower larval density.
In 1998 in Dakar, Awono-Ambéné et al [36] confirmed that predation in the "céanes" was probably mainly due to fishes (Gambusia and Tilapia) that were introduced in Dakar in the 1930s for larval control. Many works have been conducted on the impact of urban agriculture on malaria transmission [11][12][13][14][15]37,38], and most of these studies have shown that this practice was associated with a higher level of malaria transmission in the surroundings areas. Specifically, urban agricultural activities can provide breeding and resting sites for malaria

Studied areas and year
Mean number of An. gambiae s.l. per person per night during the rainy season vectors. The results of the present study tend to show that market-gardens provided resting sites for adult Anopheles, rather than increased the number of larval habitats, as was previously demonstrated in Ghana [37]. Indeed, the water collections that were located in market-gardens less frequently harboured anopheline larvae, and when anopheline larvae were present, the densities were lower than in surroundings areas (as determined by univariate analysis). This occurrence was probably due to the presence of larvivorous fishes in traditional wells and to the frequent perturbations of water due to watering; however, one cannot exclude that this effect occurred because of the use of pesticides by urban farmers [39,40].

Anopheles larvae in public and private areas
The present study has found that few larval habitats were located inside private sectors, in comparison to public areas (i.e., freely accessible areas). Regardless of the type of housing or socio-economic level, the population generally did not allow Anopheles breeding sites to appear and persist in their neighbourhoods. As a consequence, the efficacy of any larval control measures would not suffer by being focused on freely accessible public zones, and excluding private areas from this treatment could simplify the control procedures.

Spatial heterogeneity in malaria transmission in Dakar
The present results provide new evidence for malaria transmission in downtown Dakar and its nearby suburbs. The spatial heterogeneity of HBR was very marked, and it ranged from 0.1 to nearly 250 bites per person per night during the rainy season. According to the analysis undertaken in the studied areas, most of the adult Anopheles were caught between July and December, and a marked peak occurred in September-October.  been possible to compare data because the present study showed that the An. gambiae s.l. bites were concentrated in September and October, with a minimal HBR thereafter (constantly less than 5 bites per person per night from November to July). Consistency between the mosquito collections (and the HBR evaluation) was also ensured during the study period because the protocol for sampling remained unchanged. In addition, no major environmental changes have been recorded in the  studied areas because of the low evolution rate in previously urbanised areas, which allows for the comparison of adult mosquito figures between years. The ELISA method for the detection of CSP antigen in malaria vectors has been proven to detect more positive mosquitoes than dissection of salivary glands. However, the ELISA method is subject to less inter-operator variations and is less time-consuming than salivary gland dissection. For these reasons, ELISA is often considered to be the most reliable method for estimating the proportion of Plasmodium-infected malaria vectors and for analysing malaria transmission.

Anopheles species involved in malaria transmission in Dakar
The following five anopheline species were collected from the 45 studied areas of Dakar: three species from the An. gambiae s.l. complex, Anopheles pharoensis and Anopheles ziemanni. Anopheles ziemanni had been previously caught on human bait at low levels in Senegal [41], as was the case in the present study, probably because of the low affinity of this species to human blood. Anopheles pharoensis is considered to be a secondary vector for malaria and has been found in the present study at low density. The presence of An. gambiae s.s. M form was reported in 2008-2009 (and in 2005-2006 in [6]) at extremely low levels. Therefore, An. gambiae s.s., An. ziemanni and An. pharoensis could not be significantly involved in malaria transmission in Dakar.
From 2007 to 2010 in the 45 studied zones, 54 An. gambiae s.l. were infected by P. falciparum: 52 An. arabiensis and two An. melas. Both infected An. melas were caught in Golf, near the large marshy area, in 2010, and this was the first time that this species was found to be infected with P. falciparum in Dakar. From 2007 to 2010, An. melas has been reported in 18 zones, where it accounted for 1% to 40.9% of the An. gambiae complex, according to the area and the year. Because An. melas is not as good a vector as An. arabiensis [42], it could be involved in malaria transmission only in areas in Dakar where it is very abundant. None of the breeding sites were exclusive for An. melas during the entire duration of the survey, and An. arabiensis was always found at these sites on at least one date, probably because of variations in the salinity levels of the breeding habitats. The co-occurrence of both of these species indicates that An. arabiensis larvae may also have a tolerance for low salinity.

Biting behaviours of adult Anopheles
In 2007, 2008-2009 and 2009-2010, 21%, 30% and 26%, respectively, of An. gambiae s.l. were caught indoors (at one indoor versus two outdoor catching points), but it was not possible to define a specific behaviour according to the season or to the studied area. Because logistical constraints limited the number of indoor catching points to one per site and the characteristics of the rooms could be different between areas (e.g., public offices or private properties of different sizes and types), the experimental design may have not allowed for the results to be representative of the adult An. gambiae s.l. behaviour in each of the studied areas. It can only be concluded that malaria transmission can occur in Dakar both outdoors and indoors. No difference was recorded in the hours of activity of An. gambiae s.l. between the different studied sites, and most of An. gambiae s.l. bites (> 75%) were recorded to take place after midnight. This finding highlighted the interest of using insecticide-impregnated bed nets to control malaria transmission in Dakar. Nevertheless, some bites were received during the first hours of the night when people were generally not asleep and not protected by mosquito nets. The use of other mosquito-control devices, like repellents or mosquito coils, could be proposed to increase protection for people from mosquitoes during the first hours of the night in Dakar.

Temporal heterogeneity of malaria transmission in Dakar
The abundance of An. gambiae s.l. during the year was clearly dependent on rainfall because the peak of the An. gambiae s.l. biting rate followed the peak of rainfall, with a two-week lag, regardless of the year. Afterward, the densities quickly decreased after the rains ended, and only isolated specimens were caught until the beginning of the new rainy season in July. Nevertheless, in one area (Golf), some infected mosquitoes were caught during the dry season, indicating that a permanent malaria transmission is possible in this area of Dakar, which is located near the large marshy area. When comparing the human biting rates in the common areas between 2007, a year with low rainfall, and 2008, a year with heavy rainfall, the peaks of Anopheles bites were, on average, six-fold higher in 2008 than in 2007. This was consistent with the fact that rainfall is closely related to the presence of surface water, larval productivity and adult densities, the latter two being highly correlated. Indeed, several studies have highlighted the low dispersion of Anopheles in urban settings due to the high density of houses and the proximity of readily available hosts for blood meals; therefore, most An. gambiae s.l. that are caught on human bait probably hatched from eggs that have been laid in a breeding site in the same or neighbouring area [15][16][17][18][19].

Toward malaria risk mapping
Most of the factors that are associated with the presence and density of larvae in water collections are measurable or detectable, either directly or indirectly, by remote sensing (RS) and include the following: the soil temperature, vegetation and trees that are able to provide shade, the surface vegetation, the moisture level of the ground, the density of urbanization, and the persistence of water collections during the dry season. Consequently, RS and geographic information systems (GIS) can provide useful information for mapping Anopheles breeding sites, as has been shown in numerous studies [43]. Furthermore, the close relationship between larval density and adult HBR suggests that maps of Anopheles breeding sites could provide baseline data for mapping Anopheles adult densities.

Conclusion
Suitable conditions for malaria transmission in Dakar were present during the extended rainy season from July to December, and the An. gambiae s.l. HBR increased with the intensity of rainfall in most of the studied areas. Nevertheless, the probability of encountering malaria vectors was not the same between the different areas. Taking into account this highly focused risk of malaria transmission, a targeted approach of vector control must be chosen. Because they are limited, the financial, material and human resources that are available must be concentrated on areas with a higher risk of transmission. Because it would not be reasonable to conduct routine, entomological studies in the entire city of Dakar, efforts should be concentrated on drawing risk maps to help targeting those anti-vectorial activities. In conclusion, strong, spatial and temporal heterogeneity of An. gambiae s.l. larval density, HBR and malaria transmission in Dakar has been confirmed, and the environmental factors that are associated with this heterogeneity have been identified. Moreover, the results of the present study will pave the way to set up malaria risk maps using RS and GIS technologies.