Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos.

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Evidence suggesting that reindeer mothers allonurse according to the direct reciprocity and generalized reciprocity decision rules.

Allonursing is the nursing of the offspring of other mothers. Cooperation is an emergent property of evolved decision rules. Cooperation can be explained by at least three evolved decision rules: 1) direct reciprocity, i.e. help someone who previously helped you, 2) kin discrimination, i.e. preferentially direct help to kin than to non-kin, and 3) generalized reciprocity, i.e. help anyone if helped by someone. We assessed if semi-domesticated reindeer, Rangifer tarandus, mothers allonursed according to the decision rules of direct reciprocity, generalized reciprocity and kin discrimination over 2 years. To assess if reindeer mothers allonursed according to the direct reciprocity decision rule, we predicted that mothers should give more help to those who previously helped them more often. To assess if reindeer mothers allonursed according to the kin discrimination decision rule, we predicted that help given should increase as pairwise genetic relatedness increased. To assess if reindeer mothers allonursed according to the generalized reciprocity decision rule, we predicted that the overall number of help given by reindeer mothers should increase as the overall number of help received by reindeer mothers increased. The number of help given i) increased as the number of help received from the same partner increased in the 2012 group but not in both 2013 groups, ii) was not influenced by relatedness, and iii) was not influenced by an interaction between the number of help received from the same partner and relatedness. iv) The overall number of help given increased as the overall number of help received increased. The results did not support the prediction that reindeer mothers allonursed according to the kin discrimination decision rule. The results suggest that reindeer mothers may allonurse according to the direct reciprocity and generalized reciprocity decision rules.

Adaptation to the High-Arctic island environment despite long-term reduced genetic variation in Svalbard reindeer.

Typically much smaller in number than their mainland counterparts, island populations are ideal systems to investigate genetic threats to small populations. The Svalbard reindeer (Rangifer tarandus platyrhynchus) is an endemic subspecies that colonized the Svalbard archipelago ca. 6,000-8,000 years ago and now shows numerous physiological and morphological adaptations to its arctic habitat. Here, we report a de-novo chromosome-level assembly for Svalbard reindeer and analyze 133 reindeer genomes spanning Svalbard and most of the species' Holarctic range, to examine the genomic consequences of long-term isolation and small population size in this insular subspecies. Empirical data, demographic reconstructions, and forward simulations show that long-term isolation and high inbreeding levels may have facilitated the reduction of highly deleterious-and to a lesser extent, moderately deleterious-variation. Our study indicates that long-term reduced genetic diversity did not preclude local adaptation to the High Arctic, suggesting that even severely bottlenecked populations can retain evolutionary potential.

Mismatch between calf paternity and observed copulations between male and female reindeer: Multiple mating in a polygynous ungulate?

In polygynous systems, such as that exhibited by reindeer Rangifer tarandus, mate choice can be difficult to disentangle from male intrasexual competition because male behavior may constrain female choice. Multiple mating may provide an avenue for female mate choice, though it is difficult to identify using behavioral estimators alone. Molecular techniques address this issue by affording ecologists an opportunity to reassess mating systems from a genetic perspective. We assessed the frequency and possible explanations for multiple mating in reindeer using a genetic approach to determine the success of observed copulations in a semi-domesticated herd in Kaamanen, Finland. Behavioral and genetic data were synthesized with population characteristics over a 7-year period to test the hypothesis that, if present, polyandry in reindeer is driven by sexual harassment from sub-dominant males. We observed multiple mating in 42% of females, with as many as 60% exhibiting multiple mating in certain years. We found no evidence that multiple mating resulted from sexual harassment by sub-dominant males, suggesting that it is likely a deliberate strategy among females. Conversion rate of copulations into paternities varied with male size, with smaller males more likely to experience mismatch than larger males. Female preference for larger males persisted despite the occurrence of multiple mating, possibly suggesting a mechanism for cryptic post-copulatory selection. We suggest further research to delineate the possible influence of cryptic post-copulatory selection and multiple mating to defend against infertility in exhausted males.

Strong genetic structure among populations of the tick Ixodes ricinus across its range.

Ixodes ricinus is the most common and widely distributed tick species in Europe, responsible for several zoonotic diseases, including Lyme borreliosis. Population genetics of disease vectors is a useful tool for understanding the spread of pathogens and infection risks. Despite the threat to the public health due to the climate-driven distribution changes of I. ricinus, the genetic structure of tick populations, though essential for understanding epidemiology, remains unclear. Previous studies have demonstrated weak to no apparent spatial pattern of genetic differentiation between European populations. Here, we analysed the population genetic structure of 497 individuals from 28 tick populations sampled from 20 countries across Europe, the Middle-East, and northern Africa. We analysed 125 SNPs loci after quality control. We ran Bayesian and multivariate hierarchical clustering analyses to identify and describe clusters of genetically related individuals. Both clustering methods support the identification of three spatially-structured clusters. Individuals from the south and north-western parts of Eurasia form a separated cluster from northern European populations, while central European populations are a mix between the two groups. Our findings have important implications for understanding the dispersal processes that shape the spread of zoonotic diseases under anthropogenic global changes.

Temporal and structural genetic variation in reindeer (Rangifer tarandus) associated with the pastoral transition in Northwestern Siberia.

Just as the domestication of livestock is often cited as a key element in the Neolithic transition to settled, the emergence of large-scaled reindeer husbandry was a fundamental social transformation for the indigenous peoples of Arctic Eurasia. To better understand the history of reindeer domestication, and the genetic processes associated with the pastoral transition in the Eurasian Arctic, we analyzed archaeological and contemporary reindeer samples from Northwestern Siberia. The material represents Rangifer genealogies spanning from 15,000 years ago to the 18th century, as well as modern samples from the wild Taimyr population and from domestic herds managed by Nenetses. The wild and the domestic population are the largest populations of their kind in Northern Eurasia, and some Nenetses hold their domestic reindeer beside their wild cousins. Our analyses of 197 modern and 223 ancient mitochondrial DNA sequences revealed two genetic clusters, which are interpreted as representing the gene pools of contemporary domestic and past wild reindeer. Among a total of 137 different mitochondrial haplotypes identified in both the modern and archaeological samples, only 21 were detected in the modern domestic gene pool, while 11 of these were absent from the wild gene pool. The significant temporal genetic shift that we associate with the pastoral transition suggests that the emergence and spread of reindeer pastoralism in Northwestern Siberia originated with the translocation and subsequent selective breeding of a special type of animal from outside the region. The distinct and persistent domestic characteristics of the haplotype structure since the 18th century suggests little genetic exchange since then. The absence of the typical domestic clade in modern nearby wild populations suggests that the contemporary Nenets domestic breed feature an ancestry from outside its present main distribution, possibly from further South.

Genome sequence and comparative analysis of reindeer (Rangifer tarandus) in northern Eurasia.

Reindeer are semi-domesticated ruminants that have adapted to the challenging northern Eurasian environment characterized by long winters and marked annual fluctuations in daylight. We explored the genetic makeup behind their unique characteristics by de novo sequencing the genome of a male reindeer and conducted gene family analyses with nine other mammalian species. We performed a population genomics study of 23 additional reindeer representing both domestic and wild populations and several ecotypes from various geographic locations. We assembled 2.66 Gb (N50 scaffold of 5 Mb) of the estimated 2.92 Gb reindeer genome, comprising 27,332 genes. The results from the demographic history analysis suggested marked changes in the effective population size of reindeer during the Pleistocene period. We detected 160 reindeer-specific and expanded genes, of which zinc finger proteins (n = 42) and olfactory receptors (n = 13) were the most abundant. Comparative genome analyses revealed several genes that may have promoted the adaptation of reindeer, such as those involved in recombination and speciation (PRDM9), vitamin D metabolism (TRPV5, TRPV6), retinal development (PRDM1, OPN4B), circadian rhythm (GRIA1), immunity (CXCR1, CXCR2, CXCR4, IFNW1), tolerance to cold-triggered pain (SCN11A) and antler development (SILT2). The majority of these characteristic reindeer genes have been reported for the first time here. Moreover, our population genomics analysis suggested at least two independent reindeer domestication events with genetic lineages originating from different refugial regions after the Last Glacial Maximum. Taken together, our study has provided new insights into the domestication, evolution and adaptation of reindeer and has promoted novel genomic research of reindeer.

Opposing fitness consequences of habitat use in a harvested moose population.

Landscape changes are happening at an unprecedented pace, and together with high levels of wildlife harvesting humans have a large effect on wildlife populations. A thorough knowledge of their combined influence on individual fitness is important to understand factors affecting population dynamics. The goal of the study was to assess the individual consistency in the use of risky habitat types, and how habitat use was related to fitness components and life-history strategies. Using data from a closely monitored and harvested population of moose Alces alces, we examined how individual variation in offspring size, reproduction and survival was related to the use of open grasslands; a habitat type that offers high-quality forage during summer, but at the cost of being more exposed to hunters in autumn. The use of this habitat type may therefore involve a trade-off between high mortality risk and forage maximization. There was a high repeatability in habitat use, which suggests consistent behaviour within individuals. Offspring number and weight were positively related to the mothers' use of open grasslands, whereas the probability of surviving the subsequent harvest season was negatively related to the use of the same habitat type. As a consequence, we found a nonsignificant relationship between habitat use and lifetime fitness. The study suggests that harvesting, even if intended to be nonselective with regard to phenotypes, may be selective towards animals with specific behaviour and life-history strategies. As a consequence, harvesting can alter the life-history composition of the population and target life-history strategies that would be beneficial for individual fitness and population growth in the absence of hunting.

Stabilizing selection and adaptive evolution in a combination of two traits in an arctic ungulate.

Stabilizing selection is thought to be common in wild populations and act as one of the main evolutionary mechanisms, which constrain phenotypic variation. When multiple traits interact to create a combined phenotype, correlational selection may be an important process driving adaptive evolution. Here, we report on phenotypic selection and evolutionary changes in two natal traits in a semidomestic population of reindeer (Rangifer tarandus) in northern Finland. The population has been closely monitored since 1969, and detailed data have been collected on individuals since they were born. Over the length of the study period (1969-2015), we found directional and stabilizing selection toward a combination of earlier birth date and heavier birth mass with an intermediate optimum along the major axis of the selection surface. In addition, we demonstrate significant changes in mean traits toward earlier birth date and heavier birth mass, with corresponding genetic changes in breeding values during the study period. Our results demonstrate evolutionary changes in a combination of two traits, which agree closely with estimated patterns of phenotypic selection. Knowledge of the selective surface for combinations of genetically correlated traits are vital to predict how population mean phenotypes and fitness are affected when environments change.

Sea ice loss increases genetic isolation in a high Arctic ungulate metapopulation.

Sea ice loss may have dramatic consequences for population connectivity, extinction-colonization dynamics, and even the persistence of Arctic species subject to climate change. This is of particular concern in face of additional anthropogenic stressors, such as overexploitation. In this study, we assess the population-genetic implications of diminishing sea ice cover in the endemic, high Arctic Svalbard reindeer (Rangifer tarandus platyrhynchus) by analyzing the interactive effects of landscape barriers and reintroductions (following harvest-induced extirpations) on their metapopulation genetic structure. We genotyped 411 wild reindeer from 25 sampling sites throughout the entire subspecies' range at 19 microsatellite loci. Bayesian clustering analysis showed a genetic structure composed of eight populations, of which two were admixed. Overall population genetic differentiation was high (mean FST  = 0.21). Genetic diversity was low (allelic richness [AR] = 2.07-2.58; observed heterozygosity = 0.23-0.43) and declined toward the outer distribution range, where populations showed significant levels of inbreeding. Coalescent estimates of effective population sizes and migration rates revealed strong evolutionary source-sink dynamics with the central population as the main source. The population genetic structure was best explained by a landscape genetics model combining strong isolation by glaciers and open water, and high connectivity by dispersal across winter sea ice. However, the observed patterns of natural isolation were strongly modified by the signature of past harvest-induced extirpations, subsequent reintroductions, and recent lack of sea ice. These results suggest that past and current anthropogenic drivers of metapopulation dynamics may have interactive effects on large-scale ecological and evolutionary processes. Continued loss of sea ice as a dispersal corridor within and between island systems is expected to increase the genetic isolation of populations, and thus threaten the evolutionary potential and persistence of Arctic wildlife.

Chronic wasting disease associated with prion protein gene (PRNP) variation in Norwegian wild reindeer (Rangifer tarandus).

The emergence of CWD in Europe in 2016 and the first natural infection in wild reindeer warranted disease management. This led to the testing of 2424 hunted or culled reindeer during 2016-2018, from the infected subpopulation in the Nordfjella mountain range in Southern Norway. To identify any association between PRNP variation and CWD susceptibility, we characterized the open reading frame of the PRNP gene in 19 CWD positive reindeer and in 101 age category- and sex-matched CWD negative controls. Seven variant positions were identified: 6 single nucleotide variants (SNVs) and a 24 base pair (bp) deletion located between nucleotide position 238 and 272, encoding four instead of five octapeptide repeats. With a single exception, all variant positions but one were predicted to be non-synonymous. The synonymous SNV and the deletion are novel in reindeer. Various combinations of the non-synonymous variant positions resulted in the identification of five PRNP alleles (A-E) that structured into 14 genotypes. We identified an increased CWD risk in reindeer carrying two copies of the most common allele, A, coding for serine in position 225 (Ser225) and in those carrying allele A together with the 24 bp deletion.

Genomic measures of inbreeding in the Norwegian-Swedish Coldblooded Trotter and their associations with known QTL for reproduction and health traits.

BACKGROUND: Since the 1950s, the Norwegian-Swedish Coldblooded trotter (NSCT) has been intensively selected for harness racing performance. As a result, the racing performance of the NSCT has improved remarkably; however, this improved racing performance has also been accompanied by a gradual increase in inbreeding level. Inbreeding in NSCT has historically been monitored by using traditional methods that are based on pedigree analysis, but with recent advancements in genomics, the NSCT industry has shown interest in adopting molecular approaches for the selection and maintenance of this breed. Consequently, the aims of the current study were to estimate genomic-based inbreeding coefficients, i.e. the proportion of runs of homozygosity (ROH), for a sample of NSCT individuals using high-density genotyping array data, and subsequently to compare the resulting rate of genomic-based F (FROH) to that of pedigree-based F (FPED) coefficients within the breed. RESULTS: A total of 566 raced NSCT were available for analyses. Average FROH ranged from 1.78 to 13.95%. Correlations between FROH and FPED were significant (P < 0.001) and ranged from 0.27 to 0.56, with FPED and FROH from 2000 to 2009 increasing by 1.48 and 3.15%, respectively. Comparisons of ROH between individuals yielded 1403 regions that were present in at least 95% of the sampled horses. The average percentage of a single chromosome covered in ROH ranged from 9.84 to 18.82% with chromosome 31 and 18 showing, respectively, the largest and smallest amount of homozygosity. CONCLUSIONS: Genomic inbreeding coefficients were higher than pedigree inbreeding coefficients with both methods showing a gradual increase in inbreeding level in the NSCT breed between 2000 and 2009. Opportunities exist for the NSCT industry to develop programs that provide breeders with easily interpretable feedback on regions of the genome that are suboptimal from the perspective of genetic merit or that are sensitive to inbreeding within the population. The use of molecular data to identify genomic regions that may contribute to inbreeding depression in the NSCT will likely prove to be a valuable tool for the preservation of its genetic diversity in the long term.

First Detection of Chronic Wasting Disease in a Wild Red Deer (Cervus elaphus) in Europe.

Chronic wasting disease (CWD) is a fatal contagious prion disease naturally occurring in cervids in North America. In 2016, CWD was detected in wild reindeer (Rangifer tarandus) and moose (Alces alces) in Norway. Here, we report the first known naturally infected wild Norwegian red deer (Cervus elaphus).

Novel Type of Chronic Wasting Disease Detected in Moose (Alces alces), Norway.

Chronic wasting disease (CWD) persists in cervid populations of North America and in 2016 was detected for the first time in Europe in a wild reindeer in Norway. We report the detection of CWD in 3 moose (Alces alces) in Norway, identified through a large scale surveillance program. The cases occurred in 13-14-year-old female moose, and we detected an abnormal form of prion protein (PrPSc) in the brain but not in lymphoid tissues. Immunohistochemistry revealed that the moose shared the same neuropathologic phenotype, characterized by mostly intraneuronal deposition of PrPSc. This pattern differed from that observed in reindeer and has not been previously reported in CWD-infected cervids. Moreover, Western blot revealed a PrPSc type distinguishable from previous CWD cases and from known ruminant prion diseases in Europe, with the possible exception of sheep CH1641. These findings suggest that these cases in moose represent a novel type of CWD.

A genome-wide association study for harness racing success in the Norwegian-Swedish coldblooded trotter reveals genes for learning and energy metabolism.

BACKGROUND: Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array. RESULTS: Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success. CONCLUSIONS: Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Life histories and ecotype conservation in an adaptive vertebrate: Genetic constitution of piscivorous brown trout covaries with habitat stability.

Ecotype variation in species exhibiting different life history strategies may reflect heritable adaptations to optimize reproductive success, and potential for speciation. Traditionally, ecotypes have, however, been defined by morphometrics and life history characteristics, which may be confounded with individual plasticity. Here, we use the widely distributed and polytypic freshwater fish species brown trout (Salmo trutta) as a model to study piscivorous life history and its genetic characteristics in environmentally contrasting habitats; a large lake ecosystem with one major large and stable tributary, and several small tributaries. Data from 550 fish and 13 polymorphic microsatellites (He = 0.67) indicated ecotype-specific genetic differentiation (θ = 0.0170, p < .0001) among Bayesian assigned small riverine resident and large, lake migrating brown trout (>35 cm), but only in the large tributary. In contrast, large trout did not constitute a distinct genetic group in small tributaries, or across riverine sites. Whereas life history data suggest a small, river resident and a large migratory piscivorous ecotype in all studied tributaries, genetic data indicated that a genetically distinct piscivorous ecotype is more likely to evolve in the large and relatively more stable river habitat. In the smaller tributaries, ecotypes apparently resulted from individual plasticity. Whether different life histories and ecotypes result from individual plasticity or define different genetic types, have important consequence for conservation strategies.

Fitness correlates of age at primiparity in a hunted moose population.

Trade-offs between fitness-related traits are predicted from the principle of resource allocation, where increased fecundity or parental investment leads to reduced future reproduction or survival. However, fitness traits can also be positively correlated due to individual differences (e.g. body mass). Age at primiparity could potentially explain variation in individual fitness either because early primiparity is costly, or it may lead to higher lifetime reproductive success. Based on long-term monitoring and genetic parentage assignment of an island population of moose, we quantified reproductive performance and survival, and examined whether early maturing females have higher total calf production than late maturing females. We explored if harvesting of calves affected the subsequent reproductive success of their mothers, i.e. also due to a post-weaning cost of reproduction, and whether there are any intergenerational effects of female reproductive success. There was a positive relationship between current and future reproduction. The probability to reproduce was lower for females that were unsuccessful the year before, indicating a strong quality effect on productivity. Females that started to reproduce as 2-year olds had a slightly higher total calf production compared to those starting at age three or four. High-performing mothers were also correlated with daughters that performed well in terms of reproductive success. Our results suggest that the observed individual heterogeneity in fitness could be associated with differences in age at primiparity. This heterogeneity was not affected by reproductive costs associated with tending for a calf post-weaning.

Maintaining genetic integrity of coexisting wild and domestic populations: Genetic differentiation between wild and domestic Rangifer with long traditions of intentional interbreeding.

This study investigates the genetic effect of an indigenous tradition of deliberate and controlled interbreeding between wild and domestic Rangifer. The results are interpreted in the context of conservation concerns and debates on the origin of domestic animals. The study is located in Northeastern Zabaikal'e, Russia at approximately 57 degrees North latitude. Blood and skin samples, collected from wild and domestic Rangifer, are analyzed for their mtDNA and microsatellite signatures. Local husbandry traditions are documented ethnographically. The genetic data are analyzed with special reference to indigenous understandings of the distinctions between local domestic types and wild Rangifer. The genetic results demonstrate a strong differentiation between wild and domestic populations. Notably low levels of mtDNA haplotype sharing between wild and domestic reindeer, suggest mainly male-mediated gene flow between the two gene pools. The nuclear microsatellite results also point to distinct differences between regional domestic clusters. Our results indicate that the Evenki herders have an effective breeding technique which, while mixing pedigrees in the short term, guards against wholesale introgression between wild and domestic populations over the long term. They support a model of domestication where wild males and domestic females are selectively interbred, without hybridizing the two populations. Our conclusions inform a debate on the origins of domestication by documenting a situation where both wild and domestic types are in constant interaction. The study further informs a debate in conservation biology by demonstrating that certain types of controlled introgression between wild and domestic types need not reduce genetic diversity.

Lack of significant associations with early career performance suggest no link between the DMRT3 "Gait Keeper" mutation and precocity in Coldblooded trotters.

The Swedish-Norwegian Coldblooded trotter (CBT) is a local breed in Sweden and Norway mainly used for harness racing. Previous studies have shown that a mutation from cytosine (C) to adenine (A) in the doublesex and mab-3 related transcription factor 3 (DMRT3) gene has a major impact on harness racing performance of different breeds. An association of the DMRT3 mutation with early career performance has also been suggested. The aim of the current study was to investigate this proposed association in a randomly selected group of CBTs. 769 CBTs (485 raced, 284 unraced) were genotyped for the DMRT3 mutation. The association with racing performance was investigated for 13 performance traits and three different age intervals: 3 years, 3 to 6 years, and 7 to 10 years of age, using the statistical software R. Each performance trait was analyzed for association with DMRT3 using linear models. The results suggest no association of the DMRT3 mutation with precocity (i.e. performance at 3 years of age). Only two traits (race time and number of disqualifications) were significantly different between the genotypes, with AA horses having the fastest times and CC horses having the highest number of disqualifications at 3 years of age. The frequency of the AA genotype was significantly lower in the raced CBT sample compared with the unraced sample and less than 50% of the AA horses participated in a race. For the age intervals 3 to 6 and 7 to 10 years the AA horses also failed to demonstrate significantly better performance than the other genotypes. Although suggested as the most favorable genotype for racing performance in Standardbreds and Finnhorses across all ages, the AA genotype does not appear to be associated with superior performance, early or late, in the racing career of CBTs.

Phylogenetic Lineages and Postglacial Dispersal Dynamics Characterize the Genetic Structure of the Tick, Ixodes ricinus, in Northwest Europe.

Dispersal and gene flow are important mechanisms affecting the dynamics of vectors and their pathogens. Here, patterns of genetic diversity were analyzed in many North European populations of the tick, Ixodes ricinus. Population sites were selected within and between areas separated by geographical barriers in order to evaluate the importance of tick transportation by birds in producing genetic connectivity across open sea and mountain ranges. The phylogenetic analyses of the mitochondrial control region and the cytochrome b gene revealed two distinct clades with supported sub-clades, with three genetic lineages: GB and WNo associated with Great Britain and western Norway respectively, and Eu with a wider distribution across continental Europe in agreement with much lower efficiency of tick dispersal by birds than by large mammals. The results suggest different ancestry of I. ricinus colonizing Britain and the rest of northern Europe, possibly from different glacial refuges, while ticks from western Norway and continental Europe share a more recent common ancestry. Demographic history modeling suggests a period of strong increase in tick abundance coincident with progression of the European Neolithic culture, long after their post-glacial colonization of NW Europe.