Transalpine colonisation and partial phylogeographic erosion by dispersal in the common vole (Microtus arvalis)

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.

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.

Geography and past climate changes have shaped the evolution of a widespread lizard in arid Central Asia

The complex orogenic history and structure of Central Asia, coupled with Pleistocene glacial cycles have generated its stepwise aridification. Such events would have significantly influenced the evolution of many mid-latitude species in arid Central Asia (ACA). In this study, we employed two mitochondrial genes (CO1 and ND2) and genome-wide SNPs, coupled with ecological niche modeling, to investigate the lineage diversification and historical demography within a widespread lizard Phrynocepahlus helioscopus, and their associations with geography and past climate change. We obtained themtDNA dataset for 300 individuals from 96 localities within the known range of the lizard, among which 51 individuals from 27 localities were selected for generating the SNP dataset via genotyping-by-sequencing approach. Phylogenetic analyses of the concatenated mtDNA dataset revealed eight geographically correlated lineages that diverged by 4.21-10.41% for the CO1 gene, which were estimated to have coalesced ∼4.47 million years ago. However, we observed mito-nuclear discordance pattern regarding the population of Clade V (P. helioscopus sergeevi) from Tajikistan. Ancestral area estimations suggested that P. helioscopus originated from the Fergana Valley and then dispersed into the adjacent areas in ACA along with a history of multiple allopatric divergence processes, suggesting that Fergana may have been the cradle of diversification of P. helioscopus. The intensification of aridification across Central Asia during the Late Pliocene may have facilitated the rapid radiation of this arid-adapted lizard throughout this vast territory. Subsequently, the geological events (e.g., uplift of the Hissar-Alay, transgressions of the Caspian Sea) and geographic barriers (e.g., Amu Darya River, Zerarshan River) during the Pleistocene triggered the progressive diversification of P. helioscopus. Interestingly, Clade VIII (P. helioscopus varius) experienced rapid population growth coupled with range expansion while Clade IV (P. helioscopus cameranoi) underwent drastic population expansion associated with range contraction during the Last Glacial Maximum. In Clade IV, but not in Clade VIII, environmental turnover contributes more to mitochondrial genetic distinctiveness than geographic distance. Overall, the SNP dataset demonstrates that geographic distance plays a greater role than environmental distance. Both the mtDNA dataset and the SNP dataset suggest local-scale genetic differentiation in Clade IV and Clade VIII, revealing potential geographic barriers in the Ili River Valley and the Junggar Basin, respectively. Twenty-seven outlier SNPs associated with environmental factors (precipitation and temperature) were identified, which supports the signature of local adaptation to the arid desert environment. Finally, our finding suggests taxonomic implications, such as support for full species status for P. saidalievi (Clade II) and P. meridionalis (Clade I). Future analyses based on further evidence and increased taxon and geographic sampling should be carried out to corroborate our findings.

Phylogeographic analysis of Tula hantavirus highlights a single introduction to central Europe

Orthohantaviruses are zoonotic pathogens of humans, unique among the bunyaviruses in not being transmitted by an arthropod vector. Tula orthohantavirus (TULV) is an old-world hantavirus, of yet unclear human pathogenicity, with few reported cases of clinically relevant human infection. So far, phylogeographic studies exploring the global pathways of hantaviral migration are scarce and generally do not focus on a specific hantavirus species. The aim of the present study was to reconstruct the dispersal history of TULV lineages across Eurasia based on S segment sequences sampled from different geographic areas. Maximum-likelihood and Bayesian inference methods were used to perform the phylogenetic analysis and phylogeographic reconstructions. Sampling time and trapping localities were obtained for a total of 735 TULV S segment sequences available in public databases at the time of the study. The estimated substitution rate of the analyzed partial S segment alignment was 2.26 × 10-3 substitutions/site/year (95 per cent highest posterior density interval: 1.79 × 10-3 to 2.75 × 10-3). Continuous phylogeography of TULV S segment sequences placed the potential root and origin of TULV spread in the Black Sea region. In our study, we detect a single-lineage introduction of TULV to Europe, followed by local viral circulation further on.

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.

European common frog (Rana temporaria) recolonized Switzerland from multiple glacial refugia in northern Italy via trans- and circum-Alpine routes

The high mountain ranges of Western Europe had a profound effect on the biotic recolonization of Europe from glacial refugia. The Alps present a particularly interesting case because they form an absolute barrier to dispersal for most taxa, obstructing recolonization from multiple refugia in northern Italy. Here, we investigate the effect of the European Alps on the phylogeographic history of the European common frog Rana temporaria. Based on partial cytochrome b and COXI sequences from Switzerland, we find two mitochondrial lineages roughly north and south of the Alpine ridge, with contact zones between them in eastern and western Switzerland. The northern haplogroup falls within the previously identified Western European haplogroup, while the southern haplogroup is unique to Switzerland. We find that the lineages diverged ~110 kya, at approximately the onset of the last glacial glaciation; this indicates that they are from different glacial refugia. Phylogenetic analyses suggest that the northern and southern haplogroups colonized Switzerland via trans- and circum-Alpine routes from at least two separate refugia in northern Italy. Our results illustrate how a complex recolonization history of the central European Alps can arise from the semi-permeable barrier created by high mountains.

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.

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.

Secondary contact between diverged host lineages entails ecological speciation in a European hantavirus

The diversity of viruses probably exceeds biodiversity of eukaryotes, but little is known about the origin and emergence of novel virus species. Experimentation and disease outbreak investigations have allowed the characterization of rapid molecular virus adaptation. However, the processes leading to the establishment of functionally distinct virus taxa in nature remain obscure. Here, we demonstrate that incipient speciation in a natural host species has generated distinct ecological niches leading to adaptive isolation in an RNA virus. We found a very strong association between the distributions of two major phylogenetic clades in Tula orthohantavirus (TULV) and the rodent host lineages in a natural hybrid zone of the European common vole (Microtus arvalis). The spatial transition between the virus clades in replicated geographic clines is at least eight times narrower than between the hybridizing host lineages. This suggests a strong barrier for effective virus transmission despite frequent dispersal and gene flow among local host populations, and translates to a complete turnover of the adaptive background of TULV within a few hundred meters in the open, unobstructed landscape. Genetic differences between TULV clades are homogenously distributed in the genomes and mostly synonymous (93.1%), except for a cluster of nonsynonymous changes in the 5′ region of the viral envelope glycoprotein gene, potentially involved in host-driven isolation. Evolutionary relationships between TULV clades indicate an emergence of these viruses through rapid differential adaptation to the previously diverged host lineages that resulted in levels of ecological isolation exceeding the progress of speciation in their vertebrate hosts.

Like host, like virus; analysis of a natural hybrid zone in the European common vole reveals speciation processes in Tula hantaviruses triggered by evolutionary divergence in the rodent host lineages.

Natural biodiversity is driven by stochastic processes and evolutionary adaptation to ecological niches. In viruses, adaptation to specific hosts may cause diversification and eventually lead to the emergence of novel viruses. Here, we studied diversity in Tula orthohantavirus (TULV) in relation to evolutionary divergence in its natural rodent host, the European common vole (Microtus arvalis). In a geographical region in which two distinct evolutionary lineages in the common vole interact and interbreed (a hybrid zone), we found two substantially different TULV clades. Phylogenetic analyses revealed that the divergence among virus clades was likely triggered by a shift of an ancestral virus between the previously diverged host lineages in the hybrid zone. The strong association between virus clades and host lineages at a fine geographical scale results in effective separation of TULVs, despite incomplete reproductive isolation and frequent gene flow among local host populations. Virus genome sequences pointed to the amino-terminal part of the envelope protein as an important region for functional differentiation among these virus clades.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Between the Balkans and the Baltic: Phylogeography of a Common Vole Mitochondrial DNA Lineage Limited to Central Europe

The common vole (Microtus arvalis) has been a model species of small mammal for studying end-glacial colonization history. In the present study we expanded the sampling from central and eastern Europe, analyzing contemporary genetic structure to identify the role of a potential 'northern glacial refugium', i.e. a refugium at a higher latitude than the traditional Mediterranean refugia. Altogether we analyzed 786 cytochrome b (cytb) sequences (representing mitochondrial DNA; mtDNA) from the whole of Europe, adding 177 new sequences from central and eastern Europe, and we conducted analyses on eight microsatellite loci for 499 individuals (representing nuclear DNA) from central and eastern Europe, adding data on 311 new specimens. Our new data fill gaps in the vicinity of the Carpathian Mountains, the potential northern refugium, such that there is now dense sampling from the Balkans to the Baltic Sea. Here we present evidence that the Eastern mtDNA lineage of the common vole was present in the vicinity of this Carpathian refugium during the Last Glacial Maximum and the Younger Dryas. The Eastern lineage expanded from this refugium to the Baltic and shows low cytb nucleotide diversity in those most northerly parts of the distribution. Analyses of microsatellites revealed a similar pattern but also showed little differentiation between all of the populations sampled in central and eastern Europe.

High genetic structuring of Tula hantavirus

Tula virus (TULV) is a vole-associated hantavirus with low or no pathogenicity to humans. In the present study, 686 common voles (Microtus arvalis), 249 field voles (Microtus agrestis) and 30 water voles (Arvicola spec.) were collected at 79 sites in Germany, Luxembourg and France and screened by RT-PCR and TULV-IgG ELISA. TULV-specific RNA and/or antibodies were detected at 43 of the sites, demonstrating a geographically widespread distribution of the virus in the studied area. The TULV prevalence in common voles (16.7 %) was higher than that in field voles (9.2 %) and water voles (10.0 %). Time series data at ten trapping sites showed evidence of a lasting presence of TULV RNA within common vole populations for up to 34 months, although usually at low prevalence. Phylogenetic analysis demonstrated a strong genetic structuring of TULV sequences according to geography and independent of the rodent species, confirming the common vole as the preferential host, with spillover infections to co-occurring field and water voles. TULV phylogenetic clades showed a general association with evolutionary lineages in the common vole as assessed by mitochondrial DNA sequences on a large geographical scale, but with local-scale discrepancies in the contact areas.

Tracing reinforcement through asymmetrical partner preference in the European common vole Microtus arvalis

Background

The mechanistic basis of speciation and in particular the contribution of behaviour to the completion of the speciation process is often contentious. Contact zones between related taxa provide a situation where selection against hybridization might reinforce separation by behavioural mechanisms, which could ultimately fully isolate the taxa. One of the most abundant European mammals, the common vole Microtus arvalis, forms multiple natural hybrid zones where rapidly diverging evolutionary lineages meet in secondary contact. Very narrow zones of hybridization spanning only a few kilometres and sex-specific gene flow patterns indicate reduced fitness of natural hybrids and incipient speciation between some of the evolutionary lineages. In this study, we examined the contribution of behavioural mechanisms to the speciation process in these rodents by fine-mapping allopatric and parapatric populations in the hybrid zone between the Western and Central lineages and experimental testing of the partner preferences of wild, pure-bred and hybrid female common voles.

Results

Genetic analysis based on microsatellite markers revealed the presence of multiple parapatric and largely non-admixed populations at distances of about 10 km at the edge of the area of natural hybridization between the Western and Central lineages. Wild females from Western parapatric populations and lab-born F1 hybrids preferred males from the Western lineage whereas wild females of Central parapatric origin showed no measurable preference. Furthermore, wild and lab-born females from allopatric populations of the Western or Central lineages showed no detectable preference for males from either lineage.

Conclusions

The detected partner preferences are consistent with asymmetrical reinforcement of pre-mating reproductive isolation mechanisms in the European common vole and with earlier results suggesting that hybridization is more detrimental to the Western lineage. As a consequence, these differences in behaviour might contribute to a further geographical stabilization of this moving hybrid zone. Such behavioural processes could also provide a mechanistic perspective for frequently-detected asymmetrical introgression patterns in the largely allopatrically diversifying Microtus genus and other rapidly speciating rodents.

Continental-scale footprint of balancing and positive selection in a small rodent (Microtus arvalis)

Genetic adaptation to different environmental conditions is expected to lead to large differences between populations at selected loci, thus providing a signature of positive selection. Whereas balancing selection can maintain polymorphisms over long evolutionary periods and even geographic scale, thus leads to low levels of divergence between populations at selected loci. However, little is known about the relative importance of these two selective forces in shaping genomic diversity, partly due to difficulties in recognizing balancing selection in species showing low levels of differentiation. Here we address this problem by studying genomic diversity in the European common vole (Microtus arvalis) presenting high levels of differentiation between populations (average F_ST = 0.31). We studied 3,839 Amplified Fragment Length Polymorphism (AFLP) markers genotyped in 444 individuals from 21 populations distributed across the European continent and hence over different environmental conditions. Our statistical approach to detect markers under selection is based on a Bayesian method specifically developed for AFLP markers, which treats AFLPs as a nearly codominant marker system, and therefore has increased power to detect selection. The high number of screened populations allowed us to detect the signature of balancing selection across a large geographic area. We detected 33 markers potentially under balancing selection, hence strong evidence of stabilizing selection in 21 populations across Europe. However, our analyses identified four-times more markers (138) being under positive selection, and geographical patterns suggest that some of these markers are probably associated with alpine regions, which seem to have environmental conditions that favour adaptation. We conclude that despite favourable conditions in this study for the detection of balancing selection, this evolutionary force seems to play a relatively minor role in shaping the genomic diversity of the common vole, which is more influenced by positive selection and neutral processes like drift and demographic history.

The changing pace of insular life: 5000 years of microevolution in the Orkney vole (Microtus arvalis orcadensis)

Island evolution may be expected to involve fast initial morphological divergence followed by stasis. We tested this model using the dental phenotype of modern and ancient common voles (Microtus arvalis), introduced onto the Orkney archipelago (Scotland) from continental Europe some 5000 years ago. First, we investigated phenotypic divergence of Orkney and continental European populations and assessed climatic influences. Second, phenotypic differentiation among Orkney populations was tested against geography, time, and neutral genetic patterns. Finally, we examined evolutionary change along a time series for the Orkney Mainland. Molar gigantism and anterior-lobe hypertrophy evolved rapidly in Orkney voles following introduction, without any transitional forms detected. Founder events and adaptation appear to explain this initial rapid evolution. Idiosyncrasy in dental features among different island populations of Orkney voles is also likely the result of local founder events following Neolithic translocation around the archipelago. However, against our initial expectations, a second marked phenotypic shift occurred between the 4th and 12th centuries AD, associated with increased pastoral farming and introduction of competitors (mice and rats) and terrestrial predators (foxes and cats). These results indicate that human agency can generate a more complex pattern of morphological evolution than might be expected in island rodents.

Stability of the southern European border of Echinococcus multilocularis in the Alps: evidence that Microtus arvalis is a limiting factor

SUMMARY The known range of the zoonotic fox tapeworm Echinococcus multilocularis has expanded since the 1990s, and today this parasite is recorded in higher abundances throughout large parts of Europe. This phenomenon is mostly attributed to the increasing European fox populations and their invasion of urban habitats. However, these factors alone are insufficient to explain the heterogeneous distribution of the parasite in Europe. Here, we analysed the spatial interrelationship of E. multilocularis with the known distribution of seven vole species in Ticino, southern Switzerland. Among 404 necropsied foxes (1990-2006) and 79 fox faecal samples (2010-2012), E. multilocularis was consistently found in the north of the investigated area. No expansion of this endemic focus was recorded during the 22 years of the study period. This stable endemic focus is coincident with the known distribution of the vole species Microtus arvalis but not, or only partly, with the distribution of the other autochthonous vole species. Our results give evidence that this vole species plays a crucial role in the maintenance of the parasite's life cycle and that its absence could be a limiting factor for the spread of E. multilocularis in this region.

First molecular evidence for Puumala hantavirus in Poland

Puumala virus (PUUV) causes mild to moderate cases of haemorrhagic fever with renal syndrome (HFRS), and is responsible for the majority of hantavirus infections of humans in Fennoscandia, Central and Western Europe. Although there are relatively many PUUV sequences available from different European countries, little is known about the presence of this virus in Poland. During population studies in 2009 a total of 45 bank voles were trapped at three sites in north-eastern Poland, namely islands on Dejguny and Dobskie Lakes and in a forest near Mikołajki. S and M segment-specific RT-PCR assays detected PUUV RNA in three animals from the Mikołajki site. The obtained partial S and M segment sequences demonstrated the highest similarity to the corresponding segments of a PUUV strain from Latvia. Analysis of chest cavity fluid samples by IgG ELISA using a yeast-expressed PUUV nucleocapsid protein resulted in the detection of two seropositive samples, both being also RT-PCR positive. Interestingly, at the trapping site in Mikołajki PUUV-positive bank voles belong to the Carpathian and Eastern genetic lineages within this species. In conclusion, we herein present the first molecular evidence for PUUV in the rodent reservoir from Poland.

Phylogeography of Silver Pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Continuing fragmentation of a widespread species by geographical barriers as initial step in a land snail radiation on crete

The phylogeographic structure of the land snail Xerocrassa mesostena on Crete inferred from AFLP markers and mitochondrial cox1 sequences can be explained by three mechanisms: gene flow restriction, population expansion and leptokurtic dispersal. Gene flow restriction by geographic barriers caused subdivision of the gene pool into distinct clusters. Population expansion was probably facilitated by deforestation of Crete in the postglacial. Newly available areas were colonized by leptokurtic dispersal, i.e. slow active expansion resulting in isolation by distance within the clusters and occasional long distance dispersal events that resulted in departures from the isolation by distance model. Less than one percent of the AFLP markers show correlations with environmental variables. Random phylogeographic breaks in the distribution of the mitochondrial haplotype groups indicate that single locus markers, especially mitochondrial DNA, might result in a misleading picture of the phylogeographic structure of a species. Restriction of gene flow between metapopulations caused by geographical barriers can interact with sexual selection resulting in the differentiation of these metapopulations into separate species without noticeable ecological differentiation. Evidence for gene flow between parapatrically distributed evolutionary units representing different stages of the speciation process suggests that the ongoing process of fragmentation of the X. mesostena complex might be an example for parapatric speciation. The lack of ecological differentiation between these units confirms theoretical predictions that divergent selection for local adaptation is not required for rapid speciation.

Sex-specific clines support incipient speciation in a common European mammal

Hybrid zones provide excellent opportunities to study processes and mechanisms underlying reproductive isolation and speciation. Here we investigated sex-specific clines of molecular markers in hybrid zones of morphologically cryptic yet genetically highly-diverged evolutionary lineages of the European common vole (Microtus arvalis). We analyzed the position and width of four secondary contact zones along three independent transects in the region of the Alps using maternally (mitochondrial DNA) and paternally (Y-chromosome) inherited genetic markers. Given male-biased dispersal in the common vole, a selectively neutral secondary contact would show broader paternal marker clines than maternal ones. In a selective case, for example, involving a form of Haldane's rule, Y-chromosomal clines would not be expected to be broader than maternal markers because they are transmitted by the heterogametic sex and thus gene flow would be restricted. Consistent with the selective case, paternal clines were significantly narrower or at most equal in width to maternal clines in all contact zones. In addition, analyses using maximum likelihood cline-fitting detected a shift of paternal relative to maternal clines in three of four contact zones. These patterns suggest that processes at the contact zones in the common vole are not selectively neutral, and that partial reproductive isolation is already established between these evolutionary lineages. We conclude that hybrid zone movement, sexual selection and/or genetic incompatibilities are likely associated with an unusual unidirectional manifestation of Haldane's rule in this common European mammal.

Deep phylogeographic divergence and cytonuclear discordance in the grasshopper Oedaleus decorus

The grasshopper Oedaleus decorus is a thermophilic insect with a large, mostly south-Palaearctic distribution range, stretching from the Mediterranean regions in Europe to Central-Asia and China. In this study, we analyzed the extent of phylogenetic divergence and the recent evolutionary history of the species based on 274 specimens from 26 localities across the distribution range in Europe. Phylogenetic relationships were determined using sequences of two mitochondrial loci (ctr, ND2) with neighbour-joining and Bayesian methods. Additionally, genetic differentiation was analyzed based on mitochondrial DNA and 11 microsatellite markers using F-statistics, model-free multivariate and model-based Bayesian clustering approaches. Phylogenetic analyses detected consistently two highly divergent, allopatrically distributed lineages within O. decorus. The divergence among these Western and Eastern lineages meeting in the region of the Alps was similar to the divergence of each lineage to the sister species O. asiaticus. Genetic differentiation for ctr was extremely high between Western and Eastern grasshopper populations (F(ct)=0.95). Microsatellite markers detected much lower but nevertheless very significant genetic structure among population samples. The nuclear data also demonstrated a case of cytonuclear discordance because the affiliation with mitochondrial lineages was incongruent in Northern Italy. Taken together these results provide evidence of an ancient separation within Oedaleus and either historical introgression of mtDNA among lineages and/or ongoing sex-specific gene flow in this grasshopper. Our study stresses the importance of multilocus approaches for unravelling the history and status of taxa of uncertain evolutionary divergence.

Genetic structure of muskrat (Ondatra zibethicus) and its concordance with taxonomy in North America

Extrinsic factors such as physical barriers play an important role in shaping population genetic structure. A reduction in gene flow leading to population structuring may ultimately lead to population divergence. These divergent populations are often considered subspecies. Because genetic differentiation may represent differences between subspecies, patterns of genetic structure should reflect subspecies groupings. In this study, we examine the contemporary population genetic structure of muskrat (n = 331) and assess the relevance of 4 geographically distinct subspecies designations across northern North America using 9 microsatellite loci. We predicted that patterns of gene flow and genetic structure would reflect the described subspecies. We found evidence of genetic differentiation between western and eastern regions, and muskrats from Newfoundland (NF) showed significantly lower genetic diversity than central regions. A strong isolation by distance pattern was also detected within the eastern cluster. Our results did not differentiate Ondatra zibethicus spatulus (northwest) from O. z. albus (central), but they suggest a distinction between O. z. obscurus (NF) and O. z. zibethicus (east). This study highlights the need for more phylogenetic studies in order to better understand intraspecific divergence and the genetic characterization of subspecies.

Phylogenetic analysis of Puumala virus subtype Bavaria, characterization and diagnostic use of its recombinant nucleocapsid protein

Puumala virus (PUUV) is the predominant hantavirus species in Germany causing large numbers of mild to moderate cases of haemorrhagic fever with renal syndrome (HFRS). During an outbreak in South-East Germany in 2004 a novel PUUV subtype designated Bavaria was identified as the causative agent of HFRS in humans [1]. Here we present a molecular characterization of this PUUV strain by investigating novel partial and almost entire nucleocapsid (N) protein-encoding small (S-) segment sequences and partial medium (M-) segment sequences from bank voles (Myodes glareolus) trapped in Lower Bavaria during 2004 and 2005. Phylogenetic analyses confirmed their classification as subtype Bavaria, which is further subdivided into four geographical clusters. The entire N protein, harbouring an amino-terminal hexahistidine tag, of the Bavarian strain was produced in yeast Saccharomyces cerevisiae and showed a slightly different reactivity with N-specific monoclonal antibodies, compared to the yeast-expressed N protein of the PUUV strain Vranica/Hällnäs. Endpoint titration of human sera from different parts of Germany and from Finland revealed only very slight differences in the diagnostic value of the different recombinant proteins. Based on the novel N antigen indirect and monoclonal antibody capture IgG-ELISAs were established. By using serum panels from Germany and Finland their validation demonstrated a high sensitivity and specificity. In summary, our investigations demonstrated the Bavarian PUUV strain to be genetically divergent from other PUUV strains and the potential of its N protein for diagnostic applications.

Enhanced AFLP genome scans detect local adaptation in high-altitude populations of a small rodent (Microtus arvalis)

Adaptation to adverse environmental conditions such as high altitude requires physiological and/or morphological changes. Genome scans provide a means to identify the genetic basis of such adaptations without previous knowledge about the particular genetic variants or traits under selection. In this study, we scanned 3027 amplified fragment length polymorphisms (AFLP) in four populations of the common vole Microtus arvalis for loci associated with local adaptation and high altitude. We investigated voles from two populations at high elevation (~2000 m a.s.l.) representing the upper limit of the altitudinal distribution of the species and two geographically close low-altitude populations (<600 m a.s.l.). Statistical analysis incorporated a new Bayesian F(ST) outlier approach specifically developed for AFLP markers, which considers the intensity of AFLP bands instead of mere presence/absence and allows to derive population-based estimates of allele frequencies and F(IS) values. Computer simulations showed that this approach increases the statistical power of the detection of AFLP markers under selection almost to the power of single nucleotide polymorphism (SNP) data without compromising specificity. Our enhanced genome scan resulted in 20 prime candidate markers for positive selection, which show mostly extremely high allele frequency differences between the low- and high-altitude populations. The comparison of global- and pairwise-enhanced genome scans demonstrated further that very strong selective signatures may also be associated with single populations suggesting the importance of local adaptation in alpine populations of common voles.

Novel hepatitis E virus genotype in Norway rats, Germany

Human hepatitis E virus infections may be caused by zoonotic transmission of virus genotypes 3 and 4. To determine whether rodents are a reservoir, we analyzed the complete nucleotide sequence of a hepatitis E–like virus from 2 Norway rats in Germany. The sequence suggests a separate genotype for this hepatotropic virus.

Genomic scans support repetitive continental colonization events during the rapid radiation of voles (Rodentia: Microtus): the utility of AFLPs versus mitochondrial and nuclear sequence markers

Single locus studies might not resolve phylogenetic relationships and the evolutionary history of taxa. The analysis of multiple markers promises higher resolution, and congruence among loci may indicate that the phylogenies represent the underlying species history. Here, we examine the utility of a genome-wide approach based on amplified fragment length polymorphisms (AFLP) and several DNA sequence markers in resolving phylogenetic signals in the rapidly radiating rodent genus Microtus which produced about 70 vole species within the last 1.2-2 myr. The current Holarctic distribution of Microtus is assumed to have resulted from three independent colonization events out of Asia to North America, Europe, and northern Asia without subsequent colonization, which would have led to deep splits between species from different continents. We investigated this hypothesis of three single colonization events by reconstructing the phylogenetic relationships among species from all three continents based on data from the first exon of the nuclear arginine vasopressin receptor 1a gene (EXON1), an adjacent noncoding region and the mitochondrial cytochrome b gene. The phylogenetic patterns obtained from these sequence markers are contrasted to genome-wide data on more than 1800 amplified fragment length polymorphisms (AFLP) analyzed for the same samples. Our results show that the single sequence markers partially resolve the phylogenetic relationships within Microtus, but with some incongruence mostly between EXON1 and the other loci. However, deeper nodes of the radiation are only weakly supported and neither the combination of the markers nor additional nuclear sequences improved the resolution significantly. AFLPs provided much stronger support for major continent-specific clades, and show also that reciprocal monophyly of American and European voles is incomplete. Our results demonstrate that Microtus voles colonized the American and European continents each repeatedly in several independent events on similar colonization routes during their radiation. More generally, this study supports the suitability of AFLPs as an alternative to sequence markers to resolve the evolutionary history of rapidly radiating taxa.

Integrating phylogeographic patterns of microsatellite and mtDNA divergence to infer the evolutionary history of chamois (genus Rupicapra)

Background

The chamois, distributed over most of the medium to high altitude mountain ranges of southern Eurasia, provides an excellent model for exploring the effects of historical and evolutionary events on diversification. Populations have been grouped into two species, Rupicapra pyrenaica from southwestern Europe and R. rupicapra from eastern Europe. However, a previous study of cytochrome b revealed that the two proposed species were non-monophyletic. The reconstruction of phylogenetic relationships between animal species often depends on the markers studied. To further elucidate the evolutionary history of chamois, we extended earlier studies by analysing DNA sequences of four mitochondrial regions (ND1, 12S, tRNApro and Control Region) and microsatellites (20 loci) to include all subspecies and cover its entire distribution range.

Results

We found discordant microsatellite (μsat) and mitochondrial (mt) DNA phylogenies. Mitochondrial phylogenies form three clades, West, Central and East (mtW, mtC and mtE), at variance with taxonomic classification. Our divergence age estimates indicate an initial separation into branches mtW-mtC and mtE 1.7 million years ago (mya), in the late Pliocene-early Pleistocene, quickly followed by the split of clades mtW and mtC. Clade mtW contains haplotypes from the Iberian peninsula and the western Alps, Clade mtC includes haplotypes from the Apennines and the Massif of Chartreuse and Clade mtE comprises populations to the east of the Alps. Divergence among populations within these three major clades is recent (< 0.5 mya). New microsatellite multilocus genotypes added to previously published data revealed differences between every pair of subspecies, forming three well defined groups (μsatW, μsatC and μsatE) also with a strong geographic signature. Grouping does not correspond with the mitochondrial lineages but is closer to morphology and taxonomic classification. Recent drastic reductions in population size can be noted for the subspecies ornata as an extremely low diversity.

Conclusions

The phylogeographic patterns for mtDNA and microsatellites suggest an evolutionary history with limited range contractions and expansions during the Quaternary period and reflect a major effect of the Alpine barrier on west-east differentiation. The contrasting phylogenies for mtDNA and microsatellites indicate events of hybridization among highly divergent lineages in the central area of distribution. Our study points to the importance of reticulate evolution, with periods of isolation and reduction of population size followed by expansions and hybridizations, in the diversification at the level of close species or subspecies.

Recognizing the temporal distinctions between landscape genetics and phylogeography

The steadily advancing fields of landscape genetics and phylogeography share many goals. However, there are some very distinct differences between these two disciplines, including the kinds of data and analyses commonly used, the timescale over which these data are informative, and the hypotheses the data are used to examine. Recently, a number of studies appear to have confused or synonymized phylogeography and landscape genetics. The difference is not merely semantic; understanding the distinctions between these fields is important for ensuring that researchers are aware of the temporal scale over which their data are informative.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

BACKGROUND

The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations.

RESULTS

Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females.

CONCLUSION

ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Extensive host sharing of central European Tula virus

To examine the host association of Tula virus (TULV), a hantavirus present in large parts of Europe, we investigated a total of 791 rodents representing 469 Microtus arvalis and 322 Microtus agrestis animals from northeast, northwest, and southeast Germany, including geographical regions with sympatric occurrence of both vole species, for the presence of TULV infections. Based on serological investigation, reverse transcriptase PCR, and subsequent sequence analysis of partial small (S) and medium (M) segments, we herein show that TULV is carried not only by its commonly known host M. arvalis but also frequently by M. agrestis in different regions of Germany for a prolonged time period. At one trapping site, TULV was exclusively detected in M. agrestis, suggesting an isolated transmission cycle in this rodent reservoir separate from spillover infections of TULV-carrying M. arvalis. Phylogenetic analysis of the S and M segment sequences demonstrated geographical clustering of the TULV sequences irrespective of the host, M. arvalis or M. agrestis. The novel TULV lineages from northeast, northwest, and southeast Germany described here are clearly separated from each other and from other German, European, or Asian lineages, suggesting their stable geographical localization and fast sequence evolution. In conclusion, these results demonstrate that TULV represents a promiscuous hantavirus with a large panel of susceptible hosts. In addition, this may suggest an alternative evolution mode, other than a strict coevolution, for this virus in its Microtus hosts, which should be proven in further large-scale investigations on sympatric Microtus hosts.