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.

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.

Identification of Potential miRNA Biomarkers to Detect Hydrocortisone Administration in Horses.

Circulating microRNAs (miRNAs) are stable in bodily fluids and are potential biomarkers of various diseases and physiological states. Although several studies have been conducted on humans to detect drug doping by miRNAs, research on drugs and miRNAs in horses is limited. In this study, circulating miRNAs in horses after hydrocortisone administration were profiled and variations in miRNAs affected by hydrocortisone administration during endogenous hydrocortisone elevation were examined. The miRNAs were extracted from thoroughbred horse plasma before and after hydrocortisone administration and subjected to small RNA sequencing and reverse transcription quantitative PCR (RT-qPCR). RT-qPCR validation was performed for the 20 miRNAs that were most affected by hydrocortisone administration. The effects of elevated endogenous hydrocortisone levels due to exercise and adrenocorticotropic hormone administration were also confirmed. The validation results showed that approximately half of the miRNAs showed the same significant differences as those obtained using small RNA sequencing. Among the twenty miRNAs, two novel miRNAs and miR-133a were found to vary differently between exogenous hydrocortisone administration and endogenous hydrocortisone elevation. This study provides basic knowledge regarding the circulating miRNA profile of horses after hydrocortisone administration and identifies three miRNAs that could potentially be used as biomarkers to detect hydrocortisone administration.

Genetic characterization of Japanese native horse breeds by genotyping variants that are associated with phenotypic traits.

Concerns have been raised about the loss of genetic diversity in Japanese native horses because of their declining populations. In this study, we investigated the genetic variation of four genes, myostatin (MSTN), ligand-dependent nuclear receptor corepressor like (LCORL), doublesex and mab-3 related transcription factor 3 (DMRT3), and 5-hydroxytryptamine receptor 1A (HTR1A), which are associated with horse phenotypic traits, in six Japanese horse breeds (Hokkaido, Kiso, Noma, Misaki, Tokara, and Yonaguni). MSTN, LCORL, DMRT3, and HTR1A showed polymorphisms in the Kiso; Hokkaido and Noma; Hokkaido; and Kiso, Tokara, and Yonaguni breeds, respectively. The Misaki did not show polymorphisms in any of the genes. This study may serve as a basis for developing future breeding strategies focusing on traits in Japanese native horses.

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.

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.

Evaluation of genetic diversity using 31 microsatellites in Miyako horses.

The Miyako horse is a native Japanese horse breed. As with other native Japanese horses, the number of Miyako horses decreased due to mechanization and motorization, which reduced their roles, with just 14 in 1980. Although their population had increased to 55 horses by 2021, a further increase in their numbers is required to avoid extinction. Recently, their breeding has involved natural mating during group grazing; therefore, pedigree management has been difficult, and individual identification has been inconclusive. With the aim of formulating an effective breeding plan, this study used microsatellites to confirm parent-offspring relationships and evaluate the genetic diversity over time. First, the combination of microsatellite genotypes identified misunderstood parent-offspring relationships in 35.3% of the existing individuals, and a correct family tree was reconstructed. Next, the number of alleles and observed and expected values of heterozygosity were calculated separately for the populations during periods of 1998-2012 and 2013-2020. The values were 4.2, 0.705, and 0.653 and 3.9, 0.633, and 0.603, respectively, indicating that genetic diversity according to all indices decreased during period of 2013-2020. This was probably because of the bias of stallions in the 2013-2020 population. Errors in pedigree information in a small population such as Miyako horses could increase the risk of inbreeding, and confirmation of parent-offspring relationships using genotypes may be beneficial. Additionally, to maintain diversity in future breeding, it is important to avoid bias, particularly among stallions, and to ensure offspring of various individuals who are as distantly related to each other as possible.

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.

Correction: Mutations in MITF and PAX3 Cause "Splashed White" and Other White Spotting Phenotypes in Horses.

[This corrects the article DOI: 10.1371/journal.pgen.1002653.].

Genetic diversity analysis and parentage verification of Taishu horses using 31 microsatellites.

The Taishu horse in Tsushima is one of eight Japanese native breeds. The breed is on the verge of extinction due to a rapid decrease in numbers since the 1960s owing to motorization in Japan. In this study, we aimed to confirm the pedigree information of 52 horses by genotyping 31 microsatellites in order to avoid inbreeding. Parentage verification failed to identify genetic contradictions among trios (sires, dams, and foals) registered with the Japan Equine Affairs Association (JEAA). Pedigree information registered at the JEAA was obtained and adequately understood. Additionally, the genetic diversity of the Taishu horses was evaluated and compared with those of other Japanese native breeds. The average values for the number of alleles, observed heterozygosity, expected heterozygosity, and inbreeding coefficient were 4.7, 0.643, 0.632, and -0.02, respectively. Using the Structure software, the 52 horses were classified into three subgroups based on the individuals with more than 50% of specific genetic components. The phylogenetic trees created based on neighbor-joining classification tended to be consistent among the stallions. The effective population size was 27.5 and lower than that required for maintaining 90% genetic variation in the source population over a period of 100 years (47.5). Compared with the other Japanese breeds, the Taishu horse population included in the current study exhibited moderate genetic diversity. Our study will contribute to reconsideration of the breeding strategy of Taishu horses.

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.

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.

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.

Heritability estimates of fractures in Japanese Thoroughbred racehorses using a non-linear model.

Thoroughbred racehorses are produced by mating small numbers of Arabian stallions and native British mares, and have been improved by selection of horseracing performance for about 300 years. While these improvements led to good performance as racehorses, they exposed horses to numerous medical disorders, aggravated by extensive exercise. Fractures are frequent medical disorders in Thoroughbred racehorses. In this study, fracture heritability was estimated using 3,927 Japanese Thoroughbred racehorses to elucidate the risk of racehorse fractures. The heritability estimates of all examined fractures were low (h2  = 0.06), while those of fractures in carpal bone and carpus (carpal bone plus distal radius) were moderate (h2  = 0.37, 0.24, respectively). Fracture occurrence age for carpal bone and distal radius was both 3.3 years old and was younger than that for other fractures. These results indicated that a larger proportion of the variation in the studied population was due to genetic factors for carpal fractures than for other fractures, while the fractures at other bones were largely affected by environmental factors, correlated with the athlete period (number year in racing). These findings contribute to develop a management plan for suppressing racehorse fractures and improving horseracing safety.

Identification of metabolomic changes in horse plasma after racing by liquid chromatography-high resolution mass spectrometry as a strategy for doping testing.

Recently, the illegal use of novel technologies, such as gene and cell therapies, has become a great concern for the horseracing industry. As a potential way to control this, metabolomics approaches that comprehensively analyze metabolites in biological samples have been gaining attention. However, it may be difficult to identify metabolic biomarkers for doping because physiological conditions generally differ between resting and exercise states in horses. To understand the metabolic differences in horse plasma between the resting state at training centres and the sample collection stage after racing for doping test (SAD), we took plasma samples from these two stages (n=30 for each stage) and compared the metabolites present in these samples by liquid chromatography-high resolution mass spectrometry. This analysis identified 5,010 peaks, of which 1,256 peaks (approximately 25%) were annotated using KEGG analysis. Principal component analysis showed that the resting state and SAD groups had entirely different metabolite compositions. In particular, the levels of inosine, xanthosine, uric acid, and allantoin, which are induced by extensive exercise, were significantly increased in the SAD group. In addition, many metabolites not affected by extensive exercise were also identified. These results will contribute to the discovery of biomarkers for detecting doping substances that cannot be detected by conventional methods.

Genetic analysis of Taishu horses on and off Tsushima Island: Implications for conservation.

Taishu horses are a native Japanese breed, of which only 41 individuals remained on Tsushima Island in 2018. Their genetic diversity is considered lower than that of other Japanese native horse breeds; thus, it needs to be investigated for sustainable conservation of this breed. Historical records revealed that several Taishu individuals were released areas off-Tsushima Island in mid-1980s. At present, Taishu horses living outside of Tsushima Island, hereafter referred to as Non-Tsushima Taishus (NTTs), are tagged. However, the genetic structure of the NTT individuals remains unclear, and such individuals are not included in the current mating plans for Taishu horses. Herein, we examined the genetic structure of 18 NTT individuals by comparing their genomic (SNP) information with that of individuals on Tsushima Island (TT), four other native Japanese breeds, and one Anglo-Arabian breed by using ddRAD-seq. We found that all individuals related to the Taishu can be grouped in one cluster, which was separated from other horse breeds. Patterns of specific and shared SNPs in NTT individuals closely resembled those of TT individuals, suggesting very minor genetic differences. Meanwhile, the heterozygosity of NTT individuals was slightly higher than that of TT individuals, and many NTT individuals were of fertile age, suggesting that the pedigree of NTT individuals would be useful in breed conservation plans for Taishu horses. Based on their genomic information, we also reconstructed the pedigree structures of four NTT individuals with no family information. The inclusion of NTT individuals in future mating plans on Tsushima Island may be an effective and feasible method for conserving the Taishu horse breed in Japan.

Genetic relationship between Miyako and Yonaguni horses native to Okinawa based on polymorphisms of microsatellites.

The Miyako and Yonaguni horses are native horses in Okinawa. Here, we evaluated their genetic relationship using microsatellite data and Kiso horses, which have four subpopulations, as a reference population for evaluating this relationship. Microsatellite data from 35 Miyako, 78 Yonaguni, and 172 Kiso horses were evaluated using the STRUCTURE software for analyzing multilocus genotype data to investigate the population structures and their underlying relationship. The results of the STRUCTURE analysis were stable when ΔK was 2, suggesting that the Okinawan horses are different from the Kiso horses. Moreover, the results were also stable when ΔK was 6; the sample was then divided into four subpopulations of the Kiso horses and two Okinawan horse breeds. However, the diagrams from the STRUCTURE analysis were unstable when ΔK was 3. These results suggest that the genetic relationship of the Okinawan horse breeds may be close, similar to that among the subpopulations of the Kiso horses.

Genetic analyses for conservation of the traditional Tokara horse using 31 microsatellite markers.

In order to promote conservation of the traditional Tokara horse in its remaining three breeding areas in Japan (Nakanoshima, Kaimondake, and Iriki), we genotyped 123 horses using 31 microsatellite markers and determined their genetic diversity. On average, the number of alleles (NA), observed heterozygosity (HO), expected heterozygosity (HE), and inbreeding coefficient (FIS) among all horses were 3.0, 0.424, 0.481, and 0.108, respectively. Compared with other endangered horse breeds, we found that, even though the size of the Tokara horse population has recently increased, the NA, HO, and HE of Tokara horses are still notably lower than those of other breeds. Neighbor-joining tree and STRUCTURE analysis showed that the current population of Tokara horses is divided into three subpopulations, corresponding to their respective feeding and breeding areas: Nakanoshima, Kaimondake, and Iriki. This subdivision was also reflected in the NA of microsatellite DNAs, with four, three, and four different loci showing single alleles in Nakanoshima, Kaimondake, and Iriki horses, respectively. These alleles are considered to have become fixed as a consequence of breeding within the limited number of horses in each area. Since Tokara horses are currently strongly divided into subpopulations, it is vitally important to exchange several horses among their different breeding units in order to maintain the genetic diversity of the Tokara horse as a unique breed. The data obtained in this study contribute toward explaining the history of Tokara horses and also provide important information for future monitoring of population diversity and guiding conservation measures for this endangered breed.