World Association for the Advancement of Veterinary Parasitology (W.A.A.V.P.): Third edition of the guideline for evaluating efficacy of anthelmintics in ruminants (bovine, ovine, caprine)

This guideline is aimed at those who are involved in the assessment of anthelmintic efficacy in ruminant livestock species (bovine, ovine and caprine). The intent is to provide a framework that can be adopted worldwide for the testing of anthelmintics in ruminants, such that studies carried out in different countries can be compared and thereby unnecessary duplication can be reduced. Recommendations are made for the selection, housing and feeding of study animals, the type of studies required, the method used to conduct those studies, the assessment of results and the standards for defining anthelmintic efficacy.

A Chromosome-Level Genome Assembly and Annotation for the Clouded Apollo Butterfly (Parnassius mnemosyne): A Species of Global Conservation Concern

The clouded apollo (Parnassius mnemosyne) is a palearctic butterfly distributed over a large part of western Eurasia, but population declines and fragmentation have been observed in many parts of the range. The development of genomic tools can help to shed light on the genetic consequences of the decline and to make informed decisions about direct conservation actions. Here, we present a high-contiguity, chromosome-level genome assembly of a female clouded apollo butterfly and provide detailed annotations of genes and transposable elements. We find that the large genome (1.5 Gb) of the clouded apollo is extraordinarily repeat rich (73%). Despite that, the combination of sequencing techniques allowed us to assemble all chromosomes (nc = 29) to a high degree of completeness. The annotation resulted in a relatively high number of protein-coding genes (22,854) compared with other Lepidoptera, of which a large proportion (21,635) could be assigned functions based on homology with other species. A comparative analysis indicates that overall genome structure has been largely conserved, both within the genus and compared with the ancestral lepidopteran karyotype. The high-quality genome assembly and detailed annotation presented here will constitute an important tool for forthcoming efforts aimed at understanding the genetic consequences of fragmentation and decline, as well as for assessments of genetic diversity, population structure, inbreeding, and genetic load in the clouded apollo butterfly.

Batrachochytrium dendrobatidis strain affects transcriptomic response in liver but not skin in latitudinal populations of the common toad (Bufo bufo)

Batrachochytrium dendrobatidis (Bd) is a fungal pathogen that has decimated amphibian populations worldwide for several decades. We examined the changes in gene expression in response to Bd infection in two populations of the common toad, Bufo bufo, in a laboratory experiment. We collected B. bufo eggs in southern and northern Sweden, and infected the laboratory-raised metamorphs with two strains of the global panzoonotic lineage Bd-GPL. Differential expression analysis showed significant differences between infected and control individuals in both liver and skin. The skin samples showed no discernible differences in gene expression between the two strains used, while liver samples were differentiated by strain, with one of the strains eliciting no immune response from infected toads. Immune system genes were overexpressed in skin samples from surviving infected individuals, while in liver samples the pattern was more diffuse. Splitting samples by population revealed a stronger immune response in northern individuals. Differences in transcriptional regulation between populations are particularly relevant to study in Swedish amphibians, which may have experienced varying exposure to Bd. Earlier exposure to this pathogen and subsequent adaptation or selection pressure may contribute to the survival of some populations over others, while standing genetic diversity in different populations may also affect the infection outcome.

Body size mediates latitudinal population differences in the response to chytrid fungus infection in two amphibians

Factors behind intraspecific variation in sensitivity to pathogens remain poorly understood. We investigated how geographical origin in two North European amphibians affects tolerance to infection by the chytrid fungus Batrachochytrium dendrobatidis (Bd), a generalist pathogen which has caused amphibian population declines worldwide. We exposed newly metamorphosed individuals of moor frog Rana arvalis and common toad Bufo bufo from two latitudinal regions to two different BdGPL strains. We measured survival and growth as infections may cause sub-lethal effects in fitness components even in the absence of mortality. Infection loads were higher in B. bufo than in R. arvalis, and smaller individuals had generally higher infection loads. B. bufo had high mortality in response to Bd infection, whereas there was little mortality in R. arvalis. Bd-mediated mortality was size-dependent and high-latitude individuals were smaller leading to high mortality in the northern B. bufo. Bd exposure led to sub-lethal effects in terms of reduced growth suggesting that individuals surviving the infection may have reduced fitness mediated by smaller body size. In both host species, the Swedish Bd strain caused stronger sublethal effects than the British strain. We suggest that high-latitude populations can be more vulnerable to chytrids than those from lower latitudes and discuss the possible mechanisms how body size and host geographical origin contribute to the present results.

Association between the skin microbiome and MHC class II diversity in an amphibian

Microbiomes play an important role in determining the ecology and behaviour of their hosts. However, questions remain pertaining to how host genetics shape microbiomes, and how microbiome composition influences host fitness. We explored the effects of geography, evolutionary history and host genetics on the skin microbiome diversity and structure in a widespread amphibian. More specifically, we examined the association between bacterial diversity and composition and the major histocompatibility complex class II exon 2 diversity in 12 moor frog (Rana arvalis) populations belonging to two geographical clusters that show signatures of past and ongoing differential selection. We found that while bacterial alpha diversity did not differ between the two clusters, MHC alleles/supertypes and genetic diversity varied considerably depending on geography and evolutionary history. Bacterial alpha diversity was positively correlated with expected MHC heterozygosity and negatively with MHC nucleotide diversity. Furthermore, bacterial community composition showed significant variation between the two geographical clusters and between specific MHC alleles/supertypes. Our findings emphasize the importance of historical demographic events on hologenomic variation and provide new insights into how immunogenetic host variability and microbial diversity may jointly influence host fitness with consequences for disease susceptibility and population persistence.

Avian malaria parasite infections do not affect personality in the chestnut thrush (Turdus rubrocanus) on the Qinghai-Tibet Plateau

Personality traits, the consistent individual behavioral differences, are currently gaining much attention in studies of natural bird populations. However, associations between personality traits and parasite infections are not often investigated. Even less attention has been given to studies of birds in the high-elevation region such as the Tibetan plateau. This research aims to examine the relationship between avian malaria parasites and two personality traits in a population of the Chestnut Thrush (Turdus rubrocanus) breed in the Tibetan plateau. Our results revealed no evidence of sex bias in malaria parasite prevalence. Furthermore, we found no effect of infection status on two personality scores: activity and boldness. Additionally, no effects on the activity level or boldness were observed for different parasite lineages of Haemoproteus, Leucocytozoon, the sex of the birds, or their interactions. Similarly, we did not find any relationship between activity level and boldness with nestling numbers, sex, or their interactions. Notably, individuals with a larger number of offspring tended to display greater boldness. Our findings indicate that blood parasite infections are common in this population but do not significantly impact the personality of the birds.

Understanding the evolution of immune genes in jawed vertebrates

Driven by co-evolution with pathogens, host immunity continuously adapts to optimize defence against pathogens within a given environment. Recent advances in genetics, genomics and transcriptomics have enabled a more detailed investigation into how immunogenetic variation shapes the diversity of immune responses seen across domestic and wild animal species. However, a deeper understanding of the diverse molecular mechanisms that shape immunity within and among species is still needed to gain insight into-and generate evolutionary hypotheses on-the ultimate drivers of immunological differences. Here, we discuss current advances in our understanding of molecular evolution underpinning jawed vertebrate immunity. First, we introduce the immunome concept, a framework for characterizing genes involved in immune defence from a comparative perspective, then we outline how immune genes of interest can be identified. Second, we focus on how different selection modes are observed acting across groups of immune genes and propose hypotheses to explain these differences. We then provide an overview of the approaches used so far to study the evolutionary heterogeneity of immune genes on macro and microevolutionary scales. Finally, we discuss some of the current evidence as to how specific pathogens affect the evolution of different groups of immune genes. This review results from the collective discussion on the current key challenges in evolutionary immunology conducted at the ESEB 2021 Online Satellite Symposium: Molecular evolution of the vertebrate immune system, from the lab to natural populations.

A Chromosome-Level Genome Assembly for the Rock Ptarmigan (Lagopus muta)

The Rock Ptarmigan (Lagopus muta) is a cold-adapted, largely sedentary, game bird with a Holarctic distribution. The species represents an important example of an organism likely to be affected by ongoing climatic shifts across a disparate range. We provide here a high-quality reference genome and mitogenome for the Rock Ptarmigan assembled from PacBio HiFi and Hi-C sequencing of a female bird from Iceland. The total size of the genome is 1.03 Gb with a scaffold N50 of 71.23 Mb and a contig N50 of 17.91 Mb. The final scaffolds represent all 40 predicted chromosomes, and the mitochondria with a BUSCO score of 98.6%. Gene annotation resulted in 16,078 protein-coding genes out of a total 19,831 predicted (81.08% excluding pseudogenes). The genome included 21.07% repeat sequences, and the average length of genes, exons, and introns were, 33605, 394, and 4265 bp respectively. The availability of a new reference-quality genome will contribute to understanding the Rock Ptarmigan's unique evolutionary history, vulnerability to climate change, and demographic trajectories around the globe while serving as a benchmark for species in the family Phasianidae (order Galliformes).

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Major Histocompatibility Complex Variation and Haplotype Associated Survival in Response to Experimental Infection of Two Bd-GPL Strains Along a Latitudinal Gradient

While both innate and adaptive immune system mechanisms have been implicated in resistance against the chytrid fungus Batrachochytrium dendrobatidis (Bd), studies on the role of specific MHC haplotypes on Bd infection are rare. Here, we studied variation in MHC Class IIB loci in the common toad Bufo bufo along a latitudinal gradient across Sweden. In general, Swedish toad populations had few MHC Class IIB haplotypes and MHC diversity declined from south (13 haplotypes) to the north (four haplotypes). The low diversity may compromise the ability of northern populations to fight emerging disease, such as Bd. In a laboratory experiment, we infected newly metamorphosed toads with two strains of the Global Pandemic Lineage of the fungus (Bd-GPL) and compared survival with sham controls. Bd-infected toads had lower survival compared to controls. Moreover, survival was dependent on the Bd-strain and northern toads had lower Bd-mediated survival than southern individuals. MHC diversity was lower in northern toads. All northern experimental animals were monomorphic for a single MHC haplotype, whereas we found seven different haplotypes in southern experimental animals. In southern toads, survival was dependent on both Bd-strain and MHC haplotype suggesting differential infection dynamics depending on both Bd-strain and host immune system characteristics.

The era of reference genomes in conservation genomics

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics.

Demographic history and divergence of sibling grouse species inferred from whole genome sequencing reveal past effects of climate change

Background

The boreal forest is one of the largest biomes on earth, supporting thousands of species. The global climate fluctuations in the Quaternary, especially the ice ages, had a significant influence on the distribution of boreal forest, as well as the divergence and evolution of species inhabiting this biome. To understand the possible effects of on-going and future climate change it would be useful to reconstruct past population size changes and relate such to climatic events in the past. We sequenced the genomes of 32 individuals from two forest inhabiting bird species, Hazel Grouse (Tetrastes bonasia) and Chinese Grouse (T. sewerzowi) and three representatives of two outgroup species from Europe and China.

Results

We estimated the divergence time of Chinese Grouse and Hazel Grouse to 1.76 (0.46–3.37) MYA. The demographic history of different populations in these two sibling species was reconstructed, and showed that peaks and bottlenecks of effective population size occurred at different times for the two species. The northern Qilian population of Chinese Grouse became separated from the rest of the species residing in the south approximately 250,000 years ago and have since then showed consistently lower effective population size than the southern population. The Chinese Hazel Grouse population had a higher effective population size at the peak of the Last Glacial Period (approx. 300,000 years ago) than the European population. Both species have decreased recently and now have low effective population sizes.

Conclusions

Combined with the uplift history and reconstructed climate change during the Quaternary, our results support that cold-adapted grouse species diverged in response to changes in the distribution of palaeo-boreal forest and the formation of the Loess Plateau. The combined effects of climate change and an increased human pressure impose major threats to the survival and conservation of both species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01921-7.

King or royal family? Testing for species boundaries in the King Cobra, Ophiophagus hannah (Cantor, 1836), using morphology and multilocus DNA analyses

In widespread species, the diverse ecological conditions in which the populations occur, and the presence of many potential geographical barriers through their range are expected to have created ample opportunities for the evolution of distinct, often cryptic lineages. In this work, we tested for species boundaries in one such widespread species, the king cobra, Ophiophagus hannah (Cantor, 1836), a largely tropical elapid snake distributed across the Oriental realm. Based on extensive geographical sampling across most of the range of the species, we initially tested for candidate species (CS) using Maximum-Likelihood analysis of mitochondrial genes. We then tested the resulting CS using both morphological data and sequences of three single-copy nuclear genes. We used snapclust to determine the optimal number of clusters in the nuclear dataset, and Bayesian Phylogenetics and Phylogeography (BPP) to test for likely species status. We used non-metric multidimensional scaling (nMDS) analysis for discerning morphological separation. We recovered four independently evolving, geographically separated lineages that we consider Confirmed Candidate Species: (1) Western Ghats lineage; (2) Indo-Chinese lineage (3) Indo-Malayan lineage; (4) Luzon Island lineage, in the Philippine Archipelago. We discuss patterns of lineage divergence, particularly in the context of low morphological divergence, and the conservation implications of recognizing several endemic king cobra lineages.

Whole genome sequencing reveals high differentiation, low levels of genetic diversity and short runs of homozygosity among Swedish wels catfish

The use of genetic markers in the context of conservation is largely being outcompeted by whole-genome data. Comparative studies between the two are sparse, and the knowledge about potential effects of this methodology shift is limited. Here, we used whole-genome sequencing data to assess the genetic status of peripheral populations of the wels catfish (Silurus glanis), and discuss the results in light of a recent microsatellite study of the same populations. The Swedish populations of the wels catfish have suffered from severe declines during the last centuries and persists in only a few isolated water systems. Fragmented populations generally are at greater risk of extinction, for example due to loss of genetic diversity, and may thus require conservation actions. We sequenced individuals from the three remaining native populations (Båven, Emån, and Möckeln) and one reintroduced population of admixed origin (Helge å), and found that genetic diversity was highest in Emån but low overall, with strong differentiation among the populations. No signature of recent inbreeding was found, but a considerable number of short runs of homozygosity were present in all populations, likely linked to historically small population sizes and bottleneck events. Genetic substructure within any of the native populations was at best weak. Individuals from the admixed population Helge å shared most genetic ancestry with the Båven population (72%). Our results are largely in agreement with the microsatellite study, and stresses the need to protect these isolated populations at the northern edge of the distribution of the species.

A novel duplex ddPCR assay for detection and differential diagnosis of Ascaridia galli and Heterakis gallinarum eggs from chickens feces

Since the EU ban on battery cages, many studies have listed Ascaridia galli and Heterakis gallinarum as the most common roundworms in the European laying hen population. A complicating factor is that the eggs of these parasites are almost identical. Thus, lack of molecular diagnostic approaches has driven epidemiological studies to take on necropsy for species discrimination, which is labor and cost intensive. Here, we describe a novel diagnostic tool based on droplet digital PCR for simultaneous identification and absolute quantification of the eggs of both of these ascarids in chickens' droppings using two different genus-specific primer-probe sets targeting the second internal transcribed spacer region (ITS-2) in the nuclear ribosomal (rRNA) gene array. No cross-reaction was observed when different combinations of DNA and species-specific primers and probes were tested. The lowest obtained frequency threshold for the detection of H. gallinarum in the presence of a constant A. galli DNA concentration was determined to be 0.8 %. After validation, we used the assay to analyze field samples collected from several Swedish laying hen farms. Out of 134 samples, 86 (64 %) were positive for A. galli while 11 (8.3 %) samples were positive for H. gallinarum. These samples were initially analyzed with flotation technique for detection of ascarid eggs. The results of the Cohen's kappa indicated substantial agreement (85.8 %) between the two tests. In conclusion, we have validated a novel molecular-based diagnostic tool for quantification and differentiation between intestinal parasites of major importance in chickens with high precision. Although this study focuses on identification of parasites of laying hens, the findings may well have a bearing on all types of chicken production systems. The present study lays the groundwork for future research into epidemiology of these two important chicken parasite species.

The effects of drift and selection on latitudinal genetic variation in Scandinavian common toads (Bufo bufo) following postglacial recolonisation

Clinal variation is paramount for understanding the factors shaping genetic diversity in space and time. During the last glacial maximum, northern Europe was covered by glacial ice that rendered the region uninhabitable for most taxa. Different evolutionary processes during and after the recolonisation of this area from different glacial refugia have affected the genetic landscape of the present day European flora and fauna. In this study, we focus on the common toad (Bufo bufo) in Sweden and present evidence suggesting that these processes have resulted in two separate lineages of common toad, which colonised Sweden from two directions. Using ddRAD sequencing data for demographic modelling, structure analyses, and analysis of molecular variance (AMOVA), we provide evidence of a contact zone located between Uppland and Västerbotten in central Sweden. Genetic diversity was significantly higher in southern Sweden compared to the north, in accordance with a pattern of decreased genetic diversity with increasing distance from glacial refugia. Candidate genes under putative selection are identified through outlier detection and gene-environment association methods. We provide evidence of divergent selection related to stress response and developmental processes in these candidate genes. The colonisation of Sweden by two separate lineages may have implications for how future conservation efforts should be directed by identifying management units and putative local adaptations.

Effect of a targeted treatment strategy against Ascaridia galli on egg production, egg quality and bird health in a laying hen farm

Worm control is an important aspect of the successful management of the egg production industry. Of particular concern is Ascaridia galli, which at high parasite loads affect health and production in layers. Application of a targeted treatment strategy (TT) to control A. galli has shown promise. The purpose of this study was to explore the effect of such a strategy on welfare indicators and production performance of layers. Six flocks (F1-6) on a commercial farm were allocated to three treatment groups. Flocks F1 and F4 were treated (TT) with fenbendazole at 22, 27 and 36 weeks post-placement (WPP). Flocks F2 and F5 were treated at 27 WPP (conventional treatment, CT) and hens in flocks F3 and F6 served as untreated (UT) control groups. At 19, 35 and 45 WPP twenty-five hens plus thirty eggs per flock were randomly selected. Hens were weighed and their plumage conditions (PC) were assessed. The eggs were subjected to various external and interior quality analyses. Production data such as number of eggs/hen/week, egg mass and feed conversion ratio (FCR) were calculated from raw data obtained from all flocks on a weekly basis. The number of eggs/hen/week, egg mass and FCR were higher (P < 0.05) in the TT flocks and hens had better PC both at 35 and 45 WPP compared with other flocks. No differences in body weight and physical egg quality were observed between groups except for egg shell strength which was higher (P < 0.05) in the CT flocks. These data suggest that better production performance and plumage, which suggests improved health, can be achieved through the application of a TT strategy. The insights gained from this research should help to justify the extra cost and labor associated with the TT strategy.

Occurrence of Batrachochytrium dendrobatidis in Sweden: higher infection prevalence in southern species

The chytrid fungus Batrachochytrium dendrobatidis (Bd) has caused worldwide declines in amphibian populations. While Bd is widespread in southern and central Europe, its occurrence and distribution in northernmost Europe is mostly unknown. We surveyed for Bd in breeding anurans in Sweden by sampling 1917 amphibians from 101 localities and 3 regions in Sweden (southern, northern and central). We found that Bd was widespread in southern and central Sweden, occurring in all 9 investigated species and in 45.5% of the 101 localities with an overall prevalence of 13.8%. No infected individuals were found in the 4 northern sites sampled. The records from central Sweden represent the northernmost records of Bd in Europe. While the proportion of sites positive for Bd was similar between the southern and central regions, prevalence was much higher in the southern region. This was because southern species with a distribution mainly restricted to southernmost Sweden had a higher prevalence than widespread generalist species. The nationally red-listed green toad Bufotes variabilis and the fire-bellied toad Bombina bombina had the highest prevalence (61.4 and 48.9%, respectively). Across species, Bd prevalence was strongly positively, correlated with water temperature at the start of egg laying. However, no individuals showing visual signs of chytridiomycosis were found in the field. These results indicate that Bd is widespread and common in southern and central Sweden with southern species, breeding in higher temperatures and with longer breeding periods, having higher prevalence. However, the impact of Bd on amphibian populations in northernmost Europe remains unknown.

Genomic analysis of demographic history and ecological niche modeling in the endangered Chinese Grouse Tetrastes sewerzowi

Background

The Quaternary had worldwide consequences in forming the contemporary diversity of many populations, species and communities, which is characterized by marked climatic oscillations between glacial and interglacial periods. The origin and evolution of biodiversity in mountainous areas are highly dependent on historical orogenesis and associated climatic changes. The Chinese grouse Tetrastes sewerzowi is a forest-dwelling species endemic to the mountains to the east of the Qinghai–Tibet Plateau, which has been listed as Near Threatened with a decreasing trend by the IUCN because of ongoing deforestation and fragmentation of coniferous forests. It is important to place current population status into a broader ecological and evolutionary context to understand their demographic history.

Results

Analyses of the Chinese Grouse genome revealed fluctuations throughout the Pleistocene in effective population size. Populations decreased during early to middle Pleistocene but showed an expansion during late Pleistocene which was then followed by a sharp decline during the last glacial maximum (LGM). Ecological niche modeling indicated that a suitable habitat shift between high altitude regions to low altitude regions was due to a changing climate. This result parallels patterns of population size change in Chinese Grouse estimated from PSMC modelling, which suggested an expansion in population size from the last interglacial period (LIG) and then a peak and a bottleneck occurring at the last glacial maximum (LGM). Furthermore, the present-day distribution of Chinese Grouse is greatly reduced and fragmented. It will likely become even more fragmented in the future since coniferous forest cover is threatened in the region of their distribution and the availability of such habitat restricts their ecological niche.

Conclusions

The Chinese Grouse have experienced substantial population size changes from the beginning to the LIG and reached a peak before the LGM. A sharp decrease and bottleneck occurred during the LGM, when the coniferous forests were subjected to extensive loss. The results inferred from the whole genome sequencing and species distribution models both support historical population fluctuations. The distribution of the Chinese Grouse is strongly dependent on the coniferous forest cover. To protect the fragmented coniferous forests is an essential action to protect the Chinese Grouse.

Initial assessment of the economic burden of major parasitic helminth infections to the ruminant livestock industry in Europe

We report a European wide assessment of the economic burden of gastrointestinal nematodes, Fasciola hepatica (common liver fluke) and Dictyocaulus viviparus (bovine lungworm) infections to the ruminant livestock industry. The economic impact of these parasitic helminth infections was estimated by a deterministic spreadsheet model as a function of the proportion of the ruminant population exposed to grazing, the infection frequency and intensity, the effect of the infection on animal productivity and mortality and anthelmintic treatment costs. In addition, we estimated the costs of anthelmintic resistant nematode infections and collected information on public research budgets addressing helminth infections in ruminant livestock. The epidemiologic and economic input data were collected from international databases and via expert opinion of the Working Group members of the European Co-operation in Science and Technology (COST) action COMbatting Anthelmintic Resistance in ruminants (COMBAR). In order to reflect the effects of uncertainty in the input data, low and high cost estimates were obtained by varying uncertain input data arbitrarily in both directions by 20 %. The combined annual cost [low estimate-high estimate] of the three helminth infections in 18 participating countries was estimated at € 1.8 billion [€ 1.0-2.7 billion]. Eighty-one percent of this cost was due to lost production and 19 % was attributed to treatment costs. The cost of gastrointestinal nematode infections with resistance against macrocyclic lactones was estimated to be € 38 million [€ 11-87 million] annually. The annual estimated costs of helminth infections per sector were € 941 million [€ 488 - 1442 million] in dairy cattle, € 423 million [€ 205-663 million] in beef cattle, € 151million [€ 90-213 million] in dairy sheep, € 206 million [€ 132-248 million] in meat sheep and € 86 million [€ 67-107 million] in dairy goats. Important data gaps were present in all phases of the calculations which lead to large uncertainties around the estimates. Accessibility of more granular animal population datasets at EU level, deeper knowledge of the effects of infection on production, levels of infection and livestock grazing exposure across Europe would make the largest contribution to improved burden assessments. The known current public investment in research on helminth control was 0.15 % of the estimated annual costs for the considered parasitic diseases. Our data suggest that the costs of enzootic helminth infections which usually occur at high prevalence annually in ruminants, are similar or higher than reported costs of epizootic diseases. Our data can support decision making in research and policy to mitigate the negative impacts of helminth infections and anthelmintic resistance in Europe, and provide a baseline against which to measure future changes.

High variation in last male sperm precedence and genital morphology in the emerald damselfly, Lestes sponsa

In organisms in which individuals mate multiply, knowledge of the proportion of offspring sired by the last male to mate (P2) under field conditions is important for a thorough understanding of how sexual selection works in nature. In many insect groups, pronounced intraspecific variation in P2 is commonplace. Interestingly, however, in stark contrast to these observations, compilation of P2 data in dragonflies and damselflies (Odonata) indicates that a high P2, seldom below 0.95, is a feature of this taxon. Here we used double digest restriction-site associated DNA sequencing to generate a panel of single nucleotide polymorphisms (SNPs) with which we could determine paternity and estimate values of P2 in the offspring of 19 field-collected pairs of the emerald damselfly Lestes sponsa. We also estimated the relationship between P2 and male genital shape of 16 males using geometric morphometric analysis. P2 was variable (range = 0.0–1.0; mean = 0.5), and there was a marginally non-significant (P = 0.069) relationship between genital shape and P2, suggesting that males with a high P2 had an aedeagus with a broader tip. We suggest that the high P2-values reported in past studies in Odonata are partly due to the methods used to infer paternity. Use of SNPs to determine patterns of paternity and P2 in odonates is needed for a better appraisal of fitness in odonates, and would open many future avenues for use of odonates as models of sexual selection.

Genetic differentiation in Sichuan jay (Perisoreus internigrans) and its sibling species Siberian jay (P. infaustus)

The Pleistocene ice age and recent forest fragmentation have both played a significant role in shaping the population genetic variation of boreal coniferous species in the Qinghai-Tibet Plateau (QTP) and Eurasian coniferous forests. The Sichuan jay is one of the least know endemic bird species in QTP coniferous forests in western China while its sibling species, the Siberian Jay, is widespread within the coniferous forests in northern Eurasia. Here we used 11 microsatellite markers to assay genetic diversity across 58 Sichuan jay samples from China and 205 Siberian jay samples from Sweden and Russia. Results showed three distinct genetic clusters from the Sichuan jay sampling. Furthermore, the pair-wise FST values indicated high genetic differentiation not only among the two species but also between Swedish and Russian Siberian jay populations. What is more, a pattern of isolation by distance was found among the analyzed populations. Our study suggests that targeted habitat restoration in fragmented forests and more genetic work is urgently needed for conservation of the Sichuan jay.

Antimicrobial peptide and sequence variation along a latitudinal gradient in two anurans

BACKGROUND

While there is evidence of both purifying and balancing selection in immune defense genes, large-scale genetic diversity in antimicrobial peptides (AMPs), an important part of the innate immune system released from dermal glands in the skin, has remained uninvestigated. Here we describe genetic diversity at three AMP loci (Temporin, Brevinin and Palustrin) in two ranid frogs (Rana arvalis and R. temporaria) along a 2000 km latitudinal gradient. We amplified and sequenced part of the Acidic Propiece domain and the hypervariable Mature Peptide domain (~ 150-200 bp) in the three genes using Illumina Miseq and expected to find decreased AMP genetic variation towards the northern distribution limit of the species similarly to studies on MHC genetic patterns.

RESULTS

We found multiple loci for each AMP and relatively high gene diversity, but no clear pattern of geographic genetic structure along the latitudinal gradient. We found evidence of trans-specific polymorphism in the two species, indicating a common evolutionary origin of the alleles. Temporin and Brevinin did not form monophyletic clades suggesting that they belong to the same gene family. By implementing codon evolution models we found evidence of strong positive selection acting on the Mature Peptide. We also found evidence of diversifying selection as indicated by divergent allele frequencies among populations and high Theta k values.

CONCLUSION

Our results suggest that AMPs are an important source of adaptive diversity, minimizing the chance of microorganisms developing resistance to individual peptides.

Genome Assembly of the Common Pheasant Phasianus colchicus: A Model for Speciation and Ecological Genomics

The common pheasant (Phasianus colchicus) in the order Galliformes and the family Phasianidae, has 30 subspecies distributed across its native range in the Palearctic realm and has been introduced to Europe, North America, and Australia. It is an important game bird often subjected to wildlife management as well as a model species to study speciation, biogeography, and local adaptation. However, the genomic resources for the common pheasant are generally lacking. We sequenced a male individual of the subspecies torquatus of the common pheasant with the Illumina HiSeq platform. We obtained 94.88 Gb of usable sequences by filtering out low-quality reads of the raw data generated. This resulted in a 1.02 Gb final assembly, which equals the estimated genome size. BUSCO analysis using chicken as a model showed that 93.3% of genes were complete. The contig N50 and scaffold N50 sizes were 178 kb and 10.2 Mb, respectively. All these indicate that we obtained a high-quality genome assembly. We annotated 16,485 protein-coding genes and 123.3 Mb (12.05% of the genome) of repetitive sequences by ab initio and homology-based prediction. Furthermore, we applied a RAD-sequencing approach for another 45 individuals of seven representative subspecies in China and identified 4,376,351 novel single nucleotide polymorphism (SNPs) markers. Using this unprecedented data set, we uncovered the geographic population structure and genetic introgression among common pheasants in China. Our results provide the first high-quality reference genome for the common pheasant and a valuable genome-wide SNP database for studying population genomics and demographic history.

Endless forms of sexual selection

In recent years, the field of sexual selection has exploded, with advances in theoretical and empirical research complementing each other in exciting ways. This perspective piece is the product of a “stock-taking” workshop on sexual selection and sexual conflict. Our aim is to identify and deliberate on outstanding questions and to stimulate discussion rather than provide a comprehensive overview of the entire field. These questions are organized into four thematic sections we deem essential to the field. First we focus on the evolution of mate choice and mating systems. Variation in mate quality can generate both competition and choice in the opposite sex, with implications for the evolution of mating systems. Limitations on mate choice may dictate the importance of direct vs. indirect benefits in mating decisions and consequently, mating systems, especially with regard to polyandry. Second, we focus on how sender and receiver mechanisms shape signal design. Mediation of honest signal content likely depends on integration of temporally variable social and physiological costs that are challenging to measure. We view the neuroethology of sensory and cognitive receiver biases as the main key to signal form and the ‘aesthetic sense’ proposed by Darwin. Since a receiver bias is sufficient to both initiate and drive ornament or armament exaggeration, without a genetically correlated or even coevolving receiver, this may be the appropriate ‘null model’ of sexual selection. Thirdly, we focus on the genetic architecture of sexually selected traits. Despite advances in modern molecular techniques, the number and identity of genes underlying performance, display and secondary sexual traits remains largely unknown. In-depth investigations into the genetic basis of sexual dimorphism in the context of long-term field studies will reveal constraints and trajectories of sexually selected trait evolution. Finally, we focus on sexual selection and conflict as drivers of speciation. Population divergence and speciation are often influenced by an interplay between sexual and natural selection. The extent to which sexual selection promotes or counteracts population divergence may vary depending on the genetic architecture of traits as well as the covariance between mating competition and local adaptation. Additionally, post-copulatory processes, such as selection against heterospecific sperm, may influence the importance of sexual selection in speciation. We propose that efforts to resolve these four themes can catalyze conceptual progress in the field of sexual selection, and we offer potential avenues of research to advance this progress.

Genotype by environment interaction for female fertility traits under conventional and organic production systems in Danish Holsteins

Conventional and organic production systems mainly differ in feeding strategies, outdoor and pasture access, and the use of antibiotic treatments. These environmental differences could lead to a genotype by environment interaction (G × E) and a requirement for including G × E in breeding decisions. The objectives of this study were to estimate variance components and heritabilities for conventional and organic production systems and investigate G × E under these 2 production systems for female fertility traits in Danish Holsteins. The analyzed traits included the interval from calving to first insemination (ICF), the interval from first to last insemination, number of inseminations per conception (NINS), and non-return rate within 56 d after the first insemination. Records of female fertility in heifers and the first 3 lactations in cows as well as grass ratio of feed at herd level were collected during the period from 2011 to 2016. The performances of a trait in heifers and cows (lactation 1 to 3) were considered as different traits. The (co)variance components and the resulting heritabilities and genetic correlations were estimated using 2 models. One was a bivariate model treating performances of a trait under organic and conventional production systems as 2 different traits using a reduced data set, and the other was a reaction norm model with random regression on the production system and the grass ratio of feed using a full data set. The full data set comprised records of 37,836 females from 112 organic herds and 513,599 females from 1,224 conventional herds, whereas the reduced data set comprised records from all these 112 organic herds and 92,696 females from 185 convention herds extracted from the full data set with grass ratio of feed lower than 0.20. All female fertility performances of the organic production system were superior to those of the conventional production system. Besides, heterogeneities in additive genetic variances and heritabilities were observed between conventional and organic production systems for all traits. Furthermore, genetic correlations between these 2 production systems ranged from 0.607 to 1.000 estimated from bivariate models and from 0.848 to 0.999 estimated from reaction norm models. Statistically significant G × E were observed for NINS in heifers, non-return rate within 56 d after the first insemination in heifers, and ICF from the bivariate model, and for ICF and NINS in cows from the reaction norm model.

Urban environment and reservoir host species are associated with Batrachochytrium dendrobatidis infection prevalence in the common toad

Human-induced changes of the environment, including landscape alteration and habitat loss, may affect wildlife disease dynamics and have important ramifications for wildlife conservation. Amphibians are among the vertebrate taxa most threatened by anthropogenic habitat change. The emerging fungal pathogen Batrachochytrium dendrobatidis (Bd) has caused extinctions and population declines in hundreds of anuran species globally. We studied how the urban landscape is associated with the prevalence of Bd infections by sampling 655 anurans of 3 species (mainly the common toad Bufo bufo) in 42 ponds surrounded by different amounts of urban habitat (defined as towns, cities or villages). We also examined the association between Bd infections and a potential reservoir host species (the moor frog Rana arvalis). We found that 38% of the sites were positive for Bd with an infection prevalence of 4.4%. The extent of urban landscape was negatively correlated with Bd infection prevalence. However, the positive association of Bd with the presence of the possible reservoir species was substantially stronger than the urban effects. The body condition index of B. bufo was negatively associated with Bd infection. This Bd effect was stronger than the negative effect of urban landscape on body condition. Our results suggest that urban environments in Sweden have a negative impact on Bd infections, while the presence of the reservoir species has a positive impact on Bd prevalence. Our study also highlights the potential importance of Bd infection on host fitness, especially in rural landscapes.

Genomic regions of speciation and adaptation among three species of grouse

Understanding the molecular basis of adaption is one of the central goals in evolutionary biology and when investigated across sister species it can provide detailed insight into the mechanisms of speciation. Here, we sequence the genomes of 34 individuals from three closely related grouse species in order to uncover the genomic architecture of speciation and the genes involved in adaptation. We identify 6 regions, containing 7 genes that show lineage specific signs of differential selection across the species. These genes are involved in a variety of cell processes ranging from stress response to neural, gut, olfactory and limb development. Genome wide neutrality test statistics reveal a strong signal of population expansion acting across the genomes. Additionally, we uncover a 3.5 Mb region on chromosome 20 that shows considerably lower levels of differentiation across the three grouse lineages, indicating possible action of uniform selection in this region.

Effects of host species and environmental factors on the prevalence of Batrachochytrium dendrobatidis in northern Europe

The fungal pathogen Batrachochytrium dendrobatidis (Bd) poses a major threat to amphibian populations. To assist efforts to address such threats, we examined differences in Bd host infection prevalence among amphibian species and its relations to both local environmental factors in breeding habitats and landscape variables measured at three scales (500, 2000 and 5000 m radii) around breeding sites in southernmost Sweden. We sampled 947 anurans of six species in 31 ponds and assessed their infection status. We then examined correlations of infection prevalence with canopy cover, pond perimeter and pH (treated as local-scale pond characteristics), and the number of ponds, area of arable land, area of mature forest, number of resident people and presence of sea within the three radii (treated as landscape variables). The Bd infection prevalence was very low, 0.5–1.0%, in two of the six anuran species (Bufo bufo and Rana temporaria), and substantially higher (13–64%) in the other four (Bombina bombina, Bufotes variabilis, Epidalea calamita, Rana arvalis). In the latter four species Bd infection prevalence was positively associated with ponds’ pH (site range: 5.3–8.1), and negatively associated with areas of mature forest and/or wetlands in the surroundings. Our results show that the infection dynamics of Bd are complex and associated with host species, local pond characteristics and several landscape variables at larger spatial scales. Knowledge of environmental factors associated with Bd infections and differences in species’ susceptibility may help to counter further spread of the disease and guide conservation action plans, especially for the most threatened species.

Past and potential future population dynamics of three grouse species using ecological and whole genome coalescent modeling

Studying demographic history of species provides insight into how the past has shaped the current levels of overall biodiversity and genetic composition of species, but also how these species may react to future perturbations. Here we investigated the demographic history of the willow grouse (Lagopus lagopus), rock ptarmigan (Lagopus muta), and black grouse (Tetrao tetrix) through the Late Pleistocene using two complementary methods and whole genome data. Species distribution modeling (SDM) allowed us to estimate the total range size during the Last Interglacial (LIG) and Last Glacial Maximum (LGM) as well as to indicate potential population subdivisions. Pairwise Sequentially Markovian Coalescent (PSMC) allowed us to assess fluctuations in effective population size across the same period. Additionally, we used SDM to forecast the effect of future climate change on the three species over the next 50 years. We found that SDM predicts the largest range size for the cold-adapted willow grouse and rock ptarmigan during the LGM. PSMC captured intraspecific population dynamics within the last glacial period, such that the willow grouse and rock ptarmigan showed multiple bottlenecks signifying recolonization events following the termination of the LGM. We also see signals of population subdivision during the last glacial period in the black grouse, but more data are needed to strengthen this hypothesis. All three species are likely to experience range contractions under future warming, with the strongest effect on willow grouse and rock ptarmigan due to their limited potential for northward expansion. Overall, by combining these two modeling approaches, we have provided a multifaceted examination of the biogeography of these species and how they have responded to climate change in the past. These results help us understand how cold-adapted species may respond to future climate changes.

Development and worldwide use of non-lethal, and minimal population-level impact, protocols for the isolation of amphibian chytrid fungi

Parasitic chytrid fungi have emerged as a significant threat to amphibian species worldwide, necessitating the development of techniques to isolate these pathogens into culture for research purposes. However, early methods of isolating chytrids from their hosts relied on killing amphibians. We modified a pre-existing protocol for isolating chytrids from infected animals to use toe clips and biopsies from toe webbing rather than euthanizing hosts, and distributed the protocol to researchers as part of the BiodivERsA project RACE; here called the RML protocol. In tandem, we developed a lethal procedure for isolating chytrids from tadpole mouthparts. Reviewing a database of use a decade after their inception, we find that these methods have been applied across 5 continents, 23 countries and in 62 amphibian species. Isolation of chytrids by the non-lethal RML protocol occured in 18% of attempts with 207 fungal isolates and three species of chytrid being recovered. Isolation of chytrids from tadpoles occured in 43% of attempts with 334 fungal isolates of one species (Batrachochytrium dendrobatidis) being recovered. Together, these methods have resulted in a significant reduction and refinement of our use of threatened amphibian species and have improved our ability to work with this group of emerging pathogens.

Recent Asian origin of chytrid fungi causing global amphibian declines

Globalized infectious diseases are causing species declines worldwide, but their source often remains elusive. We used whole-genome sequencing to solve the spatiotemporal origins of the most devastating panzootic to date, caused by the fungus Batrachochytrium dendrobatidis, a proximate driver of global amphibian declines. We traced the source of B. dendrobatidis to the Korean peninsula, where one lineage, BdASIA-1, exhibits the genetic hallmarks of an ancestral population that seeded the panzootic. We date the emergence of this pathogen to the early 20th century, coinciding with the global expansion of commercial trade in amphibians, and we show that intercontinental transmission is ongoing. Our findings point to East Asia as a geographic hotspot for B. dendrobatidis biodiversity and the original source of these lineages that now parasitize amphibians worldwide.

Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis)

MHC genes are key components in disease resistance and an excellent system for studying selection acting on genetic variation in natural populations. Current patterns of variation in MHC genes are likely to be influenced by past and ongoing selection as well as demographic fluctuations in population size such as those imposed by post-glacial recolonization processes. Here, we investigated signatures of historical selection and demography on an MHC class II gene in 12 moor frog populations along a 1700-km latitudinal gradient. Sequences were obtained from 207 individuals and consecutively assigned into two different clusters (northern and southern clusters, respectively) in concordance with a previously described dual post-glacial colonization route. Selection analyses comparing the relative rates of non-synonymous to synonymous substitutions (dN/dS) suggested evidence of different selection patterns in the northern and the southern clusters, with divergent selection prevailing in the south but uniform positive selection predominating in the north. Also, models of codon evolution revealed considerable differences in the strength of selection: The southern cluster appeared to be under strong selection while the northern cluster showed moderate signs of selection. Our results indicate that the MHC alleles in the north diverged from southern MHC alleles as a result of differential selection patterns.

Mind the gaps in research on the control of gastrointestinal nematodes of farmed ruminants and pigs

Gastrointestinal (GI) nematode control has an important role to play in increasing livestock production from a limited natural resource base and to improve animal health and welfare. In this synthetic review, we identify key research priorities for GI nematode control in farmed ruminants and pigs, to support the development of roadmaps and strategic research agendas by governments, industry and policymakers. These priorities were derived from the DISCONTOOLS gap analysis for nematodes and follow-up discussions within the recently formed Livestock Helminth Research Alliance (LiHRA). In the face of ongoing spread of anthelmintic resistance (AR), we are increasingly faced with a failure of existing control methods against GI nematodes. Effective vaccines against GI nematodes are generally not available, and anthelmintic treatment will therefore remain a cornerstone for their effective control. At the same time, consumers and producers are increasingly concerned with environmental issues associated with chemical parasite control. To address current challenges in GI nematode control, it is crucial to deepen our insights into diverse aspects of epidemiology, AR, host immune mechanisms and the socio-psychological aspects of nematode control. This will enhance the development, and subsequent uptake, of the new diagnostics, vaccines, pharma-/nutraceuticals, control methods and decision support tools required to respond to the spread of AR and the shifting epidemiology of GI nematodes in response to climatic, land-use and farm husbandry changes. More emphasis needs to be placed on the upfront evaluation of the economic value of these innovations as well as the socio-psychological aspects to prioritize research and facilitate uptake of innovations in practice. Finally, targeted regulatory guidance is needed to create an innovation-supportive environment for industries and to accelerate the access to market of new control tools.

Phylogeography and larval spine length of the dragonfly Leucorhinia dubia in Europe

Presence or absence of predators selects for different kind of morphologies. Hence, we expect variation in traits that protect against predators to vary over geographical areas where predators vary in past and present abundance. Abdominal larval spines in dragonfly larvae provide protection against fish predators. We studied geographical variation in larval spine length of the dragonfly Leucorrhinia dubia across Western Europe using a phylogenetic approach. Larvae were raised in a common garden laboratory experiment in the absence of fish predators. Results show that larvae from northern Europe (Sweden and Finland) had significantly longer larval spines compared to larvae from western and central Europe. A phylogeny based on SNP data suggests that short larval spines is the ancestral stage in the localities sampled in this study, and that long spines have evolved in the Fenno-Scandian clade. The role of predators in shaping the morphological differences among the sampled localities is discussed.

Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian

BACKGROUND

Past events like fluctuations in population size and post-glacial colonization processes may influence the relative importance of genetic drift, migration and selection when determining the present day patterns of genetic variation. We disentangle how drift, selection and migration shape neutral and adaptive genetic variation in 12 moor frog populations along a 1700 km latitudinal gradient. We studied genetic differentiation and variation at a MHC exon II locus and a set of 18 microsatellites.

RESULTS

Using outlier analyses, we identified the MHC II exon 2 (corresponding to the β-2 domain) locus and one microsatellite locus (RCO8640) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. STRUCTURE and DAPC analyses on the neutral microsatellites assigned populations to a northern and a southern cluster, reflecting two different post-glacial colonization routes found in previous studies. Genetic variation overall was lower in the northern cluster. The signature of selection on MHC exon II was weaker in the northern cluster, possibly as a consequence of smaller and more fragmented populations.

CONCLUSION

Our results show that historical demographic processes combined with selection and drift have led to a complex pattern of differentiation along the gradient where some loci are more divergent among populations than predicted from drift expectations due to diversifying selection, while other loci are more uniform among populations due to stabilizing selection. Importantly, both overall and MHC genetic variation are lower at northern latitudes. Due to lower evolutionary potential, the low genetic variation in northern populations may increase the risk of extinction when confronted with emerging pathogens and climate change.

Blood transcriptomes and de novo identification of candidate loci for mating success in lekking great snipe (Gallinago media)

We assembled the great snipe blood transcriptome using data from fourteen lekking males, in order to de novo identify candidate genes related to sexual selection, and determined the expression profiles in relation to mating success. The three most highly transcribed genes were encoding different haemoglobin subunits. All tended to be overexpressed in males with high mating success. We also called single nucleotide polymorphisms (SNPs) from the transcriptome data and found considerable genetic variation for many genes expressed during lekking. Among these, we identified 14 polymorphic candidate SNPs that had a significant genotypic association with mating success (number of females mated with) and/or mating status (mated or not). Four of the candidate SNPs were found in HBAA (encoding the haemoglobin α-chain). Heterozygotes for one of these and one SNP in the gene PABPC1 appeared to enjoy higher mating success compared to males homozygous for either of the alleles. In a larger data set of individuals, we genotyped 38 of the identified SNPs but found low support for consistent selection as only one of the zygosities of previously identified candidate SNPs and none of their genotypes were associated with mating status. However, candidate SNPs generally showed lower levels of spatial genetic structure compared to noncandidate markers. We also scored the prevalence of avian malaria in a subsample of birds. Males infected with avian malaria parasites had lower mating success in the year of sampling than noninfected males. Parasite infection and its interaction with specific genes may thus affect performance on the lek.

Range shifts or extinction? Ancient DNA and distribution modelling reveal past and future responses to climate warming in cold-adapted birds

Global warming is predicted to cause substantial habitat rearrangements, with the most severe effects expected to occur in high-latitude biomes. However, one major uncertainty is whether species will be able to shift their ranges to keep pace with climate-driven environmental changes. Many recent studies on mammals have shown that past range contractions have been associated with local extinctions rather than survival by habitat tracking. Here, we have used an interdisciplinary approach that combines ancient DNA techniques, coalescent simulations and species distribution modelling, to investigate how two common cold-adapted bird species, willow and rock ptarmigan (Lagopus lagopus and Lagopus muta), respond to long-term climate warming. Contrary to previous findings in mammals, we demonstrate a genetic continuity in Europe over the last 20 millennia. Results from back-casted species distribution models suggest that this continuity may have been facilitated by uninterrupted habitat availability and potentially also the greater dispersal ability of birds. However, our predictions show that in the near future, some isolated regions will have little suitable habitat left, implying a future decrease in local populations at a scale unprecedented since the last glacial maximum.

Looking into the past - the reaction of three grouse species to climate change over the last million years using whole genome sequences

Tracking past population fluctuations can give insight into current levels of genetic variation present within species. Analysing population dynamics over larger timescales can be aligned to known climatic changes to determine the response of species to varying environments. Here, we applied the Pairwise Sequentially Markovian Coalescent (psmc) model to infer past population dynamics of three widespread grouse species; black grouse, willow grouse and rock ptarmigan. This allowed the tracking of the effective population size (Ne ) of all three species beyond 1 Mya, revealing that (i) early Pleistocene cooling (~2.5 Mya) caused an increase in the willow grouse and rock ptarmigan populations, (ii) the mid-Brunhes event (~430 kya) and following climatic oscillations decreased the Ne of willow grouse and rock ptarmigan, but increased the Ne of black grouse and (iii) all three species reacted differently to the last glacial maximum (LGM) - black grouse increased prior to it, rock ptarmigan experienced a severe bottleneck and willow grouse was maintained at large population size. We postulate that the varying psmc signal throughout the LGM depicts only the local history of the species. Nevertheless, the large population fluctuations in willow grouse and rock ptarmigan indicate that both species are opportunistic breeders while black grouse tracks the climatic changes more slowly and is maintained at lower Ne . Our results highlight the usefulness of the psmc approach in investigating species' reaction to climate change in the deep past, but also that caution should be taken in drawing general conclusions about the recent past.