The timid invasion: behavioural adjustments and range expansion in a non-native rodent.

Animal behaviour can moderate biological invasion processes, and the native fauna's ability to adapt. The importance and nature of behavioural traits favouring colonization success remain debated. We investigated behavioural responses associated with risk-taking and exploration, both in non-native bank voles (Myodes glareolus, N = 225) accidentally introduced to Ireland a century ago, and in native wood mice (Apodemus sylvaticus, N = 189), that decline in numbers with vole expansion. We repeatedly sampled behavioural responses in three colonization zones: established bank vole populations for greater than 80 years (2 sites), expansion edge vole populations present for 1-4 years (4) and pre-arrival (2). All zones were occupied by wood mice. Individuals of both species varied consistently in risk-taking and exploration. Mice had not adjusted their behaviour to the presence of non-native voles, as it did not differ between the zones. Male voles at the expansion edge were initially more risk-averse but habituated faster to repeated testing, compared to voles in the established population. Results thus indicate spatial sorting for risk-taking propensity along the expansion edge in the dispersing sex. In non-native prey species the ability to develop risk-averse phenotypes may thus represent a fundamental component for range expansions.

Quantitative assessment of spicule length in Heligmosomoides spp. (Nematoda, Heligmosomidae): distinction between H. bakeri, H. polygyrus and H. glareoli.

Nematode spicules play a vital role in the reproductive activity of species that possess them. Our primary objective was to compare the lengths of spicules of the laboratory mouse (Mus musculus) ­ maintained isolate H. bakeri ­ with those of H. polygyrus from naturally infected wood mice (Apodemus sylvaticus). On a more limited scale, we also included H. glareoli from bank voles (Myodes glareolus), a species reputed to possess longer spicules than either of the 2 former species. In total, we measured 1264 spicules (H. bakeri, n = 614; H. polygyrus n = 582; and H. glareoli, n = 68). There was a highly significant difference between the spicule lengths of the Nottingham-maintained H. bakeri (mean = 0.518 mm) and H. polygyrus (0.598 mm) from 11 different localities across the British Isles. A comparison of the spicules of H. bakeri maintained in 4 different laboratories in 3 continents revealed a range in the mean values from 0.518 to 0.540 mm, while those of worms from Australian wild house mice were shorter (0.480 mm). Mean values for H. polygyrus from wood mice from the British Isles ranged from 0.564 to 0.635 mm, although isolates of this species from Norway had longer spicules (0.670 mm). In agreement with the literature, the spicules of H. glareoli were considerably longer (1.098 mm). Since spicules play a vital role in the reproduction of nematode species that possess them, the difference in spicule lengths between H. bakeri and H. polygyrus adds to the growing evidence that these 2 are quite distinct species and likely reproductively isolated.

Long term patterns of association between MHC and helminth burdens in the bank vole support Red Queen dynamics.

Major histocompatibility complex (MHC) genes encode proteins crucial for adaptive immunity of vertebrates. Negative frequency-dependent selection (NFDS), resulting from adaptation of parasites to common MHC types, has been hypothesized to maintain high, functionally relevant polymorphism of MHC, but demonstration of this relationship has remained elusive. In particular, differentiation of NFDS from fluctuating selection, resulting from changes in parasite communities in time and space (FS), has proved difficult in short-term studies. Here, we used temporal data, accumulated through long-term monitoring of helminths infecting bank voles (Myodes glareolus), to test specific predictions of NFDS on MHC class II. Data were collected in three, moderately genetically differentiated subpopulations in Poland, which were characterized by some stable spatiotemporal helminth communities but also events indicating introduction of new species and loss of others. We found a complex association between individual MHC diversity and species richness, where intermediate numbers of DRB supertypes correlated with lowest species richness, but the opposite was true for DQB supertypes-arguing against universal selection for immunogenetic optimality. We also showed that particular MHC supertypes explain a portion of the variance in prevalence and abundance of helminths, but this effect was subpopulation-specific, which is consistent with both NFDS and FS. Finally, in line with NFDS, we found that certain helminths that have recently colonized or spread in a given subpopulation, more frequently or intensely infected voles with MHC supertypes that have been common in the recent past. Overall, our results highlight complex spatial and temporal patterns of MHC-parasite associations, the latter being consistent with Red Queen coevolutionary dynamics.

The effects of Borrelia infection on its wintering rodent host.

In seasonal environments, appropriate adaptations are crucial for organisms to maximize their fitness. For instance, in many species, the immune function has been noticed to decrease during winter, which is assumed to be an adaptation to the season's limited food availability. Consequences of an infection on the health and survival of the host organism could thus be more severe in winter than in summer. Here, we experimentally investigated the effect of a zoonotic, endemic pathogen, Borrelia afzelii infection on the survival and body condition in its host, the bank vole (Myodes glareolus), during late autumn-early winter under semi-natural field conditions in 11 large outdoor enclosures. To test the interaction of Borrelia infection and energetic condition, four populations received supplementary nutrition, while remaining seven populations exploited only natural food sources. Supplementary food during winter increased the body mass independent of the infection status, however, Borrelia afzelii infection did not cause severe increase in the host mortality or affect the host body condition in the late autumn-early winter. While our study suggests that no severe effects are caused by B. afzelii infection on bank vole, further studies are warranted to identify any potentially smaller effects the pathogen may cause on the host fitness over the period of whole winter.

Parasitic nematodes of the genus Syphacia Seurat, 1916 infecting Cricetidae in the British Isles: the enigmatic status of Syphacia nigeriana.

Oxyurid nematodes (Syphacia spp.) from bank (Myodes glareolus) and field/common (Microtus spp.) voles, from disparate geographical sites in the British Isles, were examined morphologically and genetically. The genetic signatures of 118 new isolates are provided, based primarily on the rDNA internal transcribed spacers (ITS1-5.8S-ITS2) region and for representative isolates also on the small subunit 18S rDNA region and cytochrome c oxidase subunit 1 (cox-1) gene locus. Genetic data on worms recovered from Microtus spp. from the European mainland and from other rodent genera from the Palaearctic, North America and West Africa are also included. We test historical hypotheses indicating that S. nigeriana is a generalist species, infecting a range of different rodent genera. Our results establish that S. nigeriana is a parasite of both bank and field voles in the British Isles. An identical genotype was also recorded from Hubert's multimammate mouse (Mastomys huberti) from Senegal, but Mastomys spp. from West Africa were additionally parasitized by a related, although genetically distinct Syphacia species. We found no evidence for S. petrusewiczi in voles from the British Isles but isolates from Russia and North America were genetically distinct and formed their own separate deep branch in maximum likelihood molecular phylogenetic trees.

Major histocompatibility complex class I diversity limits the repertoire of T cell receptors.

Major histocompatibility complex (MHC) genes encode proteins that initiate adaptive immune responses through the presentation of foreign antigens to T cells. The high polymorphism found at these genes, thought to be promoted and maintained by pathogen-mediated selection, contrasts with the limited number of MHC loci found in most vertebrates. Although expressing many diverse MHC genes should broaden the range of detectable pathogens, it has been hypothesized to also cause deletion of larger fractions of self-reactive T cells, leading to a detrimental reduction of the T cell receptor (TCR) repertoire. However, a key prediction of this TCR depletion hypothesis, that the TCR repertoire should be inversely related to the individual MHC diversity, has never been tested. Here, using high-throughput sequencing and advanced sequencing error correction, we provide evidence of such an association in a rodent species with high interindividual variation in the number of expressed MHC molecules, the bank vole (Myodes glareolus). Higher individual diversity of MHC class I, but not class II, was associated with smaller TCR repertoires. Our results thus provide partial support for the TCR depletion model, while also highlighting the complex, potentially MHC class-specific mechanisms by which autoreactivity may trade off against evolutionary expansion of the MHC gene family.

Seasonal dynamics of tick burden and associated Borrelia burgdorferi s.l. and Borrelia miyamotoi infections in rodents in a Dutch forest ecosystem.

Ixodes ricinus ticks transmit Borrelia burgdorferi sensu lato (s.l.) as well as Borrelia miyamotoi. Larvae become infected when feeding on infected rodents, with horizontal transmission of B. burgdorferi and horizontal and vertical transmission of B. miyamotoi. We studied seasonal dynamics of infection rates of I. ricinus and their rodent hosts, and hence transmission risk of these two distinctly different Borrelia species. Rodents were live-trapped and inspected for ticks from May to November in 2013 and 2014 in a forest in The Netherlands. Trapped rodents were temporarily housed in the laboratory and detached ticks were collected. Borrelia infections were determined from the trapped rodents and collected ticks. Borrelia burgdorferi s.l. and B. miyamotoi were found in ticks as well as in rodents. Rodent density was higher in 2014, whereas tick burden as well as the Borrelia infection rates in rodents were higher in 2013. The density of B. miyamotoi-infected nymphs did not differ between the years. Tick burdens were higher on Apodemus sylvaticus than on Myodes glareolus, and higher on males than on females. Borrelia-infection rate of rodents varied strongly seasonally, peaking in summer. As the larval tick burden also peaked in summer, the generation of infected nymphs was highest in summer. We conclude that the heterogeneity of environmental and host-specific factors affects the seasonal transmission of Borrelia spp., and that these effects act more strongly on horizontally transmitted B. burgdorferi spp. than on the vertically transmitted B. miyamotoi.

Does selection for behavioral and physiological performance traits alter glucocorticoid responsiveness in bank voles?

One of the key elements of an animal's Darwinian fitness is its ability to adequately respond to and cope with challenging situations. Glucocorticoid hormones, such as corticosterone, affect an organism's ability to overcome such challenges. We hypothesized that changes in the glucocorticoid response curve contribute to the evolution of increased performance during challenging conditions, and tested it on bank voles (Myodes glareolus) from a multidirectional artificial selection experiment, which involves lines selected for high aerobic exercise metabolism achieved during swimming (A - Aerobic), predatory behavior towards a cricket (P - Predatory) and ability to maintain body mass on a low-quality herbivorous diet (H - Herbivorous), as well as unselected control lines (C - Control). We elicited a glucocorticoid response either by restraining the animal or by maximum pharmacological stimulation, and measured plasma corticosterone levels at baseline, during the response and during the recovery phase. Response-level corticosterone was higher in females, and recovery from maximal level was faster than that of males. Selection did not affect baseline or stress-induced corticosterone levels, but it decreased the maximum corticosterone level in Aerobic and Predatory lines, reducing the difference between stress-induced and maximum levels. Recovery from restraint-induced corticosterone level tended to be slower in the Herbivorous than in the other lines, an effect that was stronger in females than in males. In conclusion, successful selection for increased performance in challenging conditions was not associated with changes in absolute values of the glucocorticoid response to stress, but can affect other characteristics of the glucocorticoid response curve.

Balancing selection maintains polymorphisms at neurogenetic loci in field experiments.

Most variation in behavior has a genetic basis, but the processes determining the level of diversity at behavioral loci are largely unknown for natural populations. Expression of arginine vasopressin receptor 1a (Avpr1a) and oxytocin receptor (Oxtr) in specific regions of the brain regulates diverse social and reproductive behaviors in mammals, including humans. That these genes have important fitness consequences and that natural populations contain extensive diversity at these loci implies the action of balancing selection. In Myodes glareolus, Avpr1a and Oxtr each contain a polymorphic microsatellite locus located in their 5' regulatory region (the regulatory region-associated microsatellite, RRAM) that likely regulates gene expression. To test the hypothesis that balancing selection maintains diversity at behavioral loci, we released artificially bred females and males with different RRAM allele lengths into field enclosures that differed in population density. The length of Avpr1a and Oxtr RRAMs was associated with reproductive success, but population density and the sex interacted to determine the optimal genotype. In general, longer Avpr1a RRAMs were more beneficial for males, and shorter RRAMs were more beneficial for females; the opposite was true for Oxtr RRAMs. Moreover, Avpr1a RRAM allele length is correlated with the reproductive success of the sexes during different phases of reproduction; for males, RRAM length correlated with the numbers of newborn offspring, but for females selection was evident on the number of weaned offspring. This report of density-dependence and sexual antagonism acting on loci within the arginine vasopressin-oxytocin pathway explains how genetic diversity at Avpr1a and Oxtr could be maintained in natural populations.

Zoonotic Virus Seroprevalence among Bank Voles, Poland, 2002-2010.

Bank voles in Poland are reservoirs of zoonotic viruses. To determine seroprevalence of hantavirus, arenavirus, and cowpox virus and factors affecting seroprevalence, we screened for antibodies against these viruses over 9 years. Cowpox virus was most prevalent and affected by extrinsic and intrinsic factors. Long-term and multisite surveillance is crucial.

Exposure to environmental radionuclides is associated with altered metabolic and immunity pathways in a wild rodent.

Wildlife inhabiting environments contaminated by radionuclides face putative detrimental effects of exposure to ionizing radiation, with biomarkers such as an increase in DNA damage and/or oxidative stress commonly associated with radiation exposure. To examine the effects of exposure to radiation on gene expression in wildlife, we conducted a de novo RNA sequencing study of liver and spleen tissues from a rodent, the bank vole Myodes glareolus. Bank voles were collected from the Chernobyl Exclusion Zone (CEZ), where animals were exposed to elevated levels of radionuclides, and from uncontaminated areas near Kyiv, Ukraine. Counter to expectations, we did not observe a strong DNA damage response in animals exposed to radionuclides, although some signs of oxidative stress were identified. Rather, exposure to environmental radionuclides was associated with upregulation of genes involved in lipid metabolism and fatty acid oxidation in the livers - an apparent shift in energy metabolism. Moreover, using stable isotope analysis, we identified that fur from bank voles inhabiting the CEZ had enriched isotope values of nitrogen: such an increase is consistent with increased fatty acid metabolism, but also could arise from a difference in diet or habitat between the CEZ and elsewhere. In livers and spleens, voles inhabiting the CEZ were characterized by immunosuppression, such as impaired antigen processing, and activation of leucocytes involved in inflammatory responses. In conclusion, exposure to low dose environmental radiation impacts pathways associated with immunity and lipid metabolism, potentially as a stress-induced coping mechanism.

Early life of fathers affects offspring fitness in a wild rodent.

Intergenerational fitness effects on offspring due to the early life of the parent are well studied from the standpoint of the maternal environment, but intergenerational effects owing to the paternal early life environment are often overlooked. Nonetheless, recent laboratory studies in mammals and ecologically relevant studies in invertebrates predict that paternal effects can have a major impact on the offspring's phenotype. These nongenetic, environment-dependent paternal effects provide a mechanism for fathers to transmit environmental information to their offspring and could allow rapid adaptation. We used the bank vole Myodes glareolus, a wild rodent species with no paternal care, to test the hypothesis that a high population density environment in the early life of fathers can affect traits associated with offspring fitness. We show that the protein content in the diet and/or social environment experienced during the father's early life (prenatal and weaning) influence the phenotype and survival of his offspring and may indicate adaptation to density-dependent costs. Furthermore, we show that experiencing multiple environmental factors during the paternal early life can lead to a different outcome on the offspring phenotype than stimulated by experience of a single environmental factor, highlighting the need to study developmental experiences in tandem rather than independent of each other.

Stress coping and evolution of aerobic exercise performance: corticosterone levels in voles from a selection experiment.

The locomotor performance achieved in a challenging situation depends not only on physiological limitations, such as the aerobic exercise capacity, but also on behavioral characteristics, such as adequate coping with stress. The stress response is mediated largely by the hypothalamic-pituitary-adrenal (HPA) axis, through modulated release of glucocorticoids. We used a unique experimental evolution model system to test the hypothesis that the evolution of an increased aerobic exercise performance can be facilitated by modification of the glucocorticoid-related stress-coping mechanisms. Bank voles (Myodes glareolus) from 'aerobic' (A) lines, selected for 22 generations for high maximum swim-induced rate of oxygen consumption (V̇O2,swim), achieved a 64% higher V̇O2,swim than those from unselected, control lines. The temporal pattern of exercise during the swimming trial also evolved, and the A-line voles achieved V̇O2,swim later in the course of the trial, which indicates a modification of the stress response characteristics. Both V̇O2,swim and the average metabolic rate measured during the trial tended to increase with baseline corticosterone level, and decreased with the post-exercise corticosterone level. Thus, increased baseline corticosterone level promotes high metabolic performance, but a high corticosterone response to swimming acts as an inhibitor rather than stimulator of intense activity. However, neither of the corticosterone traits differed between the A-selected and control lines. Thus, the experiment did not provide evidence that evolution of increased aerobic performance is facilitated by the modification of glucocorticoid levels. The results, however, do not exclude the possibility that other aspects of the HPA axis function evolved in response to the selection.

Experimental evolution of aerobic exercise performance and hematological traits in bank voles.

The maximum rate of aerobic exercise metabolism (VO2max) is a trait informative from both medical and evolutionary perspective, and both the physiological mechanisms limiting its level and its evolution are subject to vivid debate. Both comparative analyses and studies on the effects of training or acclimation to aerobically-demanding conditions suggest a role of oxygen transport-related properties of blood in limiting VO2max. Here we used a unique experimental evolution model - lines of bank voles selected for high rate of swim-induced aerobic metabolism (VO2swim; A lines), which evolved a 60% higher VO2swim than that observed in unselected control lines (C) - and asked how the hematological parameters evolved in response to the selection. Voles from the A lines had a decreased hemoglobin concentration in cardiac blood samples (adjusted means ±â€¯SE, A: 14.7 ±â€¯1.1, C: 15.3 ±â€¯1.1 g/dl; p = .008), hematocrit (A: 51.7 ±â€¯4.5, C: 53.8 ±â€¯4.2%; p = .042) and tended to have less erythrocytes per microliter of blood (A: 11.5 ±â€¯1.4, C: 12.7 ±â€¯1.3 mln/µl; p = .083). The effect of selection was more pronounced in males than in females. Thus, selection for high aerobic-exercise performance resulted in a paradoxical decrease in traits positively associated with oxygen-carrying capacity per unit of blood volume, and the effect was sex-dependent. However, as a decreased blood viscosity associated with the lowered hematocrit reduces the costs of blood circulation, it can be hypothesized that the change can actually boost the oxygen supply to peripheral tissues.

Borrelia afzelii alters reproductive success in a rodent host.

The impact of a pathogen on the fitness and behaviour of its natural host depends upon the host-parasite relationship in a given set of environmental conditions. Here, we experimentally investigated the effects of Borrelia afzelii, one of the aetiological agents of Lyme disease in humans, on the fitness of its natural rodent host, the bank vole (Myodes glareolus), in semi-natural conditions with two contrasting host population densities. Our results show that B. afzelii can modify the reproductive success and spacing behaviour of its rodent host, whereas host survival was not affected. Infection impaired the breeding probability of large bank voles. Reproduction was hastened in infected females without alteration of the offspring size at birth. At low density, infected males produced fewer offspring, fertilized fewer females and had lower mobility than uninfected individuals. Meanwhile, the infection did not affect the proportion of offspring produced or the proportion of mating partner in female bank voles. Our study is the first to show that B. afzelii infection alters the reproductive success of the natural host. The effects observed could reflect the sickness behaviour due to the infection or they could be a consequence of a manipulation of the host behaviour by the bacteria.

Genomics of end-Pleistocene population replacement in a small mammal.

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene-Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.

Small-scale spatial variation in infection risk shapes the evolution of a Borrelia resistance gene in wild rodents.

Spatial variation in pathogen-mediated selection is predicted to influence the evolutionary trajectory of host populations and lead to spatial variation in their immunogenetic composition. However, to date few studies have been able to directly link small-scale spatial variation in infection risk to host immune gene evolution in natural, nonhuman populations. Here, we use a natural rodent-Borrelia system to test for associations between landscape-level spatial variation in Borrelia infection risk along replicated elevational gradients in the Swiss Alps and Toll-like receptor 2 (TLR2) evolution, a candidate gene for Borrelia resistance, across bank vole (Myodes glareolus) populations. We found that Borrelia infection risk (i.e., the product of Borrelia prevalence in questing ticks and the average tick load of voles at a sampling site) was spatially variable and significantly negatively associated with elevation. Across sampling sites, Borrelia prevalence in bank voles was significantly positively associated with Borrelia infection risk along the elevational clines. We observed a significant association between naturally occurring TLR2 polymorphisms in hosts and their Borrelia infection status. The TLR2 variant associated with a reduced likelihood of Borrelia infection was most common in rodent populations at lower elevations that face a high Borrelia infection risk, and its frequency changed in accordance with the change in Borrelia infection risk along the elevational clines. These results suggest that small-scale spatial variation in parasite-mediated selection affects the immunogenetic composition of natural host populations, providing a striking example that the microbial environment shapes the evolution of the host's immune system in the wild.

Bank vole immunoheterogeneity may limit Nephropatia Epidemica emergence in a French non-endemic region.

Ecoevolutionary processes affecting hosts, vectors and pathogens are important drivers of zoonotic disease emergence. In this study, we focused on nephropathia epidemica (NE), which is caused by Puumala hantavirus (PUUV) whose natural reservoir is the bank vole, Myodes glareolus. We questioned the possibility of NE emergence in a French region that is considered to be NE-free but that is adjacent to a NE-endemic region. We first confirmed the epidemiology of these two regions and we demonstrated the absence of spatial barriers that could have limited dispersal, and consequently, the spread of PUUV into the NE-free region. We next tested whether regional immunoheterogeneity could impact PUUV chances to circulate and persist in the NE-free region. We showed that bank voles from the NE-free region were sensitive to experimental PUUV infection. We observed high levels of immunoheterogeneity between individuals and also between regions. Antiviral gene expression (Tnf and Mx2) reached higher levels in bank voles from the NE-free region. During experimental infections, anti-PUUV antibody production was higher in bank voles from the NE-endemic region. These results indicated a lower susceptibility to PUUV for bank voles from this NE-free region, which might limit PUUV persistence and therefore, the risk of NE.

Puumala Virus in Bank Voles, Lithuania.

Little is known about the presence of human pathogenic Puumala virus (PUUV) in Lithuania. We detected this virus in bank voles (Myodes glareolus) in a region of this country in which previously PUUV-seropositive humans were identified. Our results are consistent with heterogeneous distributions of PUUV in other countries in Europe.

Phylogenetic analysis of Puumala virus strains from Central Europe highlights the need for a full-genome perspective on hantavirus evolution.

Puumala virus (PUUV), carried by bank voles (Myodes glareolus), is the medically most important hantavirus in Central and Western Europe. In this study, a total of 523 bank voles (408 from Germany, 72 from Slovakia, and 43 from Czech Republic) collected between the years 2007-2012 were analyzed for the presence of hantavirus RNA. Partial PUUV genome segment sequences were obtained from 51 voles. Phylogenetic analyses of all three genome segments showed that the newfound strains cluster with other Central and Western European PUUV strains. The new sequences from Sumava (Bohemian Forest), Czech Republic, are most closely related to the strains from the neighboring Bavarian Forest, a known hantavirus disease outbreak region. Interestingly, the Slovak strains clustered with the sequences from Bohemian and Bavarian Forests only in the M but not S segment analyses. This well-supported topological incongruence suggests a segment reassortment event or, as we analyzed only partial sequences, homologous recombination. Our data highlight the necessity of sequencing all three hantavirus genome segments and of a broader bank vole screening not only in recognized endemic foci but also in regions with no reported human hantavirus disease cases.