<title>Background</title> <p>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 20<sup>th </sup>century. Malaria decreased drastically due to the draining of marshes, rearing of livestock, improvement of housing and living conditions and the use of quinine [##UREF##0##1##]. Malaria disappeared from the Camargue after World War II: the last <italic>Plasmodium vivax </italic>malaria epidemic occurred in 1943, with about 400 estimated cases [##UREF##1##2##]. Recent entomological surveys reported huge <italic>Anopheles </italic>populations in this area [##REF##17017231##3##, ####REF##17313664##4##, ##REF##18092973##5####18092973##5##], and considered <italic>Anopheles (Anopheles) hyrcanus </italic>as being the main potential malarial vector based on its anthropophilic feeding behaviour and abundance [##REF##17313664##4##,##REF##18533455##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 [##REF##17953101##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 <italic>Plasmodium </italic>introduction from endemic countries into France [##UREF##2##8##,##UREF##3##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.</p> <p>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 [##UREF##4##10##, ####REF##5306719##11##, ##UREF##5##12##, ##REF##14268587##13####14268587##13##], usually assessed at the regional scale in a semi-qualitative way (Figure ##FIG##0##1##, ##FIG##1##2##) [##UREF##6##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.</p> <p>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.</p> <p>Receptivity and infectivity were estimated for the potential vector <italic>An. hyrcanus </italic>as i) it is now considered as the main potential malaria vector and ii) other <italic>Anopheles </italic>species are rare in the Camargue area. However, the approach developed here is applicable to other mosquito species.</p>
<title>Materials, methods and data processing</title> <title>Study area</title> <p>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.</p> <p>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 [##UREF##7##15##].</p> <p>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 [##UREF##8##16##]. Livestock includes horses, cows and sheep.</p> <p>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.</p> <title>Quantitative risk assessment using a probabilistic approach</title> <p>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 [##REF##17640216##17##,##REF##9709243##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 [##REF##9709243##18##,##UREF##9##19##] (Figure ##FIG##2##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 [##UREF##10##20##,##REF##9631335##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<sup>® </sup>(Palisade Corporation) software version 4,5,3 was used.</p> <title>Meteorological data</title> <p>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 <italic>An. hyrcanus</italic>.</p> <title>Entomological data: receptivity</title> <p>In order to assess receptivity, it was necessary to evaluate the human biting rate (ma), the vector biting rate (<italic>a</italic>), the survival rate (<italic>p</italic>) and the sporogonic cycle (<italic>n</italic>) and their spatial and temporal variations (Figure ##FIG##1##2##).</p> <title>Space and time-dependency</title> <p>Although we sampled a huge amount of <italic>An. hyrcanus </italic>(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 <italic>An. hyrcanus </italic>was considered to be homogeneous in the whole area, which means that the vector biting rate (<italic>a</italic>), survival rate (<italic>p</italic>) and sporogonic cycle (<italic>n</italic>) (Figure ##FIG##1##2##) did not vary spatially. On the contrary, as <italic>An. hyrcanus </italic>presence and density depend on the biotopes and the season [##REF##17313664##4##,##REF##18092973##5##], the vector-host ratio <italic>m </italic>(and hence the human biting rate (<italic>ma</italic>)) (Figure ##FIG##1##2##) presents a strong spatial heterogeneity (Table ##TAB##0##1##). Spatial variations were assessed based on a Geographic Information System (GIS) computing data for each 30 meter-wide pixel in the Camargue [##REF##18302749##22##].</p> <p>All parameters were considered time-dependent, except the <italic>Anopheles </italic>anthropophily (h) and, as a practical approximation, the survival rate (p) (Table ##TAB##0##1##). Time-variations, which are detailed below, were described at a monthly time step.</p> <title>Human biting rate (<italic>ma</italic>)</title> <p>The presence and density of <italic>An. hyrcanus </italic>were inferredusing remote sensing, entomological adults and larvae collections [##REF##18302749##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 <italic>An. hyrcanus </italic>captured in the same pixel with CDC-light traps+CO<sub>2 </sub>showed a highly significant linear regression, allowing us to infer, using key environmental factors, the maximum number of <italic>An. hyrcanus </italic>captured during the year with CDC-light traps+CO<sub>2 </sub>for each pixel in the Camargue [##REF##18302749##22##].</p> <p>The mean annual dynamics of <italic>An. hyrcanus </italic>in the Camargue was inferred from the results of several capture campaigns conducted during several years, in several places, using several capture methods (Table ##TAB##1##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 <italic>An. hyrcanus </italic>captured during the year refers to the month of August. Hence, the spatial distribution of the annual maximum number of <italic>An. hyrcanus </italic>was combined with the mean annual dynamics in order to determine, for each pixel and for each month in the Camargue, the number of <italic>An. hyrcanus </italic>captured with CDC-light traps+CO<sub>2</sub>.</p> <p>From this, the human biting rate (<italic>ma</italic>) 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 ##TAB##2##3## and ##TAB##3##4##). For example, human landing captures carried out on the 10<sup>th </sup>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 (r<sup>2 </sup>= 0.66, p &lt; 0.01). Indeed, human biting rate (<italic>ma</italic>) was extrapolated from the following formula: <italic>ma </italic>= (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 <italic>ma </italic>in order to obtain <italic>ma </italic>for the entire night (Table ##TAB##2##3##).</p> <p>The result is the spatio-temporal distribution of <italic>ma </italic>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 ##FIG##3##4##). Variability and uncertainty were taken into account for each step leading to the assessment of variability and uncertainty of <italic>ma </italic>for each pixel.</p> <title><italic>Anopheles hyrcanus </italic>anthropophily (<italic>h</italic>)</title> <p>(Figure ##FIG##1##2##) was estimated from the comparison between human landing, light traps and horse bait trap results for <italic>An. hyrcanus </italic>and other <italic>Anopheles </italic>species [##REF##17313664##4##], and was fitted as follows: Pert distribution (0.4; 0.5; 0.8)(Figure ##FIG##4##5##). It was assumed that <italic>h </italic>did not vary throughout the year.</p> <title>The length of the trophogonic cycle (<italic>u</italic>)</title> <p>(Figure ##FIG##1##2##) was calculated using the following formula, which estimates the length of time of blood digestion: <italic>u </italic>= <italic>f</italic><sub>1</sub>/(<italic>T-g</italic><sub>1</sub>) where <italic>f</italic><sub>1 </sub>and <italic>g</italic><sub>1 </sub>are factors depending on humidity and <italic>T </italic>is the temperature. <italic>f</italic><sub>1 </sub>and <italic>g</italic><sub>1 </sub>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) [##REF##13885800##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 [##REF##13885800##23##]. Considering expert knowledge and bibliographical data [##UREF##11##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 ##TAB##4##5##.</p> <p>In order to illustrate this, the distribution of the length of the trophogonic cycle was determined for August (Figure ##FIG##5##6##).</p> <title>The daily survival rate (<italic>p</italic>)</title> <p>(Figure ##FIG##1##2##) was estimated from parity rates (P) using the following formula: <italic>p </italic>= P<sup>1/u </sup>[##REF##13214009##25##] which is relevant in the case of stable populations. Considering that the <italic>An. hyrcanus </italic>population increases and decreases progressively according to the <italic>Anopheles </italic>biology, the approximation that populations were stable during the summer was used [##REF##18092973##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) [##REF##17313664##4##]. They were fitted as indicated in Table ##TAB##5##6##, and the mean value was calculated from the five probability distributions in order to obtain the mean <italic>An. hyrcanus </italic>parity rate in the Camargue from June to September 2005.</p> <p>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 <italic>p </italic>= P<sup>1/u</sup>, a mean <italic>p </italic>value of 0.79 was obtained (Figure ##FIG##6##7##). Due to difficulty in obtaining a precise value for parity rates for the entire year, it was considered that <italic>p </italic>did not vary from March to October in this model.</p> <title>The length of the sporogonic period (<italic>n</italic>)</title> <p>(Figure ##FIG##1##2##) was calculated as follows: <italic>n </italic>= <italic>f</italic><sub>2</sub>/(<italic>T-g</italic><sub>2</sub>), where <italic>T </italic>is temperature and <italic>f</italic><sub>2 </sub>and <italic>g</italic><sub>2 </sub>are <italic>Plasmodium </italic>species-dependent factors [##UREF##12##26##]. <italic>f</italic><sub>2 </sub>and <italic>g</italic><sub>2 </sub>were experimentally evaluated at 111°C-days and 16°C, 105°C-days and 14.5°C for <italic>Plasmodium falciparum </italic>and <italic>P. vivax</italic>, respectively. However, Grassi and MacDonald estimated threshold temperatures under which the sporogonic development is not completed: 18–19°C for <italic>P. falciparum</italic>, 15–17°C for <italic>P. vivax </italic>[##UREF##11##24##,##UREF##13##27##].</p> <title>Entomological data: infectivity</title> <p><italic>Anopheles hyrcanus </italic>susceptibility was estimated based on data derived from experimental membrane-feeding experiments conducted with cultured <italic>P. falciparum </italic>in the laboratories of Radboud University Nijmegen Medical Centre (The Netherlands) following a routine protocol [##REF##7048654##28##,##REF##15463015##29##]. Among 350 <italic>An. hyrcanus </italic>tested, none were found able to transmit tropical <italic>P. falciparum </italic>strains. Hence, <italic>An. hyrcanus </italic>susceptibility to <italic>P. falciparum </italic>was fitted with a beta distribution: Beta (1; 351) (Figure ##FIG##7##8##).</p> <p><italic>Anopheles hyrcanus </italic>susceptibility to <italic>P. vivax </italic>has not yet been tested, but this species has been considered a vector of <italic>P. vivax </italic>in Afghanistan for more than 30 years [##REF##1172640##30##,##REF##18252122##31##]. Genetic comparison based on ITS2 sequences between French, Afghan, Turkish and Iranian specimens concluded that they were identical. Hence, it was considered that <italic>An. hyrcanus </italic>was susceptible to tropical <italic>P. vivax </italic>strains, and estimation of its susceptibility, based on expert knowledge, was fitted as follows: Pert distribution (0.05; 0.20; 0.70) (Figure ##FIG##8##9##).</p> <p><italic>An. hyrcanus </italic>susceptibility to either <italic>P. falciparum </italic>or <italic>P. vivax </italic>was considered to be homogeneous through time and space (Table ##TAB##0##1##).</p> <p><italic>Anopheles hyrcanus </italic>susceptibility to <italic>Plasmodium ovale </italic>and <italic>Plasmodium malariae </italic>was not inferred, as very little is known about European <italic>Anopheles </italic>susceptibility to these two species.</p> <title>Parasitological data: vulnerability</title> <p>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 [##UREF##2##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.</p> <p>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, <italic>Plasmodium </italic>species, patients' residence for the 2004–2005 period and contamination place.</p> <p>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 [##UREF##14##32##,##REF##12244913##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.</p>
<title>Results</title> <title>Vulnerability</title> <p>In 2004 and 2005, 657 imported cases were diagnosed in the region (corresponding to a total of 1251 estimated imported cases), among which <italic>P. falciparum, P. vivax, P. ovale </italic>and <italic>P. malariae </italic>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 <italic>An. hyrcanus </italic>could be present.</p> <p>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.</p> <title>Entomological risk (receptivity*infectivity)</title> <p>The mean value of the entomological risk was assessed for <italic>P. falciparum </italic>species from June to September, and for <italic>P. vivax </italic>for the month of August for each pixel of the Camargue map (Figure ##FIG##9##10##, ##FIG##10##11##, ##FIG##11##12##, ##FIG##12##13##, ##FIG##13##14##) (the entomological risk was not assessed for <italic>P. ovale </italic>and <italic>P. malariae </italic>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.</p> <p>Strong differences were observed in the entomological risk for the two <italic>Plasmodium </italic>species: <italic>P. falciparum </italic>transmission risk estimate ranges from 0 to more than 1 although <italic>P. vivax </italic>transmission risk estimate ranges from 0 to more than 100.</p> <title>Uncertainties of the risk estimate and sensitivity analysis</title> <p>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<sup>-6 </sup>to 20, with 95% of the values being between 0.016 and 4.7 for pixel having a mean value of about 1.</p> <p>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 ##FIG##14##15##).</p> <p>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 <italic>P. vivax </italic>entomological risk estimate showed the equal importance of the susceptibility, the survival rate and the human biting rate.</p>
<title>Discussion</title> <p>The entomological risk index used in this article refers to the risk of transmission (Figure ##FIG##0##1##): the risk of transmission being high when the entomological risk is high. Results clearly indicate a maximum risk of <italic>P. falciparum </italic>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 <italic>An. hyrcanus </italic>dynamics and density [##REF##18092973##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 <italic>An. hyrcanus </italic>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.</p> <p>The risk of <italic>P. vivax </italic>transmission is more than one hundred times higher than the risk of <italic>P. falciparum </italic>transmission, which is due to infectivity and the length of sporogonic cycle differences.</p> <p>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 <italic>An. hyrcanus</italic>. 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.</p> <p>Although all of France faces a large number of imported cases [##UREF##2##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.</p> <p>The entomological risk, referring to the risk of transmission in this article (Figure ##FIG##0##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 <italic>Plasmodium </italic>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 <italic>Plasmodium</italic>.</p> <p>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 [##UREF##13##27##,##REF##1678572##34##, ####UREF##15##35##, ##UREF##16##36####16##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 [##REF##9709243##18##].</p> <p>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.</p> <p>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.</p> <p>It has been shown recently the impact of anthropogenic changes on potential malaria vectors in the Camargue over the last 60 years [##REF##18258028##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.</p> <p>The sensitivity analysis underlined factors responsible for entomological risk uncertainty and variation, of which susceptibility, the survival rate and the human biting rate of <italic>An. hyrcanus </italic>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.</p> <p>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.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>The Camargue region is a former malaria endemic area, where potential <italic>Anopheles </italic>vectors are still abundant. Considering the importation of <italic>Plasmodium </italic>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.</p> <title>Methods</title> <p>Receptivity (vectorial capacity) and infectivity (vector susceptibility) were inferred using an innovative probabilistic approach and considering both <italic>Plasmodium falciparum </italic>and <italic>Plasmodium vivax</italic>. 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.</p> <title>Results</title> <p>The entomological risk presented large spatial, temporal and <italic>Plasmodium </italic>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.</p> <title>Conclusion</title> <p>The current risk of malaria re-emergence seems negligible due to the very low number of imported <italic>Plasmodium</italic>. This model demonstrated its efficiency for mosquito-borne diseases risk assessment.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>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.</p>
<title>Acknowledgements</title> <p>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 <italic>An. hyrcanus </italic>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.</p> <p>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 <ext-link ext-link-type="uri" xlink:href="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.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Risk of malaria re-emergence.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Entomological factors.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>Probability distributions.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Spatial distribution of the human biting rate in the Camargue in August</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>Distribution of <italic>An. hyrcanus </italic>anthropophily.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>Distribution of the length, in days, of the trophogonic cycle in August.</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>Distribution of the <italic>An. hyrcanus </italic>daily survival rate.</p></caption></fig> <fig position="float" id="F8"><label>Figure 8</label><caption><p>Distribution of <italic>An. hyrcanus </italic>susceptibility to <italic>P. falciparum</italic>.</p></caption></fig> <fig position="float" id="F9"><label>Figure 9</label><caption><p>Distribution of <italic>An. hyrcanus </italic>susceptibility to <italic>P. vivax</italic>.</p></caption></fig> <fig position="float" id="F10"><label>Figure 10</label><caption><p><bold>Spatial variations of <italic>P. falciparum </italic>transmission risk estimate (ranging from 0 to more than 1) in June in the Camargue</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F11"><label>Figure 11</label><caption><p><bold>Spatial variations of <italic>P. falciparum </italic>transmission risk estimate (ranging from 0 to more than 1) in July in the Camargue</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F12"><label>Figure 12</label><caption><p><bold>Spatial variations of <italic>P. falciparum </italic>transmission risk estimate (ranging from 0 to more than 1) in August in the Camargue</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F13"><label>Figure 13</label><caption><p><bold>Spatial variations of <italic>P. falciparum </italic>transmission risk estimate (ranging from 0 to more than 1) in September in the Camargue</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F14"><label>Figure 14</label><caption><p><bold>Spatial variations of the <italic>P. vivax </italic>transmission risk estimate (ranging from 0 to more than 100) in August in the Camargue</bold>. Classes were arbitrary chosen with a logarithmic scale.</p></caption></fig> <fig position="float" id="F15"><label>Figure 15</label><caption><p>Sensitivity analysis.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Space, time and <italic>Plasmodium </italic>species-dependent variations in receptivity, infectivity and vulnerability.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td/><td align="center" colspan="3">Variation</td></tr><tr><td/><td/><td colspan="3"><hr/></td></tr><tr><td align="center" colspan="2">Components</td><td align="center">Space<break/> dependent</td><td align="center">Time<break/> dependent</td><td align="center"><italic>Plasmodium</italic>-species<break/> dependent</td></tr></thead><tbody><tr><td align="left">Receptivity</td><td align="left">m (vector-host ratio)</td><td align="center">YES</td><td align="center">YES</td><td align="center">NO</td></tr><tr><td/><td align="left">a (vector-biting rate)</td><td align="center">NO</td><td align="center">YES</td><td align="center">NO</td></tr><tr><td/><td align="left">h (anthropophily)</td><td align="center">NO</td><td align="center">NO</td><td align="center">NO</td></tr><tr><td/><td align="left">u (length of the trophogonic cycle)</td><td align="center">NO</td><td align="center">YES</td><td align="center">NO</td></tr><tr><td/><td align="left">p (survival rate)</td><td align="center">NO</td><td align="center">NO</td><td align="center">NO</td></tr><tr><td/><td align="left">n (sporogonic cycle)</td><td align="center">NO</td><td align="center">YES</td><td align="center">YES</td></tr><tr><td align="left">Infectivity</td><td align="left">b (susceptibility)</td><td align="center">NO</td><td align="center">NO</td><td align="center">YES</td></tr><tr><td align="left">Vulnerability</td><td/><td align="center">YES</td><td align="center">YES</td><td align="center">YES</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Total number of <italic>An. hyrcanus </italic>collected per month and associated percentages (bracket)</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center" colspan="8">Capture places</td></tr></thead><tbody><tr><td/><td align="center">Méjannes</td><td align="center">Aimargues</td><td align="center">Mourgues</td><td align="center">Pont de gau</td><td align="center">Carbonnière</td><td align="center">Vigueirat</td><td align="center">Western<break/> Camargue<sup>1</sup></td><td align="center">Tour du Vallat</td></tr><tr><td colspan="9"><hr/></td></tr><tr><td align="center">Year</td><td align="center">2005</td><td align="center">2005</td><td align="center">2005</td><td align="center">2005</td><td align="center">2005</td><td align="center">2005</td><td align="center">1969–1994</td><td align="center">2004</td></tr><tr><td align="center">Capture method</td><td align="center">CO2 traps</td><td align="center">CO2 traps</td><td align="center">CO2 traps</td><td align="center">CO2 traps</td><td align="center">CDC-LT + CO2</td><td align="center">CDC-LT + CO2</td><td align="center">Human landing catch</td><td align="center">Horse bait trap</td></tr><tr><td align="center">Reference</td><td align="center">G. L'Ambert, unpublished data</td><td align="center">G. L'Ambert, unpublished data</td><td align="center">A. Carron, unpublished data</td><td align="center">A. Carron, unpublished data</td><td align="center">[##REF##17313664##4##]</td><td align="center">[##REF##17313664##4##]</td><td align="center">EID-Méditérannée, unpublished data<sup>2</sup></td><td align="center">[##REF##17017231##3##]</td></tr><tr><td align="center">March</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0 (0)</td><td align="center">5 (0)</td><td align="center">(0)</td><td align="center">-</td></tr><tr><td align="center">April</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0 (0)</td><td align="center">7 (0)</td><td align="center">(0)</td><td align="center">-</td></tr><tr><td align="center">May</td><td align="center">1 (0)</td><td align="center">0 (0)</td><td align="center">2 (1)</td><td align="center">2 (0)</td><td align="center">16 (0)</td><td align="center">296 (0)</td><td align="center">(0)</td><td align="center">0 (0)</td></tr><tr><td align="center">June</td><td align="center">27 (0)</td><td align="center">0 (0)</td><td align="center">18 (7)</td><td align="center">28 (3)</td><td align="center">673 (12)</td><td align="center">6737 (6)</td><td align="center">(4)</td><td align="center">3 (0)</td></tr><tr><td align="center">July</td><td align="center">357 (5)</td><td align="center">0 (0)</td><td align="center">104 (42)</td><td align="center">107 (12)</td><td align="center">374 (7)</td><td align="center">17739 (16)</td><td align="center">(11)</td><td align="center">38 (1)</td></tr><tr><td align="center">August</td><td align="center">5933 (83)</td><td align="center">102 (95)</td><td align="center">94 (38)</td><td align="center">730 (84)</td><td align="center">3257 (59)</td><td align="center">61315 (55)</td><td align="center">(43)</td><td align="center">4101 (71)</td></tr><tr><td align="center">September</td><td align="center">857 (12)</td><td align="center">5 (5)</td><td align="center">32 (13)</td><td align="center">7 (1)</td><td align="center">1228 (22)</td><td align="center">25708 (23)</td><td align="center">(42)</td><td align="center">1640 (28)</td></tr><tr><td align="center">October</td><td align="center">0 (0)</td><td align="center">0 (0)</td><td align="center">0 (0)</td><td align="center">0 (0)</td><td align="center">3 (0)</td><td align="center">124 (0)</td><td align="center">(0)</td><td align="center">4 (0)</td></tr><tr><td align="center">Total</td><td align="center">7175</td><td align="center">107</td><td align="center">250</td><td align="center">874</td><td align="center">5551</td><td align="center">111919</td><td/><td align="center">5786</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Number of <italic>An. hyrcanus </italic>captured during human landing catches sessions.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Area</td><td align="center">Pixel</td><td align="center">Adult<break/>abundance<break/>  index</td><td align="center">08/08/2005</td><td align="center">10/08/2005</td><td align="center">19/09/2005</td><td align="center">21/09/2005</td><td align="center">04/10/2005</td><td align="center">23/08/2006</td><td align="center">05/09/2006</td></tr></thead><tbody><tr><td align="center">Vigueirat</td><td align="center">Hu1</td><td align="center">8676</td><td align="center">6312</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">14</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">Hu2</td><td align="center">122</td><td align="center">26</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">Hu3</td><td align="center">308</td><td align="center">2</td><td align="center">-</td><td align="center">9</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">Hu4</td><td align="center">2430</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">67</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu5</td><td align="center">242</td><td align="center">-</td><td align="center">304</td><td align="center">-</td><td align="center">25</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu6</td><td align="center">155</td><td align="center">-</td><td align="center">2</td><td align="center">-</td><td align="center">1</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu6</td><td align="center">155</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">0</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu8</td><td align="center">687</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">39</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu5</td><td align="center">242</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">156</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu7</td><td align="center">5348</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">16</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">Hu7</td><td align="center">5348</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">4</td><td align="center">-</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Number of <italic>An. hyrcanus </italic>captured during light trap sessions. Adult abundance index relative to the pixel where the capture was conducted is indicated.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Area</td><td align="center">Pixel</td><td align="center">Adult<break/>abundance<break/>  index</td><td align="center">08/08/2005</td><td align="center">10/08/2005</td><td align="center">19/09/2005</td><td align="center">04/10/2005</td><td align="center">23/08/2006</td><td align="center">05/09/2006</td><td align="center">21/09/2005</td></tr></thead><tbody><tr><td align="center">Vigueirat</td><td align="center">LT1</td><td align="center">346</td><td align="center">218</td><td align="center">-</td><td align="center">346</td><td align="center">6</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT2</td><td align="center">208</td><td align="center">70</td><td align="center">-</td><td align="center">208</td><td align="center">1</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT3</td><td align="center">6594</td><td align="center">3546</td><td align="center">-</td><td align="center">6594</td><td align="center">20</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT4</td><td align="center">10920</td><td align="center">15971</td><td align="center">-</td><td align="center">10920</td><td align="center">15</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT5</td><td align="center">4898</td><td align="center">2016</td><td align="center">-</td><td align="center">4898</td><td align="center">32</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT6</td><td align="center">1769</td><td align="center">180</td><td align="center">-</td><td align="center">1769</td><td align="center">3</td><td align="center">-</td><td align="center">1200</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT7</td><td align="center">3732</td><td align="center">3645</td><td align="center">-</td><td align="center">3732</td><td align="center">4</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Vigueirat</td><td align="center">LT8</td><td align="center">2681</td><td align="center">1006</td><td align="center">-</td><td align="center">2681</td><td align="center">24</td><td align="center">-</td><td align="center">-</td><td align="center">-</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT1</td><td align="center">492</td><td align="center">-</td><td align="center">254</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">8</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT2</td><td align="center">1299</td><td align="center">-</td><td align="center">102</td><td align="center">-</td><td align="center">-</td><td align="center">101</td><td align="center">-</td><td align="center">10</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT3</td><td align="center">1250</td><td align="center">-</td><td align="center">39</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT4</td><td align="center">433</td><td align="center">-</td><td align="center">20</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">3</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT5</td><td align="center">176</td><td align="center">-</td><td align="center">288</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">3</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT6</td><td align="center">705</td><td align="center">-</td><td align="center">30</td><td align="center">-</td><td align="center">-</td><td align="center">78</td><td align="center">-</td><td align="center">19</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT7</td><td align="center">721</td><td align="center">-</td><td align="center">17</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">15</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT8</td><td align="center">0</td><td align="center">-</td><td align="center">1</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">0</td></tr><tr><td align="center">Carbonnière</td><td align="center">LT9</td><td align="center">5348</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">-</td><td align="center">28</td><td align="center">-</td><td align="center">-</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T5"><label>Table 5</label><caption><p>Length of the trophogonic cycle</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center">Month</td><td align="center">Mean temperature <break/>(2002–2006)</td><td align="center">Length of trophogonic cycle<break/> (from Shlenova and Detinova)</td><td align="center">Pert distribution</td></tr></thead><tbody><tr><td align="center">March</td><td align="center">10.9</td><td align="center">38.4</td><td/></tr><tr><td align="center">April</td><td align="center">14</td><td align="center">9.9</td><td align="center">(9; 10; 11)</td></tr><tr><td align="center">May</td><td align="center">17.8</td><td align="center">5.6</td><td align="center">(4.5; 5.5; 6.5)</td></tr><tr><td align="center">June</td><td align="center">23</td><td align="center">3.8</td><td align="center">(3; 4; 5)</td></tr><tr><td align="center">July</td><td align="center">24.8</td><td align="center">3.4</td><td align="center">(3; 3.5; 5)</td></tr><tr><td align="center">August</td><td align="center">24</td><td align="center">3.6</td><td align="center">(3; 3.5; 5)</td></tr><tr><td align="center">September</td><td align="center">20.3</td><td align="center">4.5</td><td align="center">(3.5; 4.5; 5.5)</td></tr><tr><td align="center">October</td><td align="center">17.1</td><td align="center">6.1</td><td align="center">(5; 6; 7)</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T6"><label>Table 6</label><caption><p>Parity rates observed in the Camargue in 2005</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="center">Carbonnière</td><td align="center">Marais du Vigueirat</td></tr></thead><tbody><tr><td align="center">June</td><td align="center">-</td><td align="center">Beta (18; 9)</td></tr><tr><td align="center">July</td><td align="center">Beta (20; 69)</td><td align="center">Beta (29; 91)</td></tr><tr><td align="center">September</td><td align="center">Beta (89; 45)</td><td align="center">Beta (40; 80)</td></tr></tbody></table></table-wrap>
<table-wrap-foot><p>Note:</p><p><sup>1 </sup>approximatively 12 points</p><p><sup>2 </sup>mean percentage</p></table-wrap-foot> <table-wrap-foot><p> Adult abundance index relative to the pixel where the capture was conducted is indicated. Human landing results are relative to the hour following sunset which represent the <italic>An. hyrcanus </italic>aggressiveness peak [##REF##17313664##4##] and approximately 80% of the total number of <italic>An. hyrcanus </italic>captured during the entire night. Each line corresponds to a human landing catch session carried out by one volunteer belonging to the research team (two people were present at Hu6 and Hu7).</p></table-wrap-foot> <table-wrap-foot><p> Light trap captures are relative to the entire night.</p></table-wrap-foot>
<graphic xlink:href="1475-2875-7-147-1"/> <graphic xlink:href="1475-2875-7-147-2"/> <graphic xlink:href="1475-2875-7-147-3"/> <graphic xlink:href="1475-2875-7-147-4"/> <graphic xlink:href="1475-2875-7-147-5"/> <graphic xlink:href="1475-2875-7-147-6"/> <graphic xlink:href="1475-2875-7-147-7"/> <graphic xlink:href="1475-2875-7-147-8"/> <graphic xlink:href="1475-2875-7-147-9"/> <graphic xlink:href="1475-2875-7-147-10"/> <graphic xlink:href="1475-2875-7-147-11"/> <graphic xlink:href="1475-2875-7-147-12"/> <graphic xlink:href="1475-2875-7-147-13"/> <graphic xlink:href="1475-2875-7-147-14"/> <graphic xlink:href="1475-2875-7-147-15"/>
[{"surname": ["Rodhain", "Charmot"], "given-names": ["F", "G"], "article-title": ["Evaluation des risques de reprise de transmission du paludisme en France"], "source": ["Med Mal Infect"], "year": ["1982"], "volume": ["12"], "fpage": ["231"], "lpage": ["236"]}, {"surname": ["Sautet"], "given-names": ["J"], "article-title": ["A propos d'une \u00e9pid\u00e9mie de paludisme en Camargue"], "source": ["Mars Med"], "year": ["1944"], "volume": ["2"], "fpage": ["53"], "lpage": ["64"]}, {"surname": ["Legros", "Arnaud", "El Minouni", "Danis"], "given-names": ["F", "A", "B", "M"], "article-title": ["Paludisme d'importation en France m\u00e9tropolitaine : donn\u00e9es \u00e9pid\u00e9miologiques 2001-2004."], "source": ["Bulletin \u00e9pid\u00e9miologique hebdomadaire"], "year": ["2006"], "volume": ["32"], "fpage": ["235"], "lpage": ["236"]}, {"collab": ["WHO-CISID"], "article-title": ["WHO - CISID-malaria"]}, {"surname": ["Ejov"], "given-names": ["M"], "article-title": ["Imported malaria and the risk of malaria reintroduction in Europe"], "source": ["CD News, Quaterly Communicable Disease Report, WHO Europe"], "year": ["2004"], "volume": ["34"], "fpage": ["3"], "lpage": ["7"]}, {"collab": ["WHO"], "article-title": ["Expert Committee on Malaria, Twelfth Report"], "source": ["Technical Report Series, N\u00b0324"], "year": ["1966"], "volume": ["324"]}, {"collab": ["WHO"], "article-title": ["Receptivity to malaria and other parasitic diseases"], "source": ["Euro Reports and Studies, 15"], "year": ["1979"], "publisher-name": ["Copenhagen "], "fpage": ["108"]}, {"surname": ["Dervieux", "Allard", "Naizot", "Auda", "Richard H, Vignot A"], "given-names": ["A", "P", "T", "Y"], "article-title": ["La Camargue est-elle un polder? Changement et stabilit\u00e9 dans le delta du Rh\u00f4ne depuis le d\u00e9but du XIX\u00e8me si\u00e8cle"], "source": ["Equilibre et ruptures dans les \u00e9cosyst\u00e8mes durant les 20 derniers mill\u00e9naires en Europe de l'Ouest"], "year": ["2002"], "volume": ["1"], "publisher-name": ["Besan\u00e7on "], "fpage": ["315"], "lpage": ["329"]}, {"surname": ["Mathevet", "Mouret", "Mesl\u00e9ard"], "given-names": ["R", "JC", "F"], "article-title": ["Du mondial au local : originalit\u00e9 et d\u00e9pendances de la riziculture camarguaise."], "source": ["Faire savoirs"], "year": ["2002"], "volume": ["2"], "fpage": ["49"], "lpage": ["56"]}, {"surname": ["Vose"], "given-names": ["DJ"], "source": ["Risk analysis. A quantitative guide"], "year": ["2000"], "publisher-name": ["Chichester, England , John Wiley and Sons, LTD"]}, {"surname": ["Anderson", "Hattis"], "given-names": ["EL", "D"], "article-title": ["When and how can you specify a probability distribution when you don't know much. Foudations: A. Uncertainty and variability."], "source": ["Risk analysis"], "year": ["1999"], "volume": ["19"], "fpage": ["47"], "lpage": ["49"]}, {"surname": ["Jetten", "Takken"], "given-names": ["TH", "W"], "source": ["Anophelism without malaria in Europe. A review of the ecology and distribution of the genus "], "italic": ["Anopheles"], "year": ["1994"], "publisher-name": ["Wageningen, The Netherlands , Landbouwuniversiteit Wageningen (Wageningen Agricultural University)"], "fpage": ["No. 94"], "lpage": ["5, 69"]}, {"surname": ["Moshkovsky", "Rashina"], "given-names": ["S", "M"], "source": ["Epid\u00e9miologie et parasitologie m\u00e9dicale \u00e0 l'usage des entomologistes"], "year": ["1951"], "publisher-name": ["Moscow "]}, {"surname": ["Macdonald"], "given-names": ["G"], "source": ["The epidemiology and control of malaria"], "year": ["1957"], "publisher-name": ["London , Oxford University Press"], "fpage": ["201"]}, {"surname": ["Legros", "Fromage", "Ancelle", "Burg", "Janot", "Maisonneuve", "Danis"], "given-names": ["F", "M", "T", "E", "C", "P", "M"], "article-title": ["Enqu\u00eate nationale de recensement des cas de paludisme d'importation en France m\u00e9tropolitaine pour l'ann\u00e9e 1997."], "source": ["Bulletin \u00e9pid\u00e9miologique hebdomadaire"], "year": ["1999"], "volume": ["11"], "fpage": ["41"], "lpage": ["42"]}, {"surname": ["Martens", "Kovats", "Nijhof", "de Vries", "Livermore", "Bradley", "Cox", "McMichael"], "given-names": ["P", "RS", "S", "P", "MTJ", "DJ", "J", "AJ"], "article-title": ["Climate change and future populations at risk of malaria"], "source": ["Glob Environ Change"], "year": ["1999"], "volume": ["9"], "fpage": ["S89"], "lpage": ["S107"], "pub-id": ["10.1016/S0959-3780(99)00020-5"]}, {"surname": ["Lindsay", "Thomas"], "given-names": ["SW", "CJ"], "article-title": ["Global warning and risk of vivax malaria in Great Britain"], "source": ["Global change and human health"], "year": ["2001"], "volume": ["2"], "fpage": ["80"], "lpage": ["84"], "pub-id": ["10.1023/A:1011995115713"]}]
Malar J. 2008 Aug 1; 7:147