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.

Graphical Abstract

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⁻⁴ 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.

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.

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.

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.

Visceral leishmaniasis is a neglected but deadly disease occurring in north-eastern India, the south-eastern lowland of Nepal, and the Ganges delta in Bangladesh; all part of the Lower Gangetic plains. Districts at higher altitude, such as those situated in the foothills of the Himalaya in Nepal, are considered non-endemic. Consequently, diagnostic, therapeutic and surveillance facilities are not available, and sporadic cases of VL occurring in residents of these districts are considered the result of Leishmania donovani infection during travel. This parasite is transmitted from man to man through the bite of a sand fly, Phlebotomus argentipes. To investigate the increasing number of cases reported from some villages in the Nepal hills, we interviewed former VL cases, as well as their healthy co-villagers, on their history of travel, we checked their blood for signs of infection, and we set traps to capture sand flies. We found that many residents had been infected, symptomatically as well as asymptomatically, regardless of travelling to endemic areas. Moreover, we managed to capture sand flies and found DNA of leishmania parasites inside. This proves that there is indeed ongoing local transmission in hill districts and that surveillance and control activities should be extended.

Excess mortality among the elderly in European countries, December 2014 to February 2015

Since December 2014 and up to February 2015, the weekly number of excess deaths from all-causes among individuals?≥?65 years of age in 14 European countries have been significantly higher than in the four previous winter seasons. The rise in unspecified excess mortality coincides with increased proportion of influenza detection in the European influenza surveillance schemes with a main predominance of influenza A(H3N2) viruses seen throughout Europe in the current season, though cold snaps and other respiratory infections may also have had an effect.

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.

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.

Population, environmental, and community effects on local bank vole (Myodes glareolus) Puumala virus infection in an area with low human incidence

In this study, the distribution of Puumala hantavirus (PUUV) infection in local bank vole Myodes glareolus populations in an area with low human PUUV infection (nephropathia epidemica [NE]) incidence in northern Belgium was monitored for 2 consecutive years. Bank voles were trapped in preferred habitat and tested for anti-PUUV IgG. Infection data were related to individual bank vole features, population demography, and environmental variables. Rare occurrence of PUUV infection was found and PUUV prevalence was low compared with data from the high NE incidence area in southern Belgium. Small-scale climatic differences seemed to play a role in PUUV occurrence, vegetation index and deciduous forest patch size both influenced PUUV prevalence and number of infected voles in a positive way. The data suggested a density threshold in vole populations below which PUUV infection does not occur. This threshold may vary between years, but the abundance of bank voles does not seem to affect the degree of PUUV seroprevalence further. We found indications for a dilution effect on PUUV prevalence, dependent on the relative proportion of nonhost wood mice Apodemus sylvaticus in a study site. In conclusion, we regard the combination of a dilution effect, a possible threshold density that depends on local conditions, and a higher fragmentation of suitable bank vole habitat in our study area as plausible explanations for the sparse occurrence of PUUV infection and low prevalence detected. Thus, beside human activity patterns, local environmental conditions and rodent community structure are also likely to play a role in determining PUUV infection risk for humans.

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.