Skip to main content

A quantitative risk assessment approach for mosquito-borne diseases: malaria re-emergence in southern France

Abstract

Background

The Camargue region is a former malaria endemic area, where potential Anopheles vectors are still abundant. Considering the importation of Plasmodium due to the high number of imported malaria cases in France, the aim of this article was to make some predictions regarding the risk of malaria re-emergence in the Camargue.

Methods

Receptivity (vectorial capacity) and infectivity (vector susceptibility) were inferred using an innovative probabilistic approach and considering both Plasmodium falciparum and Plasmodium vivax. Each parameter of receptivity (human biting rate, anthropophily, length of trophogonic cycle, survival rate, length of sporogonic cycle) and infectivity were estimated based on field survey, bibliographic data and expert knowledge and fitted with probability distributions taking into account the variability and the uncertainty of the estimation. Spatial and temporal variations of the parameters were determined using environmental factors derived from satellite imagery, meteorological data and entomological field data. The entomological risk (receptivity/infectivity) was calculated using 10,000 different randomly selected sets of values extracted from the probability distributions. The result was mapped in the Camargue area. Finally, vulnerability (number of malaria imported cases) was inferred using data collected in regional hospitals.

Results

The entomological risk presented large spatial, temporal and Plasmodium species-dependent variations. The sensitivity analysis showed that susceptibility, survival rate and human biting rate were the three most influential parameters for entomological risk. Assessment of vulnerability showed that among the imported cases in the region, only very few were imported in at-risk areas.

Conclusion

The current risk of malaria re-emergence seems negligible due to the very low number of imported Plasmodium. This model demonstrated its efficiency for mosquito-borne diseases risk assessment.

Background

In the past, malaria was endemic and constituted a major health issue in France in marshy areas, particularly the Camargue, which was an active focus until the beginning of the 20th century. Malaria decreased drastically due to the draining of marshes, rearing of livestock, improvement of housing and living conditions and the use of quinine [1]. Malaria disappeared from the Camargue after World War II: the last Plasmodium vivax malaria epidemic occurred in 1943, with about 400 estimated cases [2]. Recent entomological surveys reported huge Anopheles populations in this area [35], and considered Anopheles (Anopheles) hyrcanus as being the main potential malarial vector based on its anthropophilic feeding behaviour and abundance [4, 6]. Thus, the Camargue is currently facing an "anophelism without malaria" situation. Moreover, autochthonous transmission was recently suspected in the French Mediterranean coast in 2006 [7], supporting the idea that southern France remains suitable for malaria transmission. The number of imported malaria cases have increased dramatically since the 1970s, in parallel with increased international travels, with an average of about 6,400 cases per year for the last ten years in France, leading to a massive Plasmodium introduction from endemic countries into France [8, 9]. These observations suggest that the malaria situation needs to be re-examined, and the aim of this paper is to infer current risk of malaria re-emergence, to identify hot spots for malaria re-emergence in the Camargue and to develop a generic model for mosquito transmitted diseases.

The risk of malaria re-emergence in an area (i.e., the recurrence of malaria transmission in an area) may be estimated by three factors: receptivity, infectivity and vulnerability [1013], usually assessed at the regional scale in a semi-qualitative way (Figure 1, 2) [14]. In this article, a quantitative entomological risk, which is the product of receptivity and infectivity, is calculated and the impact of vulnerability is discussed. The main objective of this work was to estimate the risk of malaria re-emergence at the local scale, considering the temporal and spatial local variations of the three components in order to identify hot spots for malaria resurgence in the Camargue.

Figure 1
figure 1

Risk of malaria re-emergence.

Figure 2
figure 2

Entomological factors.

The main difficulty that occurs when modeling is the quantification of biological parameters, especially for entomological data, as the field and laboratory studies are very painstaking, time-consuming and only rarely permit conclusions on a precise value. In this article, a probabilistic approach, taking into account the uncertainties and variability of inputs, was applied to a vector-borne disease, which constitutes an innovative method. Even if the malaria situation needs to be re-examined in the Camargue, this disease does not constitute a major health issue. The aim of this article is not to provide a public health tool that can be used to control malaria in the Camargue, but to present an innovative approach to spatialized quantitative risk assessment applied to a vector-borne disease.

Receptivity and infectivity were estimated for the potential vector An. hyrcanus as i) it is now considered as the main potential malaria vector and ii) other Anopheles species are rare in the Camargue area. However, the approach developed here is applicable to other mosquito species.

Materials, methods and data processing

Study area

The Camargue is the main wetland area in Southern France and covers the Rhone river delta. This area has a Mediterranean climate characterized by warm, dry summers and mild, wet winters. Total annual rainfall usually ranges between 500 and 700 mm and occurs mainly in autumn, and the annual mean temperature is 14°C.

Water pools and marshes cover a large part of the Camargue. Water is provided either by rains or a very tight canal network diverted from the river Rhone used to irrigate paddies or to fill marshes. Management of water is at the level of individual field owners depending on use: grazing for horses, cows or sheep, exploitation of reeds or rice, hunting reserves for waterfowl and nature preservation. Landscapes in the Camargue are strongly affected by the duration of submersion and the salinity of the soils. They are organized roughly in a south-north gradient of salinity, with agricultural land and reed marshes in the north and natural salty ponds and salt marshes in the south [15].

Moreover, there are various forms of agriculture (including vineyard, paddies, market gardening, fruit growing and exploitation of reeds) and rice, which covers more than 18,000 hectares in the Camargue, is the main cultivation [16]. Livestock includes horses, cows and sheep.

The Camargue hosts nearly 100,000 permanent inhabitants distributed between towns, hamlets and isolated houses. Moreover, the number of people increases in summer due to tourism.

Quantitative risk assessment using a probabilistic approach

The objective of this method was to organize and analyse scientific information in order to infer the risk of malaria re-emergence taking into account the variability and uncertainty of the input components and the final risk estimate. Such analysis, using reiterated simulations, have been performed for a decade for risk assessment in food microbiology, for example [17, 18]. Information and data for the development of the entomological risk model were obtained from field surveys, literature, unpublished data and expert opinion. Biological parameters were estimated by probability distributions in a plausible way that is coherent and conceivable and they were fitted with Pert or beta distribution [18, 19] (Figure 3). An amount of 10,000 reiterated simulations generated by the Latin Hypercube method associated with the probability distributions was used to describe both variability and uncertainty within the input parameters and the model [20, 21]. The outcome is a statistical distribution of risk, as well as a mean value of the risk estimate. Sensitivity analysis was performed to point out factors responsible for the main impact on the risk estimate. The @risk® (Palisade Corporation) software version 4,5,3 was used.

Figure 3
figure 3

Probability distributions.

Meteorological data

Daily temperatures (2005) and mean monthly temperatures (from 2002 to 2006) recorded by MeteoFrance at Aigues-Mortes (western Camargue) and at Tour du Vallat (south-eastern Camargue) were used. Daily temperatures were smoothed with a moving average (running mean) of three days to filter daily variations. For both types of data (daily and monthly), we calculated the average temperature for the Camargue based on the two stations. Humidity was recorded in 2005 at Marais du Vigueirat. These meteorological data were used to estimate some of the biological parameters of An. hyrcanus.

Entomological data: receptivity

In order to assess receptivity, it was necessary to evaluate the human biting rate (ma), the vector biting rate (a), the survival rate (p) and the sporogonic cycle (n) and their spatial and temporal variations (Figure 2).

Space and time-dependency

Although we sampled a huge amount of An. hyrcanus (125,848 specimens captured, 504 females dissected), it has not been possible to estimate precisely potential spatial variations of some biological parameters. Thus and due to the small size of the Camargue, the biology of An. hyrcanus was considered to be homogeneous in the whole area, which means that the vector biting rate (a), survival rate (p) and sporogonic cycle (n) (Figure 2) did not vary spatially. On the contrary, as An. hyrcanus presence and density depend on the biotopes and the season [4, 5], the vector-host ratio m (and hence the human biting rate (ma)) (Figure 2) presents a strong spatial heterogeneity (Table 1). Spatial variations were assessed based on a Geographic Information System (GIS) computing data for each 30 meter-wide pixel in the Camargue [22].

Table 1 Space, time and Plasmodium species-dependent variations in receptivity, infectivity and vulnerability.

All parameters were considered time-dependent, except the Anopheles anthropophily (h) and, as a practical approximation, the survival rate (p) (Table 1). Time-variations, which are detailed below, were described at a monthly time step.

Human biting rate (ma)

The presence and density of An. hyrcanus were inferredusing remote sensing, entomological adults and larvae collections [22]. Analysis of larval data led to the definition of a larval index that was calculated for each pixel in the Camargue based on environmental key factors. An adult abundance index was generated from the larval index and was also calculated for each pixel in the study area. Comparison of the adult abundance index and the maximum number of An. hyrcanus captured in the same pixel with CDC-light traps+CO2 showed a highly significant linear regression, allowing us to infer, using key environmental factors, the maximum number of An. hyrcanus captured during the year with CDC-light traps+CO2 for each pixel in the Camargue [22].

The mean annual dynamics of An. hyrcanus in the Camargue was inferred from the results of several capture campaigns conducted during several years, in several places, using several capture methods (Table 2). Specimens of this species were collected from March to October, and presented huge abundance variation (none of them were collected during winter). In each capture month, the mean percentage of collected mosquitoes (among the total number of mosquitoes captured during the year) were 0, 0, 0, 4 [1; 7], 12 [2; 21], 66 [51; 80], 18 [9; 27], 0 from March to October, respectively (with the associated confidence interval in brackets), and the maximum number of An. hyrcanus captured during the year refers to the month of August. Hence, the spatial distribution of the annual maximum number of An. hyrcanus was combined with the mean annual dynamics in order to determine, for each pixel and for each month in the Camargue, the number of An. hyrcanus captured with CDC-light traps+CO2.

Table 2 Total number of An. hyrcanus collected per month and associated percentages (bracket)

From this, the human biting rate (ma) was inferred from a comparison between mean CDC-light trap captures and mean human landing captures. Captures conducted the same night, in the same area (Carbonnière vs Marais du Vigueirat) and in pixels of the same adult abundance index class were compared (the following adult abundance thresholds were chosen in order to determine the classes: 10; 100; 500; 1000; 2000; 5000 (Table 3 and 4). For example, human landing captures carried out on the 10th of August in Carbonnière in pixel Hu5 and Hu6 were compared to light trap captures carried out on the same date and in the same area, in pixel LT1, LT4 and LT5 (all of these pixels have an adult abundance index between 100 and 500). Comparison showed a highly significant linear relationship between CDC-light traps and human captures (r2 = 0.66, p < 0.01). Indeed, human biting rate (ma) was extrapolated from the following formula: ma = (0.39 LT + 0.049)*1.25. As human landing captures refer to the hour following sunset, a correction factor of 1.25 was associated with the estimation of ma in order to obtain ma for the entire night (Table 3).

Table 3 Number of An. hyrcanus captured during human landing catches sessions.
Table 4 Number of An. hyrcanus captured during light trap sessions. Adult abundance index relative to the pixel where the capture was conducted is indicated.

The result is the spatio-temporal distribution of ma in the Camargue, i.e., an estimation of the human biting rate for each 30 m × 30 m pixel and each month. The human biting rate was figured for the month of August in order to illustrate this paragraph (Figure 4). Variability and uncertainty were taken into account for each step leading to the assessment of variability and uncertainty of ma for each pixel.

Figure 4
figure 4

Spatial distribution of the human biting rate in the Camargue in August. Classes were arbitrary chosen with a logarithmic scale.

Anopheles hyrcanus anthropophily (h)

(Figure 2) was estimated from the comparison between human landing, light traps and horse bait trap results for An. hyrcanus and other Anopheles species [4], and was fitted as follows: Pert distribution (0.4; 0.5; 0.8)(Figure 5). It was assumed that h did not vary throughout the year.

Figure 5
figure 5

Distribution of An. hyrcanus anthropophily.

The length of the trophogonic cycle (u)

(Figure 2) was calculated using the following formula, which estimates the length of time of blood digestion: u = f1/(T-g1) where f1 and g1 are factors depending on humidity and T is the temperature. f1 and g1 were experimentally determined and evaluated at 36.5°C-days and 9.9°C, respectively, when the humidity reached 70–80% (from June to September 2005, mean monthly humidity varied from 63 to 90% in the Camargue) [23]. Moreover, Detinova assumed that it was reliable to add 24 hours to take into account the time necessary to find a host and the time necessary to find a pool to lay eggs in, in order to fulfil a complete gonotrophic cycle [23]. Considering expert knowledge and bibliographical data [24], a two day range of values was fixed for each month in the Camargue, and the length of the trophogonic cycle was fitted as indicated in Table 5.

Table 5 Length of the trophogonic cycle

In order to illustrate this, the distribution of the length of the trophogonic cycle was determined for August (Figure 6).

Figure 6
figure 6

Distribution of the length, in days, of the trophogonic cycle in August.

The daily survival rate (p)

(Figure 2) was estimated from parity rates (P) using the following formula: p = P1/u [25] which is relevant in the case of stable populations. Considering that the An. hyrcanus population increases and decreases progressively according to the Anopheles biology, the approximation that populations were stable during the summer was used [5]. Parity rates observed in June, July and September 2005 in Carbonnière and Marais du Vigueirat were used (due to the low number of dissected mosquitoes, the parity rate calculated in June in Carbonnière was not included) [4]. They were fitted as indicated in Table 6, and the mean value was calculated from the five probability distributions in order to obtain the mean An. hyrcanus parity rate in the Camargue from June to September 2005.

Table 6 Parity rates observed in the Camargue in 2005

During the same period (June to September 2005), the daily temperature varied from 18.5°C to 26.4°C, with a mean value of 23.2°C. The length of the trophogonic cycle, calculated using the same method as before, ranged from 3.2 days to 5.2 days, with a mean value of 3.7, and was fitted as follows: Pert distribution (3; 3.5; 5.5). Using the formula p = P1/u, a mean p value of 0.79 was obtained (Figure 7). Due to difficulty in obtaining a precise value for parity rates for the entire year, it was considered that p did not vary from March to October in this model.

Figure 7
figure 7

Distribution of the An. hyrcanus daily survival rate.

The length of the sporogonic period (n)

(Figure 2) was calculated as follows: n = f2/(T-g2), where T is temperature and f2 and g2 are Plasmodium species-dependent factors [26]. f2 and g2 were experimentally evaluated at 111°C-days and 16°C, 105°C-days and 14.5°C for Plasmodium falciparum and P. vivax, respectively. However, Grassi and MacDonald estimated threshold temperatures under which the sporogonic development is not completed: 18–19°C for P. falciparum, 15–17°C for P. vivax [24, 27].

Entomological data: infectivity

Anopheles hyrcanus susceptibility was estimated based on data derived from experimental membrane-feeding experiments conducted with cultured P. falciparum in the laboratories of Radboud University Nijmegen Medical Centre (The Netherlands) following a routine protocol [28, 29]. Among 350 An. hyrcanus tested, none were found able to transmit tropical P. falciparum strains. Hence, An. hyrcanus susceptibility to P. falciparum was fitted with a beta distribution: Beta (1; 351) (Figure 8).

Figure 8
figure 8

Distribution of An. hyrcanus susceptibility to P. falciparum .

Anopheles hyrcanus susceptibility to P. vivax has not yet been tested, but this species has been considered a vector of P. vivax in Afghanistan for more than 30 years [30, 31]. Genetic comparison based on ITS2 sequences between French, Afghan, Turkish and Iranian specimens concluded that they were identical. Hence, it was considered that An. hyrcanus was susceptible to tropical P. vivax strains, and estimation of its susceptibility, based on expert knowledge, was fitted as follows: Pert distribution (0.05; 0.20; 0.70) (Figure 9).

Figure 9
figure 9

Distribution of An. hyrcanus susceptibility to P. vivax .

An. hyrcanus susceptibility to either P. falciparum or P. vivax was considered to be homogeneous through time and space (Table 1).

Anopheles hyrcanus susceptibility to Plasmodium ovale and Plasmodium malariae was not inferred, as very little is known about European Anopheles susceptibility to these two species.

Parasitological data: vulnerability

Vulnerability is related to gametocyte carriers as the gametocyte is the stage transmissible to mosquitoes. Nevertheless, very little is known about gametocyte carriers, and any information must be considered in the context of imported malaria cases, i.e. taking account of the fact that the illness is diagnosed and treated. Since anti malarial-treatment in France is generally conducted with drugs that do not prevent gametocyte emergence [8], an assumption that every patient may develop some gametocytes was made. Hence, we used imported malaria cases to approach the gametocyte presence in the Camargue.

Data relative to imported malaria cases were obtained from five public hospitals localized in important towns in/around the Camargue (Montpellier, Nîmes, Avignon, Arles and Marseille). These data were analysed with regard to the date, Plasmodium species, patients' residence for the 2004–2005 period and contamination place.

It was estimated that imported malaria case data obtained from public hospitals represented about 50–55% of all imported malaria cases, with other cases being diagnosed by private laboratories [32, 33]. Thus, the total number of imported malaria cases was estimated using a 52.5% correcting factor assuming that epidemiological data provided by public hospitals were representative of the total number of imported malaria cases.

Results

Vulnerability

In 2004 and 2005, 657 imported cases were diagnosed in the region (corresponding to a total of 1251 estimated imported cases), among which P. falciparum, P. vivax, P. ovale and P. malariae represented 85.5%, 7%, 5.5% and 2%, respectively. Moreover, 35% of these cases occurred between August and September. Among 528 patients for whom the place of residence was known, 75% were living in Montpellier, Nîmes, or Marseille big cities, and only 7 patients were living in areas where An. hyrcanus could be present.

Countries in which patients were infected were not known for 2004–2005. For the 2001–2003 period, 96% of imported malaria cases were contracted in Africa, and the Comoros Islands represented the main place of contamination for cases diagnosed in Marseille, Arles and Nîmes. Data collected in Montpellier and Avignon showed that imported cases were mainly contracted in West Africa.

Entomological risk (receptivity*infectivity)

The mean value of the entomological risk was assessed for P. falciparum species from June to September, and for P. vivax for the month of August for each pixel of the Camargue map (Figure 10, 11, 12, 13, 14) (the entomological risk was not assessed for P. ovale and P. malariae due to the low number of imported cases and the lack of information concerning both species). It was calculated with 10,000 different randomly selected sets of values extracted from input distributions.

Figure 10
figure 10

Spatial variations of P. falciparum transmission risk estimate (ranging from 0 to more than 1) in June in the Camargue. Classes were arbitrary chosen with a logarithmic scale.

Figure 11
figure 11

Spatial variations of P. falciparum transmission risk estimate (ranging from 0 to more than 1) in July in the Camargue. Classes were arbitrary chosen with a logarithmic scale.

Figure 12
figure 12

Spatial variations of P. falciparum transmission risk estimate (ranging from 0 to more than 1) in August in the Camargue. Classes were arbitrary chosen with a logarithmic scale.

Figure 13
figure 13

Spatial variations of P. falciparum transmission risk estimate (ranging from 0 to more than 1) in September in the Camargue. Classes were arbitrary chosen with a logarithmic scale.

Figure 14
figure 14

Spatial variations of the P. vivax transmission risk estimate (ranging from 0 to more than 100) in August in the Camargue. Classes were arbitrary chosen with a logarithmic scale.

Strong differences were observed in the entomological risk for the two Plasmodium species: P. falciparum transmission risk estimate ranges from 0 to more than 1 although P. vivax transmission risk estimate ranges from 0 to more than 100.

Uncertainties of the risk estimate and sensitivity analysis

For each pixel the outcome of the model is a statistical distribution of the risk estimate, generated by variability and uncertainty within inputs. For example, the entomological risk estimate ranges from about 9.10-6 to 20, with 95% of the values being between 0.016 and 4.7 for pixel having a mean value of about 1.

The sensitivity analysis was conducted for pixels of which the mean entomological risk estimate is about 1. Results of the sensitivity analysis show the correlation coefficients between different varying inputs and the consecutive varying risk estimate (Figure 15).

Figure 15
figure 15

Sensitivity analysis.

Sensitivity analysis carried out for pixels presenting a lower risk estimate showed a predominance of susceptibility. However, correlation coefficients of the survival rate and the human biting rate were approximately equal. Sensitivity analysis carried out for the P. vivax entomological risk estimate showed the equal importance of the susceptibility, the survival rate and the human biting rate.

Discussion

The entomological risk index used in this article refers to the risk of transmission (Figure 1): the risk of transmission being high when the entomological risk is high. Results clearly indicate a maximum risk of P. falciparum transmission in August in the Camargue considering the spatial distribution and value of the entomological risk estimate, which results from length of the trophogonic cycle and sporogonic period and the human biting rate. The length of the trophogonic cycle and sporogonic development period are directly influenced by the mean temperature in our model, and the human biting rate depends directly on the An. hyrcanus dynamics and density [5]. The human biting rate being the only space-dependant factor in our model, one could be tempted to approach roughly the risk of transmission by the human biting rate. Nevertheless, the entomological risk presents also temporal variations, which depend not only on the human biting rate but also on the impact of temperature on the length of the trophogonic and sporogonic cycle. Dynamics of An. hyrcanus and mean temperature do not evolve exactly in the same way in the Camargue, which shows the necessity of estimating all the parameters of the entomological risk.

The risk of P. vivax transmission is more than one hundred times higher than the risk of P. falciparum transmission, which is due to infectivity and the length of sporogonic cycle differences.

Nevertheless, the entomological risk calculated herein is a theoretical index as the human biting rate reaches more than 10,000 bites per human per night reflecting the abundance of An. hyrcanus. This will, of course, never happen as no one can endure such a large number of bites. This suggests that, in the future, it could be necessary to combine the entomological risk with human presence and exposure to mosquito bites in order to evaluate the real human biting rate. Such analyses would require complementary geographical and sociological studies.

Although all of France faces a large number of imported cases [8], particularly in the south-east of the country, vulnerability in at risk areas is very low because most imported cases are present in large cities.

The entomological risk, referring to the risk of transmission in this article (Figure 1), has to be interpreted in a relative way due to its definition and the estimation of its parameters. Thus, this study underlined space, time and Plasmodium species-dependant of the risk of potential transmission. Of course, the risk of potential transmission is connected with the risk of malaria re-emergence when gametocytes carriers are introduced within at risk areas. Considering the entomological risk and the length in days of the infectious period of humans would allow estimation of R0 and the following absolute risk of malaria re-emergence. Finally, the current risk of malaria re-emergence seems negligible due to the very low number of imported Plasmodium.

As stated in the introduction, the aim of this study was not to build a public health tool for controlling malaria in the Camargue. The emphasis was on presenting an innovative approach of spatialized quantitative risk assessment applied to a vector-borne disease, which has not been previously conducted. The main advantages of such a probabilistic approach are the possibility of integrating the uncertainty and variability of inputs within a model and to quantify the uncertainty of the final risk estimate. The deterministic approaches used thus far have not taken uncertainty and variability into account [27, 3436] and have produced precise outcomes, which could lead to misinterpretation as the final risk estimate could vary significantly due to input variability. Integration of uncertainty and variability in deterministic models would rapidly lead to complicated models, requiring laborious mathematical developments [18].

Moreover, the approach applied in this study is based on distributions that are combined, resulting in quite sophisticated analyses that are intuitive and easily understood. The uncertainty within the risk estimate is a crucial point for decision makers, which usually apply some rough "safety margin" around a deterministic estimate to express their feeling of uncertainty. The pros of the approach developped in this article are rightly to quantify this uncertainty.

Such a method could be applied to other areas where malaria is still a threat or to emerging vector-borne diseases, such as the dengue or chikungunya virus infections. This method could be used in a controlled way, in order to identify areas and time periods that correspond to the highest risk of transmission, and to focus control measures where and when transmission is elevated.

It has been shown recently the impact of anthropogenic changes on potential malaria vectors in the Camargue over the last 60 years [37]. This method could be used also to predict the probable impact of future decisions concerning land use for example, and could be a useful tool for decision makers.

The sensitivity analysis underlined factors responsible for entomological risk uncertainty and variation, of which susceptibility, the survival rate and the human biting rate of An. hyrcanus have a major impact in this model. Variability in either of these parameters leads to variability and uncertainty in the risk estimate. In contrast to the susceptibility, survival rate and human biting rate, variability in the anthropophily range and the length of the trophogonic cycle range has only a minor impact on the risk estimate. The sensitivity analysis thus appears to be a tool particularly useful for the identification of key factors, which need to be assessed in field surveys.

In the context of emerging vector-borne diseases, emphasis is currently on developing and improving such quantitative risk assessment models integrating variability and uncertainty of biological parameters, which are usually difficult to assess, and especially parameters recorded in the field.

References

  1. Rodhain F, Charmot G: Evaluation des risques de reprise de transmission du paludisme en France. Med Mal Infect. 1982, 12 (4): 231-236.

    Article  Google Scholar 

  2. Sautet J: A propos d'une épidémie de paludisme en Camargue. Mars Med. 1944, 2 (15): 53-64.

    Google Scholar 

  3. Balenghien T, Fouque F, Sabatier P, Bicout DJ: Horse, bird, and human-seeking behavior and seasonal abundance of mosquitoes in a West Nile virus focus of southern France. J Med Entomol. 2006, 43 (5): 936-946. 10.1603/0022-2585(2006)43[936:HBAHBA]2.0.CO;2.

    Article  CAS  Google Scholar 

  4. Ponçon N, Toty C, L'Ambert G, Le Goff G, Brengues C, Schaffner F, Fontenille D: Biology and dynamics of potential malaria vectors in southern France. Malar J. 2007, 6 (1): 18-10.1186/1475-2875-6-18.

    Article  PubMed Central  Google Scholar 

  5. Ponçon N, Toty C, L'Ambert G, Le Goff G, Brengues C, Schaffner F, Fontenille D: Population dynamics of pest mosquitoes and potential malaria and West Nile virus vectors in relation to climatic factors and human activities in the Camargue, France. Med Vet Entomol. 2007, 21: 350-357.

    Article  PubMed  Google Scholar 

  6. Ponçon N, Toty C, Kengne P, Alten B, Fontenille D: Molecular evidence for similarity between Anopheles hyrcanus (Diptera: Culicidae) and Anopheles pseudopictus (Diptera: Culicidae), sympatric potential vectors of malaria in France. J Med Entomol. 2008, 45: 576-580. 10.1603/0022-2585(2008)45[576:MEFSBA]2.0.CO;2.

    Article  PubMed  Google Scholar 

  7. Doudier B, Bogreau H, DeVries A, Ponçon N, Stauffer WM, Fontenille D, Rogier C, Parola P: Autochthonous malaria from Marseille to Minneapolis ?. Emerg Infect Dis. 2007, 13 (8): 1236-1238.

    Article  PubMed Central  Google Scholar 

  8. Legros F, Arnaud A, El Minouni B, Danis M: Paludisme d'importation en France métropolitaine : données épidémiologiques 2001-2004. Bulletin épidémiologique hebdomadaire. 2006, 32: 235-236.

    Google Scholar 

  9. WHO-CISID: WHO - CISID-malaria. [http://data.euro.who.int/cisid/]

  10. Ejov M: Imported malaria and the risk of malaria reintroduction in Europe. CD News, Quaterly Communicable Disease Report, WHO Europe. 2004, 34: 3-7.

    Google Scholar 

  11. Garrett-Jones C, Shidrawi GR: Malaria vectorial capacity of a population of Anopheles gambiae: an exercise in epidemiological entomology. Bull World Health Organ. 1969, 40 (4): 531-545.

    PubMed Central  CAS  PubMed  Google Scholar 

  12. WHO: Expert Committee on Malaria, Twelfth Report. Technical Report Series, N°324. 1966, 324:

    Google Scholar 

  13. Garrett-Jones C: Prognosis for interruption of malaria transmission through assessment of the mosquito's vectorial capacity. Nature. 1964, 204: 1173-1175. 10.1038/2041173a0.

    Article  CAS  PubMed  Google Scholar 

  14. WHO: Receptivity to malaria and other parasitic diseases. Euro Reports and Studies, 15. 1979, Copenhagen , 108-

    Google Scholar 

  15. Dervieux A, Allard P, Naizot T, Auda Y: La Camargue est-elle un polder? Changement et stabilité dans le delta du Rhône depuis le début du XIXème siècle. Equilibre et ruptures dans les écosystèmes durant les 20 derniers millénaires en Europe de l'Ouest. Edited by: Richard H, Vignot A. 2002, Besançon , 1: 315-329.

    Google Scholar 

  16. Mathevet R, Mouret JC, Mesléard F: Du mondial au local : originalité et dépendances de la riziculture camarguaise. Faire savoirs. 2002, 2: 49-56.

    Google Scholar 

  17. Pouillot R, Miconnet N, Afchain AL, Delignette-Muller ML, Beaufort A, Rosso L, Denis JB, Cornu M: Quantitative risk assessment of Listeria monocytogenes in French cold-smoked salmon: I. quantitative exposure assessment. Risk analysis. 2007, 27 (3): 683-700. 10.1111/j.1539-6924.2007.00921.x.

    Article  PubMed  Google Scholar 

  18. Vose DJ: The application of quantitative risk assessment to microbial food safety. J Food Prot. 1998, 61 (5): 640-648.

    CAS  PubMed  Google Scholar 

  19. Vose DJ: Risk analysis. A quantitative guide. 2000, Chichester, England , John Wiley and Sons, LTD

    Google Scholar 

  20. Anderson EL, Hattis D: When and how can you specify a probability distribution when you don't know much. Foudations: A. Uncertainty and variability. Risk analysis. 1999, 19: 47-49.

    Google Scholar 

  21. Cassin MH, Lammerding AM, Todd EC, Ross W, McColl RS: Quantitative risk assessment for Escherichia coli O157:H7 in ground beef hamburgers. Int J Food Microbiol. 1998, 41 (1): 21-44. 10.1016/S0168-1605(98)00028-2.

    Article  CAS  PubMed  Google Scholar 

  22. Tran A, Ponçon N, Toty C, Linard C, Guis H, Ferré JB, D.L. S, Roger F, de la Rocque S, Fontenille D, Baldet T: Using remote sensing to map larval and adult populations of Anopheles hyrcanus (Diptera: Culicidae) potential malaria vector in Southern France. Int J Health Geogr. 2008, 7: 9-10.1186/1476-072X-7-9.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Detinova TS: Age-grouping methods in Diptera of medical importance with special reference to some vectors of malaria. Monogr Ser World Health Organ. 1962, 47: 13-191.

    CAS  PubMed  Google Scholar 

  24. Jetten TH, Takken W: Anophelism without malaria in Europe. A review of the ecology and distribution of the genus Anopheles in Europe. 1994, Wageningen, The Netherlands , Landbouwuniversiteit Wageningen (Wageningen Agricultural University), No. 94-5, 69.

    Google Scholar 

  25. Davidson G: Estimation of the survival rate of anopheline mosquitoes in nature. Nature. 1954, 174 (4434): 792-793. 10.1038/174792a0.

    Article  CAS  PubMed  Google Scholar 

  26. Moshkovsky S, Rashina M: Epidémiologie et parasitologie médicale à l'usage des entomologistes. 1951, Moscow

    Google Scholar 

  27. Macdonald G: The epidemiology and control of malaria. 1957, London , Oxford University Press, 201-

    Google Scholar 

  28. Ponnudurai T, Verhave JP, Meuwissen JH: Mosquito transmission of cultured Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1982, 76 (2): 278-279. 10.1016/0035-9203(82)90298-X.

    Article  CAS  Google Scholar 

  29. Ponnudurai T, Billingsley PF, Rudin W: Differential infectivity of Plasmodium for mosquitoes. Parasitol Today. 1988, 4 (11): 319-321. 10.1016/0169-4758(88)90114-7.

    Article  CAS  PubMed  Google Scholar 

  30. Onori E, Nushin MK, Cullen JE, Yakubi GH, Mohammed K, Christal FA: An epidemiological assessment of the residual effect of DDT on Anopheles hyrcanus sl and An. pulcherrimus (Theobald) in the North Eastern region of Afghanistan. Trans R Soc Trop Med Hyg. 1975, 69 (2): 236-242. 10.1016/0035-9203(75)90161-3.

    Article  CAS  PubMed  Google Scholar 

  31. Faulde MK, Hoffmann R, Fazilat KM, Hoerauft A: Malaria reemergence in Northern Afghanistan. Emerg Infect Dis. 2007, 13 (9): 1402-1404.

    Article  PubMed Central  PubMed  Google Scholar 

  32. Legros F, Fromage M, Ancelle T, Burg E, Janot C, Maisonneuve P, Danis M: Enquête nationale de recensement des cas de paludisme d'importation en France métropolitaine pour l'année 1997. Bulletin épidémiologique hebdomadaire. 1999, 11: 41-42.

    Google Scholar 

  33. Danis M, Legros F, Thellier M, Caumes E: Données actuelles sur le paludisme en France métropolitaine. Med Trop (Mars). 2002, 62 (3): 214-218.

    CAS  Google Scholar 

  34. Koella JC: On the use of mathematical models of malaria transmission. Acta Trop. 1991, 49 (1): 1-25. 10.1016/0001-706X(91)90026-G.

    Article  CAS  PubMed  Google Scholar 

  35. Martens P, Kovats RS, Nijhof S, de Vries P, Livermore MTJ, Bradley DJ, Cox J, McMichael AJ: Climate change and future populations at risk of malaria. Glob Environ Change. 1999, 9: S89-S107. 10.1016/S0959-3780(99)00020-5.

    Article  Google Scholar 

  36. Lindsay SW, Thomas CJ: Global warning and risk of vivax malaria in Great Britain. Global change and human health. 2001, 2 (1): 80-84. 10.1023/A:1011995115713.

    Article  Google Scholar 

  37. Ponçon N, Balenghien T, Toty C, Ferré JB, Thomas C, Dervieux A, L'Ambert G, Schaffner F, Bardin O, Fontenille D: Effects of local anthropogenic changes on potential malaria vector Anopheles hyrcanus and West Nile virus vector Culex modestus, Camargue, France. Emerg Infect Dis. 2007, 13: 1810-1815.

    Article  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Geert-Jan van Gemert and Marga van de Vegte-Bolmer for technical expertise in conduct of the susceptibility trials. We thank Grégory L'Ambert, Jean Baptiste Ferré, Charles Jeannin, Claire Duchet and Alexandre Carron from EID-Méditerranée for help in human landing captures and data on An. hyrcanus dynamics. We thank Philippe Parola and Marie Pierre Ranjeva from the North hospital of Marseille, Didier Basset from the hospital of Montpellier, Laurence Lachaud from the hospital of Nîmes, Muriel Roumier from the hospital of Arles, Martine Derré from the hospital of Avignon, Anzime Ali Mohamed from the "Hôpital d'Instruction des Armées Laveran" and Frabrice Legros from the CNREPIA for providing data on imported malaria cases.

This work was partially funded by the French Ministry of Agriculture and EU grant GOCE-2003-010284 EDEN. The paper is catalogued by the EDEN Steering Committee as EDEN0106 http://www.eden-fp6project.net/. The contents of this publication are the sole responsibility of the authors and can in no way be taken to reflect the views of the European Union.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Didier Fontenille.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

NP designed the study, analysed the data and drafted the manuscript. AT analysed the data and carried out spatial analysis. CT participated in data collection and analysis. AL managed the susceptibility trials. DF conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved of the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ponçon, N., Tran, A., Toty, C. et al. A quantitative risk assessment approach for mosquito-borne diseases: malaria re-emergence in southern France. Malar J 7, 147 (2008). https://doi.org/10.1186/1475-2875-7-147

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2875-7-147

Keywords