Inferences about the population history of Rangifer tarandus from Y chromosome and mtDNA phylogenies.

Reindeer, called caribou in North America, has a circumpolar distribution and all extant populations belong to the same species (Rangifer tarandus). It has survived the Holocene thanks to its immense adaptability and successful coexistence with humans in different forms of hunting and herding cultures. Here, we examine the paternal and maternal history of Rangifer based on robust Y-chromosomal and mitochondrial DNA (mtDNA) trees representing Eurasian tundra reindeer, Finnish forest reindeer, Svalbard reindeer, Alaska tundra caribou, and woodland caribou. We first assembled Y-chromosomal contigs, representing 1.3 Mb of single-copy Y regions. Based on 545 Y-chromosomal and 458 mtDNA SNPs defined in 55 males, maximum parsimony trees were created. We observed two well separated clades in both phylogenies: the "EuroBeringian clade" formed by animals from Arctic Islands, Eurasia, and a few from North America and the "North American clade" formed only by caribou from North America. The time calibrated Y tree revealed an expansion and dispersal of lineages across continents after the Last Glacial Maximum. We show for the first time unique paternal lineages in Svalbard reindeer and Finnish forest reindeer and reveal a circumscribed Y haplogroup in Fennoscandian tundra reindeer. The Y chromosome in domesticated reindeer is markedly diverse indicating that several male lineages have undergone domestication and less intensive selection on males. This study places R. tarandus onto the list of species with resolved Y and mtDNA phylogenies and builds the basis for studies of the distribution and origin of paternal and maternal lineages in the future.

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.

Chronic wasting disease in Norway-A survey of prion protein gene variation among cervids.

Susceptibility of cervids to Chronic Wasting Disease (CWD), a prion disease, can be modulated by variations in the prion protein gene (PRNP), encoding the cellular prion protein (PrPC ). In prion diseases, PrPC is conformationally converted to pathogenic conformers (PrPSc ), aggregates of which comprise infectious prions. CWD has recently been observed in its contagious form in Norwegian reindeer (Rangifer tarandus) and in novel, potentially sporadic forms, here called 'atypical CWD', in moose (Alces alces) and red deer (Cervus elaphus). To estimate relative susceptibility of different Norwegian cervid species to CWD, their non-synonymous PRNP variants were analyzed. In reindeer, seven PRNP alleles were observed and in red deer and moose two alleles were present, whereas roe deer (Capreolus capreolus) PRNP was monomorphic. One 'archetypal' PRNP allele associated with susceptibility was common to all four cervid species. The distribution of PRNP alleles differed between wild and semi-domesticated reindeer, with alleles associated with a high susceptibility occurring, on average, above 55% in wild reindeer and below 20% in semi-domesticated reindeer. This difference may reflect the diverse origins of the populations and/or selection processes during domestication and breeding. Overall, PRNP genetic data indicate considerable susceptibility to CWD among Norwegian cervids and suggest that PRNP homozygosity may be a risk factor for the atypical CWD observed in moose. The CWD isolates found in the Norwegian cervid species differ from those previously found in Canada and USA. Our study provides an overview of the PRNP genetics in populations exposed to these emerging strains that will provide a basis for understanding these strains' dynamics in relation to PRNP variability.

Mitochondrial DNA in Ixodus ricinus (Acari: Ixodidae) on birds reflects ticks' transportation routes to Lista, Norway.

Ticks are important pathogen vectors, and large mammals and birds have the greatest potential for dispersing them. To study tick dispersal by migrating birds, we have analysed genetic variations in mitochondrial DNA control region from Ixodes ricinus from northward migrating blackbird, Turdus merula, and (European) robin, Erithacus rubecula, at the Lista Bird Observatory in southwestern Norway. We compared their genetic structure with that of resident tick populations from areas covering their expected last stop (i.e. Great Britain and Jutland, Denmark) before taking off for southern Norway, and the resident tick population at Lista. The statistical analysis showed that the I. ricinus found on blackbirds differed significantly from those found on robins, which is consistent with the birds' differential migration routes. I. ricinus from robins did not differ genetically from those flagged at Jutland, suggesting that the former mainly originate in continental Europe. Bayesian analysis indicated that most of the blackbirds caught early in the spring (i.e. before or on the 1st of April) carried ticks of a mixed origin from both Great Britain and continental Europe, while blackbirds caught later in the season carried an increasing amount of ticks acquired locally.

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.

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.

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.

Decomposing demographic contributions to the effective population size with moose as a case study.

Levels of random genetic drift are influenced by demographic factors, such as mating system, sex ratio and age structure. The effective population size (Ne ) is a useful measure for quantifying genetic drift. Evaluating relative contributions of different demographic factors to Ne is therefore important to identify what makes a population vulnerable to loss of genetic variation. Until recently, models for estimating Ne have required many simplifying assumptions, making them unsuitable for this task. Here, using data from a small, harvested moose population, we demonstrate the use of a stochastic demographic framework allowing for fluctuations in both population size and age distribution to estimate and decompose the total demographic variance and hence the ratio of effective to total population size (Ne /N) into components originating from sex, age, survival and reproduction. We not only show which components contribute most to Ne /N currently, but also which components have the greatest potential for changing Ne /N. In this relatively long-lived polygynous system we show that Ne /N is most sensitive to the demographic variance of older males, and that both reproductive autocorrelations (i.e., a tendency for the same individuals to be successful several years in a row) and covariance between survival and reproduction contribute to decreasing Ne /N (increasing genetic drift). These conditions are common in nature and can be caused by common hunting strategies. Thus, the framework presented here has great potential to increase our understanding of the demographic processes that contribute to genetic drift and viability of populations, and to inform management decisions.

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.

Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus).

The reindeer is an Arctic species that exhibits distinctive biological characteristics, for which the underlying genetic basis remains largely unknown. We compared the genomes of reindeer against those of other ruminants and nonruminant mammals to reveal the genetic basis of light arrhythmicity, high vitamin D metabolic efficiency, the antler growth trait of females, and docility. We validate that two reindeer vitamin D metabolic genes (CYP27B1 and POR) show signs of positive selection and exhibit higher catalytic activity than those of other ruminants. A mutation upstream of the reindeer CCND1 gene endows an extra functional binding motif of the androgen receptor and thereby may result in female antlers. Furthermore, a mutation (proline-1172→threonine) in reindeer PER2 results in loss of binding ability with CRY1, which may explain circadian arrhythmicity in reindeer.

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.