The Thoroughbred horse is a highly valued domestic animal population under strong selection for athletic phenotypes. Here we present a high resolution genomics-based analysis of inbreeding in the ...population that may form the basis for evidence-based discussion amid concerns in the breeding industry over the increasing use of small numbers of popular sire lines, which may accelerate a loss of genetic diversity. In the most comprehensive globally representative sample of Thoroughbreds to-date (n = 10,118), including prominent stallions (n = 305) from the major bloodstock regions of the world, we show using pan-genomic SNP genotypes that there has been a highly significant decline in global genetic diversity during the last five decades (F
R
= 0.942, P = 2.19 × 10
; F
R
= 0.88, P = 1.81 × 10
) that has likely been influenced by the use of popular sire lines. Estimates of effective population size in the global and regional populations indicate that there is some level of regional variation that may be exploited to improve global genetic diversity. Inbreeding is often a consequence of selection, which in managed animal populations tends to be driven by preferences for cultural, aesthetic or economically advantageous phenotypes. Using a composite selection signals approach, we show that centuries of selection for favourable athletic traits among Thoroughbreds acts on genes with functions in behaviour, musculoskeletal conformation and metabolism. As well as classical selective sweeps at core loci, polygenic adaptation for functional modalities in cardiovascular signalling, organismal growth and development, cellular stress and injury, metabolic pathways and neurotransmitters and other nervous system signalling has shaped the Thoroughbred athletic phenotype. Our results demonstrate that genomics-based approaches to identify genetic outcrosses will add valuable objectivity to augment traditional methods of stallion selection and that genomics-based methods will be beneficial to actively monitor the population to address the marked inbreeding trend.
Abstract
The Thoroughbred breed was formed by crossing Oriental horse breeds and British native horses and is currently used in horseracing worldwide. In this study, we constructed a ...single-nucleotide variant (SNV) database using data from 101 Thoroughbred racehorses. Whole genome sequencing (WGS) revealed 11,570,312 and 602,756 SNVs in autosomal (1–31) and X chromosomes, respectively, yielding a total of 12,173,068 SNVs. About 6.9% of identified SNVs were rare variants observed only in one allele in 101 horses. The number of SNVs detected in individual horses ranged from 4.8 to 5.3 million. Individual horses had a maximum of 25,554 rare variants; several of these were functional variants, such as non-synonymous substitutions, start-gained, start-lost, stop-gained, and stop-lost variants. Therefore, these rare variants may affect differences in traits and phenotypes among individuals. When observing the distribution of rare variants among horses, one breeding stallion had a smaller number of rare variants compared to other horses, suggesting that the frequency of rare variants in the Japanese Thoroughbred population increases through breeding. In addition, our variant database may provide useful basic information for industrial applications, such as the detection of genetically modified racehorses in gene-doping control and pedigree-registration of racehorses using SNVs as markers.
Gene doping has raised concerns in human and equestrian sports and the horseracing industry. There are two possible types of gene doping in the sports and racing industry: (1) administration of a ...gene-doping substance to postnatal animals and (2) generation of genetically engineered animals by modifying eggs. In this study, we aimed to identify genetically engineered animals by whole-genome resequencing (WGR) for gene-doping control. Transgenic cell lines, in which the erythropoietin gene (EPO) cDNA form was inserted into the genome of horse fibroblasts, were constructed as a model of genetically modified horse. Genome-wide screening of non-targeted transgenes was performed to find structural variation using DELLY based on split-read and paired-end algorithms and Control-FREEC based on read-depth algorithm. We detected the EPO transgene as an intron deletion in the WGR data by the split-read algorithm of DELLY. In addition, single-nucleotide polymorphisms and insertions/deletions artificially introduced in the EPO transgene were identified by WGR. Therefore, genome-wide screening using WGR can contribute to gene-doping control even if the targets are unknown. This is the first study to detect transgenes as intron deletions for gene-doping detection.
Gene doping is banned in human sports, horseracing, and equestrian sports. One possible form of gene doping is to administer exogenous genes, called transgenes. Several transgene detection methods ...based on quantitative PCR have been developed. In this study, we investigated the robustness of digital PCR and real-time PCR in transgene detection using primers and probes that matched (P-true) or incompletely matched (P-false) the template DNA. Fluorescence intensity was significantly reduced when substituted probes were used compared to that using the matched probe in both digital and real-time PCR assays. Digital PCR yielded a similar copy number regardless of the probe (P-true: 1230.7, P-false: 1229.7), whereas real-time PCR revealed a decrease in sensitivity based on C q values (P-true: 23.5, P-false: 29.7). When substituted primers were used, the detected copy number decreased in the digital PCR assay, and the C q value in real-time PCR was much higher. Interestingly, digital PCR copy numbers improved by performing PCR at a low annealing temperature, even if a substituted probe was used. Thus, when primer and probe sequences did not completely match the template transgene, digital PCR was relatively robust, but real-time PCR was less sensitive. Although PCR specificity may be reduced, PCR sensitivity can be improved by lowering the annealing temperature. If the target sequence is substituted to escape doping detection, it may be desirable to set the annealing temperature lower and use a more robust method, such as digital PCR, to increase the detection of positive cases, which will also result in fewer false-negative results.
Gene doping, which is prohibited in horseracing and equestrian sports, can be performed by introducing exogenous genes, known as transgenes, into the bodies of postnatal animals. To detect exogenous ...genes, a method utilizing quantitative polymerase chain reaction (qPCR) with a hydrolysis probe was developed to test whole blood and plasma samples, thereby protecting the fairness of competition and the rights of stakeholders in horseracing and equestrian sports. Therefore, we aimed to develop sample storage methods suitable for A and B samples in gene doping tests using blood. For sample A, sufficient qPCR detection was demonstrated after refrigeration for 1 to 2 weeks post collection. For sample B, the following procedures were confirmed to be suitable for storage: 1) centrifugation after sample receipt, 2) frozen storage, 3) natural thawing at room temperature, and 4) centrifugation without mixing blood cell components. Our results indicated that long-term cryopreservation yielded good plasma components from frozen blood samples even though it destroyed blood cells, indicating its applicability to the gene doping test using sample B, which can be stored for later use. Sample storage procedures are as important as detection methods in doping tests. Therefore, the series of procedures that we evaluated in this study will contribute to the efficient performance of gene doping tests through qPCR using blood samples.
Individual identification and paternity testing are important for avoiding inbreeding in the management of small populations of wild and domestic animals. In horse racing industries, they are ...extremely important for identifying and registering individuals and doping control to ensure fair competition. In this study, we constructed an individual identification panel for horses by using insertion and deletion (INDEL) markers. The panel included 39 INDEL markers selected from a whole-genome INDEL database. Genotyping of 89 Thoroughbreds showed polymorphisms with minor allele frequencies (MAFs) of 0.180–0.489 in all markers. The total probability of exclusion for paternity testing, power of discrimination, and probability of identity were 0.9994271269, >0.9999999999, and 0.9999999987, respectively. The panel was applied to 13 trios (sires, dams, and foals), and no contradictions were observed in genetic inheritance among the trios. When this panel was applied to the trios (52 trios) containing false fathers, an average of 7.3 markers excluded parentage relationships. In addition, genomic DNA extracted from the urine of six horses was partially genotyped for 39 markers, and 6–28 markers were successfully genotyped. The newly constructed panel has two advantages: a low marker mutation rate compared with short tandem repeats and a genotyping procedure that is as simple as short tandem repeat typing compared with single nucleotide variant typing. This panel can be applied for individual identification, paternity determination, and urine-sample identification in Thoroughbred horses.
The domestication process of plants and animals typically involves intense inbreeding and directional selection for various traits. Here, we genotyped 370 Japanese Thoroughbred horses using the ...recently developed 670k SNP array and performed various genome-wide analysis also using genotype data of other horse breeds. We identified a number of regions showing interesting patterns of polymorphisms. For instance, the region containing the MC1R locus associated with chestnut coat color may have been targeted by selection for a different mutation much earlier on than the recent selection for chestnut color. We also identified regions that show signatures of selection specific to Thoroughbreds. In addition, we found that intense inbreeding early in the history of the Thoroughbred breed and also before the formation of the breed has a significant impact on the genomic architecture of modern Thoroughbreds. Our study demonstrates that the horse 670k array can be utilized to gain important insight into the domestication process of horses and to understand the genetic basis of the phenotypic diversity in horses.
Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and ...fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.
Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 ...individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.