Museums and other natural history collections (NHC) worldwide house millions of specimens. With the advent of molecular genetic approaches these collections have become the source of many fascinating population studies in conservation genetics that contrast historical with present-day genetic diversity. Recent developments in molecular genetics and genomics and the associated statistical tools have opened up the further possibility of studying evolutionary change directly. As we discuss here, we believe that NHC specimens provide a largely underutilized resource for such investigations. However, because DNA extracted from NHC samples is degraded, analyses of such samples are technically demanding and many potential pitfalls exist. Thus, we propose a set of guidelines that outline the steps necessary to begin genetic investigations using specimens from NHC. Museums and other natural history collections (NHC) worldwide house millions of specimens. With the advent of molecular genetic approaches these collections have become the source of many fascinating population studies in conservation genetics that contrast historical with present-day genetic diversity. Recent developments in molecular genetics and genomics and the associated statistical tools have opened up the further possibility of studying evolutionary change directly. As we discuss here, we believe that NHC specimens provide a largely underutilized resource for such investigations. However, because DNA extracted from NHC samples is degraded, analyses of such samples are technically demanding and many potential pitfalls exist. Thus, we propose a set of guidelines that outline the steps necessary to begin genetic investigations using specimens from NHC. IntroductionGiven that evolution is change over time, documenting and understanding temporal patterns has long been at the heart of evolutionary studies. In disciplines such as palaeontology, inferences about evolutionary processes are drawn from the analyses of temporal patterns in the fossil record. Similarly, our understanding of microevolutionary processes (i.e. changes in gene frequencies over time) has often involved the analyses of records taken over several years; Dobzhanksy's [1] early studies of microevolution among Drosophila used this approach, a tradition that continues among students of this model organism today [2]. However, such microevolutionary studies were often limited to certain taxa and questions because the time available to document temporal changes was often limited to a few generations. How can these limitations be overcome? Long term studies, running over several decades, are one possibility and they are yielding fascinating insights, for example, into the role of reinforcement and character displacement in adaptive radiation and speciation [3,4]. Another approach, which gives longer time series, is to extend the data back in time using well preserved fossil samples or specimens from natural history collections (NHC). Here, we review the use of specimens from NHC for the ...

The red fox (Vulpes vulpes) is one of the best-documented examples of a species that has successfully occupied cities and their suburbs during the last century. The city of Zurich (Switzerland) was colonized by red foxes 15 years ago and the number of recorded individuals has increased steadily since then. Here, we assessed the hypothesis that the fox population within the city of Zurich is isolated from adjacent rural fox populations against the alternative hypothesis that urban habitat acts as a constant sink for rural dispersers. We examined 11 microsatellite loci in 128 foxes from two urban areas, separated by the main river crossing the city, and three adjacent rural areas from the region of Zurich. Mean observed heterozygosity across individuals and the number of detected alleles were lower for foxes collected within the city as compared with their rural conspecifics. Genetic differentiation was significantly lower between rural than between rural and urban populations, and highest value of pairwise FST was recorded between the two urban areas. Our results indicate that the two urban areas were independently founded by a small number of individuals from adjacent rural areas resulting in genetic drift and genetic differentiation between rural and urban fox populations. Population admixture and immigration analysis revealed that urban-rural gene flow was higher than expected from FST statistics. In the five to seven generations since colonization, fox density has dramatically increased. Currently observed levels of migration between urban and rural populations will probably erode genetic differentiation over time.

In natural populations, quantitative trait dynamics often do not appear to follow evolutionary predictions. Despite abundant examples of natural selection acting on heritable traits, conclusive evidence for contemporary adaptive evolution remains rare for wild vertebrate populations, and phenotypic stasis seems to be the norm. This so-called “stasis paradox” highlights our inability to predict evolutionary change, which is especially concerning within the context of rapid anthropogenic environmental change. While the causes underlying the stasis paradox are hotly debated, comprehensive attempts aiming at a resolution are lacking. Here, we apply a quantitative genetic framework to individual-based long-term data for a wild rodent population and show that despite a positive association between body mass and fitness, there has been a genetic change towards lower body mass. The latter represents an adaptive response to viability selection favouring juveniles growing up to become relatively small adults, i.e., with a low potential adult mass, which presumably complete their development earlier. This selection is particularly strong towards the end of the snow-free season, and it has intensified in recent years, coinciding which a change in snowfall patterns. Importantly, neither the negative evolutionary change, nor the selective pressures that drive it, are apparent on the phenotypic level, where they are masked by phenotypic plasticity and a non causal (i.e., non genetic) positive association between body mass and fitness, respectively. Estimating selection at the genetic level enabled us to uncover adaptive evolution in action and to identify the corresponding phenotypic selective pressure. We thereby demonstrate that natural populations can show a rapid and adaptive evolutionary response to a novel selective pressure, and that explicitly (quantitative) genetic models are able to provide us with an understanding of the causes and consequences of selection that is superior to purely phenotypic estimates of selection and evolutionary change.

The amount of nuclear DNA extracted from teeth of 279 individual red fox Vulpes vulpes collected over a period spanning the last three decades was determined by quantitative polymerase chain reaction (PCR). Although teeth were autoclaved during initial collection, 73.8% of extracts contained sufficient DNA concentration (> 5 pg/ micro L) suitable for reliable microsatellite genotyping but the quantity of nuclear DNA decayed significantly over time in a nonlinear pattern. The success of PCR amplification across four examined canine microsatellites over time was dependent on fragment size. By including data from two different tests for human contamination and from frequencies of allelic dropout and false alleles, the methodological constraints of population genetic studies using microsatellite loci amplified from historic DNA are discussed.

Major histocompatibility complex (MHC) antigen-presenting genes are the most variable loci in vertebrate genomes. Host-parasite co-evolution is assumed to maintain the excessive polymorphism in the MHC loci. However, the molecular mechanisms underlying the striking diversity in the MHC remain contentious. The extent to which recombination contributes to the diversity at MHC loci in natural populations is still controversial, and there have been only few comparative studies that make quantitative estimates of recombination rates. In this study, we performed a comparative analysis for 15 different ungulates species to estimate the population recombination rate, and to quantify levels of selection. As expected for all species, we observed signatures of strong positive selection, and identified individual residues experiencing selection that were congruent with those constituting the peptide-binding region of the human DRB gene. However, in addition for each species, we also observed recombination rates that were significantly different from zero on the basis of likelihood-permutation tests, and in other non-quantitative analyses. Patterns of synonymous and non-synonymous sequence diversity were consistent with differing demographic histories between species, but recent simulation studies by other authors suggest inference of selection and recombination is likely to be robust to such deviations from standard models. If high rates of recombination are common in MHC genes of other taxa, re-evaluation of many inference-based phylogenetic analyses of MHC loci, such as estimates of the divergence time of alleles and trans-specific polymorphism, may be required.

Introgression can be an important evolutionary force but it can also lead to species extinction and as such is a crucial issue for species conservation. However, introgression is difficult to detect, morphologically as well as genetically. Hybridization with domestic cats (Felis silvestris catus) is a major concern for the conservation of European wildcats (Felis s. silvestris). The available morphologic and genetic markers for the two Felis subspecies are not sufficient to reliably detect hybrids beyond first generation. Here we present a single nucleotide polymorphism (SNP) based approach that allows the identification of introgressed individuals. Using high-throughput sequencing of reduced representation libraries we developed a diagnostic marker set containing 48 SNPs (Fst > 0.8) which allows the identification of wildcats, domestic cats, their hybrids and backcrosses. This allows assessing introgression rate in natural wildcat populations and is key for a better understanding of hybridization processes.

During the last two centuries, lynx populations have undergone severe declines and extinctions in Europe. The Alpine lynx, once distributed across the whole Alpine arc, became extinct due to direct human prosecution and deprivation of its main prey in the 1930s. Similar to the Iberian lynx Lynx pardinus, its taxonomy has been subject to several controversies. Moreover, knowing the taxonomic status of the Alpine lynx will help to define conservation units of extant lynx populations in Europe. In this study, we investigated two mitochondrial DNA regions in museum specimens (n=15) representing the autochthonous Alpine population and in samples from extant Eurasian lynx Lynx lynx populations in Europe and Asia (n=17). Phylogenetic analysis (cytochrome b, 345 bp) placed the Alpine lynx within the Eurasian lynx lineage. Among all individuals examined, seven different haplotypes (control region, 300 bp) were observed but no unique Alpine haplotype was discovered. Haplotypes of the extinct Alpine population were identical to previously described haplotypes in Scandinavian lynx signifying a recent genetic ancestry with current European populations. Moreover, our genetic data suggest two distinct glacial refugia for the Carpathian and Balkan population. Overall this study demonstrates that historical DNA from extinct populations can help to disentangle the phylogenetic relationships and historical biogeography of taxa with only a limited number of extant populations remaining.Historical DNA reveals the phylogenetic position of the extinct Alpine lynx AbstractDuring the last two centuries, lynx populations have undergone severe declines and extinctions in Europe. The Alpine lynx, once distributed across the whole Alpine arc, became extinct due to direct human prosecution and deprivation of its main prey in the 1930s. Similar to the Iberian lynx Lynx pardinus, its taxonomy has been subject to several controversies. Moreover, knowing the taxonomic status of the Alpine lynx will help to define conservation units of extant lynx populations in Europe. In this study, we investigated two mitochondrial DNA regions in museum specimens (n = 15) representing the autochthonous Alpine population and in samples from extant Eurasian lynx Lynx lynx populations in Europe and Asia (n = 17). Phylogenetic analysis (cytochrome b, 345 bp) placed the Alpine lynx within the Eurasian lynx lineage. Among all individuals examined, seven different haplotypes (control region, 300 bp) were observed but no unique Alpine haplotype was discovered. Haplotypes of the extinct Alpine population were identical to previously described haplotypes in Scandinavian lynx signifying a recent genetic ancestry with current European populations. Moreover, our genetic data suggest two distinct glacial refugia for the Carpathian and Balkan population. Overall this study demonstrates that historical DNA from extinct populations can help to disentangle the phylogenetic relationships and historical biogeography of taxa with only a limited number of extant population...

Genotyping non-invasively collected samples is challenging. Nevertheless, genetic monitoring of elusive species like the European wildcat (Felis silvestris silvestris) mainly relies on such samples. Wildcats are likely threatened through introgression with domestic cats (F. silvestris catus). To determine introgression based on single cat hairs, we developed a 96.96 Fluidigm single nucleotide polymorphism (SNP) genotyping array chip. To estimate the accuracy of this method, we compared genotypes of 17 cats called with both Sanger sequencing and Fluidigm. When Sanger sequencing genotypes were considered as a reference, the genotyping error rate with Fluidigm was 0.9 %. We subsequently compared 16 hair samples to tissue samples of the same individual. When the tissue samples were used as a reference, the genotyping error rate in hair samples was 1.6 %. This low error rate allowed reliable recognition of individuals and correct assessment of introgression levels. Thus, the genotyping method presented in this paper is suitable for non-invasively collected samples. It will help conservationists to monitor the introgression rate in wildcat populations based on non-invasive hair sampling and subsequently to conduct effective conservation measures.