Genomic signatures of climate adaptation in bank voles.

Evidence for divergent selection and adaptive variation across the landscape can provide insight into a species' ability to adapt to different environments. However, despite recent advances in genomics, it remains difficult to detect the footprints of climate-mediated selection in natural populations. Here, we analysed ddRAD sequencing data (21,892 SNPs) in conjunction with geographic climate variation to search for signatures of adaptive differentiation in twelve populations of the bank vole (Clethrionomys glareolus) distributed across Europe. To identify the loci subject to selection associated with climate variation, we applied multiple genotype-environment association methods, two univariate and one multivariate, and controlled for the effect of population structure. In total, we identified 213 candidate loci for adaptation, 74 of which were located within genes. In particular, we identified signatures of selection in candidate genes with functions related to lipid metabolism and the immune system. Using the results of redundancy analysis, we demonstrated that population history and climate have joint effects on the genetic variation in the pan-European metapopulation. Furthermore, by examining only candidate loci, we found that annual mean temperature is an important factor shaping adaptive genetic variation in the bank vole. By combining landscape genomic approaches, our study sheds light on genome-wide adaptive differentiation and the spatial distribution of variants underlying adaptive variation influenced by local climate in bank voles.

Local adaptation and future climate vulnerability in a wild rodent.

As climate change continues, species pushed outside their physiological tolerance limits must adapt or face extinction. When change is rapid, adaptation will largely harness ancestral variation, making the availability and characteristics of that variation of critical importance. Here, we used whole-genome sequencing and genetic-environment association analyses to identify adaptive variation and its significance in the context of future climates in a small Palearctic mammal, the bank vole (Clethrionomys glareolus). We found that peripheral populations of bank vole in Britain are already at the extreme bounds of potential genetic adaptation and may require an influx of adaptive variation in order to respond. Analyses of adaptive loci suggest regional differences in climate variables select for variants that influence patterns of population adaptive resilience, including genes associated with antioxidant defense, and support a pattern of thermal/hypoxic cross-adaptation. Our findings indicate that understanding potential shifts in genomic composition in response to climate change may be key to predicting species' fate under future climates.

Diverse susceptibilities and responses of human and rodent cells to orthohantavirus infection reveal different levels of cellular restriction.

Orthohantaviruses are rodent-borne emerging viruses that may cause severe diseases in humans but no apparent pathology in their small mammal reservoirs. However, the mechanisms leading to tolerance or pathogenicity in humans and persistence in rodent reservoirs are poorly understood, as is the manner in which they spread within and between organisms. Here, we used a range of cellular and molecular approaches to investigate the interactions of three different orthohantaviruses-Puumala virus (PUUV), responsible for a mild to moderate form of hemorrhagic fever with renal syndrome in humans, Tula virus (TULV) with low pathogenicity, and non-pathogenic Prospect Hill virus (PHV)-with human and rodent host cell lines. Besides the fact that cell susceptibility to virus infection was shown to depend on the cell type and virus strain, the three orthohantaviruses were able to infect Vero E6 and HuH7 human cells, but only the former secreted infectious particles. In cells derived from PUUV reservoir, the bank vole (Myodes glareolus), PUUV achieved a complete viral cycle, while TULV did not enter the cells and PHV infected them but did not produce infectious particles, reflecting differences in host specificity. A search for mature virions by electron microscopy (EM) revealed that TULV assembly occurred in part at the plasma membrane, whereas PHV particles were trapped in autophagic vacuoles in cells of the heterologous rodent host. We described differential interactions of orthohantaviruses with cellular factors, as supported by the cellular distribution of viral nucleocapsid protein with cell compartments, and proteomics identification of cellular partners. Our results also showed that interferon (IFN) dependent gene expression was regulated in a cell and virus species dependent manner. Overall, our study highlighted the complexity of the host-virus relationship and demonstrated that orthohantaviruses are restricted at different levels of the viral cycle. In addition, the study opens new avenues to further investigate how these viruses differ in their interactions with cells to evade innate immunity and how it depends on tissue type and host species.

Genic distribution modelling predicts adaptation of the bank vole to climate change.

The most likely pathway for many species to survive future climate change is by pre-existing trait variation providing a fitness advantage under the new climate. Here we evaluate the potential role of haemoglobin (Hb) variation in bank voles under future climate change. We model gene-climate relationships for two functionally distinct Hb types, HbS and HbF, which have a north-south distribution in Britain presenting an unusually tractable system linking genetic variation in physiology to geographical and temporal variation in climate. Projections to future climatic conditions suggest a change in relative climatic suitability that would result in HbS being displaced by HbF in northern Britain. This would facilitate local adaptation to future climate-without Hb displacement, populations in northern Britain would likely be suboptimally adapted because their Hb would not match local climatic conditions. Our study shows how pre-existing physiological differences can influence the adaptive capacity of species to climate change.

Niche differentiation in a postglacial colonizer, the bank vole Clethrionomys glareolus.

Species-level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species-wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid-Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.

A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole.

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present-day populations. In Europe, where both Mediterranean and extra-Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra-Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range-wide colonization. The results provide clear evidence that both extra-Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid-latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%-47%) and the Carpathian (53%-72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end-glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end-glacial colonization history of this well-studied organism and underscores the importance of genomic data in phylogeographic interpretation.

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end-glacial refugia.

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land-bridge route deriving from a "Carpathian" glacial refugium and one via a north-eastern route from an "Eastern" glacial refugium near the Ural Mountains. Clustering of genome-wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red-backed voles (Myodes rutilus) took place in Fennoscandia just after end-glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.

Playing Hide-and-Seek in Beta-Globin Genes: Gene Conversion Transferring a Beneficial Mutation between Differentially Expressed Gene Duplicates.

Increasing evidence suggests that adaptation to diverse environments often involves selection on existing variation rather than new mutations. A previous study identified a nonsynonymous single nucleotide polymorphism (SNP) in exon 2 of two paralogous ß-globin genes of the bank vole (Clethrionomys glareolus) in Britain in which the ancestral serine (Ser) and the derived cysteine (Cys) allele represent geographically partitioned functional variation affecting the erythrocyte antioxidative capacity. Here we studied the geographical pattern of the two-locus Ser/Cys polymorphism throughout Europe and tested for the geographic correlation between environmental variables and allele frequency, expected if the polymorphism was under spatially heterogeneous environment-related selection. Although bank vole population history clearly is important in shaping the dispersal of the oxidative stress protective Cys allele, analyses correcting for population structure suggest the Europe-wide pattern is affected by geographical variation in environmental conditions. The ß-globin phenotype is encoded by the major paralog HBB-T1 but we found evidence of bidirectional gene conversion of exon 2 with the low-expression paralog HBB-T2. Our data support the model where gene conversion reshuffling genotypes between high- and low- expressed paralogs enables tuning of erythrocyte thiol levels, which may help maintain intracellular redox balance under fluctuating environmental conditions. Therefore, our study suggests a possible role for gene conversion between differentially expressed gene duplicates as a mechanism of physiological adaptation of populations to new or changing environments.

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

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

Phylogenetic evidence for a new species of Barbus in the Danube River basin.

Three species of small-sized rheophilic Barbus fishes are endemic to and widely distributed throughout the mountain regions in the Danube River basin. In Hungary, barbels referred to as B. petenyi occur in streams in the foothills of the Carpathians near the borders with Slovakia, Ukraine and Romania. However, up to now, no genetic investigations were carried out on rheophilic barbels in this region. This study aims to clarify the taxonomic identity and distribution of the rheophilic barbels in the Hungarian plain based on molecular and morphological analyses. Two mitochondrial genes (cytochrome b, ATPase 6/8) and one nuclear gene (beta-actin intron 2) were sequenced and several morphometric and meristic characters were recorded. Phylogenetic and morphological analyses revealed that there are four genetically distinct lineages among the rheophilic barbels in the Carpathian Basin. The results demonstrated that North-Hungarian Barbus populations belong to B. carpathicus and that B. petenyi presumably does not occur in Hungary. As expected, B. balcanicus was only recorded in samples from the Balkans analyzed for reference. A distinct species, new to science, was discovered to be present in Sebes-Körös River (Crisul Repede) in eastern Hungary and western Romania and is formally described here as B. biharicus Antal, László, Kotlík - sp. nov.

The complete mitochondrial genome of the bank vole Clethrionomys glareolus (Rodentia: Arvicolinae).

We present the first complete sequence of the bank vole (Clethrionomys glareolus) mitochondrial genome (GenBank accession no. KF918859). The bank vole mitogenome is 16,353 base pairs long and shows the gene content, genome architecture and gene strand asymmetry typical for mammals. The sequence provides an important new genomic resource for the bank vole, which is a popular study species in ecological and evolutionary research.

Mapping 3' transcript ends in the bank vole (Clethrionomys glareolus) mitochondrial genome with RNA-Seq.

BACKGROUND: Although posttranscriptional modification of mitochondrial (mt) transcripts plays key roles in completion of the coding information and in the expression of mtDNA-encoded genes, there is little experimental evidence on the polyadenylation status and the location of mt gene poly(A) sites for non-human mammals. RESULTS: Poly(A)-enriched RNA-Seq reads collected for two wild-caught bank voles (Clethrionomys glareolus) were mapped to the complete mitochondrial genome of that species. Transcript polyadenylation was detected as unmapped adenine residues at the ends of the mapped reads. Where the tRNA punctuation model applied, there was the expected polyadenylation, except for the nad5 transcript, whose polyadenylated 3' end is at an intergenic sequence/cytochrome b boundary. As in human, two pairs of bank vole genes, nad4l/nad4 and atp8/atp6, are expressed from bicistronic transcripts. TAA stop codons of four bank vole protein-coding genes (nad1, atp6, cox3 and nad4) are incompletely encoded in the DNA and are completed by polyadenylation. This is three genes (nad2, nad3 and cob) less than in human. The bank vole nad2 gene encodes a full stop codon (TAA in one vole and TAG in the other), which is followed by a 2 bp UTR and the gene conforms to the tRNA punctuation model. In contrast, the annotations of the reference mouse and some other rodent mt genomes in GenBank include complete TAG stop codons in both nad1 and nad2, which overlap downstream trnI and trnW, respectively. Thus the RNA-Seq data of bank voles provides a model for stop codons of mt-encoded genes in mammals comparable to humans, but at odds with some of the interpretation based purely on genomic data in mouse and other rodents. CONCLUSIONS: This work demonstrates how RNA-Seq data were useful to recover mtDNA transcriptome data in a non-model rodent and to shed more light on mammalian mtDNA transcriptome and post-transcriptional modification. Even though gene content and organisation of mtDNA are strongly conserved among mammals, annotations that neglect the transcriptome may be prone to errors in relation to the stop codons.

Mitogenomic phylogenetics of the bank vole Clethrionomys glareolus, a model system for studying end-glacial colonization of Europe.

We have revisited the mtDNA phylogeny of the bank vole Clethrionomys glareolus based on Sanger and next-generation Illumina sequencing of 32 complete mitochondrial genomes. The bank vole is a key study species for understanding the response of European fauna to the climate change following the Last Glacial Maximum (LGM) and one of the most convincing examples of a woodland mammal surviving in cryptic northern glacial refugia in Europe. The genomes sequenced included multiple representatives of each of the eight bank vole clades previously described based on cytochrome b (cob) sequences. All clades with the exception of the Basque - likely a misidentified pseudogene clade - were highly supported in all phylogenetic analyses and the relationships between the clades were resolved with high confidence. Our data extend the distribution of the Carpathian clade, the marker of a northern glacial refugium in the Carpathian Mountains, to include Britain and Fennoscandia (but not adjacent areas of continental Europe). The Carpathian sub-clade that colonized Britain and Fennoscandia had a somewhat different history from the sub-clade currently found in or close to the Carpathians and may have derived from a more north-westerly refugial area. The two bank vole populations that colonized Britain at the end of the last glaciation are for the first time linked with particular continental clades, the first colonists with the Carpathian clade and the second colonists with the western clade originating in a more southerly refugium in the vicinity of the Alps. We however found no evidence that a functional divergence of proteins encoded in the mitochondrial genome promoted the partial genetic replacement of the first colonists by the second colonists detected previously in southern Britain. We did identify one codon site that changed more often and more radically in the tree than expected and where the observed amino acid change may affect the reductase activity of the cytochrome bc1 complex, but the change was not specific to a particular clade. We also found an excess of radical changes to the primary protein structure for geographically restricted clades from southern Italy and Norway, respectively, possibly related to stronger selective pressure at the latitudinal extremes of the bank vole distribution. However, overall, we find little evidence of pervasive effects of deviation from neutrality on bank vole mtDNA phylogeography.

Speciation history and widespread introgression in the European short-call tree frogs (Hyla arborea sensu lato, H. intermedia and H. sarda).

European tree frogs (Hyla) characterized by short temporal parameters of the advertisement call form six genetically differentiated but morphologically cryptic taxa, H. arborea sensu stricto, H. orientalis and H. molleri from across Europe to western Asia (together referred to as H. arborea sensu lato), two putative taxa within H. intermedia (Northern and Southern) from the Italian Peninsula and Sicily, and H. sarda from Sardinia and Corsica. Here, we assess species limits and phylogenetic relationships within these 'short-call tree frogs' based on mitochondrial DNA and nuclear protein-coding markers. The mitochondrial and nuclear genes show partly incongruent phylogeographic patterns, which point to a complex history of gene flow across taxa, particularly in the Balkans. To test the species limits in the short-call tree frogs and to infer the species tree, we used coalescent-based approaches. The monophyly of H. arborea sensu lato is supported by the mtDNA as well as by the all-gene species tree. The Northern and Southern lineages of H. intermedia have been connected by nuclear gene flow (despite their deep mtDNA divergence) and should be treated as conspecific. On the contrary, the parapatric taxa within H. arborea sensu lato should be considered distinct species (H. arborea, H. orientalis, H. molleri) based on the coalescent analysis, although signs of hybridization were detected between them (H. arborea×H. orientalis; H. arborea×H. molleri). A mitochondrial capture upon secondary contact appears to explain the close mtDNA relationship between the geographically remote Iberian H. molleri and H. orientalis from around the Black Sea. Introgressive hybridization occurred also between the Balkan H. arborea and northern Italian H. intermedia, and between the Minor Asiatic H. orientalis and Arabian H. felix arabica (the latter belonging to a different acoustic group/clade). Our results shed light on the species limits in the European short-call tree frogs and show that introgression played an important role in the evolutionary history of the short-call tree frogs and occurred even between taxa supported as distinct species.

Adaptive phylogeography: functional divergence between haemoglobins derived from different glacial refugia in the bank vole.

Over the years, researchers have used presumptively neutral molecular variation to infer the origins of current species' distributions in northern latitudes (especially Europe). However, several reported examples of genic and chromosomal replacements suggest that end-glacial colonizations of particular northern areas may have involved genetic input from different source populations at different times, coupled with competition and selection. We investigate the functional consequences of differences between two bank vole (Clethrionomys glareolus) haemoglobins deriving from different glacial refugia, one of which partially replaced the other in Britain during end-glacial climate warming. This allows us to examine their adaptive divergence and hence a possible role of selection in the replacement. We determine the amino acid substitution Ser52Cys in the major expressed ß-globin gene as the allelic difference. We use structural modelling to reveal that the protein environment renders the 52Cys thiol a highly reactive functional group and we show its reactivity in vitro. We demonstrate that possessing the reactive thiol in haemoglobin increases the resistance of bank vole erythrocytes to oxidative stress. Our study thus provides striking evidence for physiological differences between products of genic variants that spread at the expense of one another during colonization of an area from different glacial refugia.

Back from the brink: the Holocene history of the Carpathian barbel Barbus carpathicus.

As a result of specific adaptations and habitat preferences strongly rheophilic fish species may show high levels of endemism. Many temperate rheophilic fish species were subjected to a series of range contractions during the Pleistocene, and then successfully expanded during the Holocene, colonising previously abandoned areas. The Carpathian barbel (Barbus carpathicus Kotlík, Tsigenopoulos, Ráb et Berrebi 2002) occurs in the montane streams in three basins of the main Central European rivers in the northern part of the Carpathian range. We used genetic variation within 3 mitochondrial and 9 microsatellite loci to determine a pattern of postglacial expansion in B. carpathicus. We found that overall genetic variation within the species is relatively low. Estimate of time to the most recent common ancestor (tMRCA) of mitochondrial sequences falls within the Holocene. The highest levels of genetic variation found in upper reaches of the Tisa river in the Danube basin suggest that glacial refugia were located in the south-eastern part of the species range. Our data suggest that the species crossed different watersheds at least six times as three genetically distinct groups (probably established in different expansion episodes) were found in northern part of the species range. Clines of genetic variation were observed in both the Danube and Vistula basins, which probably resulted from subsequent bottlenecks while colonizing successive habitats (south eastern populations) or due to the admixture of genetically diverse individuals to a previously uniform population (Vistula basin). Therefore, B. carpathicus underwent both demographic breakdowns and expansions during the Holocene, showing its distribution and demography are sensitive to environmental change. Our findings are important in the light of the current human-induced habitats alterations.

Mitochondrial capture misleads about ecological speciation in the Daphnia pulex complex.

The North American ecological species Daphniapulicaria and Daphniapulex are thought to have diverged from a common ancestor by adaptation to sympatric but ecologically distinct lake and pond habitats respectively. Based on mtDNA relationships, European D. pulicaria is considered a different species only distantly related to its North American counterpart, but both species share a lactate dehydrogenase (Ldh) allele F supposedly involved in lake adaptation in North America, and the same allele is also carried by the related Holarctic Daphniatenebrosa. The correct inference of the species' ancestral relationships is therefore critical for understanding the origin of their adaptive divergence. Our species tree inferred from unlinked nuclear loci for D. pulicaria and D. pulex resolved the European and North American D. pulicaria as sister clades, and we argue that the discordant mtDNA gene tree is best explained by capture of D. pulex mtDNA by D. pulicaria in North America. The Ldh gene tree shows that F-class alleles in D. pulicaria and D. tenebrosa are due to common descent (as opposed to introgression), with D. tenebrosa alleles paraphyletic with respect to D. pulicaria alleles. That D. tenebrosa still segregates the ancestral and derived amino acids at the two sites distinguishing the pond and lake alleles suggests that D. pulicaria inherited the derived states from the D. tenebrosa ancestry. Our results suggest that some adaptations restricting the gene flow between D. pulicaria and D. pulex might have evolved in response to selection in ancestral environments rather than in the species' current sympatric habitats. The Arctic (D. tenebrosa) populations are likely to provide important clues about these issues.

A mitochondrial DNA phylogeny of the endangered vipers of the Vipera ursinii complex.

The last two populations of the Hungarian meadow viper Vipera ursinii rakosiensis were thought to persist in the steppe fragments of Hungary until meadow vipers were discovered in central Romania (Transylvania), suggesting a possible existence of remnant populations elsewhere. We assessed the phylogenetic position of the Transylvanian vipers using 2030 bp of mitochondrial DNA sequence. We showed that they were closely related to the Hungarian vipers, while those from northeastern Romania (Moldavia) and Danube Delta belonged to the subspecies Vipera ursinii moldavica. Montane subspecies from Europe (Vipera ursinii ursinii and Vipera ursinii macrops) formed a sister clade to the two lowland subspecies. Vipera renardi formed a sister clade to V. ursinii, with populations from the Greater Caucasus (Vipera renardi lotievi) and Tien Shan (Vipera renardi tienshanica) as the sister group to Vipera renardi renardi, and Vipera renardi eriwanensis from the Lesser Caucasus as the most basal taxon in the species. Our results illustrate that the divergence between the lowland and montane populations occurred separately in each species and several times in V. renardi. We demonstrated that the recently discovered Transylvanian population is the third surviving population of V. u. rakosiensis and the only known population outside of Hungary.

Diversity in the reproductive modes of European Daphnia pulicaria deviates from the geographical parthenogenesis.

BACKGROUND: Multiple transitions to obligate parthenogenesis have occurred in the Daphnia pulex complex in North America. These newly formed asexual lineages are differentially distributed being found predominantly at high latitudes. This conforms to the rule of geographical parthenogenesis postulating prevalence of asexuals at high latitudes and altitudes. While the reproductive mode of high-latitude populations is relatively well studied, little is known about the reproduction mode in high altitudes. This study aimed to assess the reproductive mode of Daphnia pulicaria, a species of the D. pulex complex, from high altitude lakes in Europe. METHODOLOGY/PRINCIPAL FINDINGS: Variation at eight microsatellite loci revealed that D. pulicaria from the High Tatra Mountains (HTM) had low genotype richness and showed excess of heterozygotes and significant deviations from Hardy-Weinberg expectations, and was thus congruent with reproduction by obligate parthenogenesis. By contrast, populations from the Pyrenees (Pyr) were generally in Hardy-Weinberg equilibrium and had higher genotypic richness, suggesting that they are cyclic parthenogens. Four lakes from lowland areas (LLaP) had populations with an uncertain or mixed breeding mode. All D. pulicaria had mtDNA ND5 haplotypes of the European D. pulicaria lineage. Pyr were distinct from LLaP and HTM at the ND5 gene. By contrast, HTM shared two haplotypes with LLaP and one with Pyr. Principal Coordinate Analysis of the microsatellite data revealed clear genetic differentiation into three groups. HTM isolates were intermediate to Pyr and LLaP, congruent with a hybrid origin. CONCLUSION/SIGNIFICANCE: Inferred transitions to obligate parthenogenesis have occurred only in HTM, most likely as a result of hybridizations. In contrast to North American populations, these transitions do not appear to involve meiosis suppressor genes and have not been accompanied by polyploidy. The absence of obligate parthenogenesis in Pyr, an environment highly similar to the HTM, may be due to the lack of opportunities for hybridization.

Human-aided dispersal has altered but not erased the phylogeography of the tench.

Human-aided dispersal can result in phylogeographic patterns that do not reflect natural historical processes, particularly in species prone to intentional translocations by humans. Here, we use a multiple-gene sequencing approach to assess the effects of human-aided dispersal on phylogeography of the tench Tinca tinca, a widespread Eurasian freshwater fish with a long history in aquaculture. Spatial genetic analysis applied to sequence data from four unlinked loci and 67 geographic localities (38-382 gene copies per locus) defined two groups of populations that were little structured geographically but were significantly differentiated from each other, and it identified locations of major genetic breaks, which were concordant across genes and were driven by distributions of two phylogroups. This pattern most reasonably reflects isolation in two major glacial refugia and subsequent range expansions, with the Eastern and Western phylogroups remaining largely allopatric throughout the tench range. However, this phylogeographic variation was also present in all 17 cultured breeds studied, and some populations at the western edge of the native range contained the Eastern phylogroup. Thus, natural processes have played an important role in structuring tench populations, but human-aided dispersal has also contributed significantly, with the admixed genetic composition of cultured breeds most likely contributing to the introgression.