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

Ecological determinants driving orthohantavirus prevalence in small mammals of Europe: a systematic review

Orthohantaviruses are emerging zoonotic pathogens that cause severe human disease and are considered an emerging public health threat globally. Mammalian orthohantaviruses are naturally maintained in rodent species and occasionally in other mammals. The abundance and density of natural orthohantavirus reservoir species are affected by multi annual and seasonal population cycles, community composition, ecosystem variables and climate. Horizontal transmission between host species is mostly density-driven and occurs via contact with infected host excreta, thus, fluctuations in populations and environmental variables often determine the prevalence of hantavirus in natural hosts. Given the zoonotic potential of hantaviruses, ecological factors influencing their spread and persistence in their natural reservoir and population dynamics influencing horizontal transmission require critical evaluation for human infection risk assessment. The present review paper discusses the impacts of natural host population cycles and ecosystem diversity, environmental conditions, and abiotic factors on the epidemiology of rodent-borne hantavirus infections in Europe. While significant efforts have been made to understand the drivers of hantavirus prevalence in natural hosts, we highlight key challenges in evaluating viral prevalence and assessing the role of environmental and population variables in determining hantavirus prevalence in host species.

A Bocage Landscape Restricts the Gene Flow of Pest Vole Populations

The population dynamics of most animal species inhabiting agro-ecosystems may be determined by landscape characteristics, with agricultural intensification and the reduction of natural habitats influencing dispersal and hence limiting gene flow. Increasing landscape complexity would thus benefit many endangered species by providing different ecological niches, but it could also lead to undesired effects in species that can act as crop pests and disease reservoirs. We tested the hypothesis that a highly variegated landscape influences patterns of genetic structure in agricultural pest voles. Ten populations of fossorial water vole, Arvicola scherman, located in a bocage landscape in Atlantic NW Spain were studied using DNA microsatellite markers and a graph-based model. The results showed a strong isolation-by-distance pattern with a significant genetic correlation at smaller geographic scales, while genetic differentiation at larger geographic scales indicated a hierarchical pattern of up to eight genetic clusters. A metapopulation-type structure was observed, immersed in a landscape with a low proportion of suitable habitats. Matrix scale rather than matrix heterogeneity per se may have an important effect upon gene flow, acting as a demographic sink. The identification of sub-populations, considered to be independent management units, allows the establishment of feasible population control efforts in this area. These insights support the use of agro-ecological tools aimed at recreating enclosed field systems when planning integrated managements for controlling patch-dependent species such as grassland voles.

Host genetic factors associated with the range limit of a European hantavirus

The natural host ranges of many viruses are restricted to very specific taxa. Little is known about the molecular barriers between species that lead to the establishment of this restriction or generally prevent virus emergence in new hosts. Here, we identify genomic polymorphisms in a natural rodent host associated with a strong genetic barrier to the transmission of European Tula orthohantavirus (TULV). We analysed the very abrupt spatial transition between two major phylogenetic clades in TULV across the comparatively much wider natural hybrid zone between evolutionary lineages of their reservoir host, the common vole (Microtus arvalis). Genomic scans of 79,225 single nucleotide polymorphisms (SNPs) in 323 TULV‐infected host individuals detected 30 SNPs that were consistently associated with the TULV clades CEN.S or EST.S in two replicate sampling transects. Focusing the analysis on 199 voles with evidence of genomic admixture at the individual level (0.1–0.9) supported statistical significance for all 30 loci. Host genomic variation at these SNPs explained up to 37.6% of clade‐specific TULV infections. Genes in the vicinity of associated SNPs include SAHH, ITCH and two members of the Syngr gene family, which are involved in functions related to immune response or membrane transport. This study demonstrates the relevance of natural hybrid zones as systems not only for studying processes of evolutionary divergence and speciation, but also for the detection of evolving genetic barriers for specialized parasites.

Puumala Virus Variants Circulating in Forests of Ardennes, France: Ten Years of Genetic Evolution

In Europe, Puumala virus (PUUV) transmitted by the bank vole (Myodes glareolus) is the causative agent of nephropathia epidemica (NE), a mild form of haemorrhagic fever with renal syndrome. In France, very little is known about the spatial and temporal variability of the virus circulating within bank vole populations. The present study involved monitoring of bank vole population dynamics and PUUV microdiversity over a ten-year period (2000–2009) in two forests of the Ardennes region: Elan and Croix-Scaille. Ardennes region is characterised by different environmental conditions associated with different NE epidemiology. Bank vole density and population parameters were estimated using the capture/marking/recapture method, and blood samples were collected to monitor the overall seroprevalence of PUUV in rodent populations. Phylogenetic analyses of fifty-five sequences were performed to illustrate the genetic diversity of PUUV variants between forests. The pattern of the two forests differed clearly. In the Elan forest, the rodent survival was higher, and this limited turn-over resulted in a lower seroprevalence and diversity of PUUV sequences than in the Croix-Scaille forest. Uncovering the links between host dynamics and virus microevolution is improving our understanding of PUUV distribution in rodents and the NE risk.

Impact of Landscape on Host-Parasite Genetic Diversity and Distribution Using the Puumala orthohantavirus-Bank Vole System

In nature, host specificity has a strong impact on the parasite's distribution, prevalence, and genetic diversity. The host's population dynamics is expected to shape the distribution of host-specific parasites. In turn, the parasite's genetic structure is predicted to mirror that of the host. Here, we study the tandem Puumala orthohantavirus (PUUV)-bank vole system. The genetic diversity of 310 bank voles and 33 PUUV isolates from 10 characterized localities of Northeast France was assessed. Our findings show that the genetic diversity of both PUUV and voles, was positively correlated with forest coverage and contiguity of habitats. While the genetic diversity of voles was weakly structured in space, that of PUUV was found to be strongly structured, suggesting that the dispersion of voles was not sufficient to ensure a broad PUUV dissemination. Genetic diversity of PUUV was mainly shaped by purifying selection. Genetic drift and extinction events were better reflected than local adaptation of PUUV. These contrasting patterns of microevolution have important consequences for the understanding of PUUV distribution and epidemiology.

Ecological drivers of genetic connectivity for African malaria vectors Anopheles gambiae and An. arabiensis

Anopheles gambiae and An. arabiensis are major malaria vectors in sub-Saharan Africa. Knowledge of how geographical factors drive the dispersal and gene flow of malaria vectors can help in combatting insecticide resistance spread and planning new vector control interventions. Here, we used a landscape genetics approach to investigate population relatedness and genetic connectivity of An. gambiae and An. arabiensis across Kenya and determined the changes in mosquito population genetic diversity after 20 years of intensive malaria control efforts. We found a significant reduction in genetic diversity in An. gambiae, but not in An. arabiensis as compared to prior to the 20-year period in western Kenya. Significant population structure among populations was found for both species. The most important ecological driver for dispersal and gene flow of An. gambiae and An. arabiensis was tree cover and cropland, respectively. These findings highlight that human induced environmental modifications may enhance genetic connectivity of malaria vectors.

How Bank Vole-PUUV Interactions Influence the Eco-Evolutionary Processes Driving Nephropathia Epidemica Epidemiology-An Experimental and Genomic Approach

In Europe, Puumala virus (PUUV) is responsible for nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS). Despite the presence of its reservoir, the bank vole, on most of French territory, the geographic distribution of NE cases is heterogeneous and NE endemic and non-endemic areas have been reported. In this study we analyzed whether bank vole-PUUV interactions could partly shape these epidemiological differences. We performed crossed-experimental infections using wild bank voles from French endemic (Ardennes) and non-endemic (Loiret) areas and two French PUUV strains isolated from these areas. The serological response and dynamics of PUUV infection were compared between the four cross-infection combinations. Due to logistical constraints, this study was based on a small number of animals. Based on this experimental design, we saw a stronger serological response and presence of PUUV in excretory organs (bladder) in bank voles infected with the PUUV endemic strain. Moreover, the within-host viral diversity in excretory organs seemed to be higher than in other non-excretory organs for the NE endemic cross-infection but not for the NE non-endemic cross-infection. Despite the small number of rodents included, our results showed that genetically different PUUV strains and in a lesser extent their interaction with sympatric bank voles, could affect virus replication and diversity. This could impact PUUV excretion/transmission between rodents and to humans and in turn at least partly shape NE epidemiology in France.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

Population genetics of fruit bat reservoir informs the dynamics, distribution and diversity of Nipah virus

The structure and connectivity of wildlife host populations may influence zoonotic disease dynamics, evolution and therefore spillover risk to people. Fruit bats in the genus Pteropus, or flying foxes, are the primary natural reservoir for henipaviruses-a group of emerging paramyxoviruses that threaten livestock and public health. In Bangladesh, Pteropus medius is the reservoir for Nipah virus-and viral spillover has led to human fatalities nearly every year since 2001. Here, we use mitochondrial DNA and nuclear microsatellite markers to measure the population structure, demographic history and phylogeography of P. medius in Bangladesh. We combine this with a phylogeographic analysis of all known Nipah virus sequences and strains currently available to better inform the dynamics, distribution and evolutionary history of Nipah virus. We show that P. medius is primarily panmictic, but combined analysis of microsatellite and morphological data shows evidence for differentiation of two populations in eastern Bangladesh, corresponding to a divergent strain of Nipah virus also found in bats from eastern Bangladesh. Our demographic analyses indicate that a large, expanding population of flying foxes has existed in Bangladesh since the Late Pleistocene, coinciding with human population expansion in South Asia, suggesting repeated historical spillover of Nipah virus likely occurred. We present the first evidence of mitochondrial introgression, or hybridization, between P. medius and flying fox species found in South-East Asia (P. vampyrus and P. hypomelanus), which may help to explain the distribution of Nipah virus strains across the region.

Density dependence and persistence of Morogoro arenavirus transmission in a fluctuating population of its reservoir host

A key aim in wildlife disease ecology is to understand how host and parasite characteristics influence parasite transmission and persistence. Variation in host population density can have strong impacts on transmission and outbreaks, and theory predicts particular transmission-density patterns depending on how parasites are transmitted between individuals. Here, we present the results of a study on the dynamics of Morogoro arenavirus in a population of multimammate mice (Mastomys natalensis). This widespread African rodent, which is also the reservoir host of Lassa arenavirus in West Africa, is known for its strong seasonal density fluctuations driven by food availability. We investigated to what degree virus transmission changes with host population density and how the virus might be able to persist during periods of low host density. A seven-year capture-mark-recapture study was conducted in Tanzania where rodents were trapped monthly and screened for the presence of antibodies against Morogoro virus. Observed seasonal seroprevalence patterns were compared with those generated by mathematical transmission models to test different hypotheses regarding the degree of density dependence and the role of chronically infected individuals. We observed that Morogoro virus seroprevalence correlates positively with host density with a lag of 1-4 months. Model results suggest that the observed seasonal seroprevalence dynamics can be best explained by a combination of vertical and horizontal transmission and that a small number of animals need to be infected chronically to ensure viral persistence. Transmission dynamics and viral persistence were best explained by the existence of both acutely and chronically infected individuals and by seasonally changing transmission rates. Due to the presence of chronically infected rodents, rodent control is unlikely to be a feasible approach for eliminating arenaviruses such as Lassa virus from Mastomys populations.

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.

Spatial dynamics of a zoonotic orthohantavirus disease through heterogenous data on rodents, rodent infections, and human disease

Zoonotic diseases are challenging to study from the ecological point of view as, broadly speaking, datasets tend to be either detailed on a small spatial extent, or coarse on a large spatial extent. Also, there are many ways to assess zoonotic disease transmission systems, from pathogens to hosts to humans. We explore the complementarity of datasets considering the pathogen in its host, the host and human cases in the context of Puumala orthohantavirus infection in Germany. We selected relevant environmental predictors using a conceptual framework based on resource-based habitats. This framework assesses the functions, and associated environmental resources of the pathogen and associated host. A resource-based habitat framework supports variable selection and result interpretation. Multiplying ‘keyholes’ to view a zoonotic disease transmission system is valuable, but requires a strong conceptual framework to select and interpret environmental explanatory variables. This study highlights the usefulness of a structured, ecology-based approach to study drivers of zoonotic diseases at the level of virus, host, and human - not only for PUUV but also for other zoonotic pathogens. Our results show that human disease cases are best explained by a combination of variables related to zoonotic pathogen circulation and human exposure.

Regional and local patterns of genetic variation and structure in yellow-necked mice - the roles of geographic distance, population abundance, and winter severity

The goal of this study, conducted in seven large woodlands and three areas with small woodlots in northeastern Poland in 2004-2008, was to infer genetic structure in yellow-necked mouse Apodemus flavicollis population and to evaluate the roles of environmental and population ecology variables in shaping the spatial pattern of genetic variation using 768 samples genotyped at 13 microsatellite loci. Genetic variation was very high in all studied regions. The primal genetic subdivision was observed between the northern and the southern parts of the study area, which harbored two major clusters and the intermediate area of highly admixed individuals. The probability of assignment of individual mice to the northern cluster increased significantly with lower temperatures of January and July and declined in regions with higher proportion of deciduous and mixed forests. Despite the detected structure, genetic differentiation among regions was very low. Fine-scale structure was shaped by the population density, whereas higher level structure was mainly shaped by geographic distance. Genetic similarity indices were highly influenced by mouse abundance (which positively correlated with the share of deciduous forests in the studied regions) and exhibited the greatest change between 0 and 1 km in the forests, 0 and 5 km in small woodlots. Isolation by distance pattern, calculated among regions, was highly significant but such relationship between genetic and geographic distance was much weaker, and held the linearity at very fine scale (~1.5 km), when analyses were conducted at individual level.

High Throughput Sequencing and Network Analysis Disentangle the Microbial Communities of Ticks and Hosts Within and Between Ecosystems

We aimed to develop a framework, based on graph theory, to capture the ecological meaning behind pure pair comparisons of microbiome-derived data. As a proof of concept, we applied the framework to analyze the co-occurrence of bacteria in either Ixodes ricinus ticks or the spleen of one of their main hosts, the vole Myodes glareolus. As a secondary lymphoid organ, the spleen acts as a filter of blood and represents well the exposure to microorganisms circulating in the blood; including those acquired and transmitted by ticks during feeding. The microbiome of 301 and 269 individual tick and vole samples, respectively, were analyzed using next generation sequencing (NGS) of 16S rRNA. To assess the effect of habitat on ecological communities of bacteria associated to ticks and voles, two different biotopes were included in the study, forest, and ecotone. An innovative approach of NGS data analysis combining network analysis and phylogenies of co-occuring of bacteria was used to study associations between bacteria in individual samples. Of the 126 bacterial genera found in ticks and voles, 62% were shared by both species. Communities of co-occurring bacteria were always more phylogenetically diverse in ticks than in voles. Interestingly, ~80% of bacterial phylogenetic diversity was found in ~20% of ticks. This pattern was not observed in vole-associated bacteria. Results revealed that the microbiome of I. ricinus is only slightly related to that of M. glareolus and that the biotope plays the most important role in shaping the bacterial communities of either ticks or voles. The analysis of the phylogenetic signal of the network indexes across the 16S rRNA-derived tree of bacteria suggests that the microbiome of both ticks and voles has high phylogenetic diversity and that closest bacterial genera do not co-occur. This study shows that network analysis is a promising tool to unravel complex microbial communities associated to arthropod vectors and vertebrate hosts.

High Throughput Sequencing and Network Analysis Disentangle the Microbial Communities of Ticks and Hosts Within and Between Ecosystems

We aimed to develop a framework, based on graph theory, to capture the ecological meaning behind pure pair comparisons of microbiome-derived data. As a proof of concept, we applied the framework to analyze the co-occurrence of bacteria in either Ixodes ricinus ticks or the spleen of one of their main hosts, the vole Myodes glareolus. As a secondary lymphoid organ, the spleen acts as a filter of blood and represents well the exposure to microorganisms circulating in the blood; including those acquired and transmitted by ticks during feeding. The microbiome of 301 and 269 individual tick and vole samples, respectively, were analyzed using next generation sequencing (NGS) of 16S rRNA. To assess the effect of habitat on ecological communities of bacteria associated to ticks and voles, two different biotopes were included in the study, forest, and ecotone. An innovative approach of NGS data analysis combining network analysis and phylogenies of co-occuring of bacteria was used to study associations between bacteria in individual samples. Of the 126 bacterial genera found in ticks and voles, 62% were shared by both species. Communities of co-occurring bacteria were always more phylogenetically diverse in ticks than in voles. Interestingly, ~80% of bacterial phylogenetic diversity was found in ~20% of ticks. This pattern was not observed in vole-associated bacteria. Results revealed that the microbiome of I. ricinus is only slightly related to that of M. glareolus and that the biotope plays the most important role in shaping the bacterial communities of either ticks or voles. The analysis of the phylogenetic signal of the network indexes across the 16S rRNA-derived tree of bacteria suggests that the microbiome of both ticks and voles has high phylogenetic diversity and that closest bacterial genera do not co-occur. This study shows that network analysis is a promising tool to unravel complex microbial communities associated to arthropod vectors and vertebrate hosts.

Self-clearance of Pestivirus in a Pyrenean Chamois ( Rupicapra pyrenaica) Population

Understanding the dynamics of host-pathogen interaction is key to the management of epidemics. A pestivirus belonging to the border disease virus group 4 emerged around 2001 in Pyrenean chamois ( Rupicapra pyrenaica) in Spain and France. The virus had significant demographic impact in some populations, but it was less harmful and more endemic in other places. The determinants of these local variations are still unclear. Here, we documented empirical evidence of self-clearance of the virus in a chamois population in France. This population has regularly been counted, and chamois were trapped and harvested each year, providing unique demographic and epidemiologic surveys of the population since 1984 and 1994, respectively. The virus was detected using direct (PCR) and indirect (antibody) testing. We showed that virus transmission declined in 2011-12 and likely ceased in 2013, leading to a decline in antibody prevalence since 2014. Self-clearance may be due to limited exchanges with other populations, decrease in population size after an epizootic, and herd immunity. The age structure of captured animals shifted to younger age classes after virus self-clearance, suggesting a return to a colonizing population structure. The possible consequences of virus re-entry are discussed. This observation suggests that pestivirus dynamics occurs at the scale of the metapopulation of Pyrenean chamois. Local self-clearance and re-emergence may help explain the variation of virus dynamics at the local scale.

Population genetics, community of parasites, and resistance to rodenticides in an urban brown rat (Rattus norvegicus) population

Brown rats are one of the most widespread urban species worldwide. Despite the nuisances they induce and their potential role as a zoonotic reservoir, knowledge on urban rat populations remains scarce. The main purpose of this study was to characterize an urban brown rat population from Chanteraines park (Hauts-de-Seine, France), with regards to haematology, population genetics, immunogenic diversity, resistance to anticoagulant rodenticides, and community of parasites. Haematological parameters were measured. Population genetics was investigated using 13 unlinked microsatellite loci. Immunogenic diversity was assessed for Mhc-Drb. Frequency of the Y139F mutation (conferring resistance to rodenticides) and two linked microsatellites were studied, concurrently with the presence of anticoagulant residues in the liver. Combination of microscopy and molecular methods were used to investigate the occurrence of 25 parasites. Statistical approaches were used to explore multiple parasite relationships and model parasite occurrence. Eighty-six rats were caught. The first haematological data for a wild urban R. norvegicus population was reported. Genetic results suggested high genetic diversity and connectivity between Chanteraines rats and surrounding population(s). We found a high prevalence (55.8%) of the mutation Y139F and presence of rodenticide residues in 47.7% of the sampled individuals. The parasite species richness was high (16). Seven potential zoonotic pathogens were identified, together with a surprisingly high diversity of Leptospira species (4). Chanteraines rat population is not closed, allowing gene flow and making eradication programs challenging, particularly because rodenticide resistance is highly prevalent. Parasitological results showed that co-infection is more a rule than an exception. Furthermore, the presence of several potential zoonotic pathogens, of which four Leptospira species, in this urban rat population raised its role in the maintenance and spread of these pathogens. Our findings should stimulate future discussions about the development of a long-term rat-control management program in Chanteraines urban park.

Co-infection of Ticks: The Rule Rather Than the Exception

Introduction

Ticks are the most common arthropod vectors of both human and animal diseases in Europe, and the Ixodes ricinus tick species is able to transmit a large number of bacteria, viruses and parasites. Ticks may also be co-infected with several pathogens, with a subsequent high likelihood of co-transmission to humans or animals. However few data exist regarding co-infection prevalences, and these studies only focus on certain well-known pathogens. In addition to pathogens, ticks also carry symbionts that may play important roles in tick biology, and could interfere with pathogen maintenance and transmission. In this study we evaluated the prevalence of 38 pathogens and four symbionts and their co-infection levels as well as possible interactions between pathogens, or between pathogens and symbionts.

Methodology/principal findings

A total of 267 Ixodes ricinus female specimens were collected in the French Ardennes and analyzed by high-throughput real-time PCR for the presence of 37 pathogens (bacteria and parasites), by rRT-PCR to detect the presence of Tick-Borne encephalitis virus (TBEV) and by nested PCR to detect four symbionts. Possible multipartite interactions between pathogens, or between pathogens and symbionts were statistically evaluated. Among the infected ticks, 45% were co-infected, and carried up to five different pathogens. When adding symbiont prevalences, all ticks were infected by at least one microorganism, and up to eight microorganisms were identified in the same tick. When considering possible interactions between pathogens, the results suggested a strong association between Borrelia garinii and B. afzelii, whereas there were no significant interactions between symbionts and pathogens.

Conclusion/significance

Our study reveals high pathogen co-infection rates in ticks, raising questions about possible co-transmission of these agents to humans or animals, and their consequences to human and animal health. We also demonstrated high prevalence rates of symbionts co-existing with pathogens, opening new avenues of enquiry regarding their effects on pathogen transmission and vector competence.

Ticks transmit more pathogens than any other arthropod, and one single species can transmit a large variety of bacteria and parasites. Because co-infection might be much more common than previously thought, we evaluated the prevalence of 38 known or neglected tick-borne pathogens in Ixodes ricinus ticks. Our results demonstrated that co-infection occurred in almost half of the infected ticks, and that ticks could be infected with up to five pathogens. Moreover, as it is well established that symbionts can affect pathogen transmission in arthropods, we also evaluated the prevalence of four symbiont species and demonstrated that all ticks were infected by at least one microorganism. This work highlights the co-infection phenomenon in ticks, which may have important implications for human and animal health, emphasizing the need for new diagnostic tests better adapted to tick-borne diseases. Finally, the high co-occurrence of symbionts and pathogens in ticks, reveals the necessity to also account for these interactions in the development of new alternative strategies to control ticks and tick-borne disease.

A Comparison between Transcriptome Sequencing and 16S Metagenomics for Detection of Bacterial Pathogens in Wildlife

Background

Rodents are major reservoirs of pathogens responsible for numerous zoonotic diseases in humans and livestock. Assessing their microbial diversity at both the individual and population level is crucial for monitoring endemic infections and revealing microbial association patterns within reservoirs. Recently, NGS approaches have been employed to characterize microbial communities of different ecosystems. Yet, their relative efficacy has not been assessed. Here, we compared two NGS approaches, RNA-Sequencing (RNA-Seq) and 16S-metagenomics, assessing their ability to survey neglected zoonotic bacteria in rodent populations.

Methodology/Principal Findings

We first extracted nucleic acids from the spleens of 190 voles collected in France. RNA extracts were pooled, randomly retro-transcribed, then RNA-Seq was performed using HiSeq. Assembled bacterial sequences were assigned to the closest taxon registered in GenBank. DNA extracts were analyzed via a 16S-metagenomics approach using two sequencers: the 454 GS-FLX and the MiSeq. The V4 region of the gene coding for 16S rRNA was amplified for each sample using barcoded universal primers. Amplicons were multiplexed and processed on the distinct sequencers. The resulting datasets were de-multiplexed, and each read was processed through a pipeline to be taxonomically classified using the Ribosomal Database Project. Altogether, 45 pathogenic bacterial genera were detected. The bacteria identified by RNA-Seq were comparable to those detected by 16S-metagenomics approach processed with MiSeq (16S-MiSeq). In contrast, 21 of these pathogens went unnoticed when the 16S-metagenomics approach was processed via 454-pyrosequencing (16S-454). In addition, the 16S-metagenomics approaches revealed a high level of coinfection in bank voles.

Conclusions/Significance

We concluded that RNA-Seq and 16S-MiSeq are equally sensitive in detecting bacteria. Although only the 16S-MiSeq method enabled identification of bacteria in each individual reservoir, with subsequent derivation of bacterial prevalence in host populations, and generation of intra-reservoir patterns of bacterial interactions. Lastly, the number of bacterial reads obtained with the 16S-MiSeq could be a good proxy for bacterial prevalence.

The majority of human pathogens are of animal origin, i.e. zoonoses; both domestic and wild animals act as host reservoirs. Epidemiological surveys of wildlife may help to predict, prevent and control putative episodes of emerging zoonoses. Microbial diversity and their interactions at both the individual and population level may influence epidemiological infections. Developing generic approaches able to simultaneously detect multiple pathogens without any a priori information becomes essential. Here, we assess the relative efficacy of distinct next-generation sequencing (NGS) approaches to survey neglected zoonotic bacteria in rodent populations: RNA-sequencing (RNA-Seq) and 16S-metagenomics, with the latter resolved via two sequencing techniques, 454-pyrosequencing and MiSeq. The resulting data generated a thorough inventory of zoonotic bacteria in the rodent sample without any previous knowledge of their presence. We concluded that RNA-Seq and 16S-MiSeq are equally sensitive in bacterial genus detection. Nevertheless, only the 16S approach was able to determine bacterial diversity in each individual, which then permitted the derivation of bacterial prevalence and interaction patterns within host populations. We are persuaded that NGS techniques are very affordable candidates and could become routine approaches in future large-scale epidemiological studies.

Spatiotemporal dynamics of Puumala hantavirus associated with its rodent host, Myodes glareolus

Many viruses significantly impact human and animal health. Understanding the population dynamics of these viruses and their hosts can provide important insights for epidemiology and virus evolution. Puumala virus (PUUV) is a European hantavirus that may cause regional outbreaks of hemorrhagic fever with renal syndrome in humans. Here, we analyzed the spatiotemporal dynamics of PUUV circulating in local populations of its rodent reservoir host, the bank vole (Myodes glareolus) during eight years. Phylogenetic and population genetic analyses of all three genome segments of PUUV showed strong geographical structuring at a very local scale. There was a high temporal turnover of virus strains in the local bank vole populations, but several virus strains persisted through multiple years. Phylodynamic analyses showed no significant changes in the local effective population sizes of PUUV, although vole numbers and virus prevalence fluctuated widely. Microsatellite data demonstrated also a temporally persisting subdivision between local vole populations, but these groups did not correspond to the subdivision in the virus strains. We conclude that restricted transmission between vole populations and genetic drift play important roles in shaping the genetic structure and temporal dynamics of PUUV in its natural host which has several implications for zoonotic risks of the human population.

Landscape and regional environmental analysis of the spatial distribution of hantavirus human cases in europe

Background: In Europe, the most prevalent hantavirus, Puumala virus, is transmitted by bank voles and causes nephropathia epidemica in human. The European spatial distribution of nephropathia epidemica is investigated here for the first time with a rich set of environmental variables.

Methods: The influence of variables at the landscape and regional level is studied through multilevel logistic regression, and further information on their effects across the different European ecoregions is obtained by comparing an overall niche model (boosted regression trees) with regressions by ecoregion.

Results: The presence of nephropathia epidemica is likely in populated regions with well-connected forests, more intense vegetation activity, low soil water content, mild summers, and cold winters. In these regions, landscapes with a higher proportion of built-up areas in forest ecotones and lower minimum temperature in winter are expected to be more at risk. Climate and forest connectivity have a stronger effect at the regional level. If variables are staying at their current values, the models predict that nephropathia epidemica may know intensification but should not spread (although southern Sweden, the Norwegian coast, and the Netherlands should be kept under watch).

Conclusion: Models indicate that large-scale modeling can lead to a very high predictive power. At large scale, the effect of one variable on disease may follow three response scenarios: the effect may be the same across the entire study area, the effect can change according to the variable value, and the effect can change depending on local specificities. Each of these scenarios impacts large-scale modeling differently.

Identification of factors influencing the Puumala virus seroprevalence within its reservoir in aMontane Forest Environment

Puumala virus (PUUV) is a major cause of mild to moderate haemorrhagic fever with renal syndrome and is transmitted by the bank vole (Myodes glareolus). There has been a high cumulative incidence of recorded human cases in South-eastern Germany since 2004 when the region was first recognized as being endemic for PUUV. As the area is well known for outdoor recreation and the Bavarian Forest National Park (BFNP) is located in the region, the increasing numbers of recorded cases are of concern. To understand the population and environmental effects on the seroprevalence of PUUV in bank voles we trapped small mammals at 23 sites along an elevation gradient from 317 to 1420m above sea level. Generalized linear mixed effects models(GLMEM) were used to explore associations between the seroprevalence of PUUV in bank voles and climate and biotic factors. We found that the seroprevalence of PUUV was low (6%–7%) in 2008 and 2009, and reached 29% in 2010. PUUV seroprevalence was positively associated with the local species diversity and deadwood layer, and negatively associated with mean annual temperature, mean annual solar radiation, and herb layer. Based on these findings, an illustrative risk map for PUUV seroprevalence prediction in bank voles was created for an area of the national park. The map will help when planning infrastructure in the national park (e.g., huts, shelters, and trails).

Epidemiology of leptospira transmitted by rodents in southeast Asia

Background

Leptospirosis is the most common bacterial zoonoses and has been identified as an important emerging global public health problem in Southeast Asia. Rodents are important reservoirs for human leptospirosis, but epidemiological data is lacking.

Methodology/Principal Findings

We sampled rodents living in different habitats from seven localities distributed across Southeast Asia (Thailand, Lao PDR and Cambodia), between 2009 to 2010. Human isolates were also obtained from localities close to where rodents were sampled. The prevalence of Leptospira infection was assessed by real-time PCR using DNA extracted from rodent kidneys, targeting the lipL32 gene. Sequencing rrs and secY genes, and Multi Locus Variable-number Tandem Repeat (VNTR) analyses were performed on DNA extracted from rat kidneys for Leptospira isolates molecular typing. Four species were detected in rodents, L. borgpetersenii (56% of positive samples), L. interrogans (36%), L. kirschneri (3%) and L. weilli (2%), which were identical to human isolates. Mean prevalence in rodents was approximately 7%, and largely varied across localities and habitats, but not between rodent species. The two most abundant Leptospira species displayed different habitat requirements: L. interrogans was linked to humid habitats (rice fields and forests) while L. borgpetersenii was abundant in both humid and dry habitats (non-floodable lands).

Conclusion/Significance

L. interrogans and L. borgpetersenii species are widely distributed amongst rodent populations, and strain typing confirmed rodents as reservoirs for human leptospirosis. Differences in habitat requirements for L. interrogans and L. borgpetersenii supported differential transmission modes. In Southeast Asia, human infection risk is not only restricted to activities taking place in wetlands and rice fields as is commonly accepted, but should also include tasks such as forestry work, as well as the hunting and preparation of rodents for consumption, which deserve more attention in future epidemiological studies.

Leptospirosis is the most prevalent bacterial zoonosis worldwide. Rodents are believed to be the main reservoirs of Leptospira, yet little epidemiological research has been conducted on rodents from Southeast Asia. Previous studies suggest that activities which place humans in microenvironments shared by rodents increase the probability of contracting leptospirosis. We therefore investigated the circulation of leptospiral species and strains in rodent communities and human populations in seven localities scattered throughout Southeast Asia; in Thailand, Lao PDR and Cambodia. Molecular typing assays were used to characterize leptospiral species and strains in both rodents and humans, which demonstrated common strains between humans and rodents. Additionally, we observed that the two most abundant leptospiral species; L. borgpetersenii and L. interrogans, have different habitat requirements, which supposes different modes of transmission. Lastly, in Southeast Asia, the risk of leptospiral transmission to humans is not solely limited to wetlands and rice paddy fields, but is also linked to forested areas, and activities such as the hunting and/or preparation of rodents for consumption.

Immunogenetic factors affecting susceptibility of humans and rodents to hantaviruses and the clinical course of hantaviral disease in humans

We reviewed the associations of immunity-related genes with susceptibility of humans and rodents to hantaviruses, and with severity of hantaviral diseases in humans. Several class I and class II HLA haplotypes were linked with severe or benign hantavirus infections, and these haplotypes varied among localities and hantaviruses. The polymorphism of other immunity-related genes including the C4A gene and a high-producing genotype of TNF gene associated with severe PUUV infection. Additional genes that may contribute to disease or to PUUV infection severity include non-carriage of the interleukin-1 receptor antagonist (IL-1RA) allele 2 and IL-1β (-511) allele 2, polymorphisms of plasminogen activator inhibitor (PAI-1) and platelet GP1a. In addition, immunogenetic studies have been conducted to identify mechanisms that could be linked with the persistence/clearance of hantaviruses in reservoirs. Persistence was associated during experimental infections with an upregulation of anti-inflammatory responses. Using natural rodent population samples, polymorphisms and/or expression levels of several genes have been analyzed. These genes were selected based on the literature of rodent or human/hantavirus interactions (some Mhc class II genes, Tnf promoter, and genes encoding the proteins TLR4, TLR7, Mx2 and β3 integrin). The comparison of genetic differentiation estimated between bank vole populations sampled over Europe, at neutral and candidate genes, has allowed to evidence signatures of selection for Tnf, Mx2 and the Drb Mhc class II genes. Altogether, these results corroborated the hypothesis of an evolution of tolerance strategies in rodents. We finally discuss the importance of these results from the medical and epidemiological perspectives.

Genetic characterization of the human relapsing fever spirochete Borrelia miyamotoi in vectors and animal reservoirs of Lyme disease spirochetes in France

Background

In France as elsewhere in Europe the most prevalent TBD in humans is Lyme borreliosis, caused by different bacterial species belonging to Borrelia burgdorferi sensu lato complex and transmitted by the most important tick species in France, Ixodes ricinus. However, the diagnosis of Lyme disease is not always confirmed and unexplained syndromes occurring after tick bites have become an important issue. Recently, B. miyamotoi belonging to the relapsing fever group and transmitted by the same Ixodes species has been involved in human disease in Russia, the USA and the Netherlands. In the present study, we investigate the presence of B. miyamotoi along with other Lyme Borreliosis spirochetes, in ticks and possible animal reservoirs collected in France.

Methods

We analyzed 268 ticks (Ixodes ricinus) and 72 bank voles (Myodes glareolus) collected and trapped in France for the presence of DNA from B. miyamotoi as well as from Lyme spirochetes using q-PCR and specific primers and probes. We then compared the French genotypes with those found in other European countries.

Results

We found that 3% of ticks and 5.55% of bank voles were found infected by the same B. miyamotoi genotype, while co-infection with other Lyme spirochetes (B. garinii) was identified in 12% of B. miyamotoi infected ticks. Sequencing showed that ticks and rodents carried the same genotype as those recently characterized in a sick person in the Netherlands.

Conclusions

The genotype of B. miyamotoi circulating in ticks and bank voles in France is identical to those already described in ticks from Western Europe and to the genotype isolated from a sick person in The Netherlands. This results suggests that even though no human cases have been reported in France, surveillance has to be improved. Moreover, we showed that ticks could simultaneously carry B. miyamotoi and Lyme disease spirochetes, increasing the problem of co-infection in humans.

Potential corridors and barriers for plague spread in Central Asia

Background

Plague (Yersinia pestis infection) is a vector-borne disease which caused millions of human deaths in the Middle Ages. The hosts of plague are mostly rodents, and the disease is spread by the fleas that feed on them. Currently, the disease still circulates amongst sylvatic rodent populations all over the world, including great gerbil (Rhombomys opimus) populations in Central Asia. Great gerbils are social desert rodents that live in family groups in burrows, which are visible on satellite images. In great gerbil populations an abundance threshold exists, above which plague can spread causing epizootics. The spatial distribution of the host species is thought to influence the plague dynamics, such as the direction of plague spread, however no detailed analysis exists on the possible functional or structural corridors and barriers that are present in this population and landscape. This study aims to fill that gap.

Methods

Three 20 by 20 km areas with known great gerbil burrow distributions were used to analyse the spatial distribution of the burrows. Object-based image analysis was used to map the landscape at several scales, and was linked to the burrow maps. A novel object-based method was developed – the mean neighbour absolute burrow density difference (MNABDD) – to identify the optimal scale and evaluate the efficacy of using landscape objects as opposed to square cells. Multiple regression using raster maps was used to identify the landscape-ecological variables that explain burrow density best. Functional corridors and barriers were mapped using burrow density thresholds. Cumulative resistance of the burrow distribution to potential disease spread was evaluated using cost distance analysis. A 46-year plague surveillance dataset was used to evaluate whether plague spread was radially symmetric.

Results

The burrow distribution was found to be non-random and negatively correlated with Greenness, especially in the floodplain areas. Corridors and barriers showed a mostly NWSE alignment, suggesting easier spreading along this axis. This was confirmed by the analysis of the plague data.

Conclusions

Plague spread had a predominantly NWSE direction, which is likely due to the NWSE alignment of corridors and barriers in the burrow distribution and the landscape. This finding may improve predictions of plague in the future and emphasizes the importance of including landscape analysis in wildlife disease studies.

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.

Candidatus Neoehrlichia mikurensis in bank voles, France

To further assess the geographic occurrence, possible vectors, and prevalence of Candidatus Neoehrlichia mikurensis, we analyzed spleen tissues from 276 voles trapped close to human settlements in France; 5 were infected with the organism. Sequencing showed the isolates carried the same genotype as the bacteria that caused disease in humans and animals elsewhere in Europe.

The diversity of microparasites of rodents: a comparative analysis that helps in identifying rodent-borne rich habitats in Southeast Asia

BACKGROUND

Predicting habitats prone to favor disease transmission is challenging due to confounding information on habitats, reservoirs, and diseases. Comparative analysis, which aims at investigating ecological and evolutionary patterns among species, is a tool that may help. The emergence of zoonotic pathogens is a major health concern and is closely linked to habitat modifications by human activities. Risk assessment requires a better knowledge of the interactions between hosts, parasites, and the landscape.

METHODS

We used information from a field spatial study that investigated the distribution of murid rodents, in various habitats of three countries in Southeast Asia, in combination with their status of infection by 10 taxa of microparasites obtained from the literature. Microparasite species richness was calculated by rodent species on 20,272 rodents of 13 species. Regression tree models and generalized linear models were used to explain microparasite diversity by the average distance between the trapping site and five categories of land cover: forest, steep agriculture land, flat agriculture land, water, and built-up surfaces. Another variable taken into account was the slope.

RESULTS

We found that microparasite diversity was positively associated with flat agriculture land, in this context mainly rice fields, and negatively associated with slope. Microparasite diversity decreased sharply a 100 m or less from flat agriculture land.

CONCLUSION

We conclude that there is high microparasite circulation in rodents of flooded farmlands, meaning possibly a higher risk of disease for human inhabitants.

In Search for Factors that Drive Hantavirus Epidemics

In Europe, hantaviruses (Bunyaviridae) are small mammal-associated zoonotic and emerging pathogens that can cause hemorrhagic fever with renal syndrome (HFRS). Puumala virus, the main etiological agent carried by the bank vole Myodes glareolus is responsible for a mild form of HFRS while Dobrava virus induces less frequent but more severe cases of HFRS. Since 2000 in Europe, more than 3000 cases of HFRS have been recorded, in average, each year, which is nearly double compared to the previous decade. In addition to this upside long-term trend, significant oscillations occur. Epidemic years appear, usually every 2-4 years, with an increased incidence, generally in localized hot spots. Moreover, the virus has been identified in new areas in the recent years. A great number of surveys have been carried out in order to assess the prevalence of the infection in the reservoir host and to identify links with different biotic and abiotic factors. The factors that drive the infections are related to the density and diversity of bank vole populations, prevalence of infection in the reservoir host, viral excretion in the environment, survival of the virus outside its host, and human behavior, which affect the main transmission virus route through inhalation of infected rodent excreta. At the scale of a rodent population, the prevalence of the infection increases with the age of the individuals but also other parameters, such as sex and genetic variability, interfere. The contamination of the environment may be correlated to the number of newly infected rodents, which heavily excrete the virus. The interactions between these different parameters add to the complexity of the situation and explain the absence of reliable tools to predict epidemics. In this review, the factors that drive the epidemics of hantaviruses in Middle Europe are discussed through a panorama of the epidemiological situation in Belgium, France, and Germany.

Environmental change and disease dynamics: effects of intensive forest management on Puumala hantavirus infection in boreal bank vole populations

Intensive management of Fennoscandian forests has led to a mosaic of woodlands in different stages of maturity. The main rodent host of the zoonotic Puumala hantavirus (PUUV) is the bank vole (Myodes glareolus), a species that can be found in all woodlands and especially mature forests. We investigated the influence of forest age structure on PUUV infection dynamics in bank voles. Over four years, we trapped small mammals twice a year in a forest network of different succession stages in Northern Finland. Our study sites represented four forest age classes from young (4 to 30 years) to mature (over 100 years) forests. We show that PUUV-infected bank voles occurred commonly in all forest age classes, but peaked in mature forests. The probability of an individual bank vole to be PUUV infected was positively related to concurrent host population density. However, when population density was controlled for, a relatively higher infection rate was observed in voles trapped in younger forests. Furthermore, we found evidence of a "dilution effect" in that the infection probability was negatively associated with the simultaneous density of other small mammals during the breeding season. Our results suggest that younger forests created by intensive management can reduce hantaviral load in the environment, but PUUV is common in woodlands of all ages. As such, the Fennoscandian forest landscape represents a significant reservoir and source of hantaviral infection in humans.