Genome-wide relaxation of selection and the evolution of the island syndrome in Orkney voles.

Island populations often experience different ecological and demographic conditions than their counterparts on the continent, resulting in divergent evolutionary forces affecting their genomes. Random genetic drift and selection both may leave their imprints on island populations, although the relative impact depends strongly on the specific conditions. Here we address their contributions to the island syndrome in a rodent with an unusually clear history of isolation. Common voles (Microtus arvalis) were introduced by humans on the Orkney archipelago north of Scotland >5000 years ago and rapidly evolved to exceptionally large size. Our analyses show that the genomes of Orkney voles were dominated by genetic drift, with extremely low diversity, variable Tajima's D, and very high divergence from continental conspecifics. Increased d N/d S ratios over a wide range of genes in Orkney voles indicated genome-wide relaxation of purifying selection. We found evidence of hard sweeps on key genes of the lipid metabolism pathway only in continental voles. The marked increase of body size in Orkney-a typical phenomenon of the island syndrome-may thus be associated to the relaxation of positive selection on genes related to this pathway. On the other hand, a hard sweep on immune genes of Orkney voles likely reflects the divergent ecological conditions and possibly the history of human introduction. The long-term isolated Orkney voles show that adaptive changes may still impact the evolutionary trajectories of such populations despite the pervasive consequences of genetic drift at the genome level.

Genome-wide support for incipient Tula hantavirus species within a single rodent host lineage.

Evolutionary divergence of viruses is most commonly driven by co-divergence with their hosts or through isolation of transmission after host shifts. It remains mostly unknown, however, whether divergent phylogenetic clades within named virus species represent functionally equivalent byproducts of high evolutionary rates or rather incipient virus species. Here, we test these alternatives with genomic data from two widespread phylogenetic clades in Tula orthohantavirus (TULV) within a single evolutionary lineage of their natural rodent host, the common vole Microtus arvalis. We examined voles from forty-two locations in the contact region between clades for TULV infection by reverse transcription (RT)-PCR. Sequencing yielded twenty-three TULV Central North and twenty-one TULV Central South genomes, which differed by 14.9-18.5 per cent at the nucleotide and 2.2-3.7 per cent at the amino acid (AA) level without evidence of recombination or reassortment between clades. Geographic cline analyses demonstrated an abrupt (<1 km wide) transition between the parapatric TULV clades in continuous landscape. This transition was located within the Central mitochondrial lineage of M. arvalis, and genomic single nucleotide polymorphisms showed gradual mixing of host populations across it. Genomic differentiation of hosts was much weaker across the TULV Central North to South transition than across the nearby hybrid zone between two evolutionary lineages in the host. We suggest that these parapatric TULV clades represent functionally distinct, incipient species, which are likely differently affected by genetic polymorphisms in the host. This highlights the potential of natural viral contact zones as systems for investigating the genetic and evolutionary factors enabling or restricting the transmission of RNA viruses.

Lymphocytic Choriomeningitis Virus Lineage V in Wood Mice, Germany.

We identified a novel lineage of lymphocytic choriomeningitis virus, tentatively named lineage V, in wood mice (Apodemus sylvaticus) from Germany. Wood mouse-derived lymphocytic choriomeningitis virus can be found across a substantially greater range than previously thought. Increased surveillance is needed to determine its geographic range and zoonotic potential.

Continuous presence of genetically diverse rustrela virus lineages in yellow-necked field mouse reservoir populations in northeastern Germany.

Rustrela virus (RusV; species Rubivirus strelense, family Matonaviridae) was discovered in different zoo animal species affected by fatal encephalitis. Simultaneous RusV RNA detection in multiple yellow-necked field mice (Apodemus flavicollis) suggested this rodent as a reservoir of RusV. Here, we investigated 1,264 yellow-necked field mice and sympatric other small mammals from different regions in Germany for RusV RNA using an optimized reverse transcription-quantitative polymerase chain reaction (RT-qPCR) protocol and high-throughput sequencing. The investigation resulted in the detection of RusV RNA exclusively in 50 of 396 (12.6 per cent) yellow-necked field mice but absence in other sympatric species. RT-qPCR-determined tissue distribution of RusV RNA revealed the highest viral loads in the central nervous system, with other tissues being only very rarely affected. The histopathological evaluation did not reveal any hints of encephalitis in the brains of infected animals despite the detection of viral RNA in neurons by in situ hybridization (ISH). The positive association between the body mass of yellow-necked field mice and RusV RNA detection suggests a persistent infection. Phylogenetic analysis of partial E1 and full-genome sequences showed a high diversification with at least four RusV lineages (1A-1D) in northeastern Germany. Moreover, phylogenetic and isolation-by-distance analyses indicated evolutionary processes of RusV mostly in local reservoir populations. A comparison of complete genome sequences from all detected RusV lineages demonstrated a high level of amino acid and nucleotide sequence variability within a part of the p150 peptide of the non-structural polyprotein and its coding sequence, respectively. The location of this region within the RusV genome and its genetic properties were comparable to the hypervariable region of the rubella virus. The broad range of detected RusV spillover hosts in combination with its geographical distribution in northeastern Germany requires the assessment of its zoonotic potential and further analysis of encephalitis cases in mammals. Future studies have to prove a putative co-evolution scenario for RusV in the yellow-necked field mouse reservoir.

Demographic history and genomic consequences of 10,000 generations of isolation in a wild mammal.

Increased human activities caused the isolation of populations in many species-often associated with genetic depletion and negative fitness effects. The effects of isolation are predicted by theory, but long-term data from natural populations are scarce. We show, with full genome sequences, that common voles (Microtus arvalis) in the Orkney archipelago have remained genetically isolated from conspecifics in continental Europe since their introduction by humans over 5,000 years ago. Modern Orkney vole populations are genetically highly differentiated from continental conspecifics as a result of genetic drift processes. Colonization likely started on the biggest Orkney island and vole populations on smaller islands were gradually split off, without signs of secondary admixture. Despite having large modern population sizes, Orkney voles are genetically depauperate and successive introductions to smaller islands resulted in further reduction of genetic diversity. We detected high levels of fixation of predicted deleterious variation compared with continental populations, particularly on smaller islands, yet the fitness effects realized in nature are unknown. Simulations showed that predominantly mildly deleterious mutations were fixed in populations, while highly deleterious mutations were purged early in the history of the Orkney population. Relaxation of selection overall due to benign environmental conditions on the islands and the effects of soft selection may have contributed to the repeated, successful establishment of Orkney voles despite potential fitness loss. Furthermore, the specific life history of these small mammals, resulting in relatively large population sizes, has probably been important for their long-term persistence in full isolation.

Pet Rats as the Likely Reservoir for Human Seoul Orthohantavirus Infection.

Seoul orthohantavirus (SEOV) is a rat-associated zoonotic pathogen with an almost worldwide distribution. In 2019, the first autochthonous human case of SEOV-induced hemorrhagic fever with renal syndrome was reported in Germany, and a pet rat was identified as the source of the zoonotic infection. To further investigate the SEOV reservoir, additional rats from the patient and another owner, all of which were purchased from the same vendor, were tested. SEOV RNA and anti-SEOV antibodies were found in both of the patient's rats and in two of the three rats belonging to the other owner. The complete coding sequences of the small (S), medium (M), and large (L) segments obtained from one rat per owner exhibited a high sequence similarity to SEOV strains of breeder rat or human origin from the Netherlands, France, the USA, and Great Britain. Serological screening of 490 rats from breeding facilities and 563 wild rats from Germany (2007-2020) as well as 594 wild rats from the Netherlands (2013-2021) revealed 1 and 6 seropositive individuals, respectively. However, SEOV RNA was not detected in any of these animals. Increased surveillance of pet, breeder, and wild rats is needed to identify the origin of the SEOV strain in Europe and to develop measures to prevent transmission to the human population.

Reemergence of Lymphocytic Choriomeningitis Mammarenavirus, Germany.

Lymphocytic choriomeningitis mammarenavirus (LCMV) is a globally distributed zoonotic pathogen transmitted by house mice (Mus musculus). We report the reemergence of LCMV (lineages I and II) in wild house mice (Mus musculus domesticus) and LCMV lineage I in a diseased golden lion tamarin (Leontopithecus rosalia) from a zoo in Germany.

Hantavirus Brno loanvirus is highly specific to the common noctule bat (Nyctalus noctula) and widespread in Central Europe.

Bat-associated hantaviruses have been detected in Asia, Africa and Europe. Recently, a novel hantavirus (Brno loanvirus, BRNV) was identified in common noctule bats (Nyctalus noctula) in the Czech Republic, but nothing is known about its geographical range and prevalence. The objective of this study was to evaluate the distribution and host specificity of BRNV by testing bats from neighbouring countries Germany, Austria and Poland. One thousand forty-seven bats representing 21 species from Germany, 464 bats representing 18 species from Austria and 77 bats representing 12 species from Poland were screened by L segment broad-spectrum nested reverse transcription-polymerase chain reaction (RT-PCR) or by BRNV-specific real-time RT-PCR. Three common noctules from Germany, one common noctule from Austria and three common noctules from Poland were positive in the hantavirus RNA screening. Conventional RT-PCR and primer walking resulted in the amplification of partial L segment and (almost) complete S and M segment coding sequences for samples from Germany and partial L segment sequences for samples from Poland. Phylogenetic analysis of these nucleotide sequences showed highest similarity to BRNV from Czech Republic. The exclusive detection of BRNV in common noctules from different countries suggests high host specificity. The RNA detection rate in common noctules ranged between 1 of 207 (0.5%; Austria), 3 of 245 (1.2%; Germany) and 3 of 20 (15%; Poland). In conclusion, this study demonstrates a broader distribution of BRNV in common noctules in Central Europe, but at low to moderate prevalence. Additional studies are needed to prove the zoonotic potential of this hantavirus and evaluate its transmission within bat populations.

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.

Cocirculation of Leptospira spp. and multiple orthohantaviruses in rodents, Lithuania, Northern Europe.

In Europe, zoonotic Leptospira spp. and orthohantaviruses are mainly associated with specific rodent hosts. These pathogens cause febrile human diseases with similar symptoms and disease progression. In Lithuania, the presence of Dobrava-Belgrade orthohantavirus (DOBV), Tula orthohantavirus (TULV) and Leptospira spp. in rodent reservoirs is still unknown, and Puumala orthohantavirus (PUUV) was detected in bank voles (Clethrionomys glareolus) at only one site. Therefore, we collected and screened 1617 rodents and insectivores from Lithuania for zoonotic (re-)emerging Leptospira and orthohantaviruses. We detected Leptospira DNA in six rodent species, namely striped field mouse (Apodemus agrarius), yellow-necked mouse (Apodemus flavicollis), bank vole, common vole (Microtus arvalis), field vole (Microtus agrestis) and root vole (Microtus oeconomus). Leptospira DNA was detected with an overall mean prevalence of 4.4% (range 3.7%-7.9% per rodent species). We detected DOBV RNA in 5.6% of the striped field mice, PUUV RNA in 1% of bank voles and TULV RNA in 4.6% of common voles, but no Leptospira DNA in shrews and no hantavirus-Leptospira coinfections in rodents. Based on the complete coding sequences of the three genome segments, two distant DOBV phylogenetic lineages in striped field mice, one PUUV strain in bank voles and two TULV strains in common voles were identified. The Leptospira prevalence for striped field mice and yellow-necked mice indicated a significant negative effect of the distance to water points. The detection of (re-)emerging human pathogenic Leptospira and three orthohantaviruses in rodent reservoirs in Lithuania calls for increased awareness of public health institutions and allows the improvement of molecular diagnostics for pathogen identification.

Fitness, risk taking, and spatial behavior covary with boldness in experimental vole populations.

Individuals of a population may vary along a pace-of-life syndrome from highly fecund, short-lived, bold, dispersive "fast" types at one end of the spectrum to less fecund, long-lived, shy, plastic "slow" types at the other end. Risk-taking behavior might mediate the underlying life history trade-off, but empirical evidence supporting this hypothesis is still ambiguous. Using experimentally created populations of common voles (Microtus arvalis)-a species with distinct seasonal life history trajectories-we aimed to test whether individual differences in boldness behavior covary with risk taking, space use, and fitness. We quantified risk taking, space use (via automated tracking), survival, and reproductive success (via genetic parentage analysis) in 8 to 14 experimental, mixed-sex populations of 113 common voles of known boldness type in large grassland enclosures over a significant part of their adult life span and two reproductive events. Populations were assorted to contain extreme boldness types (bold or shy) of both sexes. Bolder individuals took more risks than shyer ones, which did not affect survival. Bolder males but not females produced more offspring than shy conspecifics. Daily home range and core area sizes, based on 95% and 50% Kernel density estimates (20 ± 10 per individual, n = 54 individuals), were highly repeatable over time. Individual space use unfolded differently for sex-boldness type combinations over the course of the experiment. While day ranges decreased for shy females, they increased for bold females and all males. Space use trajectories may, hence, indicate differences in coping styles when confronted with a novel social and physical environment. Thus, interindividual differences in boldness predict risk taking under near-natural conditions and have consequences for fitness in males, which have a higher reproductive potential than females. Given extreme inter- and intra-annual fluctuations in population density in the study species and its short life span, density-dependent fluctuating selection operating differently on the sexes might maintain (co)variation in boldness, risk taking, and pace-of-life.

Genetic and biological characteristics of species A rotaviruses detected in common shrews suggest a distinct evolutionary trajectory.

Species A rotaviruses (RVAs) are important aetiological agents of severe diarrhoea in young children. They are also widely distributed in mammals and birds, and increasing evidence indicates the possibility of zoonotic transmission of RVA strains between animals and humans. Moreover, reassortment of the eleven segments of the RVA genome can result in rapid biological changes and may influence pathogenic properties. Here, the nearly complete genome of an RVA strain from a common shrew (Sorex araneus) was sequenced, which showed high nucleotide sequence similarity to additionally determined partial sequences from common shrew RVAs but only very low identity (below 68 per cent) to RVAs from other animal species and humans. New genotypes were assigned to most genome segments of the novel common shrew RVA strain KS14/269, resulting in the genome constellation G39-P[55]-I27-R26-C22-M22-A37-N26-T26-E30-H26. Phylogenetic analyses clustered the common shrew RVAs as ancestral branches of other mammalian and avian RVAs for most of the genome segments, which is in contrast to the phylogeny of the hosts. Nevertheless, conserved sequences typical for all RVAs were identified at the 5'- and 3'- non-coding segment termini. To explore whether the common shrew RVA can exchange genetic material with other mammalian RVAs by reassortment, a reverse genetics system based on the simian RVA strain SA11 was used. However, no viable reassortants could be rescued by exchanging the VP4-, VP6-, or VP7-encoding genome segment alone or in combinations. It can be concluded that highly divergent RVAs are present in common shrews, indicating an evolution of these viruses largely separated from other mammalian and avian RVAs. The zoonotic potential of the virus seems to be low but needs to be further analysed in future.

Zoonotic pathogen screening of striped field mice (Apodemus agrarius) from Austria.

The striped field mouse (Apodemus agrarius) is known to carry several zoonotic pathogens, including Leptospira spp. and Dobrava-Belgrade orthohantavirus (DOBV). Since its first detection in 1996 in south-east Austria, the striped field mouse has further expanded its range in Austria. Here, we screened 35 striped field mice collected in an Austrian region near the Hungarian border for DOBV, Leptospira spp. and seven vector-borne pathogens. Hantavirus RT-PCR screening and DOBV IgG ELISA analysis led to the detection of two DOBV-positive striped field mice. The complete coding sequences of all three genome segments of both strains were determined by a combination of target enrichment and next-generation sequencing. Both complete coding S segment sequences clustered within the DOBV genotype Kurkino clade with the highest similarity to a sequence from Hungary. In one of 35 striped field mice, Leptospira borgpetersenii sequence type (ST) 146 was detected. Bartonella spp., Borrelia miyamotoi and Neoehrlichia mikurensis DNA was detected in four, one and two of 32 mice, respectively. Babesia, Anaplasma, Ehrlichia and Rickettsia specific DNA was not detected. Future investigations will have to determine the prevalence and invasion of these pathogens with the ongoing range expansion of the striped field mouse in Austria.

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.

Seasonal migration patterns and the maintenance of evolutionary diversity in a cryptic bird radiation.

Morphological differentiation associated with evolutionary diversification is often explained with adaptive benefits but the processes and mechanisms maintaining cryptic diversity are still poorly understood. Using genome-wide data, we show here that the pale sand martin Riparia diluta in Central and East Asia consists of three genetically deeply differentiated lineages which vary only gradually in morphology but broadly reflect traditional taxonomy. We detected no signs of gene flow along the eastern edge of the Qinghai-Tibetan plateau between lowland south-eastern Chinese R. d. fohkienensis and high-altitude R. d. tibetana. Largely different breeding and migration timing between these low and high altitude populations as indicated by phenology data suggests that allochrony might act as prezygotic isolation mechanism in the area where their ranges abut. Mongolian populations of R. d. tibetana, however, displayed signs of limited mixed ancestries with Central Asian R. d. diluta. Their ranges meet in the area of a well-known avian migratory divide, where western lineages take a western migration route around the Qinghai-Tibetan plateau to winter quarters in South Asia, and eastern lineages take an eastern route to Southeast Asia. This might also be the case between western R. d. diluta and eastern R. d. tibetana as indicated by differing wintering grounds. We hypothesize that hybrids might have nonoptimal intermediate migration routes and selection against them might restrict gene flow. Although further potential isolation mechanisms might exist in the pale sand martin, our study points towards contrasting migration behaviour as an important factor in maintaining evolutionary diversity under morphological stasis.

Variation and Selection in the Putative Sperm-Binding Region of ZP3 in Muroid Rodents: A Comparison between Cricetids and Murines.

In mammals, the zona pellucida glycoprotein 3 (ZP3) is considered a primary sperm receptor of the oocyte and is hypothesized to be involved in reproductive isolation. We investigated patterns of diversity and selection in the putative sperm-binding region (pSBR) of mouse ZP3 across Cricetidae and Murinae, two hyperdiverse taxonomic groups within muroid rodents. In murines, the pSBR is fairly conserved, in particular the serine-rich stretch containing the glycosylation sites proposed as essential for sperm binding. In contrast, cricetid amino acid sequences of the pSBR were much more variable and the serine-rich motif, typical of murines, was generally substantially modified. Overall, our results suggest a general lack of species specificity of the pSBR across the two muroid families. We document statistical evidence of positive selection acting on exons 6 and 7 of ZP3 and identified several amino acid sites that are likely targets of selection, with most positively selected sites falling within or adjacent to the pSBR.

Spatial and Temporal Dynamics and Molecular Evolution of Tula orthohantavirus in German Vole Populations.

Tula orthohantavirus (TULV) is a rodent-borne hantavirus with broad geographical distribution in Europe. Its major reservoir is the common vole (Microtus arvalis), but TULV has also been detected in closely related vole species. Given the large distributional range and high amplitude population dynamics of common voles, this host-pathogen complex presents an ideal system to study the complex mechanisms of pathogen transmission in a wild rodent reservoir. We investigated the dynamics of TULV prevalence and the subsequent potential effects on the molecular evolution of TULV in common voles of the Central evolutionary lineage. Rodents were trapped for three years in four regions of Germany and samples were analyzed for the presence of TULV-reactive antibodies and TULV RNA with subsequent sequence determination. The results show that individual (sex) and population-level factors (abundance) of hosts were significant predictors of local TULV dynamics. At the large geographic scale, different phylogenetic TULV clades and an overall isolation-by-distance pattern in virus sequences were detected, while at the small scale (<4 km) this depended on the study area. In combination with an overall delayed density dependence, our results highlight that frequent, localized bottleneck events for the common vole and TULV do occur and can be offset by local recolonization dynamics.

A Putative Novel Hepatitis E Virus Genotype 3 Subtype Identified in Rabbit, Germany 2016.

Hepatitis E is an emerging viral disease that is the leading cause of viral hepatitis in the world. The vast majority of hepatitis E cases in developed countries are caused by zoonotic genotypes 3 and 4 of hepatitis E virus (HEV) for which pig and wild boar and to lesser extent rabbits are the main reservoir. According to recent reports rabbits are a source of human HEV infection and highlight the risk of zoonotic foodborne transmission. Here we report the molecular analysis of a novel HEV strain identified in a rabbit during a countrywide surveillance of rabbits and hares in Germany, 2016. The analysis of the complete genome reveals characteristics of a putative novel recombinant subtype of the species Orthohepevirus A within the clade of genotype 3 but not closely related to any known subtypes. Importantly, the genome of this strain possesses a nucleotide exchange in the overlapping region of open reading frames ORF2/ORF3 interfering with a broadly applied diagnostic real-time RT-PCR. In conclusion, a new type of HEV strain was identified in a German rabbit with atypical and novel sequence characteristics.

Geographical Distribution and Genetic Diversity of Bank Vole Hepaciviruses in Europe.

The development of new diagnostic methods resulted in the discovery of novel hepaciviruses in wild populations of the bank vole (Myodes glareolus, syn. Clethrionomys glareolus). The naturally infected voles demonstrate signs of hepatitis similar to those induced by hepatitis C virus (HCV) in humans. The aim of the present research was to investigate the geographical distribution of bank vole-associated hepaciviruses (BvHVs) and their genetic diversity in Europe. Real-time reverse transcription polymerase chain reaction (RT-qPCR) screening revealed BvHV RNA in 442 out of 1838 (24.0%) bank voles from nine European countries and in one of seven northern red-backed voles (Myodes rutilus, syn. Clethrionomys rutilus). BvHV RNA was not found in any other small mammal species (n = 23) tested here. Phylogenetic and isolation-by-distance analyses confirmed the occurrence of both BvHV species (Hepacivirus F and Hepacivirus J) and their sympatric occurrence at several trapping sites in two countries. The broad geographical distribution of BvHVs across Europe was associated with their presence in bank voles of different evolutionary lineages. The extensive geographical distribution and high levels of genetic diversity of BvHVs, as well as the high population fluctuations of bank voles and occasional commensalism in some parts of Europe warrant future studies on the zoonotic potential of BvHVs.