Thoroughbred racehorse mitochondrial DNA demonstrates closer than expected links between maternal genetic history and pedigree records

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.

Influence of horse demographics, country of training and race distance on the rating of Thoroughbreds

The aim of the research was to assess how age, sex, sire, country of foaling, country of training and race distance influenced the international racing and performance of Thoroughbreds. The research was based on performance ratings of 6216 horses assigned by the International Federation of Racing Authorities between 2004 and 2022. The most common sex was stallion (58.54 %) and more than half of the population consisted of 3- and 4-year-old horses (54.68 %). The majority of the horses had the USA as their country of foaling (25.92 %) and also as their country of training (24.87 %). The sire with the largest number of offspring in the International Federation of Horseracing Authorities (IFHA) databases was Galileo (IRE) (193 horses). Four of the 10 most frequently represented sires belonged to the Sadler's Wells (USA) paternal line. The analysis of the statistics in the database as a whole established a significant (p < 0.001 ) influence of all observed factors. Stallions achieved a significantly higher rating (117.85) compared to geldings (117.17) and mares (117.13). The horses originating in Ireland achieved a statistically higher rating (117.99) than horses from Argentina, Australia, Brazil, New Zealand, a group of other countries designated "Others" and South Africa. Statistically conclusive differences were found between horses trained in Ireland (118.80) and all other countries except Great Britain and France. Five of the 10 sires with the best offspring rating belong to the Mr. Prospector (USA) paternal line.

mtDNA CR Evidence Indicates High Genetic Diversity of Captive Forest Musk Deer in Shaanxi Province, China

Forest musk deer (Moschus berezovskii) are endangered ruminants whose adult males secrete musk. China has been breeding forest musk deer artificially since the 1950s in an effort to restore wild populations, with Shaanxi and Sichuan provinces as the two main sites for captive breeding. Genetic diversity is a significant indicator that determines the long-term viability and status of a population, particularly for species at risk of extinction. In this study, we analyzed the current genetic makeup of seven captive forest musk deer populations in the Shaanxi province, using the mitochondrial DNA (mtDNA) control region (CR) as the molecular marker. We sequenced 604 bp of mtDNA CR, with an average content of A+T higher than G+C. We observed 111 variable sites and 39 different haplotypes from 338 sequences. The nucleotide diversity (Pi) and haplotype diversity (Hd) were 0.02887 and 0.908, respectively. Genetic differentiation between these populations was not significant, and the populations might not have experienced rapid growth. By combining our sequences with previous ones, we identified 65 unique haplotypes with 26 rare haplotypes and estimated a total of 90 haplotypes in Shaanxi province captive populations. The Shaanxi province and Sichuan province obtained 88 haplotypes, the haplotypes from the two populations were mixed together, and the two populations showed moderate genetic differentiation. Our findings suggested that captive forest musk deer populations in the Shaanxi province had high genetic diversity, with a rich founder population of about 90 maternal lines. Additionally, managers could develop genetic management plans for forest musk deer based on the haplotype database. Overall, our study will provide insights and guidelines for the conservation of genetic diversity in captive forest musk deer populations in the Shaanxi province.

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.

Origin and spread of Thoroughbred racehorses inferred from complete mitochondrial genome sequences: Phylogenomic and Bayesian coalescent perspectives

The Thoroughbred horse breed was developed primarily for racing, and has a significant contribution to the qualitative improvement of many other horse breeds. Despite the importance of Thoroughbred racehorses in historical, cultural, and economical viewpoints, there was no temporal and spatial dynamics of them using the mitogenome sequences. To explore this topic, the complete mitochondrial genome sequences of 14 Thoroughbreds and two Przewalski’s horses were determined. These sequences were analyzed together along with 151 previously published horse mitochondrial genomes from a range of breeds across the globe using a Bayesian coalescent approach as well as Bayesian inference and maximum likelihood methods. The racing horses were revealed to have multiple maternal origins and to be closely related to horses from one Asian, two Middle Eastern, and five European breeds. Thoroughbred horse breed was not directly related to the Przewalski’s horse which has been regarded as the closest taxon to the all domestic horses and the only true wild horse species left in the world. Our phylogenomic analyses also supported that there was no apparent correlation between geographic origin or breed and the evolution of global horses. The most recent common ancestor of the Thoroughbreds lived approximately 8,100–111,500 years ago, which was significantly younger than the most recent common ancestor of modern horses (0.7286 My). Bayesian skyline plot revealed that the population expansion of modern horses, including Thoroughbreds, occurred approximately 5,500–11,000 years ago, which coincide with the start of domestication. This is the first phylogenomic study on the Thoroughbred racehorse in association with its spatio-temporal dynamics. The database and genetic history information of Thoroughbred mitogenomes obtained from the present study provide useful information for future horse improvement projects, as well as for the study of horse genomics, conservation, and in association with its geographical distribution.

Association of low race performance with mtDNA haplogroup L3b of Australian thoroughbred horses

Mitochondrial DNA (mtDNA) encodes the genes for respiratory chain sub-units that determine the efficiency of oxidative phosphorylation in mitochondria. The aim of this study was to determine if there were any haplogroups and variants in mtDNA that could be associated with athletic performance of Thoroughbred horses. The whole mitochondrial genomes of 53 maternally unrelated Australian Thoroughbred horses were sequenced and an association study was performed with the competition histories of 1123 horses within their maternal lineages. A horse mtDNA phylogenetic tree was constructed based on a total of 195 sequences (including 142 from previous reports). The association analysis showed that the sample groups with poor racing performance history were enriched in haplogroup L3b (p = .0003) and its sub-haplogroup L3b1a (p = .0007), while those that had elite performance appeared to be not significantly associated with haplogroups G2 and L3a1a1a (p > .05). Haplogroup L3b and L3b1a bear two and five specific variants of which variant T1458C (site 345 in 16s rRNA) is the only potential functional variant. Furthermore, secondary reconstruction of 16s RNA showed considerable differences between two types of 16s RNA molecules (with and without T1458C), indicating a potential functional effect. The results suggested that haplogroup L3b, could have a negative association with elite performance. The T1458C mutation harboured in haplogroup L3b could have a functional effect that is related to poor athletic performance.

A pedigree-based analysis of mitochondrial DNA diversity in a dog population on the example of German Hovawarts

Abstract. The purpose of the article is to illustrate the use of pedigree analysis to evaluate mtDNA diversity in a selected population of pedigree dogs, to describe the paths of mtDNA inheritance and to estimate the spread of potential pedigree errors or mutations that occurred in different generations of ancestors. Hovawart, old German breed, was used as an example. The number and frequencies of mtDNA haplotypes were calculated based on numbers of dam lines and their representatives. The scale of potential errors in calculations that can result from pedigree errors or from new mutations in ancestors from the 5th or 10th ancestral generation was evaluated. The analysis included 368 breeding bitches from four German kennel organizations. The bitches represented three dam lines, with the Ho1, Ho2 and HoU mtDNA haplotypes. Significant differences in the frequency of the haplotypes in the population, from 0.27 to 73.37 %, and among kennel organizations and regions of the country were recorded. Considerable differences in the scale of potential errors in calculations arising from mtDNA mutations or pedigree errors were noted between 0.27 and 28.69 %, depending on the number of representatives of the subline in which the error appeared and the generation taken into account in the simulations. The study revealed an interesting paradox: although the differences between the haplotypes are the result of events (mutations) from thousands of years ago, the number and the frequencies of the haplotypes in the population are the result of the modern history of the population and current breeding policy.