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.

• Klíčová slova
• Clethrionomys glareolus, climate gradient, genomic analysis, local adaptations, rodent,
• Publikační typ
• časopisecké články MeSH

BACKGROUND: Laelaps agilis C.L. Koch, 1836 is one the most abundant and widespread parasitic mite species in the Western Palearctic. It is a permanent ectoparasite associated with the Apodemus genus, which transmits Hepatozoon species via the host's blood. Phylogenetic relationships, genealogy and host specificity of the mite are uncertain in the Western Palearctic. Here, we investigated the population genetic structure of 132 individual mites across Europe from their Apodemus and Clethrionomys hosts. Phylogenetic relationships and genetic variation of the populations were analyzed using cytochrome c oxidase subunit I (COI) gene sequences. RESULTS: We recovered three main mtDNA lineages within L. agilis in the Western Palearctic, which differentiated between 1.02 and 1.79 million years ago during the Pleistocene period: (i) Lineage A, including structured populations from Western Europe and the Czech Republic, (ii) Lineage B, which included only a few individuals from Greece and the Czech Republic; and (iii) Lineage C, which comprised admixed populations from Western and Eastern Europe. Contrary to their population genetic differentiation, the lineages did not show signs of specificity to different hosts. Finally, we confirmed that the sympatric congener L. clethrionomydis is represented by a separated monophyletic lineage. CONCLUSION: Differences in the depth of population structure between L. agilis Lineages A and C, corroborated by the neutrality tests and demographic history analyses, suggested a stable population size in the structured Lineage A and a rapid range expansion for the geographically admixed Lineage C. We hypothesized that the two lineages were associated with hosts experiencing different glaciation histories. The lack of host specificity in L. agilis lineages was in contrast to the co-occurring highly host-specific lineages of Polyplax serrata lice, sharing Apodemus hosts. The incongruence was attributed to the differences in mobility between the parasites, allowing mites to switch hosts more often.

• Klíčová slova
• Demographic history, Host specificity, Laelaps agilis, Phylogenetic relationship,
• Publikační typ
• časopisecké články MeSH

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.

• Klíčová slova
• Myodes glareolus, SNP, admixture, genotyping‐by‐sequencing, postglacial colonization,
• Publikační typ
• časopisecké články MeSH

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.

• Klíčová slova
• Last Glacial Maximum, MaxEnt, Myodesglareolus, cryptic refugia, ecological niche modeling, intraspecific variation,
• Publikační typ
• časopisecké články MeSH

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.

• Klíčová slova
• Chi motif, adaptive phylogeography, antioxidative capacity, cysteine, environmental selection, gene conversion,
• Publikační typ
• časopisecké články MeSH

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.

• Klíčová slova
• Clethrionomys glareolus, Myodes glareolus, approximate Bayesian computation, genome admixture, postglacial colonization, single-nucleotide polymorphism,
• MeSH
• Arvicolinae genetika   fyziologie   MeSH
• fylogeneze MeSH
• genom * MeSH
• jednonukleotidový polymorfismus * MeSH
• migrace zvířat * MeSH
• rozšíření zvířat * MeSH
• sekvenční analýza DNA MeSH
• zvířata MeSH
• Check Tag
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Geografické názvy
• Anglie MeSH
• Skotsko MeSH
• Wales MeSH

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

• Klíčová slova
• Central Europe, Hantavirus, Homologous recombination, Myodes glareolus, Puumala orthohantavirus, Reassortment,
• MeSH
• Arvicolinae virologie   MeSH
• fylogeneze MeSH
• genotyp MeSH
• hantavirové infekce virologie   MeSH
• Hantavirus genetika   MeSH
• lidé MeSH
• molekulární evoluce MeSH
• RNA virová genetika   MeSH
• virus Puumala genetika   MeSH
• zvířata MeSH
• Check Tag
• lidé MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• Geografické názvy
• Česká republika MeSH
• Evropa MeSH
• Německo MeSH
• Slovenská republika MeSH
• Názvy látek
• RNA virová MeSH

The major histocompatibility complex (MHC) plays a central role in the adaptive immune response and is the most polymorphic gene family in vertebrates. Although high-throughput sequencing has increasingly been used for genotyping families of co-amplifying MHC genes, its potential to facilitate early steps in the characterisation of MHC variation in nonmodel organism has not been fully explored. In this study we evaluated the usefulness of de novo transcriptome assembly in characterisation of MHC sequence diversity. We found that although de novo transcriptome assembly of MHC I genes does not reconstruct sequences of individual alleles, it does allow the identification of conserved regions for PCR primer design. Using the newly designed primers, we characterised MHC I sequences in the bank vole. Phylogenetic analysis of the partial MHC I coding sequence (2-4 exons) of the bank vole revealed a lack of orthology to MHC I of other Cricetidae, consistent with the high gene turnover of this region. The diversity of expressed alleles was characterised using ultra-deep sequencing of the third exon that codes for the peptide-binding region of the MHC molecule. High allelic diversity was demonstrated, with 72 alleles found in 29 individuals. Interindividual variation in the number of expressed loci was found, with the number of alleles per individual ranging from 5 to 14. Strong signatures of positive selection were found for 8 amino acid sites, most of which are inferred to bind antigens in human MHC, indicating conservation of structure despite rapid sequence evolution.

• MeSH
• alely MeSH
• Arvicolinae genetika   MeSH
• DNA primery MeSH
• exony MeSH
• fylogeneze MeSH
• genetická variace MeSH
• genotyp MeSH
• geny MHC třídy I * MeSH
• hlavní histokompatibilní komplex genetika   MeSH
• multigenová rodina MeSH
• myši MeSH
• transkriptom * MeSH
• vysoce účinné nukleotidové sekvenování MeSH
• zvířata MeSH
• Check Tag
• myši MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• DNA primery MeSH

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.

• MeSH
• Arvicolinae genetika   MeSH
• genom mitochondriální genetika   MeSH
• mapování chromozomů MeSH
• mitochondriální DNA genetika   MeSH
• sekvenční analýza RNA MeSH
• zvířata MeSH
• Check Tag
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Názvy látek
• mitochondriální DNA MeSH