Towards integrated surveillance of zoonoses: spatiotemporal joint modeling of rodent population data and human tularemia cases in Finland.

BACKGROUND: There are an increasing number of geo-coded information streams available which could improve public health surveillance accuracy and efficiency when properly integrated. Specifically, for zoonotic diseases, knowledge of spatial and temporal patterns of animal host distribution can be used to raise awareness of human risk and enhance early prediction accuracy of human incidence. METHODS: To this end, we develop a spatiotemporal joint modeling framework to integrate human case data and animal host data to offer a modeling alternative for combining multiple surveillance data streams in a novel way. A case study is provided of spatiotemporal modeling of human tularemia incidence and rodent population data from Finnish health care districts during years 1995-2012. RESULTS: Spatial and temporal information of rodent abundance was shown to be useful in predicting human cases and in improving tularemia risk estimates in 40 and 75% of health care districts, respectively. The human relative risk estimates' standard deviation with rodent's information incorporated are smaller than those from the model that has only human incidence. CONCLUSIONS: These results support the integration of rodent population variables to reduce the uncertainty of tularemia risk estimates. However, more information on several covariates such as environmental, behavioral, and socio-economic factors can be investigated further to deeper understand the zoonotic relationship.

Serological survey in the Finnish human population implies human-to-human transmission of Ljungan virus or antigenically related viruses.

Ljungan virus (LV) is a picornavirus related to human parechoviruses (HPeV). The virus has been found in bank voles (Myodes glareolus) and several other rodent species, and suggested to have zoonotic potential. Thus far, seroepidemiological data on LV infections in humans are scarce. In this study, we aimed to characterize the demographic and geographical distribution of LV-reactive antibodies in Finland, and to investigate its occurrence in patients suspected of having a rodent-borne disease, nephropathia epidemica (NE) caused by Puumala hantavirus (PUUV). Using an immunofluorescence assay (LV strain 145SLG), we screened human sera (n = 1378) and found LV-reactive antibodies in 36% of samples. The probability of possessing LV-reactive antibodies peaked at age of 14 years, suggesting that most infections occur in childhood. The prevalence of LV-reactive antibodies was significantly higher in the urbanized area surrounding Helsinki than in more rural Central Finland. These findings are uncharacteristic of a rodent-borne pathogen, and therefore we consider human-to-human transmission of one or several Ljungan-like viruses as a likely cause for most of the observed antibody responses.

Aleutian mink disease virus in free-ranging mustelids in Finland - a cross-sectional epidemiological and phylogenetic study.

Aleutian mink disease virus (AMDV) can cause severe immune-complex-mediated disease in American mink. AMDV has also been detected in several other mustelid species with potential negative impact on their health and population. A molecular and cross-sectional epidemiological study was conducted to obtain data on the prevalence, distribution, transmission and diversity of AMDV strains in Finnish free-ranging mustelids and risk factors associated with infection. The presence of anti-AMDV antibodies and/or AMDV DNA was tested from 308 samples representing eight mustelid species and 17 administrative regions. Positive samples were detected across Finland, and in 54 % (31/57) of feral American mink, 27 % (7/26) of European badgers and 7 % (1/14) of European polecats. Samples from Eurasian otters, European pine martens, least weasels, stoat and wolverine were negative. Major risk factors for infection were the species American mink with 335 and badger with 74 times higher odds than other species, and the years 2006-2009 with five times higher odds than the years 2010-2014. No clustering according to species, geographical origin or year was evident in phylogeny, except for four divergent sequences from Estonian badgers that formed a separate phylogroup distinct from other AMDV strains. This study showed that AMDV was prevalent in certain species of Finnish free-ranging mustelids and widely distributed across Finland. Furthermore, the free-ranging mustelids carried both strains similar to those found in farmed mink, but also distinct strains that may represent novel amdoparvoviruses.

Recent increase in prevalence of antibodies to Dobrava-Belgrade virus (DOBV) in yellow-necked mice in Northern Italy.

Dobrava-Belgrade virus (DOBV) is the most pathogenic hantavirus in Europe with a case-fatality rate of up to 12%. To detect changes in risk for humans, the prevalence of antibodies to DOBV has been monitored in a population of Apodemus flavicollis in the province of Trento (northern Italy) since 2000, and a sudden increase was observed in 2010. In the 13-year period of this study, 2077 animals were live-trapped and mean hantavirus seroprevalence was 2·7% (s.e. = 0·3%), ranging from 0% (in 2000, 2002 and 2003) to 12·5% (in 2012). Climatic (temperature and precipitation) and host (rodent population density, rodent weight and sex, and larval tick burden) variables were analysed using Generalized Linear Models and multi-model inference to select the best model. Climatic changes (mean annual precipitation and maximum temperature) and individual body mass had a positive effect on hantavirus seroprevalence. Other possible drivers affecting the observed pattern need to be studied further.

Incidence and seroprevalence of tularaemia in Finland, 1995 to 2013: regional epidemics with cyclic pattern.

We studied the incidence of reported tularaemia by year and region and the prevalence of antibodies against Francisella tularensis in the adult general population in Finland. Moreover, we assessed the correlation between vole population cycles and human tularaemia outbreaks. The seroprevalence study made use of serum samples from a nationwide population-based health survey (Health 2000). The samples of 1,045 randomly selected persons, representative for the Finnish population in each region, were screened with an enzyme-linked immunosorbent assay (ELISA) for the presence of IgG antibodies against F. tularensis, and positive results were further confirmed by immunoblotting. A serological response to F. tularensis was found in 2% (95% confidence interval: 1.1­3.5) of the population. Incidence and seroprevalence were highest in the same areas, and vole population peaks clearly preceded tularaemia outbreaks one year later.

Landscape features and helminth co-infection shape bank vole immunoheterogeneity, with consequences for Puumala virus epidemiology.

Heterogeneity in environmental conditions helps to maintain genetic and phenotypic diversity in ecosystems. As such, it may explain why the capacity of animals to mount immune responses is highly variable. The quality of habitat patches, in terms of resources, parasitism, predation and habitat fragmentation may, for example, trigger trade-offs ultimately affecting the investment of individuals in various immunological pathways. We described spatial immunoheterogeneity in bank vole populations with respect to landscape features and co-infection. We focused on the consequences of this heterogeneity for the risk of Puumala hantavirus (PUUV) infection. We assessed the expression of the Tnf-α and Mx2 genes and demonstrated a negative correlation between PUUV load and the expression of these immune genes in bank voles. Habitat heterogeneity was partly associated with differences in the expression of these genes. Levels of Mx2 were lower in large forests than in fragmented forests, possibly due to differences in parasite communities. We previously highlighted the positive association between infection with Heligmosomum mixtum and infection with PUUV. We found that Tnf-α was more strongly expressed in voles infected with PUUV than in uninfected voles or in voles co-infected with the nematode H. mixtum and PUUV. H. mixtum may limit the capacity of the vole to develop proinflammatory responses. This effect may increase the risk of PUUV infection and replication in host cells. Overall, our results suggest that close interactions between landscape features, co-infection and immune gene expression may shape PUUV epidemiology.

Landscape genetics highlights the role of bank vole metapopulation dynamics in the epidemiology of Puumala hantavirus.

Rodent host dynamics and dispersal are thought to be critical for hantavirus epidemiology as they determine pathogen persistence and transmission within and between host populations. We used landscape genetics to investigate how the population dynamics of the bank vole Myodes glareolus, the host of Puumala hantavirus (PUUV), vary with forest fragmentation and influence PUUV epidemiology. We sampled vole populations within the Ardennes, a French PUUV endemic area. We inferred demographic features such as population size, isolation and migration with regard to landscape configuration. We next analysed the influence of M. glareolus population dynamics on PUUV spatial distribution. Our results revealed that the global metapopulation dynamics of bank voles were strongly shaped by landscape features, including suitable patch size and connectivity. Large effective size in forest might therefore contribute to the higher observed levels of PUUV prevalence. By contrast, populations from hedge networks highly suffered from genetic drift and appeared strongly isolated from all other populations. This might result in high probabilities of local extinction for both M. glareolus and PUUV. Besides, we detected signatures of asymmetric bank vole migration from forests to hedges. These movements were likely to sustain PUUV in fragmented landscapes. In conclusion, our study provided arguments in favour of source-sink dynamics shaping PUUV persistence and spread in heterogeneous, Western European temperate landscapes. It illustrated the potential contribution of landscape genetics to the understanding of the epidemiological processes occurring at this local scale.

Intestinal parasites of the Arctic fox in relation to the abundance and distribution of intermediate hosts.

The intestinal parasite community of Arctic foxes (Vulpes lagopus) on the Svalbard archipelago in the High Arctic was investigated in relation to the abundance and distribution of intermediate hosts. Five species of cestodes (Echinococcus multilocularis, Taenia crassiceps, Taenia polyacantha, Taenia krabbei and Diphyllobothrium sp.), ascaridoid nematodes and one unidentified acanthocephalan species were found. The cestodes E. multilocularis, T. crassiceps and T. polyacantha all showed a decreasing prevalence in the fox population with increasing distance from their spatially restricted intermediate host population of sibling voles (Microtus levis). In addition, the prevalence of E. multilocularis in a sample from the vole population was directly related to the local vole abundance. The cestode T. krabbei uses reindeer as intermediate host, and its prevalence in female foxes was positively related to the density of reindeer (Rangifer tarandus platyrhyncus). Finally, the prevalence of the ascaridoid nematodes also decreased with increasing distance from the vole population, a finding that is consistent with the idea that voles are involved in transmission, most likely as paratenic hosts. The prevalence of the remaining species (Diphyllobothrium sp. and an unidentified acanthocephalan) was very low. We conclude that the distribution and abundance of intermediate host structure the gastrointestinal parasite community of the Arctic fox on the Svalbard archipelago.

Systematic relationships of Mosgovoyia Spasskii, 1951 (Cestoda: Anoplocephalidae) and related genera inferred from mitochondrial and nuclear sequence data.

The present study evaluates the phylogenetic position and systematic relationships of two species of Mosgovoyia Spasskii, 1951 and related genera (Cestoda: Anoplocephalidae) based on sequences of 28S ribosomal RNA and mitochondrial NADH dehydrogenase subunit 1 (Nad1) genes. Both molecular data-sets show that M. pectinata (Goeze, 1782) and Schizorchis caballeroi Rausch, 1960 are sister species and that they are phylogenetically independent from M. ctenoides (Railliet, 1890). This shows unambiguously that Mosgovoyia [sensu Beveridge (1978)] is a non-monophyletic assemblage, supporting the validity of Neoctenotaenia Tenora, 1976, erected for M. ctenoides. The results also show that the morphologically related Ctenotaenia marmotae (Fröhlich, 1802) is the sister species of Andrya rhopalocephala (Riehm, 1881) and therefore represents a more derived lineage. Modified diagnoses are provided for Mosgovoyia and Neoctenotaenia.

Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics.

We investigated the factors mediating selection acting on two MHC class II genes (DQA and DRB) in water vole (Arvicola scherman) natural populations in the French Jura Mountains. Population genetics showed significant homogeneity in allelic frequencies at the DQA1 locus as opposed to neutral markers (nine microsatellites), indicating balancing selection acting on this gene. Moreover, almost exhaustive screening for parasites, including gastrointestinal helminths, brain coccidia and antibodies against viruses responsible for zoonoses, was carried out. We applied a co-inertia approach to the genetic and parasitological data sets to avoid statistical problems related to multiple testing. Two alleles, Arte-DRB-11 and Arte-DRB-15, displayed antagonistic associations with the nematode Trichuris arvicolae, revealing the potential parasite-mediated selection acting on DRB locus. Selection mechanisms acting on the two MHC class II genes thus appeared different. Moreover, overdominance as balancing selection mechanism was showed highly unlikely in this system.

A phylogeny of members of the family Taeniidae based on the mitochondrial cox1 and nad1 gene data.

The cestode family Taeniidae consists of 2 genera, Taenia and Echinococcus, which both have been the focus of intensive taxonomic and epidemiological studies because of their zoonotic importance. However, a comprehensive molecular phylogeny of this family has yet to be reconstructed. In this study, 54 isolates representing 9 Taenia species were characterized using DNA sequences in the mitochondrial cytochrome c oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (nad1) genes. Phylogenetic relationships within the family Taeniidae were inferred by combining cox1 and nad1 sequence data of the present and previous studies. In the phylogenetic analysis, the genus Echinococcus was shown to be monophyletic, but Taenia proved to be paraphyletic due to the position of T. mustelae as a probable sister taxon of Echinococcus. This indicates that T. mustelae should form a genus of its own. Taenia ovis krabbei was placed distant from T. ovis ovis, as a sister taxon of T. multiceps, supporting its recognition as a distinct species, T. krabbei. High intraspecific sequence variation within both T. polyacantha and T. taeniaeformis suggests the existence of cryptic sister species.

Human noroviruses in the faeces of wild birds and rodents-new potential transmission routes.

Human noroviruses (HuNoVs) are one of the leading global causes of diarrhoeal diseases and are transmitted mainly from person to person but also through contaminated food, water and fomites. The possible zoonotic nature of NoVs has occasionally been discussed, although the viruses are generally considered to be host-species-specific. We investigated whether wild birds and rodents could serve as carriers of HuNoVs, thereby transmitting the virus to humans directly or indirectly by contaminating foods. All samples, 115 avian and 100 rat faeces collected in springs 2009-2013 from dump sites, and 85 faeces from yellow-necked mice trapped in late autumn 2008 and 2009 after the rodents entered human settlements due to the first night frosts, were screened for HuNoV using real-time reverse transcription PCR. HuNoVs were detected in 31 (27%) faecal samples of wild birds, in two (2%) faecal samples of rats and in no samples of mice. Most (25) of the positive bird samples and both rat samples contained genogroup II, and six positive bird samples contained genogroup I HuNoV. The avian species shedding faeces containing HuNoVs were identified as gulls and crows using DNA barcoding. Our results show that wildlife, birds and rats in particular, is capable of spreading HuNoVs in the environment.

Serological survey for viral pathogens in Turkish rodents.

Wild rodents (n = 330) were trapped around the villages of Altindere and Cosandere (Maçka, Trabzon Province), Ayder, Ortan, and Yolkiyi (Camlihemsin, Rize Province), and Bozdag (Odemis, Izmir Province) in northeastern and western Turkey during April 2004. Samples were tested for arenavirus, hantavirus, and cowpox virus (family Poxviridae, genus Orthopoxvirus, CPXV) antibodies by using immunofluorescence assays (IFAs). Antibodies against arenaviruses were found in eight of 330 (2.4%) rodents. Arenavirus sero-positive animals were found from all study sites. Antibodies to Puumala virus (family Bunyaviridae, genus Hantavirus, PUUV) were detected in four of 65 Microtus voles tested. Of the PUUV-IFA-positive voles, one Microtus guentheri lydius was caught from Izmir, and one Microtus roberti and two Microtus rossiaemeridionalis were captured near Trabzon. All 264 Apodemus spp. mice tested negative for antibodies to Saaremaa virus (family Bunyaviridae, genus Hantavirus, SAAV); the single Dryomys nitedula tested negative for both PUUV and SAAV antibodies. Only one (0.3%) of the rodents, an Apodemus sylvaticus from Trabzon area, tested seropositive to CPXV. This is the first serologic survey for rodent-borne viruses in their natural hosts in Turkey. Although these preliminary results support presence of several virus groups with zoonotic potential, additional studies are needed to identify the specific viruses that are present in these populations.

Description of Paranoplocephala yoccozi n. sp. (Cestoda: Anoplocephalidae) from the snow vole Chionomys nivalis in France, with a review of anoplocephaud cestodes of snow voles in Europe.

We describe Paranoplocephala yoccozi n. sp. (Cestoda: Anoplocephalidae) from the snow vole Chionomys nivalis in Bourg-Saint-Maurice, French Alps, compare it with several related species from rodents, and review the anoplocephalid cestodes of snow voles in Europe. Paranoplocephala yoccozi n. sp. is primarily distinguished from the related species by its large scolex of characteristic shape, robust neck region, and the structure of the cirrus sac, vitellarium and vagina. We show that the anoplocephalid cestodes of snow voles in Europe, representing the genera Anoplocephaloides and Paranoplocephala, include at least seven species. This fauna consists primarily of species that snow voles share with other voles inhabiting the high-mountain areas. Some of the species, including P. yoccozi n. sp., appear to have a very localized distribution, which is assumed to be a consequence of the historical fragmentation of snow vole populations.

Genetic variation of wild Puumala viruses within the serotype, local rodent populations and individual animal.

Reverse transcriptase polymerase chain reaction cloning and sequencing were used to determine the range of S gene/N protein variability in wild Puumala virus (PUU) strains and to study phylogenetic relationships between two groups of strains which originated from Finland and from European Russia. Analyses of the nucleotide and predicted amino acid sequences showed: (1) all PUU strains shared a common ancient ancestor; and (2) the more recent ancestors were different for the Finnish branch and the Russian branch of PUU strains. A cluster of amino acid substitutions in the N protein of Finnish strains was found; this cluster was located within a highly variable region of the molecule carrying B-cell epitopes (Vapalahti et al., J. Med. Virol., 1995, in press). Different levels of S gene/N protein diversity of PUU were revealed supporting the view of geographical clustering of genetic variants. Puumala virus from individual voles was found to be a complex mixture of closely related variants-quasispecies. The ratio of non-silent to silent nucleotide mutations registered in the S genes/N proteins of PUU quasispecies was 4- to 16-fold higher than that in Puumala virus strains, resulting in a more wide range of quasispecies N protein sequence diversity.

Variability of sex ratio, mating probability and egg production in an intestinal nematode in its fluctuating host population.

In dioecious parasites, the chances of co-occurring with the opposite sex depend on the sex ratio, abundance and distribution pattern of parasites within the host population. Theory suggests that if the abundance and degree of aggregation are very low, mating probability may decrease so much that the parasite population is vulnerable to extinction. Our aim is to determine the factors affecting the mating probability and egg production in Heligmosomum mixtum (Heligmosomidae), an intestinal nematode of the bank vole, Clethrionomys glareolus, at Pallasjärvi, Finnish Lapland. We also search for factors responsible for the persistence of H. mixtum in its fluctuating host population. The results showed that during high parasite abundance practically all nematode females co-occurred with males, but during a phase of very low abundance only 15% of females had a chance to mate. Comparison of observed mating probabilities and those predicted by a theoretical model (May & Woolhouse, 1993) showed that deviation from the assumption of complete aggregation between males and females results in underestimation of the mating probability. Sex ratio and the degree of aggregation showed a minor effect on mating probability. The sex ratio (proportion of females) of H. mixtum, which was female-biased (0.58), showed a negative correlation with the mean intensity of infection in the monthly samples (decreasing female-bias at high mean intensity), but no significant relation to the intensity of infection (number of worms in a host individual). The long-term persistence of H. mixtum in its strongly varying host population seems to be due to the high transmission efficiency and long life-span of the parasite.

Disease patterns in field and bank vole populations during a cyclic decline in central Finland.

Declining field vole (Microtus agrestis) and bank vole (Clethrionomys glareolus) populations were sampled (117 field voles and 34 bank voles) in south-central Finland during the winter of 1988-89. The last surviving field voles were caught in April and bank voles in February. A subsample (16) of the April field voles were taken live to the laboratory for immunosuppression. The histopathology of the main internal organs and the presence of aerobic bacteria and certain parasites were studied. In the lungs, an increase in lymphoid tissue, probably caused by infections, was the most common finding (52% of all individuals). The prevalences in the voles, in the whole material, of Chrysosporium sp. and Pneumocystis carinii in lungs were 13 and 10% in field voles, and 9 and 0% in bank voles, respectively. Cysts of Taenia mustelae (9 and 27%) were the most common pathological changes in the liver. Enteritis was also rather common (14 and 34%). In field voles the prevalences of Frenkelia sp. in the brain and Sarcocystis sp. in leg muscles were low (both 6%). Bordetella bronchiseptica was commonly (31%) isolated from field vole lungs and Listeria monocytogenes from the intestines (34%). Salmonella spp. could not be found. The dynamics and abundance of inflammations in the lungs and intestines, as well as B. bronchiseptica isolations from the lungs, indicate that obvious epidemics took place in declining vole populations. Of the Luhanka subsample of 16 field voles brought to the laboratory in April, one died of listeriosis, two of Bordetella, and five died for unknown reasons. Even if small mustelids are the driving force in microtine cycles, it is possible that diseases also contribute to the decline.

Hantavirus infections in Finland.

Hantaviruses cause haemorrhagic fever with renal syndrome (HFRS) and hantaviruspulmonary syndrome (HPS), and occur throughout the world. The viruses aretransmitted to humans in aerosols of rodent secretions. About 200 000 people,mainly in Eurasia, develop.

Hantavirus infections in Finland.

Hantaviruses cause haemorrhagic fever with renal syndrome (HFRS) and hantaviruspulmonary syndrome (HPS), and occur throughout the world. The viruses aretransmitted to humans in aerosols of rodent secretions. About 200 000 people,mainly in Eurasia, develop

Haemogregarines of the genus Hepatozoon (Apicomplexa: Adeleina) in rodents from northern Europe.

We studied the prevalence and distribution of Hepatozoon infections in small rodents from Finland and other areas in northern Europe. Hepatozoon infections were more common in voles (Arvicolinae) than mice (Murinae) and more prevalent in voles of the genus Clethrionomys than in voles of the genus Microtus. Transmission electron microscopical examination of Hepatozoon erhardovae Krampitz, 1964 from bank voles Clethrionomys glareolus (Schreber) showed that intracellular lung meronts were located in alveolar septa. Meronts consisted of varying numbers of merozoites packed with amylopectin vacuoles inside electron-lucent parasitophorous vacuole. The size of the meronts was approximately 19 x 14 microm. Monozoic or dizoic cysts were frequent findings in the lung alveoles; the size of cysts was approximately 10 x 6 microm. Gametocytes were found inside eosinophilic granulocytes in the capillaries of lung tissue. Ultrastructurally, micronemes, microtubules, mitochondria, nuclei and lipid droplets were visible.