Assessment of the risk of malaria re-introduction in the Maremma plain (Central Italy) using a multi-factorial approach
- Roberto Romi1Email author,
- Daniela Boccolini1,
- Roberto Vallorani2, 4,
- Francesco Severini1,
- Luciano Toma1,
- Maurizio Cocchi3,
- Angelo Tamburro3,
- Gianni Messeri4,
- Antonio Crisci4,
- Luca Angeli2,
- Roberto Costantini2,
- Irene Raffaelli3,
- Giorgio Pontuale3,
- Isabelle Thiéry5,
- Annie Landier5,
- Gilbert Le Goff6,
- Anna Maria Fausto7 and
- Marco Di Luca1
© Romi et al; licensee BioMed Central Ltd. 2012
Received: 12 December 2011
Accepted: 30 March 2012
Published: 30 March 2012
KeywordsMosquito-borne diseases Residual anophelism Anopheles labranchiae Vectorial capacity Climate change Plasmodium falciparum Experimental infection
In recent years, the increase in globalization , the rise in the average temperature of the earth together with an increasing frequency and intensity of extreme weather events, as storms, floods and droughts [2, 3], and the environmental changes induced by human activities , have raised the concern about the possible introduction or reintroduction of Vector Borne Diseases in Countries where these were absent or eradicated . These considerations, coupled with the recent spread of some mosquito vector borne diseases in Europe [6, 7] and the increasing number of imported malaria cases recorded in the Continent  have renewed interest in the possible reintroduction of malaria in Southern Europe [7–9], particularly in the countries facing the Western Mediterranean Basin, where potential Anopheline vectors are still present [10–13]. Moreover, in recent years autochthonous malaria cases have been sporadically reported in Italy, France, Spain and Greece [14–20].
In 2005, a five-year study was implemented in Italy, as well as in other South European countries, with the aim to assess the status of the local potential malaria vectors and the possible re-introduction of malaria transmission [21–25]. In Italy, the selected study area was the Maremma plain, a region that was hyperendemic for malaria until 60 years ago [26–28] and that more recently was recognized as the major "at risk" area for the malaria reintroduction into Italy [14, 29, 30].
In Maremma, after the malaria eradication campaign (1947-1951), Anopheles labranchiae, the main endophilic vector of the Anopheles maculipennis complex was dramatically reduced in abundance. However, in subsequent years, the species has progressively re-colonized most of the area coming back to substantial densities [31–33]. This was mainly due to the introduction of intensive rice cultivation in the early 1970s. Since then, Maremma has been subjected to continuous entomological surveillance that was intensified after 1997, when an autochthonous Plasmodium vivax malaria case, transmitted by An. labranchiae, occurred in the Province of Grosseto . The studies carried out in the area since eradication, provides a database that allowed a follow-up the history of malaria and its vectors in Maremma over the past 60 years. Starting from the findings of the most recent entomological and environmental studies [23, 34], the present study was chosen to evaluate the malariogenic potential of the area using a multifactorial approach.
receptivity of the area, given by the presence, distribution, seasonal abundance and bionomics of the potential vector;
susceptibility of the vector, that is its ability to become infected with Plasmodium vivax and Plasmodium falciparum;
vulnerability of the territory, that is the possible introduction of malaria reservoirs, given by the number of gametocyte carriers able to infect the vector and present in the study area during months favourable to malaria transmission.
Approaches to evaluate these parameters were:
1) Field collection of further entomological data (bionomics, distribution, abundance) for mosquitoes of the An. maculipennis complex; 2) investigation of seasonal dynamics of the vector through the implementation of a weather-based statistical dynamic model; 3) production of a distribution/predictive map of An. labranchiae across the study area; 4) evaluation of the length of the possible transmission season for P. vivax and P. falciparum; 5) assessment of the vector competence of the species to P. falciparum by artificial infection; 6) evaluation of the vectorial capacity of An. labranchiae in the site where the species is most abundant; 7) risk assessment related to the possibility that the vector may feed on gametocyte carriers occasionally circulating in the study area.
Study area and collection sites
Location and features of the sites selected for assessing malaria risk in Maremma in 2005-2009
An. maculipennis s.l. range of abundance
An. labranchiae prevalence (%)
1) Principina (GR)
Farm: coastal rice-fields, 256 ha*
> 60 (100-500)
2) San Donato (GR)
Farm: coastal rice fields, 114 ha*
> 60 (100-500)
3) Val di Merse (SI)
Farm: hilly rice fields, 110 ha*
> 60 (100-500) ( An. lab . < 60)
4) The Ampio (GR)
Farm: inner-plain area
< 60 (10-20)
5) Monte Antico (GR)
Farm: inner-hilly area
< 60 (5-10)
6) Valfragida (VT)
Farm: coastal plain area
< 60 (10-40)
7) Forca di Parma (VT)
Farm: inner plain area
> 60 (20-80)
8) La Parrina (VT)
Farm: coastal plain area
> 60 (20-80)
9) Alberese Natural Park (GR)
Coastal Natural area 10,000 ha
< 60 (20-30)
10) Diaccia Botrona (GR)
Brackish water coastal marshes, 236 ha
< 60 (5-10)
11) Artificial water collections (GR;SI, VT)
Basins for agricultural, commercial and leisure purposes, 0.1-4.6 ha
Mosquito collection and laboratory processing
Collections of An. maculipennis s.l. mosquitoes were performed between 2005 and 2009 during 40 surveys. In particular, sites 1, 3, 4 and 6 were visited by 28 fortnightly surveys carried out from April 2005 to October 2006 , while in 2007-2009 these sites and the remaining ones were subject to sporadic surveys (1-4 times) in June-August. Mosquito larvae were collected by an enamel standard 500 ml dipper. The number of dipping stations and dips by station was adequate to the type and size of the breeding sites visited according to a standardized protocol [23, 29]. Adult collections, using manual or battery-powered aspirators, were mainly targeted at resting females in animal shelters; other kinds of premises, such as haylofts, woodshed, fodder and tools depots represented less than 5% of the premises inspected. A minor fraction of An. maculipennis s.l. females were collected by CDC light/CO2 traps. Three night catches on human bait were seasonally carried out in late June, mid July and late August, between 2005 and 2007, respectively in Principina, S. Donato, and Val di Merse. Catches were performed as described in Romi et al. . Mosquitoes were analyzed for species identification, blood meal source and age population structure as described in Di Luca et al. . Previous entomological data used for comparisons are from the Operative Unit of Environmental Zoology, AUSL 9, Grosseto.
Computing seasonal dynamic of An. labranchiae populations
A binomial model was developed for the less productive sites (4 and 6) categorizing the entomological data with the threshold of one adult ("A" absence category/"B" presence category), while the multinomial model, developed for the most productive site 1, was categorize with the thresholds of one and 60 adults ("A" absence category/"B" medium-low presence category/"C" medium-high presence category). The multi-logistic model input were derived from a principal component analysis that leads to an optimal dimensionality reduction of matrix predictors. Akaike Information Criterium (AIC)  was finally used for the statistical model parameters selection in order to optimize model performances and to give a better discrimination in the microclimatic variability among different sites. Concerning model selection the statistical function (i.e. STEPAIC) used is taken from MASS R package [36, 37]. Weekly model outputs consist essentially on the occurrence probability associated to each adult abundance category, hence the forecasted category was chosen as the one with highest probability. The length of development season was estimated as the number of weeks between the first and the last week with a presence category (B or C category) along the year [38–41]. A verification procedure to assess the reliability of the models was implemented for each site using skill score indexes derived from contingency Tables with observed and forecasted values: the BIAS index, the POD index (Probability Of Detection) and the FAR index (False Alarm Ratio) .
Mapping larval index and adult distribution through geospatial statistical analysis
The quantitative mapping of larvae/adult mosquito presence and abundance was implemented following the approach suggested by Tran et al. . The logistic regression model for Anopheles hyrcanus in the Camargue region, that explained the presence of larvae as a function of biotope and distance to the nearest rice field, was modified in order to obtain a larval index (probability of observing one larva in a point of a biotope at least once during the mosquito season) consistent with the collection data reported in Table 1. The logistic regression model for An. labranchiae, used in this study and implemented with Builder tool from ESRI ArcGIS™, consider also the distance from farm with livestock to the main biotopes as a new explanatory variable in addition to those used by Tran et al. Assuming that the abundance of adult mosquitoes is influenced by the presence of breeding sites in the surroundings, the adult index map was derived from the larval index map. Corine Land Cover 2000 (CLC 2000, produced by the European Environment Agency, EEA) spatial data sets for Italy (scale 1:100,000) was used to describe the environmental characteristics likely to influence the spatial distribution of An. labranchiae. Also a natural colour aerial photos (May 2007) with a spatial resolution of 1 metre was used to detect the main biotopes where An. labranchiae larvae and adults were collected, such as rice fields, reed beds, marshes, temporarily flooded rush wetland and clear water. A photo-interpretation was carried out through a workstation with ESRI ArcGIS™ software. Geographical database in shape-file format was used for localizing livestock and intensive cattle breeding farms in Grosseto Province . The highest abundance of An. labranchiae was assumed to be related to rice paddies, and thus the distance to the nearest rice field was computed for each pixel using Geographic Information System (GIS) functionality.
Evaluation of the length of the possible transmission season
Most significant climatic parameters for the 1960-1990 and 2005-2009 periods (Grosseto Airport Weather Station)
Yearly Mean Temperature (°C)
Seasonal Mean Temperature (°C)*
Yearly Mean Rainfall (mm)
Seasonal Mean Rainfall (mm)*
where GDD is growing degree-days with a base temperature of 15°C and 18°C respectively for P. vivax and P. falciparum, R is the rainfall and PET is potential evapotranspiration calculated with the empirical method of Thornthwaite 1948  as a function of mean temperature and latitude. The GMR index shows that a transmission risk exists when its value is equal or higher than 116, that is the value required for one Plasmodium spp generation.
Artificial infection assays
Field samples of An. labranchiae females, collected in site 1 (Figure 1, Table 1) were submitted to the Plateform CEPIA (Institute Pasteur, Paris, France) to artificial blood infection with gametocyte-containing cultures of the P. falciparum NF54 African strain in 2008 and 2009. Production of mature gametocytes and artificial blood infection were performed following procedures described in Mitri et al. . A laboratory colony of Anopheles gambiae (Ngousso, Cameroon) was used as a positive control. Mosquitoes were dissected on 8tand 15 days post-infection to determine prevalence and oocyst load in the midgut. For each experiment all An. labranchiae females and the control An. gambiae strain were starved 24 hours prior to blood feeding. The infected red blood cells containing P falciparum gametocytes complemented with fresh RBC and human AB serum, were deposited in a Parafilm® membrane feeder previously warmed at 37°C. After 15 minutes feeding, unfed An. labranchiae were offered a second blood meal on the next day, when possible. Engorged females were kept at 26 ± 1°C inside small cages and were provided with 10% sucrose until dissection 8 days or 15 days post infection. In the 2009, the detection of sporozoites was carried out at the Institute pour la Recherche et le Développement in Montpellier (IRD, France), using the cut head-thorax from the 15th day survived mosquitoes. DNA extraction was performed by a single-round, multiplex PCR, according to Padley et al. . Legs of all infected females were used for species identification by Multiplex PCR .
Assessment of the vectorial capacity and host feeding preference of An. labranchiae
Vectorial capacity was assessed according to the Macdonald formula  revised by Garret-Jones . The experimental variables needed for estimating it were evaluated as follows: the human biting rate (ma) by night catches on human bait, the human blood index (HBI) by the origin of the blood meal of the engorged females collected early in the morning in different premises and the parity rate by ovarian dissection . The factors temperature-dependent, such as the length of the sporogonic cycle (n) of P. falciparum and P. vivax and the duration of gonotrophic cycle (gc) were calculated according to Macdonald . The host feeding preference An. labranchiae was estimated by considering different feeding preference indices in addition to the HBI, i.e. the forage ratio (FR) and the feeding index (FI). FR quantifies vector selection of a particular vertebrate host rather than other available hosts. It was calculated by dividing the percentage of females fed on a given host by the percentage which that host represented in the total census of available animals and humans at the collecting site . FRs significantly > 1.0 indicate a selective bias and values < 1.0 indicate avoidance in favour of other hosts; FRs ≈ 1.0 show neither preference nor avoidance. FI is defined as the observed proportion of females fed on a certain animal host with respect to another one divided by the expected comparative proportion of feeds on these two hosts . This crude index was adjusted by taking in account factors that affect feeding, such as host abundance, their size and their temporal and spatial concurrence with the mosquito species. FI = 1 indicates equal feeding on the two hosts, while smaller or larger values indicate a decrease or increase in feeding on the first host relative to the second, respectively. FI were calculated for each pair of hosts.
The possible relationship between global female abundance and size of the fraction biting man during night catches was also evaluated by the Pearson's statistical test, comparing the data recorded in the same area of site 1 with those from site 3 over a period of 14 years (1995-2008 - our own unpublished data). A coefficient of endophagy of An. labranchiae (i.e. the ratio of the number of specimens caught biting indoors versus those caught biting outdoors) was also assessed by the analysis of the retrospective data from human bait catches, performed both outdoors and indoors dwellings in 1994-1996.
Evaluation of the presence of potential reservoirs of infection
The number of gametocyte carriers that may have been circulating in the territory during the period favourable to malaria transmission (June-October) was obtained by the analysis of the cases of imported malaria in Italy in 2000-2009 (cases confirmed by the Malaria Reference Centre at Istituto Superiore di Sanità), selecting those reported from hospitals located into the study area.
Results are grouped and showed by parameter adopted for assessing the malariogenic potential of the study area.
Out of a total of 8,006 females belonging to the An. maculipennis complex considered in this study, 1,772 (22.1%) were morphologically and molecularly identified at species level. Although at different levels of prevalence and abundance (Table 1), An. labranchiae occurred in all the study sites where it represents the dominant species of the maculipennis complex, with the exception of site 3 and 5, the rice fields of Val di Merse and the farm of Monte Antico, both located in an hilly area over 300 m a.s.l., where its prevalence accounted for 1-3% and 16% respectively, being predominant An. maculipennis s.s. because of the different climate conditions [23, 34]. The rice fields of the coastal plain (sites 1-2) remained the most productive areas of An. labranchiae (100-500 females/shelter), where it represents 96-98% of the species belonging to the complex. A high prevalence of An labranchiae (90-98%), but with lower levels of abundance (range 5-80 females/shelter) was recorded in the remaining study sites, where changes in land use occurred during the last three decades, have contributed to make the territory less favourable to the development of anopheline mosquitoes. A comparison of recent findings with those available for the previous decade showed a reduction of the abundance of An. maculipennis, s.l. resting females in the study area of about 75-80%.
Seasonality of Anopheles labranchiae
Differences between predicted and observed starting and ending week of development season of Anopheles labranchiae
Anopheles labranchiae adult population predictive map
Length of the potential transmission season
Results of the experimental infections of Anopheles maculipennis s.l. females from Principina (Grosseto, Italy)
Positive females/total dissected
Mean oocysts/female (int. limit)
25/96 (after 8 days)
9/130 (after 15 days)
In both 2008 and 2009, some infection attempts of F1 first batch of An. labranchiae were also carried out but, no infected mosquitoes were detected in these samples. On the whole, fewer than 5% of the field collected females took an infected blood meal, and only 20% of these mosquitoes survived 15 days post infection. In all infection experiments, the An. gambiae became infected with oocyst prevalence ranging from 65 to 100%.
Vectorial capacity (VC) of Anopheles labranchiae in the rice fields of Principina (Grosseto, Italy), 2005-2006
Potential transmission season
Average Temperature (°C)
Average Relative Humidity (%)
Length of gonotrophic cycle (n. days)
P. falciparum sporogonic cycle (n) Pf
P. vivax sporogonic cycle (n) Pv
Human Blood index (HBI)
Human biting rate (ma)
Parity rate % (P)
Vector daily survival probability (p)
Expectancy of infective life (n. days) (p n ) Pf
Expectancy of infective life (n. days) (p n ) Pv
Expectancy of infective life (1/-log e p)
Longevity factor (p n /-log e p) Pf
Longevity factor (p n /-log e p) Pv
VC P. falciparum
VC P. vivax
Feeding preference and anthropophily of An. labranchiae
Blood meal sources and forage ratio for Anopheles labranchiae in the rice-fields of Principina (Grosseto)
Mean weight (kg)
No. of fed An. Labranchia e females
Forage ratio (FR)
Gametocyte carrier introduction
Number of gametocyte carriers circulating in June-October in the study area (2000-2009)
No. of gametocytes carriers/out of national cases
Gametocyte carriers circulating in the study area
1 W. Africa
The presence of potential vectors, the progressive climate increase and the possible introduction of parasite reservoirs raises the concern about the possibility of malaria re-emerging in Italy. The results of previous and present studies and the analysis of historic data showed a continuing receptivity in the Maremma, due to the presence of An. labranchiae at all selected sites, despite a marked reduction of the abundance of the vector with respect to the previous 3 decades. Rice fields (site 1, in particular) remained the most productive areas for An. maculipennis s.l. However, while in sites1-2 in the coastal plain An. labranchiae represents 96-98% of the species belonging to the complex, in the hilly area of site 3 the prevalence of An. labranchiae is only 2%. It is noteworthy that in site 3 (where the dominant species is An. maculipennis s.s.) An. labranchiae historically absent in that area, was first recorded in 2005, possibly indicating a north-eastern expansion of the range of this potential vector.
The vectorial capacity values assessed in site 1 were very low for both P. falciparum and P. vivax, because a very low size of the HBI (Table 4), in all cases below 0.5, commonly considered as the threshold that characterize a situation of instability or even below 0.02, that should represents the threshold below which the malaria transmission may be interrupted . Nevertheless, it should be considered that the competence of a malaria vector may be strongly affected by environmental factors (temperature and land cover) and by some other traits, related to the trophic activity, such as host feeding preference, which have genetic components [55, 56]. There are two critical points in determining vectorial capacity: "ma", that represents human exposure to mosquito bites, may lead to strongly overestimated values , and the HBI, that may underestimate vector-human contacts because collections of resting females in the human dwellings are not considered in the current protocols of the entomological surveys in Europe. For these reasons even a small change in accessibility to humans will have a marked impact on the VC values.
Moreover, it should be noted that the daily survival rate pn for the sporogonic development of the parasite in the vector, a factor that strongly affects the Macdonald formula, is a function of temperature.
Hence the rise in temperature appears to represent the most important factor that may influence the receptivity of Maremma. The climatic analysis (Table 2 and Figure 3) outlines a sharp increase of the mean temperature in the study period. From 2005 to 2009 an increase of 0.8°C and 1.2°C respectively for the yearly and the seasonal (May/August) mean temperature or even an increase of 1.2°C and 1.6°C respectively for the yearly and the seasonal maximum temperature was assessed (Table 2). These results are confirmed by recent climatic studies [57–62] which highlight a marked warming and an increase in extreme temperature events in Tuscany, and more generally in Italy, as well as a warming trend elsewhere in the Mediterranean area [63–67]. The potential transmission risk analysis for P. vivax and P. falciparum evaluated through the GMR index calculation showed in 2005-2009 a larger favourable transmission period during the year than the climatic reference period 1961-1990 (Figures 4, 5 and 6). Climate change scenarios are typically oriented towards higher temperatures but there is greater uncertainty about climate influences on rainfall . However, these uncertainties are irrelevant to the malaria-risk evaluation in the sites of Maremma were the most productive breeding sites are rice-fields that are independent of rainfall. GMR index results obtained with constant R/PET ratio suggest a need for vigilance in the future through surveillance and monitoring activities.
Although the malariogenic potential of Maremma seemed to be very low, it is worth considering the following points:
i) The study carried out in site 1 showed a high abundance of An. labranchiae, and even if the species may appear to be opportunistic in its behaviour, past and present data confirm its ability to bite humans in presence or absence of alternative hosts, indoors as well as outdoors. Moreover, despite the very difficult access to humans, and the very low VC values, it is quite remarkable that to have been found also some An. labranchiae gorged on human indicate a little (but not zero) risk of contact human-vector in the area.
ii) Tourism development and changes in land use have resulted often unfavourable to the development of mosquitoes, but in some cases, such as the extension of resorts and holiday farms close to cultivated fields, may be promoting the promiscuity between mosquitoes, the increased availability of non immune humans and of gametocyte carriers.
iii) The very low vulnerability of the study area, due to the scarce presence of gametocyte carriers circulating in the area during the favourable months of the summer may be increased by an unknown number of non regular immigrants entering Italy seasonally as farm labourers (most of them coming from French-speaking West African countries, where malaria is highly endemic).
iv) The general rise in average temperature during the late spring and summer could favour the parasite development, shortening the extrinsic cycle of Plasmodium spp, as well as the gonothrophic cycle of the vector and increasing the length of the transmission season.
In conclusion, Maremma, as well as other "at risk" areas recently investigated in the Mediterranean coastal countries, is excluded for the eventual return to a situation of endemic malaria [21, 25, 29, 69–71], while the occurrence of sporadic, isolated cases of introduced P. vivax malaria may be considered possible.
We thank for technical assistance M.T. Lecoq for producing Anopheles gambiae control strain used for Plasmodium falciparum infection, A. Massa and F. Mancini Barbieri for help in field and lab work. This research was supported by the Istituto Superiore di Sanità and funded under the EU 6th Framework Program (GOCE-CT-2003-010284 EDEN). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission.
- Saker L, Lee K, Cannito B, Gilmore A, Campbell-Lendrum D: Globalization and infectious diseases: a review of the linkages. 2004, Geneva: World Health Organization; TDR/STR/SEB/ST/04.2, [Series: Special topics in social, economic, and behavioural (SEB) research Special Topics No. 3]Google Scholar
- Sutherst RW: Global change and human vulnerability to vector borne diseases. Clin Microbiol Rev. 2004, 17: 136-173. 10.1128/CMR.17.1.136-173.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Gould EA, Higgs S: Impact of climate change and other factors on emerging arbovirus diseases. Trans R Soc Trop Med Hyg. 2009, 103: 109-121. 10.1016/j.trstmh.2008.07.025.PubMed CentralView ArticlePubMedGoogle Scholar
- Vora N: Impact of anthropogenic environmental alterations on vector-borne diseases. Medscape J Med. 2008, 10: 238-244.PubMed CentralPubMedGoogle Scholar
- Harrus S, Baneth G: Drivers for the emergence and re-emergence of vector-borne protozoal and bacterial diseases. Int J Parasitol. 2005, 35: 1309-1318. 10.1016/j.ijpara.2005.06.005.View ArticlePubMedGoogle Scholar
- Hubálek Z: Mosquito-borne viruses in Europe. Parasitol Res. 2008, 103 (Suppl 1): S29-S43.View ArticlePubMedGoogle Scholar
- Takken W, Knols BGJ: Emerging pests and vector-borne diseases in Europe. Ecology and control of vector-borne diseases. Volume 1. 2007, Wageningen: Wageningen Academic PublishersView ArticleGoogle Scholar
- WHO Regional Office for Europe: CISID malaria. [http://data.euro.who.int/cisid]
- Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 8: 1763-1766.Google Scholar
- Kuhn KG, Campbell-Lendrum DH, Davies CR: A continental risk map for malaria mosquito (Diptera: Culicidae) vectors in Europe. J Med Entomol. 2002, 39: 621-630. 10.1603/0022-2585-39.4.621.View ArticlePubMedGoogle Scholar
- Zahar AR: Vector bionomics in the epidemiology and control of malaria. Part II: The WHO European region and the WHO Eastern Mediterranean region. Volume II: Applied field studies. Section III: Vector bionomics, malaria epidemiology and control by geographical areas (a) the Mediterranean basin. 1990, Geneva: World Health Organization, VBC/90.2-MAL/90.2Google Scholar
- Romi R: Anopheles labranchiae, an important vector in Italy, and other potential malaria vectors in Southern Europe. Eur Mosq Bull. 1999, 4: 8-10.Google Scholar
- Snow K, Ramsdale C: Distribution chart for European Mosquitoes. Eur Mosq Bull. 1999, 3: 14-31.Google Scholar
- Baldari M, Tamburro A, Sabatinelli G, Romi R, Severini C, Cuccagna P, Fiorilli G, Allegri MP, Buriani C, Toti M: Malaria in Maremma, Italy. Lancet. 1998, 351: 1246-1247. 10.1016/S0140-6736(97)10312-9.View ArticlePubMedGoogle Scholar
- Kampen H, Proft J, Etti S, Maltezos E, Pagonaki M, Maier WA, Seitz HM: Individual cases of autochthonous malaria in Evros Province, northern Greece: entomological aspects. Parasitol Res. 2003, 89: 252-258.PubMedGoogle Scholar
- Cuadros J, Calvente MJ, Benito A, Arévalo J, Calero MA, Segura J, Rubio JM: Plasmodium ovale malaria acquired in central Spain. Emerg Infect Dis. 2002, 8: 1506-1508.PubMed CentralView ArticlePubMedGoogle Scholar
- Santa-Olalla-Peralta P, Vazquez-Torres MC, Latorre-Fandós E, Mairal-Claver P, Cortina-Solano P, Puy-Azón A, Diego-Sancho B, Leitmeyer K, Lucientes-Curdi J, Sierra-Moros MJ: First autochthonous malaria case due to Plasmodium vivax since eradication, Spain, October 2010. Euro Surveill. 2010, 15: pii = 19684Google Scholar
- Danis K, Baka A, Lenglet A, Van-Bortel W, Terzaki I, Tseroni M, Detsis M, Papanikolaou E, Balaska A, Gewehr S, Dougas G, Sideroglou T, Economopoulou A, Vakalis Tsiodras N, Bonovas N, Kremastinou J: Autochthonous Plasmodium vivax malaria in Greece, 2011. Euro Surveill. 2011, 16: pii = 19993Google Scholar
- Armengaud A, Legros F, D'Ortenzio E, Quatresous I, Barre H, Houze S, Valayer P, Fanton Y, Schaffner F: A case of autochthonous Plasmodium vivax malaria, Corsica, August 2006. Travel Med Infect Dis. 2008, 6: 36-40. 10.1016/j.tmaid.2007.09.042.View ArticlePubMedGoogle Scholar
- Doudier B, Bogreau H, DeVries A, Ponçon N, Stauffer WM, Fontenille D, Rogier C, Parola P: Possible autochthonous malaria from Marseille to Minneapolis. Emerg Infect Dis. 2007, 13: 1236-1638. 10.3201/eid1308.070143.PubMed CentralView ArticlePubMedGoogle Scholar
- Ponçon N, Tran A, Toty C, Luty AJ, Fontenille D: A quantitative risk assessment approach for mosquito-borne diseases: malaria re-emergence in southern France. Malar J. 2008, 7: 147-159. 10.1186/1475-2875-7-147.PubMed CentralView ArticlePubMedGoogle Scholar
- Faraj C, Ouahabi S, Adlaoui E, Boccolini D, Romi R, El Aouad R: Assessment of malaria resurgence in Morocco. Study of the vectorial capacity of Anopheles labranchiae in a rice cultivation area in the north of the country. Parasite. 2008, 15: 605-610.View ArticlePubMedGoogle Scholar
- Di Luca M, Boccolini D, Severini F, Toma L, Mancini Barbieri F, Massa A, Romi R: A 2-year entomological study of potential malaria vectors in Central Italy. Vector Borne Zoonotic Dis. 2009, 9: 703-711. 10.1089/vbz.2008.0129.View ArticlePubMedGoogle Scholar
- Vicente JL, Sousa CA, Alten B, Caglar S, Falcuta E, Latorre JM, Toty C, Barrè H, Demirci B, Di Luca M, Toma L, Alves R, Salgueiro P, Silva TL, Bargues MA, Mas-Coma S, Boccolini D, Romi R, Nicolescu G, do Rosario VE, Ozer N, Fontenille D, Pinto J: Genetic and phenotypic variation of the malaria vector Anopheles atroparvus in Southern Europe. Malar J. 2011, 10: 5-13. 10.1186/1475-2875-10-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Sainz-Elipe S, Latorre JM, Escosa R, Montserrat M, Fuentes MV, Mas-Coma S, Bargues MD: Malaria resurgence risk in Southern Europe: climate assessment in an historically endmic area of rice fields at the Mediterranean shore of Spain. Malar J. 2010, 9: 221-236. 10.1186/1475-2875-9-221.PubMed CentralView ArticlePubMedGoogle Scholar
- Buonomini G, Gori D: Anofelismo e malaria in alcune zone della provincia di Siena. Riv Malariol. 1937, 16: 142-159.Google Scholar
- Bettini S: L'anofelismo e la malaria in Maremma. Quaderni ARA 97/1. 1997, Orbetello (Grosseto): Associazione Ricerche AmbienteGoogle Scholar
- Bueti S, Corti M: La malaria nella Maremma grossetana: una malattia dagli aspetti sociali. Med Secoli. 1998, 10: 541-555.PubMedGoogle Scholar
- Romi R, Pierdominici G, Severini C, Tamburro A, Cocchi M, Menichetti D, Pili E, Marchi A: Status of malaria vectors in Italy. J Med Entomol. 1997, 34: 263-271.View ArticlePubMedGoogle Scholar
- Romi R, Sabatinelli G, Majori G: Could Malaria Reappear in Italy?. Emerg Infect Dis. 2001, 7: 915-919. 10.3201/eid0706.010601.PubMed CentralView ArticlePubMedGoogle Scholar
- Coluzzi A, Finizio E: Il problema degli insetti ematofagi nella zona litoranea della provincia di Grosseto. Riv Malariol. 1966, 45: 137-145.PubMedGoogle Scholar
- Bettini S, Gradoni L, Cocchi M, Tamburro A: Rice culture and Anopheles labranchiae in Central Italy. WHO/Mal. 1978, 78: 1-6.Google Scholar
- Majori G, Maroli M, Bettini S, Pierdominici G: Osservazioni sull'anofelismo residuo nel Grossetano. Riv Parassitol. 1970, 31: 147-154.PubMedGoogle Scholar
- Boccolini D, Toma L, Di Luca M, Severini F, Cocchi M, Bella A, Massa A, Mancini Barbieri F, Bongiorno G, Pontuale G, Raffaelli I, Fausto AM, Tamburro A, Romi R: Impact of environmental changes and human-related factors on the potential malaria vector, Anopheles labranchiae (Diptera: Culicidae), in Maremma, Central Italy. J Med Entomology.Google Scholar
- Akaike H: Information theory and an extension of the maximum likelihood principle. Proceedings of the 2nd International Symposium on Information Theory: 1973. Edited by: Petrov BN, Csaki F. 1973, Budapest: Akadémiai Kiadó, 267-281.Google Scholar
- Yee TW, Hastie TJ: Reduced-rank vector generalized linear models. Stat Modelling. 2003, 3: 15-41. 10.1191/1471082X03st045oa.View ArticleGoogle Scholar
- Venables WN, Ripley BD: Modern Applied Statistics Fourth edition. 2002, New York: Springer Science + Business MediaGoogle Scholar
- Malone JB: Biology-based mapping of vector-borne parasites by Geographic Information Systems and Remote Sensing. Parassitologia. 2005, 47: 27-50.PubMedGoogle Scholar
- De Gaetano AT: Meteorological effects on adult mosquito (Culex) populations in metropolitan New Jersey. Int J Biometeorol. 2005, 49: 345-353. 10.1007/s00484-004-0242-2.View ArticleGoogle Scholar
- Atchley FO, Traylor WR, Weathersbee AA: Effects of variations in reservoir levels, rainfall, and temperature on anopheline densities in a coastal plains area in South Carolina. J Parasitol. 1995, 41: 273-280.View ArticleGoogle Scholar
- Ribbands CR: The influence of rainfall, tides and periodic fluctuations on a population of Anopheles melas. Bull Entomol Res. 1944, 35: 271-295.View ArticleGoogle Scholar
- WWRP/WGNE Joint Working Group on Forecast Verification Research. [http://www.cawcr.gov.au/projects/verification]
- Tran A, Ponçon N, Toty C, Linard C, Guis H, Ferré JB, Seen DL, Roger F, de la Rocque S, Fontenille D, Baldet T: Using remote sensing to map larval and adult populations of Anophels hyrcanus (Diptera: Culicidae) a potential malaria vector in Southern France. Int J Health Geog. 2008, 7: 9-20. 10.1186/1476-072X-7-9.View ArticleGoogle Scholar
- Regione Toscana. Agricoltura e Foreste. Annuario Statistico, 1998-2008. [http://www.regione.toscana.it/agricoltura.html]
- Bruce-Chwatt LJ: Essential malariology. 1985, London: William Heinemann Medical Books, SecondGoogle Scholar
- Thornthwaite CW: An approach toward a rational classification of climate. Geogr Rev. 1948, 38: 55-94. 10.2307/210739.View ArticleGoogle Scholar
- Mitri C, Thiery I, Bourgouin C, Paul R: Density-dependent impact of the human malaria parasite Plasmodium falciparum gametocyte sex ratio on mosquito infection rates. Proc Biol Sci. 2009, 276: 3721-3726. 10.1098/rspb.2009.0962.PubMed CentralView ArticlePubMedGoogle Scholar
- Padley D, Moody AH, Chiodini PL, Saldanha J: Use of a rapid, single-round, multiplex PCR to detect malarial parasites and identify the species present. Ann Trop Med Parasitol. 2003, 29: 131-137.View ArticleGoogle Scholar
- Proft J, Majer WA, Jampen H: Identification of six sibling species of the Anopheles maculipennis complex (Diptera: Culicidae) by a polymerase chain reaction assay. Parasitol Res. 1999, 85: 637-643.View ArticleGoogle Scholar
- MacDonald G: The epidemiology and control of malaria. 1957, London: Oxford University PressGoogle Scholar
- Garret-Jones C: The Human Blood Index of Malaria Vectors in Relation to Epidemiological Assessment. Bull World Health Organ. 1964, 30: 241-261.Google Scholar
- Davidson G: Estimation of the survival rate of anopheline mosquitoes in nature. Nature. 1954, 174: 792-793. 10.1038/174792a0.View ArticlePubMedGoogle Scholar
- Hess AD, Hayes RO, Tempelis CH: Use of forage ratio technique in mosquito host preference studies. Mosq News. 1968, 28: 386-387.Google Scholar
- Kay BH, Boreham PFL, Williams GM: Host preferences and feeding patterns of mosquitoes at Kowanyama, Cape York Peninsula, Northern Quensland. Bull Entomol Res. 1979, 69: 411-417.Google Scholar
- Beerntsen BT, James AA, Christensen BM: Genetic of mosquito vector competence. Microbiol Mol Biol Rev. 2000, 64: 115-137. 10.1128/MMBR.64.1.115-137.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Chaves LF, Harington LC, Keogh C, Nguyen AM, Kitron UD: Blood feeding patterns of mosquitoes: random or structured?. Front Zool. 2010, 7: 1-16. 10.1186/1742-9994-7-1.View ArticleGoogle Scholar
- Menne B, Ebi KL: Climate change and adaptation strategies for human health. 2006, Darmstadt: Steinkopff-VerlagGoogle Scholar
- Bartolini G, Morabito M, Crisci A, Grifoni D, Torrigiani T, Petralli M, Maracchi G, Orlandini S: Recent trends in Tuscany (Italy) summer temperature and indices of extremes. Int J Climatol. 2008, 26: 1751-1760.View ArticleGoogle Scholar
- Bartolini G, Di Stefano V, Maracchi G, Orlandini S: Recent Mediterranean warming is especially due to summer season, Evidences from Tuscany (central Italy). Theor Appl Climatol. 2012, 107: 279-295. 10.1007/s00704-011-0481-1.View ArticleGoogle Scholar
- Brunetti M, Maugeri M, Monti F, Nanni T: Temperature and precipitation variability in Italy in the last two centuries from homogenised instrumental time series. Int J Climatol. 2006, 26: 345-381. 10.1002/joc.1251.View ArticleGoogle Scholar
- Toreti A, Desiato F: Changes in temperature extremes over Italy in the last 44 years. Int J Climatol. 2007, 28: 733-745.View ArticleGoogle Scholar
- Toreti A, Desiato F: Temperature trend over Italy from 1961 to 2004. Theor Appl Climatol. 2006, 91: 51-58.View ArticleGoogle Scholar
- Giannakopoulos C, Le Sager P, Bindi M, Moriondo M, Kostopoulou E, Goodess CM: Climatic changes and associated impacts in the Mediterranean resulting from a 2°C global warming. Glob Planet Change. 2009, 68: 209-224. 10.1016/j.gloplacha.2009.06.001.View ArticleGoogle Scholar
- Giorgi F, Lionello P: Climate change projections for the Mediterranean region. Glob Planet Change. 2008, 63: 90-104. 10.1016/j.gloplacha.2007.09.005.View ArticleGoogle Scholar
- Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL: Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. 2007, Cambridge: Cambridge University PressGoogle Scholar
- Schär C, Vidale PL, Lüthi D, Frei C, Häberli C, Liniger MA, Appenzeller C: The role of increasing temperature variability in European summer heatwaves. Nature. 2004, 443: 332-336.View ArticleGoogle Scholar
- Somot S, Sevault F, Deque M, Crepon M: 1st century climate change scenario for the Mediterranean using a coupled atmosphere-ocean regional climate model. Glob Planet Change. 2008, 63: 112-126. 10.1016/j.gloplacha.2007.10.003.View ArticleGoogle Scholar
- Kyselý J, Beranova R: Climate-change effects on extreme precipitation in central Europe: uncertainties of scenarios based on regional climate models. Theor Appl Climatol. 2009, 95: 361-374. 10.1007/s00704-008-0014-8.View ArticleGoogle Scholar
- Capinha C, Gomes E, Reis E, Rocha J, Sousa CA, do Rosario VE, Ameida AP: Present habitat suitability for anopheles atroparvus (Diptera:Culicidae) and its coincidence with former malaria vector areas in main Portugal. Geospat Health. 2009, 3: 177-187.View ArticlePubMedGoogle Scholar
- Linard C, Ponçon N, Fontenille D, Lambin EF: Risk of malaria reemergence in southern France: testing scenarios with a multiagent simulation model. EcoHealth. 2009, 6: 135-147. 10.1007/s10393-009-0236-y.View ArticlePubMedGoogle Scholar
- Toty C, Barré H, Le Goff G, Larget-Thiéry I, Rahola N, Couret D, Fontenille D: Malaria risk in Corsica, former hot spot of malaria in France. Malaria J. 2010, 9: 231-238. 10.1186/1475-2875-9-231.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.