Research | Open | Published:
A quantitative risk assessment approach for mosquito-borne diseases: malaria re-emergence in southern France
Malaria Journalvolume 7, Article number: 147 (2008)
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.
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.
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.
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.
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 . Malaria disappeared from the Camargue after World War II: the last Plasmodium vivax malaria epidemic occurred in 1943, with about 400 estimated cases . Recent entomological surveys reported huge Anopheles populations in this area [3–5], 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 , 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 [10–13], usually assessed at the regional scale in a semi-qualitative way (Figure 1, 2) . 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.
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
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 .
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 . 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.
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 .
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 . 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 .
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.
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).
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.
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 , 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.
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) . 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 . Considering expert knowledge and bibliographical data , 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.
In order to illustrate this, the distribution of the length of the trophogonic cycle was determined for August (Figure 6).
The daily survival rate (p)
(Figure 2) was estimated from parity rates (P) using the following formula: p = P1/u  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 . 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) . 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.
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.
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 . 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).
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).
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 , 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.
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.
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).
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.
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 . 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 , 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, 34–36] 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 .
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 . 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.
Rodhain F, Charmot G: Evaluation des risques de reprise de transmission du paludisme en France. Med Mal Infect. 1982, 12 (4): 231-236.
Sautet J: A propos d'une épidémie de paludisme en Camargue. Mars Med. 1944, 2 (15): 53-64.
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.
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.
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.
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.
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.
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.
WHO-CISID: WHO - CISID-malaria. [http://data.euro.who.int/cisid/]
Ejov M: Imported malaria and the risk of malaria reintroduction in Europe. CD News, Quaterly Communicable Disease Report, WHO Europe. 2004, 34: 3-7.
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.
WHO: Expert Committee on Malaria, Twelfth Report. Technical Report Series, N°324. 1966, 324:
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.
WHO: Receptivity to malaria and other parasitic diseases. Euro Reports and Studies, 15. 1979, Copenhagen , 108-
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.
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.
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.
Vose DJ: The application of quantitative risk assessment to microbial food safety. J Food Prot. 1998, 61 (5): 640-648.
Vose DJ: Risk analysis. A quantitative guide. 2000, Chichester, England , John Wiley and Sons, LTD
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.
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.
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.
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.
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.
Davidson G: Estimation of the survival rate of anopheline mosquitoes in nature. Nature. 1954, 174 (4434): 792-793. 10.1038/174792a0.
Moshkovsky S, Rashina M: Epidémiologie et parasitologie médicale à l'usage des entomologistes. 1951, Moscow
Macdonald G: The epidemiology and control of malaria. 1957, London , Oxford University Press, 201-
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.
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.
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.
Faulde MK, Hoffmann R, Fazilat KM, Hoerauft A: Malaria reemergence in Northern Afghanistan. Emerg Infect Dis. 2007, 13 (9): 1402-1404.
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.
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.
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.
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.
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.
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.
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.
The authors declare that they have no competing interests.
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.