Search

Culex pipiens mosquito is a species widely spread across Europe and represents a competent vector for many arboviruses such as West Nile virus (WNV), which has been recently circulating in many European countries, causing hundreds of human cases. In order to identify the main determinants of the high heterogeneity in Cx. pipiens abundance observed in Piedmont region (Northwestern Italy) among different seasons, we developed a density-dependent stochastic model that takes explicitly into account the role played by temperature, which affects both developmental and mortality rates of different life stages. The model was calibrated with a Markov chain Monte Carlo approach exploring the likelihood of recorded capture data gathered in the study area from 2000 to 2011; in this way, we disentangled the role played by different seasonal eco-climatic factors in shaping the vector abundance. Illustrative simulations have been performed to forecast likely changes if temperature or density-dependent inputs would change. Our analysis suggests that inter-seasonal differences in the mosquito dynamics are largely driven by different temporal patterns of temperature and seasonal-specific larval carrying capacities. Specifically, high temperatures during early spring hasten the onset of the breeding season and increase population abundance in that period, while, high temperatures during the summer can decrease population size by increasing adult mortality. Higher densities of adult mosquitoes are associated with higher larval carrying capacities, which are positively correlated with spring precipitations. Finally, an increase in larval carrying capacity is expected to proportionally increase adult mosquito abundance.

The spread of tick-borne pathogens represents an important threat to human and animal health in many parts of Eurasia. Here, we analysed a 9-year time series of Ixodes ricinus ticks feeding on Apodemus flavicollis mice (main reservoir-competent host for tick-borne encephalitis, TBE) sampled in Trentino (Northern Italy). The tail of the distribution of the number of ticks per host was fitted by three theoretical distributions: Negative Binomial (NB), Poisson-LogNormal (PoiLN), and Power-Law (PL). The fit with theoretical distributions indicated that the tail of the tick infestation pattern on mice is better described by the PL distribution. Moreover, we found that the tail of the distribution significantly changes with seasonal variations in host abundance. In order to investigate the effect of different tails of tick distribution on the invasion of a non-systemically transmitted pathogen, we simulated the transmission of a TBE-like virus between susceptible and infective ticks using a stochastic model. Model simulations indicated different outcomes of disease spreading when considering different distribution laws of ticks among hosts. Specifically, we found that the epidemic threshold and the prevalence equilibria obtained in epidemiological simulations with PL distribution are a good approximation of those observed in simulations feed by the empirical distribution. Moreover, we also found that the epidemic threshold for disease invasion was lower when considering the seasonal variation of tick aggregation.

The structure of a network dramatically affects the spreading phenomena unfolding upon it. The contact distribution of the nodes has long been recognized as the key ingredient in influencing the outbreak events. However, limited knowledge is currently available on the role of the weight of the edges on the persistence of a pathogen. At the same time, recent works showed a strong influence of temporal network dynamics on disease spreading. In this work we provide an analytical understanding, corroborated by numerical simulations, about the conditions for infected stable state in weighted networks. In particular, we reveal the role of heterogeneity of edge weights and of the dynamic assignment of weights on the ties in the network in driving the spread of the epidemic. In this context we show that when weights are dynamically assigned to ties in the network, a heterogeneous distribution is able to hamper the diffusion of the disease, contrary to what happens when weights are fixed in time.

Background: Flavivirus is a highly heterogeneous viral genus that includes important human pathogens and several viral strains with unknown zoonotic potential. Mosquito-borne flaviviruses have been isolated and characterized in Northern Italy: West Nile virus and Usutu virus were detected in mosquitoes and in different host species and recent studies provided evidence about the circulation of "insect Flavivirus" strains., Methods: In order to clarify the diffusion and the distribution of the mosquito-transmitted flaviviruses in Italy, we analyzed Culex and Ochlerotatus mosquitoes collected in 2009 and 2010 in an area divided evenly between hills and plains and where the landscape is dominated by mixed agricultural patches, rice fields, deciduous tree forests, and urban environments. Each mosquito pool was tested for the presence of Flavivirus strains and we characterized positive samples by genetic sequencing., Results: Positive mosquito pools revealed low infection prevalence, but suggested a continuous circulation of both Usutu virus and insect Flavivirus. Interestingly, phylogenetic analyses based on NS5 gene partial sequences showed a closer relationship among new Usutu virus strains from Piedmont and the reference sequences from the Eastern Europe, with respect to Italian samples characterized so far. Moreover, NS5 gene phylogeny suggested that mosquito flaviviruses found in Italy could belong to different lineages., Conclusions: Our results contribute to a wider point of view on the heterogeneity of viruses infecting mosquitoes suggesting a taxonomical revision of the Mosquito-borne Flavivirus group.

Background: West Nile Virus (WNV) transmission in Italy was first reported in 1998 as an equine outbreak near the swamps of Padule di Fucecchio, Tuscany. No other cases were identified during the following decade until 2008, when horse and human outbreaks were reported in Emilia Romagna, North Italy. Since then, WNV outbreaks have occurred annually, spreading from their initial northern foci throughout the country. Following the outbreak in 1998 the Italian public health authority defined a surveillance plan to detect WNV circulation in birds, horses and mosquitoes. By applying spatial statistical analysis (spatial point pattern analysis) and models (Bayesian GLMM models) to a longitudinal dataset on the abundance of the three putative WNV vectors [Ochlerotatus caspius (Pallas 1771), Culex pipiens (Linnaeus 1758) and Culex modestus (Ficalbi 1890)] in eastern Piedmont, we quantified their abundance and distribution in space and time and generated prediction maps outlining the areas with the highest vector productivity and potential for WNV introduction and amplification., Results: The highest abundance and significant spatial clusters of Oc. caspius and Cx. modestus were in proximity to rice fields, and for Cx. pipiens, in proximity to highly populated urban areas. The GLMM model showed the importance of weather conditions and environmental factors in predicting mosquito abundance. Distance from the preferential breeding sites and elevation were negatively associated with the number of collected mosquitoes. The Normalized Difference Vegetation Index (NDVI) was positively correlated with mosquito abundance in rice fields (Oc. caspius and Cx. modestus). Based on the best models, we developed prediction maps for the year 2010 outlining the areas where high abundance of vectors could favour the introduction and amplification of WNV., Conclusions: Our findings provide useful information for surveillance activities aiming to identify locations where the potential for WNV introduction and local transmission are highest. Such information can be used by vector control offices to stratify control interventions in areas prone to the invasion of WNV and other mosquito-transmitted pathogens.

To estimate the basic reproduction number (R(0)) of Borrelia lusitaniae and Borrelia afzelii, we formulated a mathematical model considering the interactions among the tick vector, vertebrate hosts, and pathogens in a 500-ha enclosed natural reserve on Le Cerbaie hills, Tuscany, central Italy. In the study area, Ixodes ricinus were abundant and were found infected by B. lusitaniae and B. afzelii. Lizards (Podarcis spp.) and mice (Apodemus spp.), respectively, are the reservoir hosts of these two Borrelia burgdorferi sensu lato (s.l.) genospecies and compete for immature ticks. B. lusitaniae R(0) estimation is in agreement with field observations, indicating the maintenance and diffusion of this genospecies in the study area, where lizards are abundant and highly infested by I. ricinus immature stages. In fact, B. lusitaniae shows a focal distribution in areas where the tick vector and the vertebrate reservoir coexist. Mouse population dynamics and their relatively low suitability as hosts for nymphs seem to determine, on the other hand, a less efficient transmission of B. afzelii, whose R(0) differs between scenarios in the study area. Considering host population dynamics, the proposed model suggests that, given a certain combination of the two host population sizes, both spirochete genospecies can coexist in our study area. Additional incompetent hosts for B. burgdorferi s.l. have a negative effect on B. afzelii maintenance, whose R(0) results > 1 only with high mouse population densities and/or low lizards abundance, but they do not seem to influence B. lusitaniae transmission cycle on Le Cerbaie. Secondly, our model confirms the importance of nymphs' infestation, of host population density and diversity, and spirochetes host association for the maintenance of the transmission cycle of B. burgdorferi s.l.

We explore the Hawk-Dove game on networks with topologies ranging from regular lattices to random graphs with small-world networks in between. This is done by means of computer simulations using several update rules for the population evolutionary dynamics. We find the overall result that cooperation is sometimes inhibited and sometimes enhanced in those network structures, with respect to the mixing population case. The differences are due to different update rules and depend on the gain-to-cost ratio. We analyze and qualitatively explain this behavior by using local topological arguments.