INTRODUCTION
The first volume of Thoroughbred Stud Book was published in England in 1791 [
1], and the first volume was published in Korea in 1998. Domestic foals can be listed in the Korean Stud Book (KSB) through appearance screening and paternity verification. In Korea, the paternity of Thoroughbred horses was tested via blood type until the late 1990s but was changed to microsatellite DNA analysis in 2002 [
2].
Microsatellite polymorphism, i.e., short tandem repeat (STR) in a horse, was first reported by Ellegren et al [
3] and Marklund et al [
4]. Microsatellites have been used for genetic diversity assessment, genetic maps constructions, quantitative trait loci mapping and parentage testing [
5,
6]. The STR became a valuable tool for horse parentage verification and pedigree registration [
7]. In 1998, at the Equine Genetic and Thoroughbred parentage Testing workshop hosted by the International Society for Animal Genetics (ISAG), an international panel with nine STRs for horse DNA typing was created, in 2011, three STRs were added, with 12 STRs currently recommended by ISAG. The International Stud Book Committee (ISBC) recommends using these 12 STRs for registration of Thoroughbred horses. It also requires that the total probability of exclusion (PE) value of parentage testing for Thoroughbred horses exceeds 0.9995 [
8].
In Korea, Thoroughbred horse parentage testing started in 2002 using 14 STRs, including the nine (AHT4, AHT5, ASB2, HMS3, HMS6, HMS7, HTG4, HTG10, VHL20) STRs recommended by ISAG. In total, 15 STRs, including 12 (9 STRs plus ASB17, ASB23, HMS2) recommended by ISAG, have been routinely used since 2013.
At present, in South Korea, ~1,000 donkeys and ~27,000 domestic horses exist, including 12,000 Thoroughbreds, 1,000 individuals from other horse breeds (e.g., Warmblood, Quarter horse), and 14,000 native horses (the Jeju Halla horse), of which ~5,000 Jeju horses which was designated as natural monument No. 347 by the government [
9].
In Korea, ≥1,500 Thoroughbred foals are born annually. In order for domestic-born Thoroughbred horses to be used as racehorses, they must be identified through appearance tests or genetic tests to be listed in the KSB. In addition, the horse's microsatellite DNA marker used in Korea was recognized for the registration of Thoroughbred horses from ISAG.
Stallions in Korea, have been replaced a lot so far; hence genetic diversity is also presumed to have changed. Therefore, in this study, we investigated the recent changes such as allele frequencies and PE in Thoroughbred horses in Korea using STR parentage panels between 2006 and 2016.
MATERIALS AND METHODS
Ethical statement
This study was a routine test for the registration of Thoroughbred horses in Korea and not approved by Research Ethical Committee, but sampling of this study was performed according to the guidelines of the Animal Ethics Committee of Korea Racing Authority in Korea for the care and use of experimental animals.
Sample collection and DNA extraction
A total of 5,988 Thoroughbred horses were examined. Genomic DNA was extracted from hair root samples using a MagExtractor System MFX-6100 (Toyobo, Osaka, Japan) according to the manufacturer's instructions [
8].
Microsatellite markers and analysis
A total of 15 microsatellite loci (AHT4, AHT5, ASB2, ASB17, ASB23, CA425, HMS1, HMS2, HMS3, HMS6, HMS7, HTG4, HTG10, LEX3, and VHL20) were used for the analysis of Thoroughbred horses. Detailed characteristics of 15 STRs is as shown in
Table 1. The fluorescent dye attached to the markers used in this study and the marker sizes (DNA base pair) shown in
Figure 1. Polymerase chain reaction (PCR) was performed according to the manufacturer’s protocols (Applied Biosystems, Waltham, MA, USA). Of the 15 markers, HTG10 was examined by a single PCR.
The PCR was performed with a total volume of 15 μL of the following mixture: 40 ng of genomic DNA, 4.0 μL of each primer, 1.25 mM of dNTPs, 2.5 μL of 10× reaction buffer, and 5 U of Taq polymerase (Applied Biosystems, USA). For the single PCR, template DNA 2 μL, 3 pmol forward and reverse primers 2 μL each, and sterile distilled water 6.5 μL were mixed with PCR Premix buffer (Qiagen, Hilden, Germany), and adjusted to 25 μL.
Multiplex PCR amplification was as follows: initial denaturation for 10 min at 95°C, followed by 30 cycles at 95°C for 30 s, 60°C for 30 s and 72°C for 1 min. An extension step at 72°C for 60 min was added after the final cycle [
10]. The single PCR was as follows: initial denaturation for 3 min at 95°C, followed by 30 cycles at 95°C for 45 s, 56°C for 30 s and 72°C for 1 min. An extension step at 72°C for 60 min was added after the final cycle.
The PCR was performed using a GeneAmp PCR System 9700 (Applied Biosystems, USA). The PCR products were electrophoresed and then analyzed using an ABI 3130 xl Genetic Analyzer (Applied Biosystems, USA). The genotype of each of the 15 STRs was determined using Genotyper/Genemapper software Ver.3.7 (Applied Biosystems, USA). Then, with the peak row data, the size of allelic genes (base pair) for each marker was determined based on the result of horse comparison test for STR conducted by ISAG.
Statistical analysis
Allelic frequencies, the number of alleles per locus were estimated by direct counting from observed genotype, and the observed heterozygosity (H
O), expected heterozygosity (H
E), polymorphism information content (PIC), and PE were computed using the Cervus ver. 3.0.3 [
11].
RESULTS
Allele frequencies and expected heterozygosities
In our study, 15 microsatellites were used to identify the genetic diversity of Thoroughbred horses. Allele frequency at each locus is shown in
Table 2. The observed number of alleles per locus ranged from 3 (HMS1) to 9 (ASB17) in 2006 and 4 (HMS1) to 9 (ASB2) in 2016, with a mean value of 6.28 and 6.40, respectively. Of the 15 markers, HMS2, HTG4, and CA425 loci had relatively low PIC (<0.5000) in Thoroughbred horses. Mean levels of genetic variation in 2006 and 2016 were H
O = 0.708 and H
E = 0.685 as well as H
O = 0.699 and H
E = 0.682, respectively (
Tables 3 and
4).
Exclusion probability
The total PE was calculated for each group based on the allele frequencies of 14 or 15 STRs. The 2006 survey analyzed that PE was 0.9998, but it increased to 0.9999 in 2016 after the HMS2 marker was introduced in 2011.
DISCUSSION
Several microsatellite markers have been isolated from the horse genome [
12], and the microsatellites showed multiple alleles and high heterozygosity among European horse breeds such as Thoroughbred horse. Several types of microsatellites are informative because of their high polymorphism, and they are useful in paternity testing of animals. In cattle, pig, horse, dog, and other populations, pedigree control has been performed on a routine basis in most countries [
12–
18]. These controls rely on microsatellite typing standardized through regular comparison tests with aid from the ISAG.
Blood typing tests of horses require fresh whole blood and include the costs of sampling by veterinarians, and transportation due to rapid transportation; the entire test process is manually conducted by humans. The results were also slightly inaccurate, with 97% to 99.5% in probability of parentage verification. However, since 2002, microsatellite DNA typing was introduced and used for the paternity assessment of horses, most of these deficiencies in existing blood typing tests have been resolved. In Korea, a comparative analysis of the number of alleles, frequency of alleles, Ho, He, and PIC with microsatellite markers used in 2006 and 2016 showed that many stallions and mares had been replaced, but genetic diversity did not change significantly because of the small numbers of mares and stallions raised in the country. The total PE of the STR panel has changed significantly. This trend continues after adding HMS2 in 2011. Nevertheless, the total PE values for 14 or 15 STRs in Thoroughbred horses in Korea have been greater than the ISBC recommended values (>0.9995) over the past decade.
In Korea, 15 STR markers are used to conduct parentage verification for approximately 1,500 Thoroughbred horses in a year. Parentage verification is determined using Mendel's genetic law. If a discrepancy occurs in one marker after the test judgment, it is assumed to be a mutation and additionally examined with three to five TKY marker; if two or more markers deviate from Mendel's genetic law, it is finally determined that they are not paternity. In general, approximately 10 to 12 heads (0.7% to 0.8%) per year are not paternity.
On the basis of the PIC of each marker, the validity and reliability of the marker can be estimated. If the PIC value is ≥0.5000, the reliability of the marker is valid for pedigree analysis. If it is ≥0.7000, it has universal validity for analysis and can get a highly reliable result can be obtained.
It is currently under discussion and research to introduce single nucleotide polymorphism (SNP) instead of microsatellite DNA markers for parentage verification and individual identification of Thoroughbred horses. However, the results of this study are expected to be valid until the introduction of the SNP in the Thoroughbred population in Korea.