Rare and common variant discovery by whole-genome sequencing of 101 Thoroughbred racehorses

Analytical advances in horseracing medication and doping control from 2018 to 2023

The analytical approaches taken by laboratories to implement robust and efficient regulation of horseracing medication and doping control are complex and constantly evolving. Each laboratory's approach will be dictated by differences in regulatory, economic and scientific drivers specific to their local environment. However, in general, laboratories will all be undertaking developments and improvements to their screening strategies in order to meet new and emerging threats as well as provide improved service to their customers. In this paper, the published analytical advances in horseracing medication and doping control since the 22nd International Conference of Racing Analysts and Veterinarians will be reviewed. Due to the unprecedented impact of COVID-19 on the worldwide economy, the normal 2-year period of this review was extended to over 5 years. As such, there was considerable ground to cover, resulting in an increase in the number of relevant publications included from 107 to 307. Major trends in publications will be summarised and possible future directions highlighted. This will cover developments in the detection of 'small' and 'large' molecule drugs, sample preparation procedures and the use of alternative matrices, instrumental advances/applications, drug metabolism and pharmacokinetics, the detection and prevalence of 'endogenous' compounds and biomarker and OMICs approaches. Particular emphasis will be given to research into the potential threat of gene doping, which is a significant area of new and continued research for many laboratories. Furthermore, developments in analytical instrumentation relevant to equine medication and doping control will be discussed.

Genetic influence of a STAU2 frameshift mutation and RELN regulatory elements on performance in Icelandic horses

Selection for performance in horse breeding benefits from precise genetic insights at a molecular level, but knowledge remains limited. This study used whole-genome sequences of 39 elite and non-elite Icelandic horses to identify candidate causal variants linked to previously identified haplotypes in the STAU2 and RELN genes affecting pace and other gaits. A frameshift variant in linkage disequilibrium with the previously identified haplotypes in the STAU2 gene (r2 = 0.85) was identified within a predicted STAU2 transcript. This variant alters the amino acid sequence and introduces a premature stop codon but does not appear harmful or disease-causing and is potentially unique to equine biology. A large portion of the RELN haplotype overlapped with an H3K27me3 modification mark, suggesting a regulatory role of this region. Despite the small sample size, the RELN haplotype's effects were validated for tölt, trot, and canter/gallop. Additionally, the RELN haplotype significantly influenced the age at which horses were presented for breeding field tests, indicating a potential role of the region in precocity and trainability. Functional experiments are needed to further investigate the regions' influences on biological processes and their potential impact on horse performance.

Genomic regions and candidate genes associated with forehead whorl positioning in Thoroughbred horses

Previous studies have examined the relationship between hair characteristics and temperament traits in various animals. A partial genetic association has been suggested in humans because whorl formation and neurological development occur simultaneously during the fetal period. In the horse industry, anecdotal evidence suggests a link between the forehead whorl position and temperament. In our previous study, the heritability of forehead whorl positioning was h2=0.653, indicating significant genetic contributions in thoroughbreds. Therefore, in this study, we designed a genome-wide association study using 192 Thoroughbred horses to detect candidate genes associated with forehead whorl positioning. The results revealed 11 suggestive markers on chromosomes 2, 14, 15, 19, and 26, although no significant markers were found. Five genes, namely PTTG1, CCNH, RASA1, COX7C, and CLDN1, located near these markers had functions related to skin cell or hair follicle development and temperament. Therefore, these genes may be involved in the positioning of the forehead whorl in Thoroughbreds. The candidate genes identified in this study are expected to contribute to elucidating the genetic factors behind forehead whorl in horses and to help interpret the common process of whorl and temperament development.

A missense mutation in the KCNE4 gene is not predictive of equine anhidrosis

Anhidrosis is defined as a decreased or absent ability to sweat in response to heat and exercise. In horses, this condition can increase the risk of life-threatening hyperthermia. A prior study has suggested that equine anhidrosis is associated with a missense variant (rs68643109) in the Potassium Voltage-Gated Channel Subfamily E Regulatory Subunit 4 (KCNE4) gene. This project aimed to validate this association in a population of well-phenotyped horses and to determine the allele frequency of this variant in publicly available whole-genome sequence data. Fifty horses within the University of California Davis Center for Equine Health herd were evaluated for anhidrosis using a series of intradermal terbutaline injections. From existing whole-genome sequence data, the rs68643109 genotype of each horse was identified. When stimulated with terbutaline, all 50 horses produced sweat. All three genotypes at rs68643109 were present in this population of horses; the allele previously associated with anhidrosis (G) was present at a frequency of 0.72. No statistical difference in total sweat score was found (p = 0.31). In whole-genome sequences from 820 other horses reported across three prior studies, the alternative (candidate) allele frequency was similarly high, ranging from 0.52 to 0.68. Since all 50 horses tested in our population produced sweat regardless of genotype, and the previously associated allele is present at a high frequency across datasets, these data fail to validate the missense variant within the KCNE4 gene as causative of or contributing to equine anhidrosis.

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.

Non-Synonymous Substitutions in Cadherin 13, Solute Carrier Family 6 Member 4, and Monoamine Oxidase A Genes are Associated with Personality Traits in Thoroughbred Horses

Retraining retired racehorses for various purposes can help correct behavioral issues. However, ensuring efficiency and preventing accidents present global challenges. Based on the hypothesis that a simple personality assessment could help address these challenges, the present study aimed to identify genetic markers associated with personality. Eight genes were selected from 18 personality-related candidate genes that are orthologs of human personality genes, and their association with personality was verified based on actual behavior. A total of 169 Thoroughbred horses were assessed for their tractability (questionnaire concerning tractability in 14 types of situations and 3 types of impressions) during the training process. Personality factors were extracted from the data using principal component analysis and analyzed for their association with single nucleotide variants as non-synonymous substitutions in the target genes. Three genes, CDH13, SLC6A4, and MAOA, demonstrated significant associations based on simple linear regression, marking the identification of these genes for the first time as contributors to temperament in Thoroughbred horses. All these genes, as well as the previously identified HTR1A, are involved in the serotonin neurotransmitter system, suggesting that the tractability of horses may be correlated with their social personality. Assessing the genotypes of these genes before retraining is expected to prevent problems in the development of a racehorse's second career and shorten the training period through individual customization of training methods, thereby improving racehorse welfare.

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.

A Functional Single-Nucleotide Polymorphism Upstream of the Collagen Type III Gene Is Associated with Catastrophic Fracture Risk in Thoroughbred Horses

Fractures caused by bone overloading are a leading cause of euthanasia in Thoroughbred racehorses. The risk of fatal fracture has been shown to be influenced by both environmental and genetic factors but, to date, no specific genetic mechanisms underpinning fractures have been identified. In this study, we utilised a genome-wide polygenic risk score to establish an in vitro cell system to study bone gene regulation in horses at high and low genetic risk of fracture. Candidate gene expression analysis revealed differential expression of COL3A1 and STAT1 genes in osteoblasts derived from high- and low-risk horses. Whole-genome sequencing of two fracture cases and two control horses revealed a single-nucleotide polymorphism (SNP) upstream of COL3A1 that was confirmed in a larger cohort to be significantly associated with fractures. Bioinformatics tools predicted that this SNP may impact the binding of the transcription factor SOX11. Gene modulation demonstrated SOX11 is upstream of COL3A1, and the region binds to nuclear proteins. Furthermore, luciferase assays demonstrated that the region containing the SNP has promoter activity. However, the specific effect of the SNP depends on the broader genetic background of the cells and suggests other factors may also be involved in regulating COL3A1 expression. In conclusion, we have identified a novel SNP that is significantly associated with fracture risk and provide new insights into the regulation of the COL3A1 gene.

Construction of an individual identification panel for horses using insertion and deletion markers

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.

Short Insertion and Deletion Discoveries via Whole-Genome Sequencing of 101 Thoroughbred Racehorses

Thoroughbreds are some of the most famous racehorses worldwide and are currently animals of high economic value. To understand genomic variability in Thoroughbreds, we identified genome-wide insertions and deletions (INDELs) and obtained their allele frequencies in this study. INDELs were obtained from whole-genome sequencing data of 101 Thoroughbred racehorses by mapping sequence reads to the horse reference genome. By integrating individual data, 1,453,349 and 113,047 INDELs were identified in the autosomal (1–31) and X chromosomes, respectively, while 18 INDELs were identified on the mitochondrial genome, totaling 1,566,414 INDELs. Of those, 779,457 loci (49.8%) were novel INDELs, while 786,957 loci (50.2%) were already registered in Ensembl. The sizes of diallelic INDELs ranged from −286 to +476, and the majority, 717,736 (52.14%) and 220,672 (16.03%), were 1-bp and 2-bp variants, respectively. Numerous INDELs were found to have lower frequencies of alternative (Alt) alleles. Many rare variants with low Alt allele frequencies (<0.5%) were also detected. In addition, 5955 loci were genotyped as having a minor allele frequency of 0.5 and being heterogeneous genotypes in all the horses. While short-read sequencing and its mapping to reference genome is a simple way of detecting variants, fake variants may be detected. Therefore, our data help to identify true variants in Thoroughbred horses. The INDEL database we constructed will provide useful information for genetic studies and industrial applications in Thoroughbred horses, including a gene-editing test for gene-doping control and a parentage test using INDELs for horse registration and identification.

Identification of Personality-Related Candidate Genes in Thoroughbred Racehorses Using a Bioinformatics-Based Approach Involving Functionally Annotated Human Genes

Considering the personality traits of racehorses (e.g., flightiness, anxiety, and affability) is considered essential to improve training efficiency and decrease accident frequency, especially when retraining for a second career that may involve contact with inexperienced personnel after retiring from racing. Studies on human personality-related genes are frequently conducted; however, such studies are rare in horses because a consistent methodology for personality evaluation is lacking. Using the recently published whole genome variant database of 101 Thoroughbred horses, we compared horse genes orthologous to human genes related to the Big Five personality traits, and identified 18 personality-related candidate genes in horses. These genes include 55 variants that involve non-synonymous substitutions that highly impact the encoded protein. Moreover, we evaluated the allele frequencies and functional impact on the proteins in terms of the difference in molecular weights and hydrophobicity levels between reference and altered amino acids. We identified 15 newly discovered genes that may affect equine personality, but their associations with personality are still unclear. Although more studies are required to compare genetic and behavioral information to validate this approach, it may be useful under limited conditions for personality evaluation.

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

Detection of Indiscriminate Genetic Manipulation in Thoroughbred Racehorses by Targeted Resequencing for Gene-Doping Control

The creation of genetically modified horses is prohibited in horse racing as it falls under the banner of gene doping. In this study, we developed a test to detect gene editing based on amplicon sequencing using next-generation sequencing (NGS). We designed 1012 amplicons to target 52 genes (481 exons) and 147 single-nucleotide variants (SNVs). NGS analyses showed that 97.7% of the targeted exons were sequenced to sufficient coverage (depth > 50) for calling variants. The targets of artificial editing were defined as homozygous alternative (HomoALT) and compound heterozygous alternative (ALT1/ALT2) insertion/deletion (INDEL) mutations in this study. Four models of gene editing (three homoALT with 1-bp insertions, one REF/ALT with 77-bp deletion) were constructed by editing the myostatin gene in horse fibroblasts using CRISPR/Cas9. The edited cells and 101 samples from thoroughbred horses were screened using the developed test, which was capable of identifying the three homoALT cells containing 1-bp insertions. Furthermore, 147 SNVs were investigated for their utility in confirming biological parentage. Of these, 120 SNVs were amenable to consistent and accurate genotyping. Surrogate (nonbiological) dams were excluded by 9.8 SNVs on average, indicating that the 120 SNV could be used to detect foals that have been produced by somatic cloning or embryo transfer, two practices that are prohibited in thoroughbred racing and breeding. These results indicate that gene-editing tests that include variant calling and SNV genotyping are useful to identify genetically modified racehorses.

Identification of processed pseudogenes in the genome of Thoroughbred horses: Possibility of gene-doping detection considering the presence of pseudogenes

Processed pseudogenes, also known as retrocopy genes, are copies of messenger RNAs that have been reverse transcribed into DNA and inserted into the genome. In this study, we identified 62 processed pseudogene candidates as intron-less genes from whole-genome sequencing (WGS) data of Thoroughbred horses using delly structural variation software. The 62 processed pseudogene candidates were confirmed by PCR amplification of intron-less products. A total of 11 processed pseudogenes were confirmed in the genome of all 23 analysed horses, whereas three processed pseudogenes with structures of ATP11B, DPH3 and RPL17 were detected in only one of 115 horses by PCR amplification of intron-less products. Currently, most of the gene doping tests proposed in human and horse sports are adapted PCR-based methods using hydrolysis probes to detect exon/exon junctions in transgenes because the operation is simple and economical. However, when the pseudogene is present in the host genome, the PCR-based methods may have a potential risk of detecting false positives. In this study, because processed pseudogenes that exist less frequently in the horse genome may affect PCR-based transgene detection in gene-doping tests, we propose and demonstrate that PCR amplification and sequencing using primers designed on transgene and promotors and/or polyadenylation signal for gene expression are useful for gene-doping detection as an additional confirmatory test to prevent false positives.

Design and storage stability of reference materials for microfluidic quantitative PCR-based equine gene doping tests

One method of gene doping in horseracing is administering of exogenous genetic materials, known as transgenes. Several polymerase chain reaction (PCR)-based methods have been developed for detecting transgenes with high sensitivity and specificity. However, novel designs for reference materials (RMs) and/or positive template controls (PTCs) are necessary for simultaneous analysis of multiple transgene targets. In this study, we designed and developed a novel RM for simultaneously detecting multiple targets via microfluidic quantitative PCR (MFQPCR). Twelve equine genes were selected as targets in this study. A sequence region including primers and probes for quantitative PCR was designed, and a 10 bp sequence was inserted to allow the RM to be distinguished from the original transgene sequences. The sequences of individual detection sites were then connected for 12 genes and cloned into a single plasmid vector. We performed fragment size analysis to distinguish between the PCR products of the original transgene sequence and those of the RM, enabling identification of RM contamination. PTCs diluted to 10,000, 1,000, 100, and 10 copies/µl with horse genomic DNA from RM were stably stored at 4°C for 1 year. As digital PCR enabled absolute quantification, the designed substances can serve as an RM. These findings indicate that the RM design and storage conditions were suitable for gene doping tests using MFQPCR.