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.

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.

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.

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.

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.

Genomic insight into diet adaptation in the biological control agent Cryptolaemus montrouzieri.

BACKGROUND: The ladybird beetle Cryptolaemus montrouzieri Mulsant, 1853 (Coleoptera, Coccinellidae) is used worldwide as a biological control agent. It is a predator of various mealybug pests, but it also feeds on alternative prey and can be reared on artificial diets. Relatively little is known about the underlying genetic adaptations of its feeding habits. RESULTS: We report the first high-quality genome sequence for C. montrouzieri. We found that the gene families encoding chemosensors and digestive and detoxifying enzymes among others were significantly expanded or contracted in C. montrouzieri in comparison to published genomes of other beetles. Comparisons of diet-specific larval development, survival and transcriptome profiling demonstrated that differentially expressed genes on unnatural diets as compared to natural prey were enriched in pathways of nutrient metabolism, indicating that the lower performance on the tested diets was caused by nutritional deficiencies. Remarkably, the C. montrouzieri genome also showed a significant expansion in an immune effector gene family. Some of the immune effector genes were dramatically downregulated when larvae were fed unnatural diets. CONCLUSION: We suggest that the evolution of genes related to chemosensing, digestion, and detoxification but also immunity might be associated with diet adaptation of an insect predator. These findings help explain why this predatory ladybird has become a successful biological control agent and will enable the optimization of its mass rearing and use in biological control programs.

Genetic, phenotypic and ecological differentiation suggests incipient speciation in two Charadrius plovers along the Chinese coast.

BACKGROUND: Speciation with gene flow is an alternative to the nascence of new taxa in strict allopatric separation. Indeed, many taxa have parapatric distributions at present. It is often unclear if these are secondary contacts, e.g. caused by past glaciation cycles or the manifestation of speciation with gene flow, which hampers our understanding of how different forces drive diversification. Here we studied genetic, phenotypic and ecological aspects of divergence in a pair of incipient shorebird species, the Kentish (Charadrius alexandrinus) and the White-faced Plovers (C. dealbatus), shorebirds with parapatric breeding ranges along the Chinese coast. We assessed divergence based on molecular markers with different modes of inheritance and quantified phenotypic and ecological divergence in aspects of morphometric, dietary and climatic niches. RESULTS: Our integrative analyses revealed small to moderate levels of genetic and phenotypic distinctiveness with symmetric gene flow across the contact area at the Chinese coast. The two species diverged approximately half a million years ago in dynamic isolation with secondary contact occurring due to cycling sea level changes between the Eastern and Southern China Sea in the mid-late Pleistocene. We found evidence of character displacement and ecological niche differentiation between the two species, invoking the role of selection in facilitating divergence despite gene flow. CONCLUSION: These findings imply that ecology can indeed counter gene flow through divergent selection and thus contributes to incipient speciation in these plovers. Furthermore, our study highlights the importance of using integrative datasets to reveal the evolutionary history and assist the inference of mechanisms of speciation.

Genomic and spatial variability of a European common vole hepevirus.

Rodents host different orthohepeviruses, namely orthohepevirus C genotype HEV-C1 (rat hepatitis E virus, HEV) and the additional putative genotypes HEV-C3 and HEV-C4. Here, we screened 2,961 rodents from Central Europe by reverse transcription polymerase chain reaction (RT-PCR) and identified HEV RNA in 13 common voles (Microtus arvalis) and one bank vole (Myodes glareolus) with detection rates of 2% (95% confidence interval [CI]: 1-3.4) and 0.08% (95% CI: 0.002-0.46), respectively. Sequencing of a 279-nucleotide RT-PCR amplicon corresponding to a region within open reading frame (ORF) 1 showed a high degree of similarity to recently described common vole-associated HEV (cvHEV) sequences from Hungary. Five novel complete cvHEV genome sequences from Central Europe showed the typical HEV genome organization with ORF1, ORF2 and ORF3 and RNA secondary structure. Uncommon features included a noncanonical start codon in ORF3, multiple insertions and deletions within ORF1 and ORF2/ORF3, and the absence of a putative ORF4. Phylogenetic analysis showed all of the novel cvHEV sequences to be monophyletic, clustering most closely with an unassigned bird-derived sequence and other sequences of the species Orthohepevirus C. The nucleotide and amino acid sequence divergence of the common vole-derived sequences was significantly correlated with the spatial distance between the trapping sites, indicating mostly local evolutionary processes. Detection of closely related HEV sequences in common voles in multiple localities over a distance of 800 kilometers suggested that common voles are infected by cvHEV across broad geographic distances. The common vole-associated HEV strain is clearly divergent from HEV sequences recently found in narrow-headed voles (Microtus gregalis) and other cricetid rodents.

Field vole-associated Traemmersee hantavirus from Germany represents a novel hantavirus species.

Vole-associated hantaviruses occur in the Old and New World. Tula orthohantavirus (TULV) is widely distributed throughout the European continent in its reservoir, the common vole (Microtus arvalis), but the virus was also frequently detected in field voles (Microtus agrestis) and other vole species. TULV and common voles are absent from Great Britain. However, field voles there harbor Tatenale and Kielder hantaviruses. Here we screened 126 field voles and 13 common voles from Brandenburg, Germany, for hantavirus infections. One common vole and four field voles were anti-TULV antibody and/or TULV RNA positive. In one additional, seropositive field vole a novel hantavirus sequence was detected. The partial S and L segment nucleotide sequences were only 61.1% and 75.6% identical to sympatrically occurring TULV sequences, but showed highest similarity of approximately 80% to British Tatenale and Kielder hantaviruses. Subsequent determination of the entire nucleocapsid (N), glycoprotein (GPC), and RNA-dependent RNA polymerase encoding sequences and determination of the pairwise evolutionary distance (PED) value for the concatenated N and GPC amino acid sequences confirmed a novel orthohantavirus species, tentatively named Traemmersee orthohantavirus. The identification of this novel hantavirus in a field vole from eastern Germany underlines the necessity of a large-scale, broad geographical hantavirus screening of voles to understand evolutionary processes of virus-host associations and host switches.

Host-Associated Absence of Human Puumala Virus Infections in Northern and Eastern Germany.

Human hantavirus disease cases, caused by Puumala virus (PUUV), are mainly recorded in western and southern areas of Germany. This bank vole reservoir survey confirmed PUUV presence in these regions but its absence in northern and eastern regions. PUUV occurrence is associated with the presence of the Western bank vole phylogroup.

Revised time scales of RNA virus evolution based on spatial information.

The time scales of pathogen evolution are of major concern in the context of public and veterinary health, epidemiology and evolutionary biology. Dating the emergence of a pathogen often relies on estimates of evolutionary rates derived from nucleotide sequence data. For many viruses, this has yielded estimates of evolutionary origins only a few hundred years in the past. Here we demonstrate through the incorporation of geographical information from virus sampling that evolutionary age estimates of two European hantaviruses are severely underestimated because of pervasive mutational saturation of nucleotide sequences. We detected very strong relationships between spatial distance and genetic divergence for both Puumala and Tula hantavirus-irrespective of whether nucleotide or derived amino acid sequences were analysed. Extrapolations from these relationships dated the emergence of these viruses most conservatively to at least 3700 and 2500 years ago, respectively. Our minimum estimates for the age of these hantaviruses are ten to a hundred times older than results from current non-spatial methods, and in much better accordance with the biogeography of these viruses and their respective hosts. Spatial information can thus provide valuable insights on the deeper time scales of pathogen evolution and improve our understanding of disease emergence.

A highly divergent Puumala virus lineage in southern Poland.

Puumala virus (PUUV) represents one of the most important hantaviruses in Central Europe. Phylogenetic analyses of PUUV strains indicate a strong genetic structuring of this hantavirus. Recently, PUUV sequences were identified in the natural reservoir, the bank vole (Myodes glareolus), collected in the northern part of Poland. The objective of this study was to evaluate the presence of PUUV in bank voles from southern Poland. A total of 72 bank voles were trapped in 2009 at six sites in this part of Poland. RT-PCR and IgG-ELISA analyses detected three PUUV positive voles at one trapping site. The PUUV-infected animals were identified by cytochrome b gene analysis to belong to the Carpathian and Eastern evolutionary lineages of bank vole. The novel PUUV S, M and L segment nucleotide sequences showed the closest similarity to sequences of the Russian PUUV lineage from Latvia, but were highly divergent to those previously found in northern Poland, Slovakia and Austria. In conclusion, the detection of a highly divergent PUUV lineage in southern Poland indicates the necessity of further bank vole monitoring in this region allowing rational public health measures to prevent human infections.

Phylogeographic past and invasive presence of Arion pest slugs in Europe.

Arion vulgaris (syn. A. lusitanicus) is the most destructive pest slug in Europe. The species has been regarded a classic case of an ongoing biological invasion with negative economic and ecological impact in many European countries, but this status has recently been contested. In this study, we assessed mitochondrial and autosomal genetic diversity in populations of A. vulgaris across the entire distribution range in order to characterize its evolutionary history. Mitochondrial diversity in A. vulgaris was strongly reduced compared with the closely related and largely codistributed noninvasive species A. rufus and A. ater, indicating a very rapid spread of A. vulgaris through Europe. Autosomal diversity assessed in 632 individuals from 32 populations decreased towards eastern and northern Europe which is consistent with the reported expansion of the species towards these regions in the last decades. Demographic simulations supported very recent population founding events in most of the European range. The short periods between the first detection of A. vulgaris in different countries and only a very weak association of genetic structuring among populations with geographical distances suggest a human contribution in the ongoing expansion of the slug. We propose that this contribution may ultimately prevent the exact localization of the debated region of origin of A. vulgaris. However, the reclassification of A. vulgaris as noninvasive would be premature. Without counter measures, the Eastern and Northern European countries can expect to see this biological invasion continued in the future.

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.

Ignoring heterozygous sites biases phylogenomic estimates of divergence times: implications for the evolutionary history of microtus voles.

Phylogenetic reconstruction of the evolutionary history of closely related organisms may be difficult because of the presence of unsorted lineages and of a relatively high proportion of heterozygous sites that are usually not handled well by phylogenetic programs. Genomic data may provide enough fixed polymorphisms to resolve phylogenetic trees, but the diploid nature of sequence data remains analytically challenging. Here, we performed a phylogenomic reconstruction of the evolutionary history of the common vole (Microtus arvalis) with a focus on the influence of heterozygosity on the estimation of intraspecific divergence times. We used genome-wide sequence information from 15 voles distributed across the European range. We provide a novel approach to integrate heterozygous information in existing phylogenetic programs by repeated random haplotype sampling from sequences with multiple unphased heterozygous sites. We evaluated the impact of the use of full, partial, or no heterozygous information for tree reconstructions on divergence time estimates. All results consistently showed four deep and strongly supported evolutionary lineages in the vole data. These lineages undergoing divergence processes split only at the end or after the last glacial maximum based on calibration with radiocarbon-dated paleontological material. However, the incorporation of information from heterozygous sites had a significant impact on absolute and relative branch length estimations. Ignoring heterozygous information led to an overestimation of divergence times between the evolutionary lineages of M. arvalis. We conclude that the exclusion of heterozygous sites from evolutionary analyses may cause biased and misleading divergence time estimates in closely related taxa.

Complete genome of a Puumala virus strain from Central Europe.

Puumala virus (PUUV) is one of the predominant hantavirus species in Europe causing mild to moderate cases of haemorrhagic fever with renal syndrome. Parts of Lower Saxony in north-western Germany are endemic for PUUV infections. In this study, the complete PUUV genome sequence of a bank vole-derived tissue sample from the 2007 outbreak was determined by a combined primer-walking and RNA ligation strategy. The S, M and L genome segments were 1,828, 3,680 and 6,550 nucleotides in length, respectively. Sliding-window analyses of the nucleotide sequences of all available complete PUUV genomes indicated a non-homogenous distribution of variability with hypervariable regions located at the 3'-ends of the S and M segments. The overall similarity of the coding genome regions to the other PUUV strains ranged between 80.1 and 84.7 % at the level of the nucleotide sequence and between 89.5 and 98.1 % for the deduced amino acid sequences. In comparison to the phylogenetic trees of the complete coding sequences, trees based on partial segments revealed a general drop in phylogenetic support and a lower resolution. The Astrup strain S and M segment sequences showed the highest similarity to sequences of strains from geographically close sites in the Osnabrück Hills region. In conclusion, a primer-walking-mediated strategy resulted in the determination of the first complete nucleotide sequence of a PUUV strain from Central Europe. Different levels of variability along the genome provide the opportunity to choose regions for analyses according to the particular research question, e.g., large-scale phylogenetics or within-host evolution.