Lyme borreliosis in Belgium: a cost-of-illness analysis.

BACKGROUND: Lyme borreliosis (LB) is the most common tick-borne disease in Europe and North America, yet its economic burden remains largely unknown. This study aimed to estimate the economic cost associated with the different clinical manifestations of LB in Belgium. METHODS: An incidence approach and societal perspective were used to estimate the total cost-of-illness for LB in Belgium. Costs were calculated for patients with erythema migrans (EM) or disseminated/late LB, including patients who developed post-treatment Lyme disease syndrome (PTLDS). Direct medical, direct non-medical (transportation & paid help) and indirect non-medical costs (productivity losses) were included in the analysis. Ambulatory cost data were collected through a prospective cohort study from June 2016 to March 2020, in which patients with LB were followed up 6 to 12 months after diagnosis. Hospitalization costs were retrieved from the Minimal Clinical Data registry, a mandatory registry for all Belgian hospitals, linked to the Minimal Financial Data registry. Costs were expressed in 2019 euros. RESULTS: The total annual cost associated with clinical manifestations of LB in Belgium was estimated at €5.59 million (95% UI 3.82-7.98). Of these, €3.44 million (95% UI 2.05-5.48) or 62% was related to disseminated/late LB diagnoses and €2.15 million (95% UI 1.30-3.26) to EM. In general, direct medical costs and productivity losses accounted for 49.8% and 46.4% of the total costs, respectively, while direct non-medical costs accounted for only 3.8%. The estimated mean costs were €193 per EM patient and €5,148 per disseminated/late LB patient. While patients with PTLDS seemed to have somewhat higher costs compared to patients without PTLDS, the number of patients was too small to have representative estimates. CONCLUSIONS: We estimate the total annual direct medical costs, direct non-medical and indirect non-medical costs associated with LB to exceed €5.5 million per year, almost evenly distributed between EM (40%) and disseminated/late LB (60%). EM costs 26 times less per patient but occurs also 16 times more frequently than disseminated/late LB. The cost burden remains limited by comparison to other infectious diseases due to the relative lower incidence.

Non-specific symptoms and post-treatment Lyme disease syndrome in patients with Lyme borreliosis: a prospective cohort study in Belgium (2016-2020).

BACKGROUND: Patients with Lyme borreliosis (LB) may report persisting non-specific symptoms such as fatigue, widespread musculoskeletal pain or cognitive difficulties. When present for more than 6 months and causing a reduction in daily activities, this is often referred to as post-treatment Lyme disease syndrome (PTLDS). This study aimed to compare the occurrence of symptoms between LB patients and controls, to estimate the proportion of LB patients developing PTLDS and to identify risk factors. METHODS: A prospective cohort study was set up including three subpopulations: patients with an erythema migrans (EM) (i) or disseminated/late LB (ii) and a non-LB control group (iii). At 6- and 12-months follow-up, the occurrence of several symptoms, including six symptoms used to define PTLDS, i.e. muscle pain, joint pain, fatigue, memory problems, difficulties concentrating and problems finding words, and impact on daily activities, was compared between LB patients and controls. Finally, the proportion of LB patients developing PTLDS as defined by the Infectious Disease Society of America was estimated, including a time frame for symptoms to be present. RESULTS: Although the risk of presenting PTLDS-related symptoms was significantly higher in EM patients (n = 120) compared to controls (n = 128) at 6 months follow-up, the risk of presenting at least one of these symptoms combined with impact on daily activities was not significantly higher in EM patients, at either 6- or 12-months follow-up. A significant association was found between disseminated/late LB (n = 15) and the occurrence of any PTLDS-symptom with an impact on daily activities at both time points. The proportion of patients with PTLDS was estimated at 5.9% (95% CI 2.7-12.9) in EM patients and 20.9% (95% CI 6.8-64.4) in patients with disseminated/late LB (RR = 3.53, 95% CI 0.98-12.68, p = 0.053). No significant risk factors were identified, which may be explained by small sample sizes. CONCLUSIONS: In our study, PTLDS was present in both LB cohorts, yet with a higher percentage in disseminated/late LB patients. Additional research is needed into risk factors for and causes of this syndrome. In addition, development and validation of standardized methods to assess the PTLDS case definition, easily applicable in practice, is of great importance.

No molecular detection of tick-borne pathogens in the blood of patients with erythema migrans in Belgium.

BACKGROUND: A number of tick-borne pathogens circulate in the Belgian tick population in addition to the causative agent of Lyme borreliosis. However, so far, only a few patients with tick-borne diseases other than Lyme borreliosis have been reported in Belgium. The aim of this study was to investigate the occurrence of other human tick-borne infections in Belgium and their possible clinical manifestation. METHODS: Patients with fever (> 37.5 °C) after a tick bite or those with erythema migrans (EM) were included in the study. EDTA-blood samples were screened for the presence of DNA from Borrelia burgdorferi sensu lato, Borrelia miyamotoi, Anaplasma phagocytophilum, Neoehrlichia mikurensis, spotted fever group rickettsiae (genus Rickettsia), Babesia spp., Bartonella spp., Spiroplasma ixodetis and tick-borne encephalitis virus, using multiplex PCR methods. A questionnaire on, among others, demographics and clinical symptoms, was also filled in. RESULTS: Over a period of 3 years, 119 patients with EM and 14 patients with fever after a recent tick bite were enrolled in the study. Three samples initially tested positive for N. mikurensis by quantitative PCR (qPCR), but the results could not be confirmed by other PCR methods, and repetition of the DNA extraction procedure and qPCR test was not successful. The qPCR test results for the other tick-borne pathogens were negative. CONCLUSIONS: In general, only a few patients with fever after a tick bite could be identified. Although no tick-borne pathogens were detected, their occurrence cannot be excluded based on the limited number of patients and the limitations inherent to current methodologies. This study underscores the possibility of false-positive PCR results and the necessity for the development of multiple independent tools for the sensitive and specific detection of emerging tick-borne pathogens.

Why Hantavirus Prevalence Does Not Always Increase With Host Density: Modeling the Role of Host Spatial Behavior and Maternal Antibodies.

For wildlife diseases, one often relies on host density to predict host infection prevalence and the subsequent force of infection to humans in the case of zoonoses. Indeed, if transmission is mainly indirect, i.e., by way of the environment, the force of infection is expected to increase with host density, yet the laborious field data supporting this theoretical claim are often absent. Hantaviruses are among those zoonoses that have been studied extensively over the past decades, as they pose a significant threat to humans. In Europe, the most widespread hantavirus is the Puumala virus (PUUV), which is carried by the bank vole and causes nephropathia epidemica (NE) in humans. Extensive field campaigns have been carried out in Central Finland to shed light on this supposed relationship between bank vole density and PUUV prevalence and to identify other drivers for the infection dynamics. This resulted in the surprising observation that the relationship between bank vole density and PUUV prevalence is not purely monotonic on an annual basis, contrary to what previous models predicted: a higher vole density does not necessary result in a higher infection prevalence, nor in an increased number of humans reported having NE. Here, we advance a novel individual-based spatially-explicit model which takes into account the immunity provided by maternal antibodies and which simulates the spatial behavior of the host, both possible causes for this discrepancy that were not accounted for in previous models. We show that the reduced prevalence in peak years can be attributed to transient immunity, and that the density-dependent spatial vole behavior, i.e., the fact that home ranges are smaller in high density years, plays only a minor role. The applicability of the model is not limited to the study and prediction of PUUV (and NE) occurrence in Europe, as it could be easily adapted to model other rodent-borne diseases, either with indirect or direct transmission.

Prevalence of pathogens in ticks collected from humans through citizen science in Belgium.

BACKGROUND: In order to evaluate the risk of human exposure to tick-borne pathogens in Belgium, a study on the prevalence of several pathogens was conducted on feeding ticks removed from humans in 2017. METHODS: Using a citizen science approach based on an existing notification tool for tick bites, a sample of ticks was collected across the country. Collected ticks were screened by PCR for the presence of the following pathogens: Anaplasma phagocytophilum, Babesia spp., Borrelia burgdorferi (sensu lato), Borrelia miyamotoi, Neoehrlichia mikurensis, Rickettsia helvetica and tick-borne encephalitis virus (TBEV). RESULTS: In total, 1599 ticks were included in the sample. The great majority of ticks belonged to Ixodes ricinus (99%); other tick species were identified as Ixodes hexagonus (0.7%) and Dermacentor reticulatus (0.3%). Borrelia burgdorferi (s.l.) was detected in 14% of nymphs and adult ticks. Adult ticks (20%) were more likely to be infected than nymphs (12%). The most common genospecies were B. afzelii (52%) and B. garinii (21%). Except for TBEV, the other tick-borne pathogens studied were all detected in the tick sample, although at a lower prevalence: 1.5% for Babesia spp.; 1.8% for A. phagocytophilum; 2.4% for B. miyamotoi; 2.8% for N. mikurensis; and 6.8% for R. helvetica. Rickettsia raoultii, the causative agent of tick-borne lymphadenopathy, was identified for the first time in Belgium, in two out of five D. reticulatus ticks. Co-infections were found in 3.9% of the examined ticks. The most common co-infection was B. burgdorferi (s.l.) + N. mikurensis. CONCLUSIONS: Although for most of the tick-borne diseases in Belgium, other than Lyme borreliosis, no or few cases of human infection are reported, the pathogens causing these diseases were all (except for TBEV) detected in the tick study sample. Their confirmed presence can help raise awareness among citizens and health professionals in Belgium on possible diseases other than Lyme borreliosis in patients presenting fever or other non-characteristic symptoms after a tick bite.

Identifying the patterns and drivers of Puumala hantavirus enzootic dynamics using reservoir sampling.

Hantaviruses are zoonotic hemorrhagic fever viruses for which prevention of human spillover remains the first priority in disease management. Tailored intervention measures require an understanding of the drivers of enzootic dynamics, commonly inferred from distorted human incidence data. Here, we use longitudinal sampling of approximately three decades of Puumala orthohantavirus (PUUV) evolution in isolated reservoir populations to estimate PUUV evolutionary rates, and apply these to study the impact of environmental factors on viral spread. We find that PUUV accumulates genetic changes at a rate of ∼10-4 substitutions per site per year and that land cover type defines the dispersal dynamics of PUUV, with forests facilitating and croplands impeding virus spread. By providing reliable short-term PUUV evolutionary rate estimates, this work facilitates the evaluation of spatial risk heterogeneity starting from timed phylogeographic reconstructions based on virus sampling in its animal reservoir, thereby side-stepping the need for difficult-to-collect human disease incidence data.

Combining primary care surveillance and a meta-analysis to estimate the incidence of the clinical manifestations of Lyme borreliosis in Belgium, 2015-2017.

Lyme borreliosis (LB) is an important tick-borne disease which can cause a broad range of symptoms mainly affecting the skin, the nervous system and the joints. This study aims to estimate the incidence of the different clinical manifestations of LB in Belgium. The incidence of erythema migrans (EM) was estimated through the network of sentinel general practices at 97.6/100,000 inhabitants (uncertainty interval [UI] 82.0-113.0) for the period 2015-2017. This result was used to estimate the incidence of other LB manifestations based on their proportional distribution (ratios) to EM reported in the neighboring countries of Belgium. To estimate these ratios, we performed a systematic review of studies published between February 1, 2008 and January 31, 2018 and pooled the results using a random effects meta-analysis. Six studies were retained in the systematic review, and the meta-analysis estimated the occurrence ratios for Lyme neuroborreliosis/EM, Lyme arthritis/EM and other manifestations/EM at 0.024 (95% confidence interval [CI] 0.016-0.037), 0.022 (95% CI 0.020-0.024) and 0.014 (95% CI 0.012-0.016) respectively. Applying these ratios to the EM incidence in Belgium resulted in an incidence estimation of 2.4/100,000 inhabitants (95% UI 1.5-3.7) for Lyme neuroborreliosis, 2.1/100,000 (95% UI 1.7-2.6) for Lyme arthritis and 1.4/100,000 (95% UI 1.1-1.7) for other less frequent manifestations. Some of these LB manifestations, other than EM, are more severe, hence these estimates are essential to assess the health burden and economic cost of LB which would be highly relevant for patients, healthcare providers and policymakers. As both over- and underestimation of different clinical LB manifestations remain possible due to characteristics of the primary surveillance systems and the disease itself, future studies to validate these estimates would be of great value.

The HUMTICK study: protocol for a prospective cohort study on post-treatment Lyme disease syndrome and the disease and cost burden of Lyme borreliosis in Belgium.

BACKGROUND: In Belgium, different routine surveillance systems are in place to follow-up Lyme borreliosis trends. However, accurate data on the disease and monetary burden for the different clinical manifestations are lacking. Despite recommended antibiotic treatment, a proportion of Lyme patients report persisting aspecific symptoms for six months or more (e.g. fatigue, widespread musculoskeletal pain, cognitive difficulties), a syndrome now named "post-treatment Lyme disease syndrome" (PTLDS). Controversy exists on the cause, incidence and severity of PTLDS. This study aims to estimate the incidence of PTLDS in patients with Lyme borreliosis and to quantify the disease burden and economic costs associated with the different clinical manifestations of Lyme borreliosis in Belgium. METHODS: The project is a prospective cohort study in which about 600 patients with an erythema migrans and 100 patients with disseminated Lyme borreliosis will be followed up. Questionnaires, including the SF-36 vitality and pain subscale, the Cognitive Failure Questionnaire and the EQ-5D-5L, will be used to collect information on acute and persisting symptoms and the impact on quality of life. Symptom frequency and severity will be compared with self-reported pre-Lyme health status, a control group and existing Belgian population norms. Additionally, information on the associated costs and possible risk factors for the development of PTLDS will be collected. DISCUSSION: A study of the health burden will allow evaluation of the relative importance of Lyme borreliosis in Belgium and information on the economic cost will help to formulate cost-effective measures. There are only few prospective studies conducted estimating the incidence of PTLDS and even though discussion exists about the prevalence of subjective symptoms in the general population, a control group of non-Lyme borreliosis participants has often not been included.

Excess all-cause and influenza-attributable mortality in Europe, December 2016 to February 2017.

Since December 2016, excess all-cause mortality was observed in many European countries, especially among people aged ≥ 65 years. We estimated all-cause and influenza-attributable mortality in 19 European countries/regions. Excess mortality was primarily explained by circulation of influenza virus A(H3N2). Cold weather snaps contributed in some countries. The pattern was similar to the last major influenza A(H3N2) season in 2014/15 in Europe, although starting earlier in line with the early influenza season start.

Transmission of Leishmania donovani in the Hills of Eastern Nepal, an Outbreak Investigation in Okhaldhunga and Bhojpur Districts.

BACKGROUND: In the Indian subcontinent, Visceral leishmaniasis is endemic in a geographical area coinciding with the Lower Gangetic Plain, at low altitude. VL occurring in residents of hill districts is therefore often considered the result of Leishmania donovani infection during travel. Early 2014 we conducted an outbreak investigation in Okhaldhunga and Bhojpur districts in the Nepal hills where increasing number of VL cases have been reported. METHODOLOGY/PRINCIPAL FINDINGS: A house-to-house survey in six villages documented retrospectively 35 cases of Visceral Leishmaniasis (VL). Anti-Leishmania antibodies were found in 22/23 past-VL cases, in 40/416 (9.6%) persons without VL and in 12/155 (7.7%) domestic animals. An age- and sex- matched case-control study showed that exposure to known VL-endemic regions was no risk factor for VL, but having a VL case in the neighbourhood was. SSU-rDNA PCR for Leishmania sp. was positive in 24 (5%) of the human, in 18 (12%) of the animal samples and in 16 (14%) bloodfed female Phlebotomus argentipes sand flies. L. donovani was confirmed in two asymptomatic individuals and in one sand fly through hsp70-based sequencing. CONCLUSIONS/SIGNIFICANCE: This is epidemiological and entomological evidence for ongoing local transmission of L. donovani in villages at an altitude above 600 meters in Nepal, in districts considered hitherto non-endemic for VL. The VL Elimination Initiative in Nepal should therefore consider extending its surveillance and control activities in order to assure VL elimination, and the risk map for VL should be redesigned.

Puumala hantavirus infection alters the odour attractiveness of its reservoir host.

The random-mixing assumptions of many parasite-transmission models are challenged if healthy individuals can alter their behaviour to reduce their risk of infection. Some pathogens reduce the attractiveness of their hosts' excretions, for example, potentially altering contact rates and thus the predicted force of infection for pathogens transmissible by contact with excretions. For bank voles (Myodes glareolus), contact with contaminated urine is an important route of transmission for Puumala hantavirus (PUUV); however, it is not known whether PUUV infection changes the voles' urinary odours or their attractiveness. Here, we use a Y-maze to test whether PUUV infection alters the attractiveness of male bank voles' urine. We presented wild-caught PUUV-free male and female bank voles with PUUV-infected conspecific urine, uninfected urine and a water control, and measured the relative and absolute latency to first visit, number of visits, and total time bank voles spent investigating each treatment over 30 min. PUUV infection significantly altered the bank voles' initial response to conspecific urine, with fewer visits and less time spent close to infected urine relative to uninfected urine, and less total time spent near the infected urine than the uninfected urine or control. These strong preferences weakened over the 30-min trial, however, partly due to a general decline in male activity, and there were no absolute differences between the treatments overall. This suggests that PUUV infection does change the attractiveness of bank vole urine to conspecifics, and we discuss the implications of these results for random-mixing assumptions.

Ecological niche modelling of bank voles in Western Europe.

The bank vole (Myodes glareolus) is the natural host of Puumala virus (PUUV) in vast areas of Europe. PUUV is one of the hantaviruses which are transmitted to humans by infected rodents. PUUV causes a general mild form of hemorrhagic fever with renal syndrome (HFRS) called nephropathia epidemica (NE). Vector-borne and zoonotic diseases generally display clear spatial patterns due to different space-dependent factors. Land cover influences disease transmission by controlling both the spatial distribution of vectors or hosts, as well as by facilitating the human contact with them. In this study the use of ecological niche modelling (ENM) for predicting the geographical distribution of bank vole population on the basis of spatial climate information is tested. The Genetic Algorithm for Rule-set Prediction (GARP) is used to model the ecological niche of bank voles in Western Europe. The meteorological data, land cover types and geo-referenced points representing the locations of the bank voles (latitude/longitude) in the study area are used as the primary model input value. The predictive accuracy of the bank vole ecologic niche model was significant (training accuracy of 86%). The output of the GARP models based on the 50% subsets of points used for testing the model showed an accuracy of 75%. Compared with random models, the probability of such high predictivity was low (χ(2) tests, p < 10(-6)). As such, the GARP models were predictive and the used ecologic niche model indeed indicates the ecologic requirements of bank voles. This approach successfully identified the areas of infection risk across the study area. The result suggests that the niche modelling approach can be implemented in a next step towards the development of new tools for monitoring the bank vole's population.

Impact of Puumala virus infection on maturation and survival in bank voles: a capture-mark-recapture analysis.

Many zoonotic diseases are caused by rodent-borne viruses. Major fluctuations in the transmission of these viruses have been related to large changes in reservoir host population numbers due to external factors. However, the impact of the pathogen itself on the demography of its reservoir host is often overlooked. We investigated the impact of Puumala virus (PUUV) on survival and reproductive maturation probability of its reservoir host, the bank vole (Myodes glareolus). Three years (2004-06) of data from nine independent sites in southern Belgium were collected and analyzed with a capture-mark-recapture (CMR) method that includes statistical correction for the variation in capture probability of voles. A multistate model based on four states of reproductive activity and PUUV immunoglobulin G (IgG) antibody status was used to estimate survival and probability of transition from one reproductive or infection state to another. Although survival estimates for reproductively active voles were similar between infected and noninfected individuals, PUUV infection in reproductively inactive voles decreased mean monthly survival by 14%. PUUV infection was associated with a threefold increase in the probability of reproductive maturation in bank voles. Moreover, the probability of PUUV IgG seroconversion was three times higher for reproductively active voles compared to reproductively inactive voles. Our model indicates that PUUV infection may alter bank vole population dynamics by affecting both survival and maturation in its host. Additional studies, using CMR methodology with shorter time intervals between trapping sessions and possibly a longer duration, are needed to confirm these findings.

Temporal variation in individual factors associated with hantavirus infection in bank voles during an epizootic: implications for Puumala virus transmission dynamics.

Puumala virus (PUUV), the causal agent of nephropathia epidemica in humans, is one of the many hantaviruses included in the list of emerging pathogens. Hantavirus infection is not distributed evenly among PUUV reservoir hosts (i.e., bank voles [Myodes glareolus]). Besides environmental factors and local population features, individual characteristics play an important role in vole PUUV infection risk. Identifying the relative importance of these individual characteristics can provide crucial information on PUUV transmission processes. In the present study, bank voles were monitored during the nephropathia epidemica outbreak of 2005 in Belgium. Vole sera were tested for presence of immunoglobulin G against PUUV, and a logistic mixed model was built to investigate the temporal variation in individual characteristics and their relative importance to PUUV infection risk in bank voles. Relative risk calculations for individual vole characteristics related to PUUV infection in the reservoir host show that reproductive activity dominates infection risk. The gender effect is only found in reproductively active voles, where reproductively active males have the highest infection risk. Results also revealed a clear seasonal variation in the importance of reproductive activity linked to PUUV infection. In contrast to the main effect found in other trapping sessions, no difference in infection risk ratio was found between reproductively active and nonactive voles in the spring period. Combined with increased infection risk for the reproductively nonactive group at that time, these results indicate a shift in the transmission process due to changes in bank vole behavior, physiology, or climate conditions. Hence, our results suggest that mathematical models should take into account seasonal shifts in transmission mechanisms. When these results are combined with the seasonal changes in population structure during the epizootic period, we identify vole reproductive activity and length of the breeding season as potential drivers of PUUV epizootics in west-central European regions.

Associations between MHC genes and Puumala virus infection in Myodes glareolus are detected in wild populations, but not from experimental infection data.

We analysed the influence of MHC class II Dqa and Drb genes on Puumala virus (PUUV) infection in bank voles (Myodes glareolus). We considered voles sampled in five European localities or derived from a previous experiment that showed variable infection success of PUUV. The genetic variation observed in the Dqa and Drb genes was assessed by using single-strand conformation polymorphism and pyrosequencing methods, respectively. Patterns were compared with those obtained from 13 microsatellites. We revealed significant genetic differentiation between PUUV-seronegative and -seropositive bank voles sampled in wild populations, at the Drb gene only. The absence of genetic differentiation observed at neutral microsatellites confirmed the important role of selective pressures in shaping these Drb patterns. Also, we found no significant associations between infection success and MHC alleles among laboratory-colonized bank voles, which is explained by a loss of genetic variability that occurred during the captivity of these voles.

Environmental conditions and Puumala virus transmission in Belgium.

BACKGROUND: Non-vector-borne zoonoses such as Puumala hantavirus (PUUV) can be transmitted directly, by physical contact between infected and susceptible hosts, or indirectly, with the environment as an intermediate. The objective of this study is to better understand the causal link between environmental features and PUUV prevalence in bank vole population in Belgium, and hence with transmission risk to humans. Our hypothesis was that environmental conditions controlling the direct and indirect transmission paths differ, such that the risk of transmission to humans is not only determined by host abundance. We explored the relationship between, on one hand, environmental variables and, on the other hand, host abundance, PUUV prevalence in the host, and human cases of nephropathia epidemica (NE). Statistical analyses were carried out on 17 field sites situated in Belgian broadleaf forests. RESULTS: Linear regressions showed that landscape attributes, particularly landscape configuration, influence the abundance of hosts in broadleaf forests. Based on logistic regressions, we show that PUUV prevalence among bank voles is more linked to variables favouring the survival of the virus in the environment, and thus the indirect transmission: low winter temperatures are strongly linked to prevalence among bank voles, and high soil moisture is linked to the number of NE cases among humans. The transmission risk to humans therefore depends on the efficiency of the indirect transmission path. Human risk behaviours, such as the propensity for people to go in forest areas that best support the virus, also influence the number of human cases. CONCLUSION: The transmission risk to humans of non-vector-borne zoonoses such as PUUV depends on a combination of various environmental factors. To understand the complex causal pathways between the environment and disease risk, one should distinguish between environmental factors related to the abundance of hosts such as land-surface attributes, landscape configuration, and climate - i.e., host ecology, - and environmental factors related to PUUV prevalence, mainly winter temperatures and soil moisture - i.e., virus ecology. Beyond a threshold abundance of hosts, environmental factors favouring the indirect transmission path (soil and climate) can better predict the number of NE cases among humans than factors influencing the abundance of hosts.

Determinants of the geographic distribution of Puumala virus and Lyme borreliosis infections in Belgium.

BACKGROUND: Vector-borne and zoonotic diseases generally display clear spatial patterns due to different space-dependent factors. Land cover and land use influence disease transmission by controlling both the spatial distribution of vectors or hosts, and the probability of contact with susceptible human populations. The objective of this study was to combine environmental and socio-economic factors to explain the spatial distribution of two emerging human diseases in Belgium, Puumala virus (PUUV) and Lyme borreliosis. Municipalities were taken as units of analysis. RESULTS: Negative binomial regressions including a correction for spatial endogeneity show that the spatial distribution of PUUV and Lyme borreliosis infections are associated with a combination of factors linked to the vector and host populations, to human behaviours, and to landscape attributes. Both diseases are associated with the presence of forests, which are the preferred habitat for vector or host populations. The PUUV infection risk is higher in remote forest areas, where the level of urbanisation is low, and among low-income populations. The Lyme borreliosis transmission risk is higher in mixed landscapes with forests and spatially dispersed houses, mostly in wealthy peri-urban areas. The spatial dependence resulting from a combination of endogenous and exogenous processes could be accounted for in the model on PUUV but not for Lyme borreliosis. CONCLUSION: A large part of the spatial variation in disease risk can be explained by environmental and socio-economic factors. The two diseases not only are most prevalent in different regions but also affect different groups of people. Combining these two criteria may increase the efficiency of information campaigns through appropriate targeting.

Immunotoxicology in wood mice along a heavy metal pollution gradient.

We carried out an immunotoxicological field study of wood mice in three populations along a heavy metal pollution gradient. Heavy metal concentrations in liver tissue indicated that exposure to silver, arsenic, cadmium, cobalt and lead decreased with increasing distance from a non-ferrous smelter. Host resistance to the endoparasite Heligmosomoides polygyrus decreased with increasing exposure, while the abundance of tick larvae and the nematode Syphacia stroma was unrelated to heavy metal exposure. Spleen mass was increased at the intermediate and the most polluted sites and was positively correlated with the number of H. polygyrus and tick larvae. Proportion of early apoptotic leukocytes increased towards the smelter and was positively related to cadmium exposure. Red and white blood cell counts and lysozyme activity showed no relationship with metal exposure. All together, our observations suggest negative effects of heavy metal exposure on the immune function of wood mice under field conditions.