Temporal patterns in Ixodes ricinus microbial communities: an insight into tick-borne microbe interactions.

BACKGROUND: Ticks transmit pathogens of medical and veterinary importance and are an increasing threat to human and animal health. Assessing disease risk and developing new control strategies requires identifying members of the tick-borne microbiota as well as their temporal dynamics and interactions. METHODS: Using high-throughput sequencing, we studied the Ixodes ricinus microbiota and its temporal dynamics. 371 nymphs were monthly collected during three consecutive years in a peri-urban forest. After a Poisson lognormal model was adjusted to our data set, a principal component analysis, sparse network reconstruction, and differential analysis allowed us to assess seasonal and monthly variability of I. ricinus microbiota and interactions within this community. RESULTS: Around 75% of the detected sequences belonged to five genera known to be maternally inherited bacteria in arthropods and to potentially circulate in ticks: Candidatus Midichloria, Rickettsia, Spiroplasma, Arsenophonus and Wolbachia. The structure of the I. ricinus microbiota varied over time with interannual recurrence and seemed to be mainly driven by OTUs commonly found in the environment. Total network analysis revealed a majority of positive partial correlations. We identified strong relationships between OTUs belonging to Wolbachia and Arsenophonus, evidence for the presence of the parasitoid wasp Ixodiphagus hookeri in ticks. Other associations were observed between the tick symbiont Candidatus Midichloria and pathogens belonging to Rickettsia. Finally, more specific network analyses were performed on TBP-infected samples and suggested that the presence of pathogens belonging to the genera Borrelia, Anaplasma and Rickettsia may disrupt microbial interactions in I. ricinus. CONCLUSIONS: We identified the I. ricinus microbiota and documented marked shifts in tick microbiota dynamics over time. Statistically, we showed strong relationships between the presence of specific pathogens and the structure of the I. ricinus microbiota. We detected close links between some tick symbionts and the potential presence of either pathogenic Rickettsia or a parasitoid in ticks. These new findings pave the way for the development of new strategies for the control of ticks and tick-borne diseases. Video abstract.

Consequences of organ choice in describing bacterial pathogen assemblages in a rodent population.

High-throughput sequencing technologies now allow for rapid cost-effective surveys of multiple pathogens in many host species including rodents, but it is currently unclear if the organ chosen for screening influences the number and identity of bacteria detected. We used 16S rRNA amplicon sequencing to identify bacterial pathogens in the heart, liver, lungs, kidneys and spleen of 13 water voles (Arvicola terrestris) collected in Franche-Comté, France. We asked if bacterial pathogen assemblages within organs are similar and if all five organs are necessary to detect all of the bacteria present in an individual animal. We identified 24 bacteria representing 17 genera; average bacterial richness for each organ ranged from 1·5 ± 0·4 (mean ± standard error) to 2·5 ± 0·4 bacteria/organ and did not differ significantly between organs. The average bacterial richness when organ assemblages were pooled within animals was 4·7 ± 0·6 bacteria/animal; Operational Taxonomic Unit accumulation analysis indicates that all five organs are required to obtain this. Organ type influences bacterial assemblage composition in a systematic way (PERMANOVA, 999 permutations, pseudo-F 4,51 = 1·37, P = 0·001). Our results demonstrate that the number of organs sampled influences the ability to detect bacterial pathogens, which can inform sampling decisions in public health and wildlife ecology.

[RATS OF THE GENUS RATTUS AS HOSTS FOR NATURAL FOCAL INFECTIOUS AGENTS].

The paper reviews the significance of rats of the genus Rattus as hosts for zoogenous infections in a genus formation area (Southeast Asia) as compared to the invasion part of the genus area. The rats of the genus Rattus and their related disease agents are shown to be a uniqe model for the formation and development of a host-pathogen system. In the modern period of urbanization growth, the rats are among few species of warm-blooded vectors that can maintain the anthropurgic foci of feral nidal infections in the cities and towns and transmit their pathogens to the urban population. There are all prerequisites for the high activity of these foci in the native area of rats. By having settled, the rats have carried infectious agents outside this area along all continents in historical times. During invasions, the rats have become carriers of many other infections.

Environmental factors and public health policy associated with human and rodent infection by leptospirosis: a land cover-based study in Nan province, Thailand.

Leptospirosis incidence has increased markedly since 1995 in Thailand, with the eastern and northern parts being the most affected regions, particularly during flooding events. Here, we attempt to overview the evolution of human prevalence during the past decade and identify the environmental factors that correlate with the incidence of leptospirosis and the clinical incidence in humans. We used an extensive survey of Leptospira infection in rodents conducted in 2008 and 2009 and the human incidence of the disease from 2003 to 2012 in 168 villages of two districts of Nan province in Northern Thailand. Using an ad-hoc developed land-use cover implemented in a geographical information system we showed that humans and rodents were not infected in the same environment/habitat in the land-use cover. High village prevalence was observed in open habitat near rivers for the whole decade, or in 2008-2009 mostly in rice fields prone to flooding, whereas infected rodents (2008-2009) were observed in patchy habitat with high forest cover, mostly situated on sloping ground areas. We also investigated the potential effects of public health campaigns conducted after the dramatic flood event of 2006. We showed that, before 2006, human incidence in villages was explained by the population size of the village according to the environmental source of infection of this disease, while as a result of the campaigns, human incidence in villages after 2006 appeared independent of their population size. This study confirms the role of the environment and particularly land use, in the transmission of bacteria, emphasized by the effects of the provincial public health campaigns on the epidemiological pattern of incidence, and questions the role of rodents as reservoirs.

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.

Coumatetralyl resistance of Rattus tanezumi infesting oil palm plantations in Indonesia.

Rodent control is an important issue in human health and agriculture. Oil palm plantations are rapidly expanding in Indonesia and this is having a major economic and ecological impact. Rodent control in oil palm plantations is based principally on the use of anti-vitamin K (AVK), the main anticoagulant used being coumatetralyl, a first-generation AVK. We conducted a comparative study in two well established oil palm plantations in Indonesia: (1) one without chemical control in Riau and (2) another with intensive coumatetralyl use on Bangka Island. Rat species were identified by the molecular barcoding method. Susceptibility to coumatetralyl was then assessed within the two populations and we screened for mutations in vkorc1, which encodes the molecular target of AVK. Different species were found in the two areas: Rattus tiomanicus in Riau, and a mix of R. tanezumi and a close relative one in Bangka. The rats in Riau were much more susceptible to coumatetralyl than those in Bangka. This study is the first to demonstrate physiological tolerance to AVK in these species. vkorc1 displayed low levels of polymorphism, and no SNP was associated with the high-tolerance phenotypes of R. tanezumi clade, even those exposed to very high concentrations (32 × the effective dose of 0.36 mg kg(-1)). The biochemical basis of this tolerance remains unknown, but may involve the vkorc1 promoter and/or cytochrome P450 metabolism. We discuss our results and the selective role of anticoagulant use in the occurrence of phenotypic tolerance.

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.

Sex-biased dispersal patterns depend on the spatial scale in a social rodent.

Dispersal is a fundamental process in ecology because it influences the dynamics, genetic structure and persistence of populations. Furthermore, understanding the evolutionary causes of dispersal pattern, particularly when they differ between genders, is still a major question in evolutionary ecology. Using a panel of 10 microsatellite loci, we investigated at different spatial scales the genetic structure and the sex-specific dispersal patterns in the common vole Microtus arvalis, a small colonial mammal. This study was conducted in an intensive agricultural area of western France. Hierarchical F(ST) analyses, relatedness and assignment tests suggested (i) that females are strongly kin-clustered within colonies; (ii) that dispersal is strongly male-biased at a local scale; and (iii) long-distance dispersal is not rare and more balanced between genders. We conclude that males migrate continuously from colony to colony to reproduce, whereas females may disperse just once and would be mainly involved in new colony foundation.

Characterization and PCR multiplexing of polymorphic microsatellite loci in the whipworm Trichuris arvicolae, parasite of arvicoline rodents and their cross-species utilization in T. muris, parasite of murines.

Trichuris arvicolae and T. muris are gastro-intestinal nematodes of respectively arvicoline and murine rodents. We aim to investigate the ecology of these Trichuris species using population genetics. We sampled nematodes from rodents trapped in the East of France. After confirming the species identification of the nematodes using ITS1-5.8S-ITS2 ribosomal DNA sequences, we isolated and characterized twelve dinucleotide microsatellite loci in T. arvicolae. A multiplex panel was developed. Application to a set of 30 individuals allowed clear and easy characterization of allele sizes. The number of alleles ranged from 2 to 6 per locus with observed heterozygosities ranging from 0 to 0.93. A test on eleven T. muris revealed that eight loci among twelve amplified, and five were polymorphic. These sets of microsatellite loci provide high throughput capacity for population genetic studies.

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.

Kinship, dispersal and hantavirus transmission in bank and common voles.

Hantaviruses are among the main emerging infectious agents in Europe. Their mode of transmission in natura is still not well known. In particular, social features and behaviours could be crucial for understanding the persistence and the spread of hantaviruses in rodent populations. Here, we investigated the importance of kinclustering and dispersal in hantavirus transmission by combining a fine-scale spatiotemporal survey (4 km2) and a population genetics approach. Two specific host-hantavirus systems were identified and monitored: the bank vole Myodes, earlier Clethrionomys glareolus--Puumala virus and the common vole Microtus arvalis--Tula virus. Sex, age and landscape characteristics significantly influenced the spatial distribution of infections in voles. The absence of temporal stability in the spatial distributions of viruses suggested that dispersal is likely to play a role in virus propagation. Analysing vole kinship from microsatellite markers, we found that infected voles were more closely related to each other than non-infected ones. Winter kin-clustering, shared colonies within matrilineages or delayed dispersal could explain this pattern. These two last results hold, whatever the host-hantavirus system considered. This supports the roles of relatedness and dispersal as general features for hantavirus transmission.

Spatial genetic structure of a small rodent in a heterogeneous landscape.

Gene flow in natural populations may be strongly influenced by landscape features. The integration of landscape characteristics in population genetic studies may thus improve our understanding of population functioning. In this study, we investigated the population genetic structure and gene flow pattern for the common vole, Microtus arvalis, in a heterogeneous landscape characterised by strong spatial and temporal variation. The studied area is an intensive agricultural zone of approximately 500 km(2) crossed by a motorway. We used individual-based Bayesian methods to define the number of population units and their spatial borders without prior delimitation of such units. Unexpectedly, we determined a single genetic unit that covered the entire area studied. In particular, the motorway considered as a likely barrier to dispersal was not associated with any spatial genetic discontinuity. Using computer simulations, we demonstrated that recent anthropogenic barriers to effective dispersal are difficult to detect through analysis of genetic variation for species with large effective population sizes. We observed a slight, but significant, pattern of isolation by distance over the whole study site. Spatial autocorrelation analyses detected genetic structuring on a local scale, most probably due to the social organisation of the study species. Overall, our analysis suggests intense small-scale dispersal associated with a large effective population size. High dispersal rates may be imposed by the strong spatio-temporal heterogeneity of habitat quality, which characterises intensive agroecosystems.

Dispersal in Mastomys natalensis mice: use of fine-scale genetic analyses for pest management.

Mastomys natalensis is the major pest rodent in sub-Saharan Africa. In this study, population genetic techniques were used to gain new insights into its dispersal behaviour, a critical parameter in pest management. Using 11 microsatellites, 272 individuals from a 300 ha area in Tanzania were genotyped. Genetic diversity was high, with no isolation by distance and little differentiation between field plots far apart, indicating a large effective population size and high dispersal rates in agreement with ecological observations. On the other hand, genetic differentiation between nearby field plots, isolation by distance within a single field plot and kin clustering were also observed. This apparent contradiction may be explained by yearly founder effects of a small number of breeding individuals per square area, which is consistent with the presence of linkage disequilibrium. An alternative, not mutually exclusive explanation is that there are both dispersing and sedentary animals in the population. The low-density field plots were characterized by low relatedness and small genetic distances to other field plots, indicating a high turnover rate and negative density-dependent dispersal. In one field plot female-biased dispersal was observed, which may be related to inbreeding avoidance or female competition for resources. Most juveniles appeared to be local recruits, but they did not seem to stay in their native area for more than two months. Finally, possible implications for pest management are discussed.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles.

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.

Molecular evidence of Pleistocene bidirectional faunal exchange between Europe and the Near East: the case of the bicoloured shrew (Crocidura leucodon, Soricidae).

We sequenced 1077 bp of the mitochondrial cytochrome b gene and 511 bp of the nuclear Apolipoprotein B gene in bicoloured shrew (Crocidura leucodon, Soricidae) populations ranging from France to Georgia. The aims of the study were to identify the main genetic clades within this species and the influence of Pleistocene climatic variations on the respective clades. The mitochondrial analyses revealed a European clade distributed from France eastwards to north-western Turkey and a Near East clade distributed from Georgia to Romania; the two clades separated during the Middle Pleistocene. We clearly identified a population expansion after a bottleneck for the European clade based on mitochondrial and nuclear sequencing data; this expansion was not observed for the eastern clade. We hypothesize that the western population was confined to a small Italo-Balkanic refugium, whereas the eastern population subsisted in several refugia along the southern coast of the Black Sea.

Mediterranean populations of the lesser white-toothed shrew (Crocidura suaveolens group): an unexpected puzzle of Pleistocene survivors and prehistoric introductions.

An earlier study revealed the strong phylogeographical structure of the lesser white-toothed shrew (Crocidura suaveolens group) within the northern Palaearctic. Here, we aim to reconstruct the colonization history of Mediterranean islands and to clarify the biogeography and phylogeographical relationships of the poorly documented Middle East region with the northern Palaearctic. We performed analyses on 998-bp-long haplotypes of the mitochondrial cytochrome b gene of 143 samples collected around the Mediterranean basin, including islands and the Middle East. The analyses suggest that the Cypriot shrew belongs to the rare group of relict insular Pleistocene mammal taxa that have survived to the present day. In contrast, the Cretan, Corsican and Menorcan populations were independently introduced from the Middle East during the Holocene. The phylogeographical structure of this temperate Palaearctic species within the Middle East appears to be complex and rich in diversity, probably reflecting fragmentation of the area by numerous mountain chains. Four deeply divergent clades of the C. suaveolens group occur in the area, meaning that a hypothetical contact zone remains to be located in central western Iran.

The intestinal nematode Trichuris arvicolae affects the fecundity of its host, the common vole Microtus arvalis.

Parasites have detrimental effects on host fitness. Consequently, they play a major role for host population dynamics. In this study, we investigated experimentally the impact of the nematode Trichuris arvicolae on the reproduction of its host, the common vole Microtus arvalis. Wild common voles were trapped in east of France and reared in standardized conditions before being experimentally infected. Infection with Trichuris arvicolae did not affect host consumption of food or water. Parasitized females gave birth to slightly less pups (mean 3.36 +/- 0.38) than unparasitized females (mean 3.60 +/- 0.40). Controlling for natal litter size using analysis of covariance (ANCOVA), T. arvicolae infection had a significant effect on the individual mass at birth, with pups from parasitized females having significantly lower mass (2.11 g +/- 0.01) than pups from unparasitized females (2.20 g +/- 0.01). Other measures of host reproductive outputs (time to first reproduction, mass of pups at weaning, litter survival) were not affected by maternal parasite infection. We discuss how these changes in M. arvalis reproductive investments associated with T. arvicolae infection must now be investigated in the context of physiological trade-offs.

Migration and recovery of the genetic diversity during the increasing density phase in cyclic vole populations.

In cyclic populations, high genetic diversity is currently reported despite the periodic low numbers experienced by the populations during the low phases. Here, we report spatio-temporal monitoring at a very fine scale of cyclic populations of the fossorial water vole (Arvicola terrestris) during the increasing density phase. This phase marks the transition from a patchy structure (demes) during low density to a continuous population in high density. We found that the genetic diversity was effectively high but also that it displayed a local increase within demes over the increasing phase. The genetic diversity remained relatively constant when considering all demes together. The increase in vole abundance was also correlated with a decrease of genetic differentiation among demes. Such results suggest that at the end of the low phase, demes are affected by genetic drift as the result of being small and geographically isolated. This leads to a loss of local genetic diversity and a spatial differentiation among demes. This situation is counterbalanced during the increasing phase by the spatial expansion of demes and the increase of the effective migration among differentiated demes. We provide evidences that in cyclic populations of the fossorial water voles, the relative influence of drift operating during low density populations and migration occurring principally while population size increases interacts closely to maintain high genetic diversity.

Genetic structure is influenced by landscape features: empirical evidence from a roe deer population.

The delimitation of population units is of primary importance in population management and conservation biology. Moreover, when coupled with landscape data, the description of population genetic structure can provide valuable knowledge about the permeability of landscape features, which is often difficult to assess by direct methods (e.g. telemetry). In this study, we investigated the genetic structuring of a roe deer population which recently recolonized a fragmented landscape. We sampled 1148 individuals from a 40 x 55-km area containing several putative barriers to deer movements, and hence to gene flow, namely a highway, rivers and several canals. In order to assess the effect of these landscape features on genetic structure, we implemented a spatial statistical model known as geneland which analyses genetic structure, explicitly taking into account the spatial nature of the problem. Two genetic units were inferred, exhibiting a very low level of differentiation (F(ST) = 0.008). The location of their boundaries suggested that there are no absolute barriers in this study area, but that the combination of several landscape features with low permeability can lead to population differentiation. Our analysis hence suggests that the landscape has a significant influence on the structuring of the population under study. It also illustrates the use of geneland as a powerful method to infer population structure, even in situations of young populations exhibiting low genetic differentiation.

Dispersal is not female biased in a resource-defence mating ungulate, the European roe deer.

Dispersal is frequently more prevalent in one sex compared to the other. Greenwood proposed that patterns of sex-biased dispersal among birds and mammals are linked to their mating strategies. For species where males defend resources rather than females, he predicted female-biased dispersal, because males should remain at their birth site where they are familiar with the distribution of the resources that they must defend. Greenwood's hypothesis has been extensively supported among birds, where most species exhibit a resource-defence mating strategy. However, almost no equivalent information is available for mammals as males generally defend mates in this group. An exception is the European roe deer, a resource-defence mating ungulate. We thus tested Greenwood's hypothesis on this atypical mammalian model, looking for female-biased dispersal using sex-specific inter-individual genetic distances. We conclusively show that gene flow is not higher among females compared to males in the studied roe deer population, and hence that dispersal is not female-biased, suggesting that male mating strategy is not the primary selective force driving the evolution of dispersal in roe deer. We discuss the role of female mate choice and intra-sexual competition as possible alternative selective pressures involved.