Genomic inbreeding trends, influential sire lines and selection in the global Thoroughbred horse population

Pedigree tracing to determine the origins of the modern Palomino Horse Breeders of America stock-type horse

Registration numbers for the Palomino Horse Breeders of America (PHBA) have risen with the assistance of new registration guidelines. These changes in guidelines allowed for registration of any horse with a palomino body color, including those with unknown pedigrees. While this promoted registry growth, documenting the breed's orgins from those known pedigrees making up the breed's foundation is important with the introduction of unknown pedigrees. This can be accomplished using pedigree tracing, which is of particular value for color registries as the study of coat color patterns is plagued with misunderstandings due to pedigrees with inaccurate color designations. Thus, the objective of this study was to develop an understanding of the historical origins of the modern PHBA stock-type horse through pedigree tracing. Horses (n = 460) were randomly sampled from published entries from 10 years (2012-2022) of the PHBA World Show Championship Horse Show. The All-Breed Database was utilized to obtain documentation of extended pedigrees that included uninterrupted generation-to-generation inheritance patterns of a cream dilution allele-influenced coat color. Pedigrees were traced back to the earliest recorded relative. We observed six foundation sires contributed to the modern PHBA stock-type horse with the majority (58%) tracing back to a sire named Old Fred. Two identified sires, Yellow Jacket (4%) and Cotton Eyed Joe (3%), have conflicts concerning color designation. While these results assist in documenting the breed's origins, genetic testing is required to accurately identify the cream dilution allele within the horses registered to ensure a phenotype-based solely on the cream dilution allele.

Analyses of whole-genome sequences from 185 North American Thoroughbred horses, spanning 5 generations

Whole genome sequences (WGS) of 185 North American Thoroughbred horses were compared to quantify the number and frequency of variants, diversity of mitotypes, and autosomal runs of homozygosity (ROH). Of the samples, 82 horses were born between 1965 and 1986 (Group 1); the remaining 103, selected to maximize pedigree diversity, were born between 2000 and 2020 (Group 2). Over 14.3 million autosomal variants were identified with 4.5-5.0 million found per horse. Mitochondrial sequences associated the North American Thoroughbreds with 9 of 17 clades previously identified among diverse breeds. Individual coefficients of inbreeding, estimated from ROH, averaged 0.266 (Group 1) and 0.283 (Group 2). When SNP arrays were simulated using subsets of WGS markers, the arrays over-estimated lengths of ROH. WGS-based estimates of inbreeding were highly correlated (r > 0.98) with SNP array-based estimates, but only moderately correlated (r = 0.40) with inbreeding based on 5-generation pedigrees. On average, Group 1 horses had more heterozygous variants (P < 0.001), more total variants (P < 0.001), and lower individual inbreeding (FROH; P < 0.001) than horses in Group 2. However, the distribution of numbers of variants, allele frequency, and extent of ROH overlapped among all horses such that it was not possible to identify the group of origin of any single horse using these measures. Consequently, the Thoroughbred population would be better monitored by investigating changes in specific variants, rather than relying on broad measures of diversity. The WGS for these 185 horses is publicly available for comparison to other populations and as a foundation for modeling changes in population structure, breeding practices, or the appearance of deleterious variants.

Pedigree analysis of Korean native chickens: unraveling inbreeding and genetic diversity

This study assessed the trends in inbreeding, effective population size, and genetic diversity across six Korean native chicken lines using pedigree records from 54,383 chickens. Understanding these genetic parameters is significantly important for maintaining healthy and viable chicken populations. The primary objective was to analyze the pedigree data to assess the levels of inbreeding and genetic diversity and to evaluate the effective population size across the different lines. Pedigree analysis revealed that pedigree completeness peaked in the first generation and declined in subsequent generations for all lines. Line A exhibited a mean inbreeding coefficient of 0.0201, whereas the other lines displayed lower mean values ranging from 0.0009 to 0.0098, indicating that inbreeding levels were within an acceptable range and considered safe from extinction. Average relatedness consistently increased with time. Individual increases in inbreeding were the highest in Line A (0.62%), with smaller increases in the other lines ranging from 0.02 to 0.23%. Effective population sizes varied from 81 to 2500, with average coancestry within parental populations ranging from 0.0032 to 0.0290. The fe/fa ratio between 1.00 and 1.69 in the 6 lines suggested a moderate impact during bottleneck events, with subsequent populations recovering well. The genetic diversity loss due to genetic drift and unequal founder contributions ranged from 0.66-3.15%, indicating that considerable genetic variability remains within the populations. The results of this study have practical applications in the management and conservation of genetic resources in poultry breeding programs. By highlighting the importance of monitoring inbreeding and maintaining genetic diversity, the findings can help develop strategies to ensure the long-term sustainability of these chicken lines. This study provides valuable insights into the genetic management of Korean native chicken lines, emphasizing the need for strategic breeding practices to preserve genetic health and diversity.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Does inbreeding contribute to pregnancy loss in Thoroughbred horses?

BACKGROUND

Excessive inbreeding increases the probability of uncovering homozygous recessive genotypes and has been associated with an increased risk of retained placenta and lower semen quality. No genomic analysis has investigated the association between inbreeding levels and pregnancy loss.

OBJECTIVES

To compare genetic inbreeding coefficients (F) of naturally occurring Thoroughbred Early Pregnancy Loss (EPLs), Mid and Late term Pregnancy Loss (MLPL) and Controls. The F value was hypothesised to be higher in cases of pregnancy loss (EPLs and MLPLs) than Controls.

STUDY DESIGN

Observational case-control study.

METHODS

Allantochorion and fetal DNA from EPL (n = 37, gestation age 14-65 days), MLPL (n = 94, gestational age 70 days-24 h post parturition) and Controls (n = 58) were genotyped on the Axiom Equine 670K SNP Genotyping Array. Inbreeding coefficients using Runs of Homozygosity (FROH) were calculated using PLINK software. ROHs were split into size categories to investigate the recency of inbreeding.

RESULTS

MLPLs had significantly higher median number of ROH (188 interquartile range [IQR], 180.8-197.3), length of ROH (3.10, IQR 2.93-3.33), and total number of ROH (590.8, IQR 537.3-632.3), and FROH (0.26, IQR 0.24-0.28) when compared with the Controls and the EPLs (p < 0.05). There was no significant difference in any of the inbreeding indices between the EPLs and Controls. The MLPLs had a significantly higher proportion of long (>10 Mb) ROH (2.5%, IQR 1.6-3.6) than the Controls (1.7%, IQR 0.6-2.5), p = 0.001. No unique ROHs were found in the EPL or MLPL populations.

MAIN LIMITATIONS

SNP-array data does not allow analysis of every base in the sequence.

CONCLUSIONS

This first study of the effect of genomic inbreeding levels on pregnancy loss showed that inbreeding is a contributor to MLPL, but not EPL in the UK Thoroughbred population. Mating choices remain critical, because inbreeding may predispose to MLPL by increasing the risk of homozygosity for specific lethal allele(s).

Effect of an endothelial regulatory module on plasma proteomics in exercising horses

Elite performing exercise requires an intricate modulation of the blood pressure to support the working muscles with oxygen. We have previously identified a genomic regulatory module that associates with differences in blood pressures of importance for elite performance in racehorses. This study aimed to determine the effect of the regulatory module on the protein repertoire. We sampled plasma from 12 Coldblooded trotters divided into two endothelial regulatory module haplotype groups, a sub-elite performing haplotype (SPH) and an elite performing haplotype (EPH), each at rest and exercise. The haplotype groups and their interaction were interrogated in two analyses, i) individual paired ratio analysis for identifying differentially abundant proteins of exercise (DAPE) and interaction (DAPI) between haplotype and exercise, and ii) unpaired ratio analysis for identifying differentially abundant protein of haplotype (DAPH). The proteomics analyses revealed a widespread change in plasma protein content during exercise, with a decreased tendency in protein abundance that is mainly related to lung function, tissue fluids, metabolism, calcium ion pathway and cellular energy metabolism. Furthermore, we provide the first investigation of the proteome variation due to the interaction between exercise and related blood pressure haplotypes, which this difference was related to a faster switch to the lipoprotein and lipid metabolism during exercise for EPH. The molecular signatures identified in the present study contribute to an improved understanding of exercise-related blood pressure regulation.

An endothelial regulatory module links blood pressure regulation with elite athletic performance

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Historic Horse Family Displaying Malformations of the Cervicothoracic Junction and Their Connection to Modern German Warmblood Horses

Malformations of the equine cervicothoracic junction affect the C6 and C7 cervical vertebrae, the T1 thoracic vertebra and in variable extent the first and second sternal ribs. To date, the clinical impact of this malformation, its prevalence and mode of inheritance in equine populations are not yet determined. We examined five skeletons for signs of malformation of the cervicothoracic junction, including three skeletons from widely used Thoroughbred stallions affected with the malformation and two skeletons serving as a comparison. The three affected historical horses were the Thoroughbred stallions Der Loewe XX, Birkhahn XX and their common great grandsire Dark Ronald XX. Malformations of C6 and C7 showed a large variation between the three stallions, as Dark Ronald XX, Der Loewe XX and Birkhahn XX were affected uni-laterally at C6 and C7, uni-laterally at C6 and bi-laterally at C6 and C7, respectively, with varying grades. In order to evaluate whether or not these malformations are incidental, we took a random sample of 20 living German Warmblood horses, which are distant descendants of these stallions. This sample consisted of ten controls and ten horses with malformations of C6/C7. Blood proportions of the historical sires in the modern Warmblood horses ranged from 0.10 to 6.25%. The contribution to inbreeding in each individual horse of our selected horse group by those sires was expressed as a percentage of the total inbreeding coefficient and ranged from 0.01 to 17.96%, demonstrating their influence on the modern Warmblood. In the present study, we were able to describe the variability of the malformation of C6/C7 within a horse family including historic and modern horses. Additionally, we detected variations appearing in connection with malformations of the cervicothoracic junction that have not been described in the literature yet. This is the first time that the malformations of C6 and C7 have been described within a familial context, providing hints on inheritance in Thoroughbreds and Warmbloods. It is worthwhile to carry out further studies in a larger setting to gain more comprehensive insights into the inheritance of the malformation and the role of important ancestors.

Influence of horse demographics, country of training and race distance on the rating of Thoroughbreds

The aim of the research was to assess how age, sex, sire, country of foaling, country of training and race distance influenced the international racing and performance of Thoroughbreds. The research was based on performance ratings of 6216 horses assigned by the International Federation of Racing Authorities between 2004 and 2022. The most common sex was stallion (58.54 %) and more than half of the population consisted of 3- and 4-year-old horses (54.68 %). The majority of the horses had the USA as their country of foaling (25.92 %) and also as their country of training (24.87 %). The sire with the largest number of offspring in the International Federation of Horseracing Authorities (IFHA) databases was Galileo (IRE) (193 horses). Four of the 10 most frequently represented sires belonged to the Sadler's Wells (USA) paternal line. The analysis of the statistics in the database as a whole established a significant (p < 0.001 ) influence of all observed factors. Stallions achieved a significantly higher rating (117.85) compared to geldings (117.17) and mares (117.13). The horses originating in Ireland achieved a statistically higher rating (117.99) than horses from Argentina, Australia, Brazil, New Zealand, a group of other countries designated "Others" and South Africa. Statistically conclusive differences were found between horses trained in Ireland (118.80) and all other countries except Great Britain and France. Five of the 10 sires with the best offspring rating belong to the Mr. Prospector (USA) paternal line.

Novel polymorphisms in the prion protein gene (PRNP) and stability of the resultant prion protein in different horse breeds

Prion diseases are fatal neurodegenerative disorders in which the main pathogenic event is the conversion of the cellular prion protein (PrPC) into an abnormal and misfolded isoform known as PrPSc. Most prion diseases and their susceptibility and pathogenesis are mainly modulated by the PRNP gene that codes for PrP. Mutations and polymorphisms in the PRNP gene can alter PrPC amino acid sequence, leading to a change in transmission efficiency depending on the place where it occurs. Horses are animals that are considered to be highly resistant to prions. Several studies have attempted to identify polymorphisms in the PRNP gene that explain the reason for this high resistance. In this study, we have analysed 207 horses from 20 different breeds, discovering 3 novel PRNP polymorphisms. By using computer programmes such as PolyPhen-2, PROVEAN, PANTHER, Meta-SNP and PredictSNP, we have predicted the possible impact that these new polymorphisms would have on the horse prion protein. In addition, we measured the propensity for amyloid aggregation using AMYCO and analysed the lack of hydrogen bridges that these changes would entail together with their electrostatic potentials using Swiss-PdbViewer software, showing that an increased amyloid propensity could be due to changes at the level of electrostatic potentials.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Inbreeding depression and durability in the North American Thoroughbred horse

The proportion of the genome containing runs of homozygosity (ROH) affects production traits in livestock populations. In European and Australasian Thoroughbreds inbreeding, quantified using ROH (F_ROH ), is associated with the probability of ever racing. Here, we measured F_ROH using 333 K SNP genotypes from 768 Thoroughbred horses born in North America to evaluate the effect of inbreeding on racing traits in that region. Among North American horses, F_ROH was not associated (p = 0.518) with the probability of ever racing but was significantly associated with the number of race starts (p = 0.002). Among raced horses, those with a 10% higher F_ROH than the mean inbreeding coefficient were predicted to have 3.5 fewer race starts compared to horses with a mean inbreeding coefficient. Considering the trend of increasing inbreeding and a decline in the average number of race starts per runner in North America, mitigating inbreeding in the population could positively influence racing durability.

Prevalence of Sex-Related Chromosomal Abnormalities in a Large Cohort of Spanish Purebred Horses

Chromosomal abnormalities are largely associated with fertility impairments in the domestic horse. To date, over 600 cases of individuals carrying abnormal chromosome complements have been reported, making the domestic horse the species with the highest prevalence. However, studies analyzing the prevalence of chromosomal diseases in whole populations are scarce. We, therefore, employed a two-step molecular tool to screen and diagnose chromosomal abnormalities in a large population of 25,237 Pura Raza Español horses. Individuals were first screened using short tandem repeats parentage testing results and phenotypic evaluations. Those animals showing results suggesting chromosomal abnormalities were re-tested using a single nucleotide polymorphism (SNP)-based diagnostic methodology to accurately determine the chromosomal complements. Thirteen individuals showed a positive screening, all of which were diagnosed as chromosomally abnormal, including five 64,XY mares with sex development disorders (DSD) and four cases of blood chimerism (two male/female and two female/female cases). In addition, we detected one Turner and one Klinefelter syndrome and two individuals carrying complex karyotypes. The overall prevalence in the entire population was ~0.05%, with the prevalence of 64,XY DSD and blood chimerism ~0.02% and ~0.016%, respectively. However, the overall results should be taken with caution since the individuals carrying Turner syndrome (in full (63,X) or mosaic (mos 63,X/64,XX) forms) cannot be detected due to limitations in the methodology employed. Finally, the lack of agreement between populational studies performed using karyotyping or molecular methods is discussed. To our knowledge, this is the largest populational study performed evaluating the prevalence of the most common chromosomal abnormalities in the domestic horse.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Unlocking Horse Y Chromosome Diversity

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Inbreeding depression and the probability of racing in the Thoroughbred horse

Small effective population sizes and active inbreeding can lead to inbreeding depression due to deleterious recessive mutations exposed in the homozygous state. The Thoroughbred racehorse has low levels of population genetic diversity, but the effects of genomic inbreeding in the population are unknown. Here, we quantified inbreeding based on runs of homozygosity (ROH) using 297 K SNP genotypes from 6128 horses born in Europe and Australia, of which 13.2% were unraced. We show that a 10% increase in inbreeding (F_ROH) is associated with a 7% lower probability of ever racing. Moreover, a ROH-based genome-wide association study identified a haplotype on ECA14 which, in its homozygous state, is linked to a 32.1% lower predicted probability of ever racing, independent of F_ROH. The haplotype overlaps a candidate gene, EFNA5, that is highly expressed in cartilage tissue, which when damaged is one of the most common causes of catastrophic musculoskeletal injury in racehorses. Genomics-informed breeding aiming to reduce inbreeding depression and avoid damaging haplotype carrier matings will improve population health and racehorse welfare.

Fine-Scale Analysis of Runs of Homozygosity Islands Affecting Fertility in Mares

The loss of genetic variability in livestock populations bred under strict selection processes is a growing concern, as it may lead to increased inbreeding values and lower fertility, as a consequence of the “inbreeding depression” effect. This is particularly important in horses, where inbreeding levels tend to rise as individuals become more and more closely related. In this study, we evaluated the effect of increased inbreeding levels on mare fertility by combining an SNP-based genomic approach using runs of homozygosity and the estimation of genetic breeding values for reproductive traits in a large population of Pura Raza Española mares. Our results showed a negative correlation between whole-genome homozygosity and fertility estimated breeding values (EBVs) at the genome level (ρ = −0.144). However, the analysis at chromosome level revealed a wide variability, with some chromosomes showing higher correlations than others. Interestingly, the correlation was stronger (−0.241) when we repeated the analysis in a reduced dataset including the 10% most and least fertile individuals, where the latter showed an increase in average inbreeding values (F_ROH) of around 30%. We also found 41 genomic regions (ROHi, runs of homozygosity islands) where homozygosity increased 100-fold, 13 of which were significantly associated with fertility after cross-validation. These regions encompassed 17 candidate genes previously related to oocyte and embryo development in several species. Overall, we demonstrated the relationship between increased homozygosis at the genomic level and fertility in mares. Our findings may help to deal with the occurrence of inbreeding depression, as well as further our understanding of the mechanisms underlying fertility in mares.

Genomic legacy of migration in endangered caribou

Wide-ranging animals, including migratory species, are significantly threatened by the effects of habitat fragmentation and habitat loss. In the case of terrestrial mammals, this results in nearly a quarter of species being at risk of extinction. Caribou are one such example of a wide-ranging, migratory, terrestrial, and endangered mammal. In populations of caribou, the proportion of individuals considered as "migrants" can vary dramatically. There is therefore a possibility that, under the condition that migratory behavior is genetically determined, those individuals or populations that are migratory will be further impacted by humans, and this impact could result in the permanent loss of the migratory trait in some populations. However, genetic determination of migration has not previously been studied in an endangered terrestrial mammal. We examined migratory behavior of 139 GPS-collared endangered caribou in western North America and carried out genomic scans for the same individuals. Here we determine a genetic subdivision of caribou into a Northern and a Southern genetic cluster. We also detect >50 SNPs associated with migratory behavior, which are in genes with hypothesized roles in determining migration in other organisms. Furthermore, we determine that propensity to migrate depends upon the proportion of ancestry in individual caribou, and thus on the evolutionary history of its migratory and sedentary subspecies. If, as we report, migratory behavior is influenced by genes, caribou could be further impacted by the loss of the migratory trait in some isolated populations already at low numbers. Our results indicating an ancestral genetic component also suggest that the migratory trait and their associated genetic mutations could not be easily re-established when lost in a population.

Y-Chromosomal Insights into Breeding History and Sire Line Genealogies of Arabian Horses

The Y chromosome is a valuable genetic marker for studying the origin and influence of paternal lineages in populations. In this study, we conducted Y-chromosomal lineage-tracing in Arabian horses. First, we resolved a Y haplotype phylogeny based on the next generation sequencing data of 157 males from several breeds. Y-chromosomal haplotypes specific for Arabian horses were inferred by genotyping a collection of 145 males representing most Arabian sire lines that are active around the globe. These lines formed three discrete haplogroups, and the same haplogroups were detected in Arabian populations native to the Middle East. The Arabian haplotypes were clearly distinct from the ones detected in Akhal Tekes, Turkoman horses, and the progeny of two Thoroughbred foundation sires. However, a haplotype introduced into the English Thoroughbred by the stallion Byerley Turk (1680), was shared among Arabians, Turkomans, and Akhal Tekes, which opens a discussion about the historic connections between Oriental horse types. Furthermore, we genetically traced Arabian sire line breeding in the Western World over the past 200 years. This confirmed a strong selection for relatively few male lineages and uncovered incongruences to written pedigree records. Overall, we demonstrate how fine-scaled Y-analysis contributes to a better understanding of the historical development of horse breeds.

Analysis of Genetic Diversity in the American Standardbred Horse Utilizing Short Tandem Repeats and Single Nucleotide Polymorphisms

American Standardbreds were developed as a harness racing horse breed. The United States Trotting Association closed the studbook in 1973 and implemented a book size cap in 2009. This study aimed to investigate genetic diversity in the American Standardbred after the studbook cap was introduced using short tandem repeats (STRs) and single-nucleotide polymorphisms (SNPs). Sixteen STRs from horses foaled from 2010 to 2015 and their sires and dams (n = 50 621) were utilized to examine allelic richness (A_r), expected heterozygosity (H_E), observed heterozygosity (H_O), unbiased heterozygosity (H_U), inbreeding coefficient (F_IS), and fixation index (F_ST). These analyses found that trotting and pacing sires were less genetically diverse than dams (H_EP_Bonferroni = 0.029 and 6.3 × 10⁻⁵, respectively) and their offspring (A_rP_Bonferroni = 0.034 and 6.9 × 10^-6, respectively), and pacing offspring were significantly less diverse than their dams (H_EP_Bonferroni = 2 × 10^-3). Inbreeding coefficients for trotters (F_IS = −0.014) and pacers (F_IS = −0.012) suggest that breeding practices have maintained diversity. Moderate levels of genetic differentiation (0.066 < F_ST < 0.11) were found between pacing and trotting groups. Additionally, 10 of the most prolific trotting sires and their male offspring (n = 84) were genotyped on the 670K Axiom Equine HD Array. H_O values higher than H_E (P < 0.001), low inbreeding coefficients (mean F = −0.064), and mean F_ROH = 21% indicate relatively high levels of diversity in this cohort, further supporting the STR data. However, in contrast, H_O values were higher for trotting sires (0.41) than their offspring (0.36). This observation warrants further monitoring of diversity over time. These data provide an updated foundation of diversity indices for further, long-term analysis in the breed.

16 Years of breed management brings substantial improvement in population genetics of the endangered Cleveland Bay Horse

The consequences of poor breed management and inbreeding can range from gradual declines in individual productivity to more serious fertility and mortality concerns. However, many small and closed groups, as well as larger unmanaged populations, are plagued by genetic regression, often due to a dearth in breeding support tools which are accessible and easy to use in supporting decision-making. To address this, we have developed a population management tool (BCAS, Breed Conservation and Management System) based on individual relatedness assessed using pedigree-based kinship, which offers breeding recommendations for such populations. Moreover, we demonstrate the success of this tool in 16 years of employment in a closed equine population native to the UK, most notably, the rate of inbreeding reducing from more than 3% per generation, to less than 0.5%, or that attributed to genetic drift, as assessed over the last 16 years of implementation. Furthermore, with adherence to this program, the long-term impact of poor management has been reversed and the genetic resource within the breed has grown from an effective population size of 20 in 1994 to more than 140 in 2020. The development and availability of our BCAS for breed management and selection establish a new paradigm for the successful maintenance of genetic resources in animal populations.

Thoroughbred stallion fertility is significantly associated with FKBP6 genotype but not with inbreeding or the contribution of a leading sire

This is a follow-up study to validate the previously detected association of the FKBP6 gene with stallion subfertility. Using a select cohort of 150 Thoroughbred stallions with detailed breeding records, we confirm significant association (P < 0.0001) between low per-cycle pregnancy rates (≤50%) and a combined A/A-A/A genotype of SNPs chr13:11 353 372G>A and chr13:11 353 436A>C in FKBP6 exon 5. We also show that stallion subfertility and the combined genotype A/A-A/A are not associated with the level of genetic diversity based on 12 autosomal microsatellite markers, or with pedigree-based inbreeding rate, or the extent of contribution of a leading Thoroughbred sire, Northern Dancer, in a stallion's pedigree. We develop a TaqMan allelic discrimination assay for the two SNPs to facilitate accurate and high-throughput genotyping. We determine allele, genotype and combined genotype frequencies of FKBP6 exon 5 SNPs in a global cohort of 518 Thoroughbreds (76% stallions or geldings and 24% mares) and show that the frequency of the A/A-A/A genotype is 4%. Because there is no similar association between the FKBP6 exon 5 genotype and stallion subfertility in Hanoverians, we suggest that the two SNPs are not causative but rather tagging a breed-specific haplotype with genetic variants unique to Thoroughbreds. Further WGS-based research is needed to identify the molecular causes underlying the observed genotype-phenotype association in Thoroughbred stallions.

The First Report of Genetic Polymorphisms of the Equine SPRN Gene in Outbred Horses, Jeju and Halla Horses

Simple Summary

Prion disease is a fatal neurodegenerative disease caused by the accumulation of pathogenic prion protein (PrP^Sc) in various mammalian hosts. However, to date, prion disease has not been reported in horses. Since the Sho protein encoded by the shadow of the prion protein gene (SPRN) plays an essential role in the progression of prion diseases, we investigated the genetic characteristics of the equine SPRN gene in horses. We found four single nucleotide polymorphisms (SNPs) of the equine SPRN gene and significant different distributions among three horse breeds including Jeju, Halla and Thoroughbred horses. Although the polymorphisms affect the property of mRNA of the equine SPRN gene, it did not affect the sequence and structure of Sho protein. Since several non-synonymous SNPs of the SPRN gene have been reported in prion diseases-susceptible animals, the absence of non-synonymous SNP of the equine SPRN gene in the horses is noticeable.

Abstract

Prion disease is a fatal infectious disease caused by the accumulation of pathogenic prion protein (PrP^Sc) in several mammals. However, to date, prion disease has not been reported in horses. The Sho protein encoded by the shadow of the prion protein gene (SPRN) plays an essential role in the pathomechanism of prion diseases. To date, the only genetic study of the equine SPRN gene has been reported in the inbred horse, Thoroughbred horse. We first discovered four SPRN single nucleotide polymorphisms (SNPs) in 141 Jeju and 88 Halla horses by direct DNA sequencing. In addition, we found that the genotype, allele and haplotype frequencies of these SNPs of Jeju horses were significantly different from those of Halla and Thoroughbred horses, this latter breed is also included in this study. Furthermore, we observed that the minimum free energy and mRNA secondary structure were significantly different according to haplotypes of equine SPRN polymorphisms by the RNAsnp program. Finally, we compared the SNPs in the coding sequence (CDS) of the SPRN gene between horses and prion disease-susceptible species. Notably, prion disease-susceptible animals had polymorphisms that cause amino acid changes in the open reading frame (ORF) of the SPRN gene, while these polymorphisms were not found in horses.

Screening and detection of chromosomal copy number alterations in the domestic horse using SNP-array genotyping data

Chromosomal abnormalities are a common cause of infertility in horses. However, they are difficult to detect using automated methods. Here, we propose a simple methodology based on single nucleotide polymorphism (SNP)-array data that allows us to detect the main chromosomal abnormalities in horses in a single procedure. As proof of concept, we were able to detect chromosomal abnormalities in 33 out of 268 individuals, including monosomies, chimerisms, and male and female sex-reversions, by analyzing the raw signal intensity produced by an SNP array-based genotyping platform. We also demonstrated that the procedure is not affected by the SNP density of the array employed or by the inbreeding level of the individuals. Finally, the methodology proposed in this study could be performed in an open bioinformatic environment, thus permitting its integration as a flexible screening tool in diagnostic laboratories and genomic breeding programs.

The value of genomic relationship matrices to estimate levels of inbreeding

BACKGROUND

Genomic relationship matrices are used to obtain genomic inbreeding coefficients. However, there are several methodologies to compute these matrices and there is still an unresolved debate on which one provides the best estimate of inbreeding. In this study, we investigated measures of inbreeding obtained from five genomic matrices, including the Nejati-Javaremi allelic relationship matrix (FNEJ), the Li and Horvitz matrix based on excess of homozygosity (FL&H), and the VanRaden (methods 1, FVR1, and 2, FVR2) and Yang (FYAN) genomic relationship matrices. We derived expectations for each inbreeding coefficient, assuming a single locus model, and used these expectations to explain the patterns of the coefficients that were computed from thousands of single nucleotide polymorphism genotypes in a population of Iberian pigs.

RESULTS

Except for FNEJ, the evaluated measures of inbreeding do not match with the original definitions of inbreeding coefficient of Wright (correlation) or Malécot (probability). When inbreeding coefficients are interpreted as indicators of variability (heterozygosity) that was gained or lost relative to a base population, both FNEJ and FL&H led to sensible results but this was not the case for FVR1, FVR2 and FYAN. When variability has increased relative to the base, FVR1, FVR2 and FYAN can indicate that it decreased. In fact, based on FYAN, variability is not expected to increase. When variability has decreased, FVR1 and FVR2 can indicate that it has increased. Finally, these three coefficients can indicate that more variability than that present in the base population can be lost, which is also unreasonable. The patterns for these coefficients observed in the pig population were very different, following the derived expectations. As a consequence, the rate of inbreeding depression estimated based on these inbreeding coefficients differed not only in magnitude but also in sign.

CONCLUSIONS

Genomic inbreeding coefficients obtained from the diagonal elements of genomic matrices can lead to inconsistent results in terms of gain and loss of genetic variability and inbreeding depression estimates, and thus to misleading interpretations. Although these matrices have proven to be very efficient in increasing the accuracy of genomic predictions, they do not always provide a useful measure of inbreeding.

Impaired Reproductive Function in Equines: From Genetics to Genomics

Fertility is one of the key factors in the economic and productive success of the equine industry. Despite this, studies on the genetic causes affecting reproductive performance are scarce, especially in mares, where the genetic architecture of the reproductive traits is extremely complex. Today, with the increasing availability of new genomic methodologies for this species, we are presented with an interesting opportunity to understand the genetic basis of equine reproductive disorders. These include, among others, novel techniques for detecting chromosomal abnormalities, whose association with infertility in horses was established over 50 years ago; new sequencing technologies permitting an accurate detection of point mutations influencing fertility, as well as the study of inbreeding and molecular homozygosity, which has been widely suggested as one of the main causes of low reproductive performance in horses. Finally, over the last few years, reproductive performance has also been associated with copy number variants and candidate genes detected by genome-wide association studies on fertility traits. However, such studies are still scarce, probably because they depend on the existence of large and accurate phenotypic datasets of reproductive and/or fertility traits, which are still difficult to obtain in equines.

Genomics and the Evolutionary History of Equids

The equid family contains only one single extant genus, Equus, including seven living species grouped into horses on the one hand and zebras and asses on the other. In contrast, the equine fossil record shows that an extraordinarily richer diversity existed in the past and provides multiple examples of a highly dynamic evolution punctuated by several waves of explosive radiations and extinctions, cross-continental migrations, and local adaptations. In recent years, genomic technologies have provided new analytical solutions that have enhanced our understanding of equine evolution, including the species radiation within Equus; the extinction dynamics of several lineages; and the domestication history of two individual species, the horse and the donkey. Here, we provide an overview of these recent developments and suggest areas for further research.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.