A genome-wide association study for body weight in Japanese Thoroughbred racehorses clarifies candidate regions on chromosomes 3, 9, 15, and 18

Unveiling Genetic Potential for Equine Meat Production: A Bioinformatics Approach

In view of the predicted significant increase in global meat production, alternative sources such as horsemeat are becoming increasingly important due to their lower environmental impact and high nutritional value. This study aimed to identify SNP markers on the GeneSeek® Genomic Profiler™ Equine (Neogen, Lansing, MI, USA) that are important for horsemeat production traits. First, orthologous genes related to meat yield in cattle and common genes between horses and cattle within QTLs for body size and weight were identified. Markers for these genes were then evaluated based on predicted variant consequences, GERP scores, and positions within constrained elements and orthologous regulatory regions in pigs. A total of 268 markers in 57 genes related to meat production were analyzed. This resulted in 27 prioritized SNP markers in 22 genes, including notable markers in LCORL, LASP1, IGF1R, and MSTN. These results will benefit smallholder farmers by providing genetic insights for selective breeding that could improve meat yield. This study also supports future large-scale genetic analyses such as GWAS and Genomic Best Linear Unbiased Prediction (GBLUP). The results of this study may be helpful in improving the accuracy of genomic breeding values. However, limitations include reliance on bioinformatics without experimental validation. Future research can validate these markers and consider a wider range of traits to ensure accuracy in equine breeding.

A genome-wide association study on growth traits of Korean commercial pig breeds using Bayesian methods

OBJECTIVE

This study aims to identify the significant regions and candidate genes of growth-related traits (adjusted backfat thickness [ABF], average daily gain [ADG], and days to 90 kg [DAYS90]) in Korean commercial GGP pig (Duroc, Landrace, and Yorkshire) populations.

METHODS

A genome-wide association study (GWAS) was performed using single-nucleotide polymorphism (SNP) markers for imputation to Illumina PorcineSNP60. The BayesB method was applied to calculate thresholds for the significance of SNP markers. The identified windows were considered significant if they explained ≥1% genetic variance.

RESULTS

A total of 28 window regions were related to genetic growth effects. Bayesian GWAS revealed 28 significant genetic regions including 52 informative SNPs associated with growth traits (ABF, ADG, DAYS90) in Duroc, Landrace, and Yorkshire pigs, with genetic variance ranging from 1.00% to 5.46%. Additionally, 14 candidate genes with previous functional validation were identified for these traits.

CONCLUSION

The identified SNPs within these regions hold potential value for future markerassisted or genomic selection in pig breeding programs. Consequently, they contribute to an improved understanding of genetic architecture and our ability to genetically enhance pigs. SNPs within the identified regions could prove valuable for future marker-assisted or genomic selection in pig breeding programs.

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.

Genetic structure and genome-wide association study of the traditional Kazakh horses

Horses are traditionally used in Kazakhstan as a source of food and as working and saddle animals as well. Here, for the first time, microarray-based medium-density single nucleotide polymorphism (SNP) genotyping of six traditionally defined types and breeds of indigenous Kazakh horses was conducted to reveal their genetic structure and find markers associated with animal size and weight. The results showed that the predefined separation between breeds and sampled populations was not supported by the molecular data. The lack of genetic variation between breeds and populations was revealed by the principal component analysis, ADMIXTURE, and distance-based analyses, as well as the general population parameters expected and observed heterozygosity (He and Ho) and between-group fixation index (Fst). The analysis revealed that the studied types and breeds should be considered as a single breed, namely the 'Kazakh horse'. The comparison with previously published data on global horse breed diversity revealed the relatively high level of individual diversity of Kazakh horses in comparison with the well-known foreign breeds. The Mongolian and Tuva breeds were identified as the closest horse landraces, demonstrating similar patterns of internal variability. The genome-wide association analysis was performed for animal size and weight as the traits directly related with the meat productivity of horses. The analysis identified a set of 60 SNPs linked with horse genes involved in the regulation of processes of development of connective tissues and the bone system, neural system, immune system regulation, and other processes. The present study is novel and introduces Kazakh horses as a promising genetic source for horse breeding and selection both on the domestic and international levels.

Common protein-coding variants influence the racing phenotype in galloping racehorse breeds

Selection for system-wide morphological, physiological, and metabolic adaptations has led to extreme athletic phenotypes among geographically diverse horse breeds. Here, we identify genes contributing to exercise adaptation in racehorses by applying genomics approaches for racing performance, an end-point athletic phenotype. Using an integrative genomics strategy to first combine population genomics results with skeletal muscle exercise and training transcriptomic data, followed by whole-genome resequencing of Asian horses, we identify protein-coding variants in genes of interest in galloping racehorse breeds (Arabian, Mongolian and Thoroughbred). A core set of genes, G6PC2, HDAC9, KTN1, MYLK2, NTM, SLC16A1 and SYNDIG1, with central roles in muscle, metabolism, and neurobiology, are key drivers of the racing phenotype. Although racing potential is a multifactorial trait, the genomic architecture shaping the common athletic phenotype in horse populations bred for racing provides evidence for the influence of protein-coding variants in fundamental exercise-relevant genes. Variation in these genes may therefore be exploited for genetic improvement of horse populations towards specific types of racing.

Population structure and genomic footprints of selection in five major Iranian horse breeds

The genetic structure and characteristics of Iranian native breeds are yet to be comprehensibly investigated and studied. Therefore, we employed genomic information of 364 Iranian native horses representing the Asil (n = 109), Caspian (n = 40), Dareshuri (n = 44), Kurdish (n = 95), and Turkoman (n = 76) breeds to reveal the genetic structure and characteristics. For these and 19 other horse breeds, principal component analysis, Bayesian model-based, Neighbor-Net, and bootstrap-based TreeMix approaches were applied to investigate and compare their genetic structure. Additionally, three haplotype-based methods including haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length were applied to trace genomic footprints of selection of Asil, Caspian, Dareshuri, Kurdish, and Turkoman groups. Then, the Mahalanobis distance based on the negative-log₁₀ rank-based P-values was estimated based on the haplotype homozygosity pooled, integrated haplotype score, and number of segregating sites by length values. Asil, Caspian, Dareshuri, Kurdish, and Turkoman can be categorized into five different genetic clusters. Based on the top 1% of Mahalanobis distance based on the negative-log₁₀ rank-based P-values of SNPs, we identified 24 SNPs formerly reported to be associated with different traits and >100 genes undergoing selection pressures in Asil, Caspian, Dareshuri, Kurdish, and Turkoman. The detected QTL undergoing selection pressures were associated with withers height, equine metabolic syndrome, overall body size, insect bite hypersensitivity, guttural pouch tympany, white markings, Rhodococcus equi infection, jumping test score, alternate gaits, and body weight traits. Our findings will aid to have a better perspective of the genetic characteristics and population structure of Asil, Caspian, Dareshuri, Kurdish, and Turkoman horses as Iranian native horse breeds.

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.

Genetic characterization of Kushum horses in Kazakhstan based on haplotypes of mtDNA and Y chromosome, and genes associated with important traits of the horses

The Kushum is a relatively new breed of horses in Kazakhstan that was established in the middle of the 20th century through a cross between mares of Kazakhstan local horses and stallions of Thoroughbred, Trotter, and Russian Don breeds to supply military horses. To reveal the genetic characteristics of this breed, we investigated haplotypes of mitochondrial DNA (mtDNA) and single-nucleotide polymorphisms of the Y chromosome, as well as genotypes of five functional genes associated with coat color, body composition, and locomotion traits. We detected 10 mtDNA haplotypes that fell into 8 of the 17 major haplogroups of horse mtDNA, indicating a unique haplotype composition with high genetic diversity. We also found two Y-chromosomal haplotypes in Kushum horses, which likely originated from Trotter and/or Don breeds. The findings regarding the mtDNA and Y-chromosomal haplotypes are concordant with the documented maternal and paternal origins of the Kushum horses. The allele frequencies of ASIP, MC1R, and MATP associated with coat color were consistent with the coat color variations of Kushum horses. The allele frequencies of MSTN associated with endurance performance and those of DMRT3 associated with gait suggested that the observed allele frequencies of these genes were the result of selective breeding for these traits. As a result of this study, we were able to obtain useful information for a better understanding of the origin and breeding history of the Kushum horse breed using molecular markers.

Satellite cells and their regulation in livestock

Satellite cells are the myogenic stem and progenitor population found in skeletal muscle. These cells typically reside in a quiescent state until called upon to support repair, regeneration, or muscle growth. The activities of satellite cells are orchestrated by systemic hormones, autocrine and paracrine growth factors, and the composition of the basal lamina of the muscle fiber. Several key intracellular signaling events are initiated in response to changes in the local environment causing exit from quiescence, proliferation, and differentiation. Signals emanating from Notch, wingless-type mouse mammary tumor virus integration site family members, and transforming growth factor-β proteins mediate the reversible exit from growth 0 phase while those initiated by members of the fibroblast growth factor and insulin-like growth factor families direct proliferation and differentiation. Many of these pathways impinge upon the myogenic regulatory factors (MRF), myogenic factor 5, myogenic differentiation factor D, myogenin and MRF4, and the lineage determinate, Paired box 7, to alter transcription and subsequent satellite cell decisions. In the recent past, insight into mouse transgenic models has led to a firm understanding of regulatory events that control satellite cell metabolism and myogenesis. Many of these niche-regulated functions offer subtle differences from their counterparts in livestock pointing to the existence of species-specific controls. The purpose of this review is to examine the mechanisms that mediate large animal satellite cell activity and their relationship to those present in rodents.

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.

Microfluidic Quantitative PCR Detection of 12 Transgenes from Horse Plasma for Gene Doping Control

Gene doping, an activity which abuses and misuses gene therapy, is a major concern in sports and horseracing industries. Effective methods capable of detecting and monitoring gene doping are urgently needed. Although several PCR-based methods that detect transgenes have been developed, many of them focus only on a single transgene. However, numerous genes associated with athletic ability may be potential gene-doping material. Here, we developed a detection method that targets multiple transgenes. We targeted 12 genes that may be associated with athletic performance and designed two TaqMan probe/primer sets for each one. A panel of 24 assays was prepared and detected via a microfluidic quantitative PCR (MFQPCR) system using integrated fluidic circuits (IFCs). The limit of detection of the panel was 6.25 copy/μL. Amplification-specificity was validated using several concentrations of reference materials and animal genomic DNA, leading to specific detection. In addition, target-specific detection was successfully achieved in a horse administered 20 mg of the EPO transgene via MFQPCR. Therefore, MFQPCR may be considered a suitable method for multiple-target detection in gene-doping control. To our knowledge, this is the first application of microfluidic qPCR (MFQPCR) for gene-doping control in horseracing.

A candidate-SNP retrospective cohort study for fracture risk in Japanese Thoroughbred racehorses

Fractures are medical conditions that compromise the athletic potential of horses and/or the safety of jockeys. Therefore, the reduction of fracture risk is an important horse and human welfare issue. The present study used molecular genetic approaches to determine the effect of genetic risk for fracture at four candidate SNPs spanning the myostatin (MSTN) gene on horse chromosome 18. Among the 3706 Japanese Thoroughbred racehorses, 1089 (29.4%) had experienced fractures in their athletic life, indicating the common occurrence of this injury in Thoroughbreds. In the case/control association study, fractures of the carpus (carpal bones and distal radius) were statistically associated with g.65809482T/C (P = 1.17 x 10^-8 ), g.65868604G/T (P = 2.66 x 10^-9 ), and g.66493737C/T (P = 6.41 x 10^-8 ). In the retrospective cohort study using 1710 racehorses born in 2000, the relative risk (RR) was highest for male horses at g.65868604G/T, based on the dominant allele risk model (RR = 2.251, 95% confidence interval 1.407-3.604, P = 0.00041), and for female horses at g.65868604G/T, based on the recessive allele risk model (RR = 2.313, 95% confidence interval 1.380-3.877, P = 0.00163). Considering the association of these SNPs with racing performance traits such as speed, these genotypes may affect the occurrence of carpus fractures in Japanese Thoroughbred racehorses as a consequence of the non-genetic influence of the genotype on the distance and/or intensity of racing and training. The genetic information presented here may contribute to the development of strategic training programs and racing plans for racehorses that improve their health and welfare.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

Genetic variants influencing obesity-related traits in Japanese population

Context: Adipose tissue is the main organ that stores energy and participates in adaptive thermogenesis of the human body. The adipose tissue content in an individual is determined by a combination of genetic factors and lifestyle related factors. While Japanese people, along with the closely related East Asians, are generally thinner than individuals of European ancestry, they are prone to accumulating visceral adipose tissues. Genome-wide discovery of loci influencing obesity-related traits, and application of the genome sequence data to assess natural selection, provides evidence that the obesity-related traits in East Asians might be shaped by natural selection. Objective: This review aims to summarise health and evolutionary implications of genetic variants influencing obesity-related traits in Japanese. Methods: This study gathered recently published papers of medical, genetic and evolutionary studies regarding obesity-related traits in the Japanese and closely related East Asians. Results and conclusion: A high susceptibility to central obesity of Japanese and closely related East Asians might have been shaped by natural selection favouring thrifty genotypes. Moreover, natural selection favouring higher thermogenic activity of brown adipose tissues would contribute to increased non-thrifty alleles in ancestors of East Asians.

Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse

The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.

Identification of loci affecting sexually dimorphic patterns for height and recurrent laryngeal neuropathy risk in American Belgian Draft Horses

Equine recurrent laryngeal neuropathy (RLN) is a bilateral mononeuropathy with an unknown etiology. In Thoroughbreds (TB), we previously demonstrated that the haplotype association for height (LCORL/NCAPG locus on ECA3, which affects body size) and RLN was coincident. In the present study, we performed a genome-wide association scan (GWAS) for RLN in 458 American Belgian Draft Horses, a breed fixed for the LCORL/NCAPG risk alelle. In this breed, RLN risk is associated with sexually dimorphic differences in height, and we identified a novel locus contributing to height in a sex-specific manner: MYPN (ECA1). Yet this specific locus contributes little to RLN risk, suggesting that other growth traits correlated to height may underlie the correlation to this disease. Controlling for height, we identified a locus on ECA15 contributing to RLN risk specifically in males. These results suggest that loci with sex-specific gene expression play an important role in altering growth traits impacting RLN etiology, but not necessarily adult height. These newly identified genes are promising targets for novel preventative and treatment strategies.

Digital PCR detection of plasmid DNA administered to the skeletal muscle of a microminipig: a model case study for gene doping detection

Objective

Doping control is an important and indispensable aspect of fair horse racing; genetic doping has been recently included to this. In this study, we aimed to develop a detection method of gene doping. A plasmid cloned with human erythropoietin gene (p.hEPO, 250 μg/head) was intramuscularly injected into a microminipig. Subsequently, p.hEPO was extracted from 1 mL of plasma and detected by droplet digital polymerase chain reaction.

Results

The results confirmed that the maximum amount of plasmid was detected at 15 min after administration and the majority of the plasmid was degraded in the bloodstream within 1–2 days after administration. In contrast, low amounts of p.hEPO were detected at 2–3 weeks after administration. These results suggest that the proposed method to detect gene doping can help obtain information for experiments using horses.