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.

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.

Evaluation of parentage testing using single nucleotide polymorphism markers for draft horses in Japan.

We evaluated the utility of single nucleotide polymorphism (SNP) markers for parentage testing in Breton (BR) and Percheron (PR) horses in Japan using the proposed International Society for Animal Genetics (P-ISAG) 147 SNP panel and 414 autosomal SNPs. Genomic DNA was extracted from 98 horses of two breeds, BR (n = 47) and PR (n = 51), and sequenced using next-generation sequencing. The average minor allele frequencies for the P-ISAG panel for BR and PR were 0.306 and 0.301, respectively. The combined probabilities of exclusion (PEs) given two parents and one offspring: exclude a relationship (PE01) and given one parent and one offspring: exclude their relationship (PE02) were over 0.9999 for both breeds. Using the P-ISAG panel, no exclusion or doubtful cases were identified in 35 valid parent-offspring pairs, suggesting that the P-ISAG panel is helpful for parentage verification in both breeds. In contrast, as 0.18% of falsely accepted parentages were observed in the parentage discovery cases, additional markers such as the combination of the P-ISAG panel and 414 autosomal SNPs (561-SNP set) presented here should be used to identify valid parent-offspring pairs of horses with unknown parentage relationships.

Investigation of optimal procedures for storage and use of plasma samples suitable for gene doping tests.

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.

Multiplex Detection of Transgenes Using πCode Technology for Gene Doping Control.

To ensure fair competition and sports integrity, gene doping is prohibited in horseracing and equine sports. One gene doping method is by administering exogenous genes, called transgenes, to postnatal animals. Although several transgene detection methods have been developed for horses, many are unsuitable for multiplex detection. In this proof-of-concept study, we developed a highly sensitive and multiplex transgene detection method using multiple πCode with identification patterns printed on the surface. The following steps were employed: (1) multiplex polymerase chain reaction amplification of 12 targeted transgenes in a single tube, (2) detection using a mixture of 12 probes labeled with different πCodes, and (3) median fluorescence intensity measurement of fluorescent πCodes. Twelve transgenes cloned into plasmid vectors were targeted, and 1500 copies of each plasmid were spiked into 1.5 mL of horse plasma. Subsequently, a novel method using πCode succeeded in detecting all the transgenes using their DNA extracts. Additionally, we detected the erythropoietin (EPO) transgene in blood samples from a horse administered solely with the EPO transgene using this method. Therefore, the πCode detection method is considered suitable for multitarget gene detection in gene doping tests.

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.

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.

Use of mitochondrial sequencing to detect gene doping in horses via gene editing and somatic cell nuclear transfer.

Gene editing and subsequent cloning techniques offer great potential not only in genetic disease correction in domestic animals but also in livestock production by enhancement of desirable traits. The existence of the technology, however, leaves it open to potential misuse in performance-led sports such as horseracing and other equestrian events. Recent advances in equine gene editing, regarding the generation of gene-edited embryos using CRISPR/Cas9 technology and somatic cell nuclear transfer, have highlighted the need to develop tools to detect potential prohibited use of the technology. One possible method involves the characterisation of the mitochondrial genome (which is not routinely preserved during cloning) and comparing it with the sequence of the registered dam. We present here our approach to whole-mitochondrial sequencing using tiled long-range PCR and next-generation sequencing. To determine whether the background mutation rate in the mitochondrial genome could potentially confound results, we sequenced 10 sets of dam and foal duos. We found variation between duos but none within duos, indicating that this method is feasible for future screening systems. Analysis of WGS data from over 100 Thoroughbred horses revealed wide variation in the mitochondria sequence within the breed, further displaying the utility of this approach.

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.

Low-copy transgene detection using nested digital polymerase chain reaction for gene-doping control.

Gene doping is prohibited for fair competition in human and horse sports. One style of gene doping is the administration of an exogeneous gene, called a transgene, to postnatal humans and horses. Although many transgene detection methods based on quantitative polymerase chain reaction (PCR), including real-time PCR and digital PCR, have been recently developed, it remains difficult to reliably detect low-copy transgenes. In this study, we developed and validated a nested digital PCR method to specifically detect low-copy transgenes. The nested digital PCR consists of (1) preamplification using conventional PCR and (2) droplet digital PCR detection using a hydrolysis probe. Using 5, 10, 20, 60 and 120 transgene copies as template, 496.0, 1089.7, 1820.7, 4313.3 and 7840.0 copies per microlitre, respectively, were detected using our nested digital PCR. Although high concentrations of phenol, proteinase K, ethanol, EDTA, heparin and genomic DNA all inhibited preamplification, their effects on the digital PCR detection were limited. Once preamplification was successful, even substitution of bases within the primers and probes had minimal effects on transgene detection. The nested digital PCR developed in this study successfully detected low-copy transgenes and can be used to perform a qualitative test, indicating its usefulness in the prevention of false positives and false negatives in gene-doping detection.

Sequence determination of phosphorothioated oligonucleotides using MALDI-TOF mass spectrometry for controlling gene doping in equestrian sports.

In human and equestrian sporting events, one method of gene doping is the illegal use of therapeutic oligonucleotides to alter gene expression. In this study, we aimed to identify therapeutic oligonucleotides via sequencing using matrix-assisted laser desorption/ionisation-time-of-flight mass spectrometry (MALDI-TOF MS). As a model of therapeutic oligonucleotides, 22 bp-long phosphorothioated oligonucleotides (PSOs) were used. By using a Clarity OTX kit for extracting short-length oligonucleotides, a spectrum of singly charged PSO with a mean intensity of 6.08 × 104 (standard deviation: 4.34 × 103 ) was detected from 500 pmol PSO in 1 ml horse plasma using the linear negative mode of MALDI-TOF MS. In addition, a 17 bp sequence was determined using in-source decay (ISD) mode, indicating that 500 pmol of a PSO in 1 ml plasma is the detection limit for sequencing. Using the determined sequences (17 bp), a targeted gene for PSO was singly identified on the horse reference genome, EquCab2.0, via a GGGenome search. These procedures can be potentially used to identify therapeutic oligonucleotides, whose nucleotides are unknown, for gene doping control.

Two secured FACT recruitment mechanisms are essential for heterochromatin maintenance.

FACT (facilitate chromatin transcription) is involved in heterochromatic silencing, but its mechanisms and function remain unclear. We reveal that the Spt16 recruitment mechanism operates in two distinct ways in heterochromatin. First, Pob3 mediates Spt16 recruitment onto the heterochromatin through its Spt16 dimerization and tandem PH domains. Without Pob3, Spt16 recruitment is partially reduced, exhibiting a silencing defect and impaired H2A/H2B organization. Second, heterochromatin protein 1 (HP1)/Swi6 mediates Spt16 recruitment onto the heterochromatin by physical interaction of the Swi6 chromo-shadow domain (CSD) and Spt16 peptidase-like domains. Several CSD mutants are tested for Spt16 binding activity, and the charged loop connecting ß1 and ß2 is critical for Spt16 binding and heterochromatic silencing. Loss of these pathways causes a severe defect in H3K9 methylation and HP1/Swi6 localization in the pericentromeric region, exhibiting transcriptional silencing defects and disordered heterochromatin. Our findings suggest that FACT and HP1/Swi6 work intimately to regulate heterochromatin organization.

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

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.

Robustness of digital PCR and real-time PCR against inhibitors in transgene detection for gene doping control in equestrian sports.

Gene doping is a threat to fair competition in sports, both human and equestrian. One method of gene doping is to administer exogenous genetic materials, called transgenes, into the bodies of postnatal humans and horses. Polymerase chain reaction (PCR)-based transgene detection methods such as digital PCR and real-time PCR have been developed for gene doping testing in humans and horses. However, the significance of PCR inhibitors in gene doping testing has not been well evaluated. In this study, we evaluated the effects of PCR inhibitors on transgene detection using digital PCR and real-time PCR against gene doping. Digital PCR amplification was significantly inhibited by high concentrations of proteinase K (more than 0.1 µg/µl), ethylenediaminetetraacetic acid (more than 5 nmol/µl), and heparin (more than 0.05 unit/µl) but not by ethanol or genomic DNA. In addition, phenol affected droplet formation in the digital PCR amplification process. Real-time PCR amplification was inhibited by high concentrations of phenol (more than 1% v/v), proteinase K (more than 0.001 µg/µl), ethylenediaminetetraacetic acid (more than 1 nmol/µl), heparin (more than 0.005 unit/µl), and genomic DNA (more than 51.9 ng/µl) but not by ethanol. Although both PCR systems were inhibited by nearly the same substances, digital PCR was more robust than real-time PCR against the inhibitors. We believe that our findings are important for the development of better methods for transgene detection and prevention of false negative results in gene doping testing.

Robustness of Digital PCR and Real-Time PCR in Transgene Detection for Gene-Doping Control.

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 Cq 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 Cq 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.

Detection of non-targeted transgenes by whole-genome resequencing for gene-doping control.

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.

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.

Whole-genome resequencing using genomic DNA extracted from horsehair roots for gene-doping control in horse sports.

Gene doping is prohibited in horseracing and equestrian sports. In previous studies, we developed non-targeted transgene and genome editing detection methods based on whole genome resequencing (WGR) using genomic DNA extracted from whole blood. In this study, we aimed to develop a WGR method using DNA extracts from hair roots. Hair roots are a preferred substrate because their collection is less invasive than blood collection. Hair is also easier to store for long periods of time. Although almost all genomic DNA extracted from hair root samples stored for years at room temperature was degraded, the quality of genomic DNA from samples stored for years at refrigerated temperatures (4-8°C) was maintained. High-molecular-weight genomic DNA was isolated from hair roots using a magnetic silica beads method of extraction, enabling WGR from horsehair root extracts. Nucleotide sequencing results and numbers of single-nucleotide polymorphisms and insertions/deletions concurred with those previously reported for WGR of DNA extracted from whole blood. Therefore, we consider that storing hair samples at refrigerated temperatures prevents degradation of DNA, allowing the detection of gene doping in these samples based on WGR. It is likely this finding will also have a deterrent effect, as it is now possible to test horses with archived samples even if they or their parents are deceased. To our knowledge, this is the first report employing WGR on horsehair roots stored for a long term.