A genome-wide association study reveals loci influencing height and other conformation traits in horses

Genetic composition of Kazakh horses of Zhabe type evaluated by SNP genotyping

Horses are animals traditionally playing prominent role as both food source and working animals for Kazakh people. Zhabe horses are traditional type of indigenous Kazakh horses characterized by versatility and adaptation to conditions of Central Asia. The present work focuses on examination of genetic structure of Zhabe horses using SNP genotyping with addition of previously published data. Total 1038 individuals including 403 new samples of Zhabe horses and 42 sample of white horses 'Zhetysu Asyly' have been considered. DNA was extracted from hair roots using commercial DNA isolation kit and further used for analysis of SNP by Illumina iScan system with Equine80k SNP array. The analysis of population genetic parameters (expected and observed heterozygosity, linkage disequilibrium, Wright's F st ) and genetic structure (PCA, ADMIXTURE) in comparison with publicly available data on selected foreign cultivars demonstrated low between population differentiations and lack of selection factors. Genome wide association study performed for body size and weight have revealed low occurrence of SNPs with significant associations, total 57 SNPs linked to various genes with low density across all genome. The obtained results highlight difference between traditional horse breeding practices of Kazakh people and stable based breeding of foreign breeds. In contrast, the 'Zhetysu Asyly' horse breed derived from Kazakh horses demonstrate the effect of intense breeding process on the same landrace. The results provide new data on the traditional Kazakh horses of type Zhabe and will assist further studies of this original landrace.

A genome-wide association study for recurrent laryngeal neuropathy in the Thoroughbred horse identifies a candidate gene that regulates myelin structure

BACKGROUND

Equine recurrent laryngeal neuropathy (RLN) is an economically important upper respiratory tract (URT) disease with a genetic contribution to risk, but genetic variants independent of height have not been identified for Thoroughbreds. The method of clinical assessment for RLN is critical to accurately phenotype groups for genetic studies.

OBJECTIVES

To identify genetic risk loci for RLN in Thoroughbreds in a genome-wide association study (GWAS) following high-resolution phenotyping.

STUDY DESIGN

Case-control.

METHODS

Thoroughbred horses were characterised as RLN cases and controls using resting and exercising URT endoscopic examinations and laryngeal ultrasonography, with the case-cohort supplemented using a questionnaire. Genotypes for 43 831 autosomal single-nucleotide polymorphisms (SNPs) from n = 235 horses (n = 110 cases; n = 125 controls) were used to estimate trait heritability and identify significantly associated SNPs in a GWAS. Haplotypes were examined in cases and controls and risk allele frequencies were examined in a population cohort (n = 3126).

RESULTS

Heritability was h2 = 0.30 including sex and 5PCs as covariates. A SNP on ECA20 located between candidate genes, DAAM2 and LRFN2, was significantly associated with RLN. Six index SNPs with allelic effect sizes OR = 1.5-2.9 were identified on ECA1, ECA14, and ECA20 close to candidate genes ATPA10, KCNN2, and TFAP2A. Eleven ECA20 SNPs defined seven haplotypes with homozygous H2/H2 horses having a 3.1× higher risk of RLN. Risk alleles segregate in the population, and stallions are carriers.

MAIN LIMITATIONS

The main study population was young. Horses in the control group had no evidence of RLN as 2- or 3-year olds but may have developed RLN later.

CONCLUSIONS

Genetic markers for RLN were identified which may be useful for the development of a polygenic risk score. Candidate genes with functions in neuropathies may further the understanding of RLN pathobiology.

Single-Step Genome-Wide Association Study of Factors for Evaluated and Linearly Scored Traits in Swedish Warmblood Horses

Swedish Warmblood horses (SWB) are bred for show jumping and/or dressage with young horse test scores as indicator traits. This study aimed to investigate possible candidate genes and regions of importance for evaluated and linearly scored young horse test traits. A single-step genome-wide association study (ssGWAS) was done using the BLUPF90 suite of programs for factors scores from factor analysis of traits assessed at young horse tests together with height at withers. The ssGWAS included 20,814 SWB with factors scores for four factors for evaluated traits. A total of 6436 of these horses also had factor scores for 13 factors for linearly scored traits. Genotypes from a 670K SNP array were available for 380 of the horses in this study. All genotyped horses had factor scores for evaluated traits, and 379 also had factors scores for linearly scored traits. Significant SNPs associated with three factors related to size were located on ECA3 within or nearby a well-known region, including the genes ligand dependent nuclear receptor corepressor like (LCORL), non-SMC condensin I complex subunit G (NCAPG), DDB1 and CUL4 Associated Factor 16 (DCAF16), and the Family with Sequence Similarity 184 Member B (FAM184B). Significant SNPs were also detected for two factors for evaluated traits representing conformation and jumping, and four factors for linearly scored traits related to body length, neck conformation, walk and trot (hindleg position and activity), respectively. Among nearby genes, calcium/calmodulin-dependent protein kinase type 1D (CAMK1D) for the factor for linearly scored traits related to neck conformation and GLI Family Zinc Finger 2 (GLI2) for the factor for evaluated jumping traits, were most promising. For these, top associated SNPs were detected within the genes, and the known gene functions seems to be related to the phenotypes. In conclusion, ssGWAS is beneficial to detect plausible candidate genes/regions for desired traits in warmblood horses.

Comparative Whole-Genome Analysis of Production Traits and Genetic Structure in Baiyu and Chuanzhong Black Goats

Natural selection and artificial breeding are crucial methods for developing new animal groups. The Baiyu black goats and Chuanzhong black goats are indigenous goat breeds from distinct ecological regions in Sichuan Province, with dramatically different growth and reproductivity. This study aimed to systematically elucidate the differences in production performance and genetic traits between Baiyu black goats and Chuanzhong black goats. We quantified growth and reproductive attributes for both breeds. Furthermore, we conducted a comprehensive analysis of genetic diversity, population structure, and selection signatures using whole-genome resequencing data. This dataset included 30 individuals from the Baiyu black goat breed, 41 from the Chuanzhong black goat breed, and an additional 59 individuals representing Chengdu grey goats, Tibetan cashmere goats, and Jianchang black goats, totaling 130 individuals across five goat breeds. The comparative analysis of production performance revealed that the weight and body size of Chuanzhong black goats were significantly higher than those of Baiyu black goats (p < 0.01). At the same time, the average kidding rate and kid-weaning survival rate of Chuanzhong black goats were also notably superior to those of Baiyu black goats (p < 0.01). The Baiyu black goats exhibited a more abundant genetic diversity and distinct genetic differences compared to the Chuanzhong black goat, according to an analysis grounded on genomic variation. The Baiyu black goats are more closely related to Tibetan cashmere goats, whereas Chuanzhong black goats share a closer genetic relationship with Chengdu grey goats. Additionally, we employed the π, Fst, and XP-EHH methodologies to identify genes related to immunity (TRIM10, TRIM15, TRIM26, and TRIM5), neurodevelopment (FOXD4L1, PCDHB14, PCDHB4, PCDHB5, PCDHB6, and PCDHB7), reproduction (BTNL2 and GABBR1), body size (NCAPG, IBSP, and MKNK1), and meat quality traits (SUCLG2 and PGM5). These results provide a theoretical basis for further resource conservation and breeding improvement of the Baiyu black goat and Chuanzhong black goat.

Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat

Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.

Transcriptional Signatures of Domestication Revealed through Meta-Analysis of Pig, Chicken, Wild Boar, and Red Junglefowl Gene Expression Data

Domesticated animals have undergone significant changes in their behavior, morphology, and physiological functions during domestication. To identify the changes in gene expression associated with domestication, we collected the RNA-seq data of pigs, chickens, wild boars, and red junglefowl from public databases and performed a meta-analysis. Gene expression was quantified, and the expression ratio between domesticated animals and their wild ancestors (DW-ratio) was calculated. Genes were classified as "upregulated", "downregulated", or "unchanged" based on their DW-ratio, and the DW-score was calculated for each gene. Gene set enrichment analysis revealed that genes upregulated in pigs were related to defense from viral infection, whereas those upregulated in chickens were associated with aminoglycan and carbohydrate derivative catabolic processes. Genes commonly upregulated in pigs and chickens are involved in the immune response, olfactory learning, epigenetic regulation, cell division, and extracellular matrix. In contrast, genes upregulated in wild boar and red junglefowl are related to stress response, cell proliferation, cardiovascular function, neural regulation, and energy metabolism. These findings provide valuable insights into the genetic basis of the domestication process and highlight potential candidate genes for breeding applications.

Genomic tools for early selection among Thoroughbreds and Polo Argentino horses for practicing polo

The Polo Argentino (PA) horse is a recognized breed, developed originally by mixing crossbred and Thoroughbred (TB) horses to play polo. Early PA selection is difficult due to unreliable performance estimations. This study investigated the usefulness of genomic markers previously linked to morphological and functional traits as a tool for the early selection of PA. To this, we genotyped 520 PA and 30 TB horses using the Equine GGPArray (Illumina, n = 71,778 SNPs). Analyses included a genetic characterization of six genetic markers associated with behavioral (DRD4), muscular development (MSTN), and body size (LCORL, HMGA6, ZFAT, and LASP1) genes. Genetic differences in the DRD4, MSTN, and LCORL SNP were found between the two breeds, in the last two FST index between breeds was 0.13 and 0.6, respectively (p < 0.01). In DRD4, G allele was the more prevalent in PA (0.56 vs 0.45 in TB, p < 0.05), but no differences were observed between the genotypes associated with phenotypes. In MSTN, heterozygous genotypes were the most common in PA (48 %), with a significant decrease in AA (Hardy-Weinberg p < 0.05), suggesting a negative selection against it in polo horses. In body size, HMGA2 was monomorphic in all horses, while ZFAT and LASP1 SNP showed higher variability. Interestingly, 99 % of PA showed a TT genotype in LCORL (only 66 % in TB), demonstrating selection for smaller horses. Our results suggest that empirical selection in PA has generated an incipient genomic differentiation in discrete traits which could be used as a marker-assisted selection tool for early selection of polo horses.

Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection

The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.

Genomic analyses of withers height and linear conformation traits in German Warmblood horses using imputed sequence-level genotypes

BACKGROUND

Body conformation, including withers height, is a major selection criterion in horse breeding and is associated with other important traits, such as health and performance. However, little is known about the genomic background of equine conformation. Therefore, the aim of this study was to use imputed sequence-level genotypes from up to 4891 German Warmblood horses to identify genomic regions associated with withers height and linear conformation traits. Furthermore, the traits were genetically characterised and putative causal variants for withers height were detected.

RESULTS

A genome-wide association study (GWAS) for withers height confirmed the presence of a previously known quantitative trait locus (QTL) on Equus caballus (ECA) chromosome 3 close to the LCORL/NCAPG locus, which explained 16% of the phenotypic variance for withers height. An additional significant association signal was detected on ECA1. Further investigations of the region on ECA3 identified a few promising candidate causal variants for withers height, including a nonsense mutation in the coding sequence of the LCORL gene. The estimated heritability for withers height was 0.53 and ranged from 0 to 0.34 for the conformation traits. GWAS identified significantly associated variants for more than half of the investigated conformation traits, among which 13 showed a peak on ECA3 in the same region as withers height. Genetic parameter estimation revealed high genetic correlations between these traits and withers height for the QTL on ECA3.

CONCLUSIONS

The use of imputed sequence-level genotypes from a large study cohort led to the discovery of novel QTL associated with conformation traits in German Warmblood horses. The results indicate the high relevance of the QTL on ECA3 for various conformation traits, including withers height, and contribute to deciphering causal mutations for body size in horses.

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

Defining a Haplotype Encompassing the LCORL-NCAPG Locus Associated with Increased Lean Growth in Beef Cattle

Numerous studies have shown genetic variation at the LCORL-NCAPG locus is strongly associated with growth traits in beef cattle. However, a causative molecular variant has yet to be identified. To define all possible candidate variants, 34 Charolais-sired calves were whole-genome sequenced, including 17 homozygous for a long-range haplotype associated with increased growth (QQ) and 17 homozygous for potential ancestral haplotypes for this region (qq). The Q haplotype was refined to an 814 kb region between chr6:37,199,897-38,014,080 and contained 218 variants not found in qq individuals. These variants include an insertion in an intron of NCAPG, a previously documented mutation in NCAPG (rs109570900), two coding sequence mutations in LCORL (rs109696064 and rs384548488), and 15 variants located within ATAC peaks that were predicted to affect transcription factor binding. Notably, rs384548488 is a frameshift variant likely resulting in loss of function for long isoforms of LCORL. To test the association of the coding sequence variants of LCORL with phenotype, 405 cattle from five populations were genotyped. The two variants were in complete linkage disequilibrium. Statistical analysis of the three populations that contained QQ animals revealed significant (p < 0.05) associations with genotype and birth weight, live weight, carcass weight, hip height, and average daily gain. These findings affirm the link between this locus and growth in beef cattle and describe DNA variants that define the haplotype. However, further studies will be required to define the true causative mutation.

Unlocking the Transcriptional Control of NCAPG in Bovine Myoblasts: CREB1 and MYOD1 as Key Players

Muscle formation directly determines meat production and quality. The non-SMC condensin I complex subunit G (NCAPG) is strongly linked to the growth features of domestic animals because it is essential in controlling muscle growth and development. This study aims to elucidate the tissue expression level of the bovine NCAPG gene, and determine the key transcription factors for regulating the bovine NCAPG gene. In this study, we observed that the bovine NCAPG gene exhibited high expression levels in longissimus dorsi and spleen tissues. Subsequently, we cloned and characterized the promoter region of the bovine NCAPG gene, consisting of a 2039 bp sequence, through constructing the deletion fragment double-luciferase reporter vector and site-directed mutation-identifying core promoter region with its key transcription factor binding site. In addition, the key transcription factors of the core promoter sequence of the bovine NCAPG gene were analyzed and predicted using online software. Furthermore, by integrating overexpression experiments and the electrophoretic mobility shift assay (EMSA), we have shown that cAMP response element binding protein 1 (CREB1) and myogenic differentiation 1 (MYOD1) bind to the core promoter region (-598/+87), activating transcription activity in the bovine NCAPG gene. In conclusion, these findings shed important light on the regulatory network mechanism that underlies the expression of the NCAPG gene throughout the development of the muscles in beef cattle.

Determination of breeding criteria for gait proficiency in leisure riding and racing dromedary camels: a stepwise multivariate analysis of factors predicting overall biomechanical performance

To date, the biomechanical dynamics in camelids have not been addressed, although it might be a factor that can affect selection and breeding in this species. Therefore, the aim of this article is to conduct curve fitting and discriminant canonical analysis to identify the mathematical function that best captures the dynamics of camel locomotion and to study the impact of kinematic, morphometric, physiological, and phaneroptic variables on gait performance in leisure riding and racing activities in dromedaries, respectively. The cubic function emerged as the most suitable mathematical model to represent the locomotive behavior of camels. Various factors were found to play a pivotal role in the athletic performance of leisure riding and racing dromedary camels. Concretely, angular measurements at the distal fore and rear extremity areas, pelvis inclination, relative volume of the hump, impact forces of the front limbs, post-neutering effects, and the kinematic behavior of the scapula, shoulder, carpus, hip, and foot are the factors that greatly impact gait performance in leisure riding and racing camels. The biomechanical performance at these specific body regions has a profound impact on weight absorption and minimization of mechanic impact during camel locomotion, static/dynamic balance, force distribution, energy of propulsion, movement direction and amplitude, and storage of elastic strain in leisure riding and racing dromedaries. In contrast, other animal- and environment-dependent factors do not exert significant influence on camel gait performance, which can be attributed to species-specific, inherited adaptations developed in response to desert conditions, including the pacing gait, broad foot pads, and energy-efficient movements. The outcomes of our functional data analysis can provide valuable insights for making informed breeding decisions aimed at enhancing animal functional performance in camel riding and racing activities. Furthermore, these findings can open avenues for exploring alternative applications, such as camel-assisted therapy.

Sequenced-based GWAS for linear classification traits in Belgian Blue beef cattle reveals new coding variants in genes regulating body size in mammals

BACKGROUND

Cohorts of individuals that have been genotyped and phenotyped for genomic selection programs offer the opportunity to better understand genetic variation associated with complex traits. Here, we performed an association study for traits related to body size and muscular development in intensively selected beef cattle. We leveraged multiple trait information to refine and interpret the significant associations.

RESULTS

After a multiple-step genotype imputation to the sequence-level for 14,762 Belgian Blue beef (BBB) cows, we performed a genome-wide association study (GWAS) for 11 traits related to muscular development and body size. The 37 identified genome-wide significant quantitative trait loci (QTL) could be condensed in 11 unique QTL regions based on their position. Evidence for pleiotropic effects was found in most of these regions (e.g., correlated association signals, overlap between credible sets (CS) of candidate variants). Thus, we applied a multiple-trait approach to combine information from different traits to refine the CS. In several QTL regions, we identified strong candidate genes known to be related to growth and height in other species such as LCORL-NCAPG or CCND2. For some of these genes, relevant candidate variants were identified in the CS, including three new missense variants in EZH2, PAPPA2 and ADAM12, possibly two additional coding variants in LCORL, and candidate regulatory variants linked to CCND2 and ARMC12. Strikingly, four other QTL regions associated with dimension or muscular development traits were related to five (recessive) deleterious coding variants previously identified.

CONCLUSIONS

Our study further supports that a set of common genes controls body size across mammalian species. In particular, we added new genes to the list of those associated with height in both humans and cattle. We also identified new strong candidate causal variants in some of these genes, strengthening the evidence of their causality. Several breed-specific recessive deleterious variants were identified in our QTL regions, probably as a result of the extreme selection for muscular development in BBB cattle.

Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds

BACKGROUND

The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type.

RESULTS

We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis.

CONCLUSIONS

Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.

Population genomic structures and signatures of selection define the genetic uniqueness of several fancy and meat rabbit breeds

Following the recent domestication process of the European rabbit (Oryctolagus cuniculus), many different breeds and lines, distinguished primarily by exterior traits such as coat colour, fur structure and body size and shape, have been constituted. In this study, we genotyped, with a high-density single-nucleotide polymorphism panel, a total of 645 rabbits from 10 fancy breeds (Belgian Hare, Champagne d'Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex and Rhinelander) and three meat breeds (Italian White, Italian Spotted and Italian Silver). ADMIXTURE analysis indicated that breeds with similar phenotypic traits (e.g. coat colour and body size) shared common ancestries. Signatures of selection using two haplotype-based approaches (iHS and XP-EHH), combined with the results obtained with other methods previously reported that we applied to the same breeds, we identified a total of 5079 independent genomic regions with some signatures of selection, covering about 1777 Mb of the rabbit genome. These regions consistently encompassed many genes involved in pigmentation processes (ASIP, EDNRA, EDNRB, KIT, KITLG, MITF, OCA2, TYR and TYRP1), coat structure (LIPH) and body size, including two major genes (LCORL and HMGA2) among many others. This study revealed novel genomic regions under signatures of selection and further demonstrated that population structures and signatures of selection, left into the genome of these rabbit breeds, may contribute to understanding the genetic events that led to their constitution and the complex genetic mechanisms determining the broad phenotypic variability present in these untapped rabbit genetic resources.

Comparative population genomics analysis uncovers genomic footprints and genes influencing body weight trait in Chinese indigenous chicken

Body weight of chicken is a typical quantitative trait, which shows phenotypic variations due to selective breeding. Despite some QTL loci have been obtained, the body weight of native chicken breeds in different geographic regions varies greatly, its genetic basis remains unresolved questions. To address this issue, we analyzed 117 Chinese indigenous chickens from 10 breeds (Huiyang Bearded, Xinhua, Hotan Black, Baicheng You, Liyang, Yunyang Da, Jining Bairi, Lindian, Beijing You, Tibetan). We applied fixation index (FST) analysis to find selected genomic regions and genes associated with body weight traits. Our study suggests that NELL1, XYLT1, and NCAPG/LCORL genes are strongly selected in the body weight trait of Chinese indigenous chicken breeds. In addition, the IL1RAPL1 gene was strongly selected in large body weight chickens, while the PCDH17 and CADM2 genes were strongly selected in small body weight chickens. This result suggests that the patterns of genetic variation of native chicken and commercial chicken, and/or distinct local chicken breeds may follow different evolutionary mechanisms.

Patterns of enrichment and acceleration in evolutionary rates of promoters suggest a role of regulatory regions in cetacean gigantism

BACKGROUND

Cetaceans (whales, porpoises, and dolphins) are a lineage of aquatic mammals from which some species became giants. Only recently, gigantism has been investigated from the molecular point of view. Studies focused mainly on coding regions, and no data on the influence of regulatory regions on gigantism in this group was available. Accordingly, we investigated the molecular evolution of non-coding regulatory regions of genes already described in the literature for association with size in mammals, focusing mainly on the promoter regions. For this, we used Ciiider and phyloP tools. Ciiider identifies significantly enriched transcription factor binding sites, and phyloP estimates the molecular evolution rate of the promoter.

RESULTS

We found evidence of enrichment of transcription binding factors related to large body size, with distinct patterns between giant and non-giant cetaceans in the IGFBP7 and NCAPG promoters, in which repressive agents are present in small cetaceans and those that stimulate transcription, in giant cetaceans. In addition, we found evidence of acceleration in the IGF2, IGFBP2, IGFBP7, and ZFAT promoters.

CONCLUSION

Our results indicate that regulatory regions may also influence cetaceans' body size, providing candidate genes for future research to understand the molecular basis of the largest living animals.

Genome-Wide Detection of Copy Number Variations Associated with Miniature Features in Horses

Copy number variations (CNVs) are crucial structural genomic variants affecting complex traits in humans and livestock animals. The current study was designed to conduct a comprehensive comparative copy number variation analysis among three breeds, Debao (DB), Baise (BS), and Warmblood (WB), with a specific focus on identifying genomic regions associated with miniature features in horses. Using whole-genome next-generation resequencing data, we identified 18,974 CNVs across 31 autosomes. Among the breeds, we found 4279 breed-specific CNV regions (CNVRs). Baise, Debao, and Warmblood displayed 2978, 986, and 895 distinct CNVRs, respectively, with 202 CNVRs shared across all three breeds. After removing duplicates, we obtained 1545 CNVRs from 26 horse genomes. Functional annotation reveals enrichment in biological functions, including antigen processing, cell metabolism, olfactory conduction, and nervous system development. Debao horses have 970 genes overlapping with CNVRs, possibly causing their small size and mountainous adaptations. We also found that the genes GHR, SOX9, and SOX11 may be responsible for the miniature features of the Debao horse by analyzing their overlapping CNVRs. Overall, this study offers valuable insights into the widespread presence of CNVs in the horse genome. The findings contribute to mapping horse CNVs and advance research on unique miniature traits observed in the Debao horse.

Genetic analysis of geometric morphometric 3D visuals of French jumping horses

BACKGROUND

For centuries, morphology has been the most commonly selected trait in horses. A 3D video recording enabled us to obtain the coordinates of 43 anatomical landmarks of 2089 jumping horses. Generalized Procrustes analysis provided centered and scaled coordinates that were independent of volume, i.e., centroid size. Genetic analysis of these coordinates (mixed model; 17,994 horses in the pedigree) allowed us to estimate a variance-covariance matrix. New phenotypes were then defined: the "summarized shapes". They were obtained by linear combinations of Procrustes coordinates with, as coefficients, the eigenvectors of the genetic variance-covariance matrix. These new phenotypes were used in genome-wide association analyses (GWAS) and multitrait genetic analysis that included judges' scores and competition results of the horses.

RESULTS

We defined ten shapes that represented 86% of the variance, with heritabilities ranging from 0.14 to 0.42. Only one of the shapes was found to be genetically correlated with competition success (rg = - 0.12, standard error = 0.07). Positive and negative genetic correlations between judges' scores and shapes were found. This means that the breeding objective defined by judges involves improvement of anatomical parts of the body that are negatively correlated with each other. Known single nucleotide polymorphisms (SNPs) on chromosomes 1 and 3 for height at withers were significant for centroid size but not for any of the shapes. As these SNPs were not associated with the shape that distinguished rectangular horses from square horses (with height at withers greater than body length), we hypothesize that these SNPs play a role in the overall development of horses, i.e. in height, width, and length but not in height at withers when standardized to unit centroid size. Several other SNPs were found significant for other shapes.

CONCLUSIONS

The main application of 3D morphometric analysis is the ability to define the estimated breeding value (EBV) of a sire based on the shape of its potential progeny, which is easier for breeders to visualize in a single synthetic image than a full description based on linear profiling. However, the acceptance of these new phenotypes by breeders and the complex nature of summarized shapes may be challenging. Due to the low genetic correlations of the summarized shapes with jumping performance, the methodology did not allow indirect performance selection criteria to be defined.

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.

New genomic insights into the conformation of Lipizzan horses

Conformation traits are important selection criteria in equine breeding, as they describe the exterior aspects of the horse (height, joint angles, shape). However, the genetic architecture of conformation is not well understood, as data of these traits mainly consist of subjective evaluation scores. Here, we performed genome-wide association studies on two-dimensional shape data of Lipizzan horses. Based on this data, we identified significant quantitative trait loci (QTL) associated with cresty neck on equine chromosome (ECA)16 within the MAGI1 gene, and with type, hereby differentiating heavy from light horses on ECA5 within the POU2F1 gene. Both genes were previously described to affect growth, muscling and fatty deposits in sheep, cattle and pigs. Furthermore, we pin-pointed another suggestive QTL on ECA21, near the PTGER4 gene, associated with human ankylosing spondylitis, for shape differences in the back and pelvis (roach back vs sway back). Further differences in the shape of the back and abdomen were suggestively associated with the RYR1 gene, involved in core muscle weakness in humans. Therefore, we demonstrated that horse shape space data enhance the genomic investigations of horse conformation.

Integrating plant and fungal quantitative genetics to improve the ecological and agricultural applications of mycorrhizal symbioses

Finding and targeting genes that quantitatively contribute to agricultural and ecological processes progresses food production and conservation efforts. Typically, quantitative genetic approaches link variants in a single organism's genome with a trait of interest. Recently, genome-to-genome mapping has found genome variants interacting between species to produce the result of a multiorganism (including multikingdom) interaction. These were plant and bacterial pathogen genome interactions; plant-fungal coquantitative genetics have not yet been applied. Plant-mycorrhizae symbioses exist across most biomes, for a majority of land plants, including crop plants, and manipulate many traits from single organisms to ecosystems for which knowing the genetic basis would be useful. The availability of Rhizophagus irregularis mycorrhizal isolates, with genomic information, makes dual-genome methods with beneficial mutualists accessible and imminent.

Current Analytical Methods and Research Trends Are Used to Identify Domestic Pig and Wild Boar DNA in Meat and Meat Products

The pig, one of the most important livestock species, is a meaningful source of global meat production. It is necessary, however, to prove whether a food product that a discerning customer selects in a store is actually made from pork or venison, or does not contain it at all. The problem of food authenticity is widespread worldwide, and cases of meat adulteration have accelerated the development of food and the identification methods of feed species. It is worth noting that several different molecular biology techniques can identify a porcine component. However, the precise differentiation between wild boar and a domestic pig in meat products is still challenging. This paper presents the current state of knowledge concerning the species identification of the domestic pig and wild boar DNA in meat and its products.

Genome-wide association analysis of milk production, somatic cell score, and body conformation traits in Holstein cows

Milk production and body conformation traits are critical economic traits for dairy cows. To understand the basic genetic structure for those traits, a genome wide association study was performed on milk yield, milk fat yield, milk fat percentage, milk protein yield, milk protein percentage, somatic cell score, body form composite index, daily capacity composite index, feed, and leg conformation traits, based on the Illumina Bovine HD100k BeadChip. A total of 57, 12 and 26 SNPs were found to be related to the milk production, somatic cell score and body conformation traits in the Holstein cattle. Genes with pleiotropic effect were also found in this study. Seven significant SNPs were associated with multi-traits and were located on the PLEC, PLEKHA5, TONSL, PTGER4, and LCORL genes. In addition, some important candidate genes, like GPAT3, CEBPB, AGO2, SLC37A1, and FNDC3B, were found to participate in fat metabolism or mammary gland development. These results can be used as candidate genes for milk production, somatic cell score, and body conformation traits of Holstein cows, and are helpful for further gene function analysis to improve milk production and quality.

Along the Bos taurus genome, uncover candidate imprinting control regions

BACKGROUND

In mammals, Imprinting Control Regions (ICRs) regulate a subset of genes in a parent-of-origin-specific manner. In both human and mouse, previous studies identified a set of CpG-rich motifs occurring as clusters in ICRs and germline Differentially Methylated Regions (gDMRs). These motifs consist of the ZFP57 binding site (ZFBS) overlapping a subset of MLL binding units known as MLL morphemes. MLL or MLL1 (Mixed Lineage Leukemia 1) is a relatively large multidomain protein that plays a central role in the regulation of transcription. The structures of both MLL1 and MLL2 include a domain (MT) that binds CpG-rich DNA and a conserved domain (SET) that methylates lysine 4 in histone H3 producing H3K4me3 marks in chromatin.

RESULTS

Since genomic imprinting impacts many developmental and key physiological processes, we followed a previous bioinformatics strategy to pinpoint ICR positions in the Bos taurus genome. Initial genome-wide analyses involved finding the positions of ZFP57 binding sites, and the CpG-rich motifs (ZFBS-morph overlaps) along cattle chromosomal DNA. By creating plots displaying the density of ZFBS-morph overlaps, we removed background noise and thus improved signal detection. With the density-plots, we could view the positions of peaks locating known and candidate ICRs in cattle DNA. Our evaluations revealed the correspondence of peaks in plots to reported known and inferred ICRs/DMRs in cattle. Beside peaks pinpointing such ICRs, the density-plots also revealed additional peaks. Since evaluations validated the robustness of our approach, we inferred that the additional peaks may correspond to candidate ICRs for imprinted gene expression.

CONCLUSION

Our bioinformatics strategy offers the first genome-wide approach for systematically localizing candidate ICRs. Furthermore, we have tailored our datasets for upload onto the UCSC genome browser so that researchers could find known and candidate ICRs with respect to a wide variety of annotations at all scales: from the positions of Single Nucleotide Polymorphisms (SNPs), to positions of genes, transcripts, and repeated DNA elements. Furthermore, the UCSC genome browser offers tools to produce enlarged views: to uncover the genes in the vicinity of candidate ICRs and thus discover potential imprinted genes for experimental validations.

Detection of selection signatures in South African Mutton Merino sheep using whole-genome sequencing data

The South African Mutton Merino (SAMM), a dual-purpose (meat and wool) sheep breed, is characterized by its excellent performance on growth, carcass traits and meat quality compared to other fine-wool Merino breeds. Nowadays, the SAMM breed has been widely used to cross with commercial and indigenous fine-wool or coarse-wool breeds to improve the growth and meat performance in many countries. To date, however, little is known about the genetic basis for its prominent characteristics. In this study, whole-genome sequences of 10 SAMM were sequenced and the selection signatures were analyzed together with those of 39 Australian Merino and Chinese Merino (wool-type Merino) by F_ST , iHS, and XP-EHH methods. In total, 313 genes in 277 regions were identified by at least 2 methods with the signal of selection and 21 of them were identified by all three methods. We highlighted a list of interesting genes, including GHR, LCORL, SMO, NCAPG, DCC, IBSP, PPARGC1A, PACRGL, PRDM5, XYLB, AHCYL2, TEFM, AFG1L, and FAM184B, which have been shown to be involved in growth, carcass traits, and meat quality by previous studies. Herein, GHR, encoding a transmembrane receptor for growth hormone, is the most notable one. We report the first study on selection signatures analysis of SAMM at whole-genome sequence level. These results provide new insights into the genetic mechanisms underlying the growth and carcass traits in SAMM.

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.

Genomic diversity and signatures of selection in meat and fancy rabbit breeds based on high-density marker data

Background

Domestication of the rabbit (Oryctolagus cuniculus) has led to a multi-purpose species that includes many breeds and lines with a broad phenotypic diversity, mainly for external traits (e.g. coat colours and patterns, fur structure, and morphometric traits) that are valued by fancy rabbit breeders. As a consequence of this human-driven selection, distinct signatures are expected to be present in the rabbit genome, defined as signatures of selection or selective sweeps. Here, we investigated the genome of three Italian commercial meat rabbit breeds (Italian Silver, Italian Spotted and Italian White) and 12 fancy rabbit breeds (Belgian Hare, Burgundy Fawn, Champagne d’Argent, Checkered Giant, Coloured Dwarf, Dwarf Lop, Ermine, Giant Grey, Giant White, Rex, Rhinelander and Thuringian) by using high-density single nucleotide polymorphism data. Signatures of selection were identified based on the fixation index (F_ST) statistic with different approaches, including single-breed and group-based methods, the latter comparing breeds that are grouped based on external traits (different coat colours and body sizes) and types (i.e. meat vs. fancy breeds).

Results

We identified 309 genomic regions that contained signatures of selection and that included genes that are known to affect coat colour (ASIP, MC1R and TYR), coat structure (LIPH), and body size (LCORL/NCAPG, COL11A1 and HOXD) in rabbits and that characterize the investigated breeds. Their identification proves the suitability of the applied methodologies for capturing recent selection events. Other regions included novel candidate genes that might contribute to the phenotypic variation among the analyzed breeds, including genes for pigmentation-related traits (EDNRA, EDNRB, MITF and OCA2) and body size, with a strong candidate for dwarfism in rabbit (COL2A1).

Conclusions

We report a genome-wide view of genetic loci that underlie the main phenotypic differences in the analyzed rabbit breeds, which can be useful to understand the shift from the domestication process to the development of breeds in O. cuniculus. These results enhance our knowledge about the major genetic loci involved in rabbit external traits and add novel information to understand the complexity of the genetic architecture underlying body size in mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00696-9.

Genetic and genomic characterization followed by single-step genomic evaluation of withers height in German Warmblood horses

Reliability of genomic predictions is influenced by the size and genetic composition of the reference population. For German Warmblood horses, compilation of a reference population has been enabled through the cooperation of five German breeding associations. In this study, preliminary data from this joint reference population were used to genetically and genomically characterize withers height and to apply single-step methodology for estimating genomic breeding values for withers height. Using data on 2113 mares and their genomic information considering about 62,000 single nucleotide polymorphisms (SNPs), analysis of the genomic relationship revealed substructures reflecting breed origin and different breeding goals of the contributing breeding associations. A genome-wide association study confirmed a known quantitative trait locus (QTL) for withers height on equine chromosome (ECA) 3 close to LCORL and identified a further significant peak on ECA 1. Using a single-step approach with a combined relationship matrix, the estimated heritability for withers height was 0.31 (SE = 0.08) and the corresponding genomic breeding values ranged from − 2.94 to 2.96 cm. A mean reliability of 0.38 was realized for these breeding values. The analyses of withers height showed that compiling a reference population across breeds is a suitable strategy for German Warmblood horses. The single-step method is an appealing approach for practical genomic prediction in horses, because not many genotypes are available yet and animals without genotypes can by this way directly contribute to the estimation system.

Insight into the Candidate Genes and Enriched Pathways Associated with Height, Length, Length to Height Ratio and Body-Weight of Korean Indigenous Breed, Jindo Dog Using Gene Set Enrichment-Based GWAS Analysis

As a companion and hunting dog, height, length, length to height ratio (LHR) and body-weight are the vital economic traits for Jindo dog. Human selection and targeted breeding have produced an extraordinary diversity in these traits. Therefore, the identification of causative markers, genes and pathways that help us to understand the genetic basis of this variability is essential for their selection purposes. Here, we performed a genome-wide association study (GWAS) combined with enrichment analysis on 757 dogs using 118,879 SNPs. The genomic heritability (h²) was 0.33 for height and 0.28 for weight trait in Jindo. At p-value < 5 × 10^-5, ten, six, thirteen and eleven SNPs on different chromosomes were significantly associated with height, length, LHR and body-weight traits, respectively. Based on our results, HHIP, LCORL and NCAPG for height, IGFI and FGFR3 for length, DLK1 and EFEMP1 for LHR and PTPN2, IGFI and RASAL2 for weight can be the potential candidate genes because of the significant SNPs located in their intronic or upstream regions. The gene-set enrichment analysis highlighted here nine and seven overlapping significant (p < 0.05) gene ontology (GO) terms and pathways among traits. Interestingly, the highlighted pathways were related to hormone synthesis, secretion and signalling were generally involved in the metabolism, growth and development process. Our data provide an insight into the significant genes and pathways if verified further, which will have a significant effect on the breeding of the Jindo dog's population.

Genomic analyses reveal distinct genetic architectures and selective pressures in Chinese donkeys

Donkey (Equus asinus) is an important livestock animal in China because of its draft and medicinal value. After a long period of natural and artificial selection, the variety and phenotype of donkeys have become abundant. We clarified the genetic and demographic characteristics of Chinese domestic donkeys and the selection pressures by analyzing 78 whole genomes from 12 breeds. According to population structure, most Chinese domestic donkeys showed a dominant ancestral type. However, the Chinese donkeys still represented a significant geographical distribution trend. In the selective sweep, gene annotation, functional enrichment, and differential expression analyses between large and small donkey groups, we identified selective signals, including NCAPG and LCORL, which are related to rapid growth and large body size. Our findings elucidate the evolutionary history and formation of different donkey breeds and provide theoretical insights into the genetic mechanism underlying breed characteristics and molecular breeding programs of donkey clades.

From the Eurasian Steppes to the Roman Circuses: A Review of Early Development of Horse Breeding and Management

Simple Summary

Horses were domesticated later than any other major livestock species. Their role in shaping ancient civilizations cannot be overestimated. As a primary means of transportation, an essential asset in warfare, and later one of the key elements of circus entertainment, horses quickly became luxurious goods. Vast amounts of money were invested in the horse industry resulted resulting in the rapid development of horse breeding and husbandry. This review examines paleogenetic, archeological, and classical studies on managing horses in antiquity. Many ancient approaches and practices in horse management are still relevant today and some of them, now abandoned, are worth re-examination.

Abstract

The domestication of the horse began about 5500 years ago in the Eurasian steppes. In the following millennia horses spread across the ancient world, and their role in transportation and warfare affected every ancient culture. Ownership of horses became an indicator of wealth and social status. The importance of horses led to a growing interest in their breeding and management. Many phenotypic traits, such as height, behavior, and speed potential, have been proven to be a subject of selection; however, the details of ancient breeding practices remain mostly unknown. From the fourth millennium BP, through the Iron Age, many literature sources thoroughly describe horse training systems, as well as various aspects of husbandry, many of which are still in use today. The striking resemblance of ancient and modern equine practices leaves us wondering how much was accomplished through four thousand years of horse breeding.

Single nucleotide polymorphisms affect miRNA target prediction in bovine

Single nucleotide polymorphisms (SNPs) can have significant effects on phenotypic characteristics in cattle. MicroRNAs (miRNAs) are small, non-coding RNAs that act as post-transcriptional regulators by binding them to target mRNAs. In the present study, we scanned ~56 million SNPs against 1,064 bovine miRNA sequences and analyzed, in silico, their possible effects on target binding prediction, primary miRNA formation, association with QTL regions and the evolutionary conservation for each SNP locus. Following target prediction, we show that 71.6% of miRNA predicted targets were altered as a consequence of SNPs located within the seed region of the mature miRNAs. Next, we identified variations in the Minimum Free Energy (MFE), which represents the capacity to alter molecule stability and, consequently, miRNA maturation. A total of 48.6% of the sequences analyzed showed values within those previously reported as sufficient to alter miRNA maturation. We have also found 131 SNPs in 46 miRNAs, with altered target prediction, occurring in QTL regions. Lastly, analysis of evolutionary conservation scores for each SNP locus suggested that they have a conserved biological function through the evolutionary process. Our results suggest that SNPs in microRNAs have the potential to affect bovine phenotypes and could be of great value for genetic improvement studies, as well as production.

A QTL for conformation of back and croup influences lateral gait quality in Icelandic horses

BACKGROUND

The back plays a vital role in horse locomotion, where the spine functions as a spring during the stride cycle. A complex interaction between the spine and the muscles of the back contribute to locomotion soundness, gait ability, and performance of riding and racehorses. Conformation is commonly used to select horses for breeding and performance in multiple horse breeds, where the back and croup conformation plays a significant role. The conformation of back and croup plays an important role on riding ability in Icelandic horses. However, the genes behind this trait are still unknown. Therefore, the aim of this study was to identify genomic regions associated with conformation of back and croup in Icelandic horses and to investigate their effects on riding ability. One hundred seventy-seven assessed Icelandic horses were included in the study. A genome-wide association analysis was performed using the 670 K+ Axiom Equine Genotyping Array, and the effects of different haplotypes in the top associated region were estimated for riding ability and additional conformation traits assessed during breeding field tests.

RESULTS

A suggestive quantitative trait loci (QTL) for the score of back and croup was detected on Equus caballus (ECA) 22 (p-value = 2.67 × 10- 7). Haplotype analysis revealed two opposite haplotypes, which resulted in higher and lower scores of the back and croup, respectively (p-value < 0.001). Horses with the favorable haplotype were more inclined to have a well-balanced backline with an uphill conformation and had, on average, higher scores for the lateral gaits tölt (p-value = 0.02) and pace (p-value = 0.004). This genomic region harbors three genes: C20orf85, ANKRD60 and LOC100056167. ANKRD60 is associated with body height in humans. C20orf85 and ANKRD60 are potentially linked to adolescent idiopathic scoliosis in humans.

CONCLUSIONS

Our results show that the detected QTL for conformation of back and croup is of importance for quality of lateral gaits in Icelandic horses. These findings could result in a genetic test to aid in the selection of breeding horses, thus they are of major interest for horse breeders. The results may also offer a gateway to comparative functional genomics by potentially linking both motor laterality and back inclination in horses with scoliosis in humans.

Pathways involved in pony body size development

BACKGROUND

The mechanism of body growth in mammals is poorly understood. Here, we investigated the regulatory networks involved in body growth through transcriptomic analysis of pituitary and epiphyseal tissues of smaller sized Debao ponies and Mongolian horses at the juvenile and adult stages.

RESULTS

We found that growth hormone receptor (GHR) was expressed at low levels in long bones, although growth hormone (GH) was highly expressed in Debao ponies compared with Mongolian horses. Moreover, significant downregulated of the GHR pathway components m-RAS and ATF3 was found in juvenile ponies, which slowed the proliferation of bone osteocytes. However, WNT2 and PLCβ2 were obviously upregulated in juvenile Debao ponies, which led to premature mineralization of the bone extracellular matrix. Furthermore, we found that the WNT/Ca²⁺ pathway may be responsible for regulating body growth. GHR was demonstrated by q-PCR and Western blot analyses to be expressed at low levels in long bones of Debao ponies. Treatment with WNT antagonistI decreased the expression of WNT pathway components (P < 0.05) in vitro. Transduction of ATDC5 cells with a GHR-RNAi lentiviral vector decreased the expression of the GHR pathway components (P < 0.05). Additionally, the expression of the IGF-1 gene in the liver was lower in Debao ponies than in Mongolian horses at the juvenile and adult stages. Detection of plasma hormone concentrations showed that Debao ponies expressed higher levels of IGF-1 as juveniles and higher levels of GH as adults than Mongolian horses, indicating that the hormone regulation in Debao ponies differs from that in Mongolian horses.

CONCLUSION

Our work provides insights into the genetic regulation of short stature growth in mammals and can provide useful information for the development of therapeutic strategies for small size.

Genetic consistency between gait analysis by accelerometry and evaluation scores at breeding shows for the selection of jumping competition horses

The aim was to assess the efficiency of gaits characteristics in improving jumping performance of sport horses and confront accelerometers and judge scores for this purpose. A sample of 1,477 young jumping horses were measured using accelerometers for walk, trot, and canter. Of these, 702 were genotyped with 541,175 SNPs after quality control. Dataset of 26,914 horses scored by judges in breeding shows for gaits and dataset of 142,682 horses that performed in jumping competitions were used. Analysis of accelerometric data defined three principal components from 64% to 89% of variability explained for each gait. Animal mixed models were used to estimate genetic parameters with the inclusion to up 308,105 ancestors for the relationship matrix. Fixed effects for the accelerometric variables included velocity, gender, age, and event. A GWAS was performed on residuals with the fixed effect of each SNP. The GWAS did not reveal other QTLs for gait traits than the one related to the height at withers. The accelerometric principal components were highly heritable for the one linked to stride frequency and dorsoventral displacement at trot (0.53) and canter (0.41) and moderately for the one linked to longitudinal activities (0.33 for trot, 0.19 for canter). Low heritabilities were found for the walk traits. The genetic correlations of the accelerometric principal components with the jumping competition were essentially nil, except for a negative correlation with longitudinal activity at canter (-0.19). The genetic correlation between the judges’ scores and the jumping competition reached 0.45 for canter (0.31 for trot and 0.17 for walk). But these correlations turned negative when the scores were corrected for the known parental breeding value for competition at the time of the judging. In conclusion, gait traits were not helpful to select for jumping performances. Different gaits may be suitable for a good jumping horse.

Selection in the Finnhorse, a native all-around horse breed

Selection by breeders modifies the morphology, behaviour and performance of domesticated species. Here, we examined signs of selection in Finnhorse, the only native horse breed in Finland. We first searched divergent genomic regions between Finnhorses and other breeds, as well as between different breeding sections of the Finnhorse with data from Illumina Equine SNP70 BeadChip, and then studied several of the detected regions in more detail. We found altogether 35 common outlier SNPs between Finnhorses and other breeds using two different selection tests. Many of the SNPs were located close to genes affecting coat colour, performance, size, sugar metabolism, immune response and olfaction. We selected genes affecting coat colour (KIT, MITF, PMEL), performance (MSTN) and locomotion (DMRT3) for a more detailed examination. In addition, we looked for, and found, associations with height at withers and SNPs located close to gene LCORL. Among the four breeding sections of Finnhorses (harness trotters, riding horses, draught horses and pony‐sized horses), a single SNP located close to the DMRT3 gene was significantly differentiated and only between harness trotters and pony‐sized horses.

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.

Comparative population genomic analysis uncovers novel genomic footprints and genes associated with small body size in Chinese pony

Background

Body size is considered as one of the most fundamental properties of an organism. Due to intensive breeding and artificial selection throughout the domestication history, horses exhibit striking variations for heights at withers and body sizes. Debao pony (DBP), a famous Chinese horse, is known for its small body size and lives in Guangxi mountains of southern China. In this study, we employed comparative population genomics to study the genetic basis underlying the small body size of DBP breed based on the whole genome sequencing data. To detect genomic signatures of positive selection, we applied three methods based on population comparison, fixation index (F_ST), cross population composite likelihood ratio (XP-CLR) and nucleotide diversity (θπ), and further analyzed the results to find genomic regions under selection for body size-related traits.

Results

A number of protein-coding genes in windows with the top 1% values of F_ST (367 genes), XP-CLR (681 genes), and log₂ (θπ ratio) (332 genes) were identified. The most significant signal of positive selection was mapped to the NELL1 gene, probably underlies the body size and development traits, and may also have been selected for short stature in the DBP population. In addition, some other loci on different chromosomes were identified to be potentially involved in the development of body size.

Conclusions

Results of our study identified some positively selected genes across the horse genome, which are possibly involved in body size traits. These novel candidate genes may be useful targets for clarifying our understanding of the molecular basis of body size and as such they should be of great interest for future research into the genetic architecture of relevant traits in horse breeding program.

Whole-genome sequencing of European autochthonous and commercial pig breeds allows the detection of signatures of selection for adaptation of genetic resources to different breeding and production systems

Background

Natural and artificial directional selection in cosmopolitan and autochthonous pig breeds and wild boars have shaped their genomes and resulted in a reservoir of animal genetic diversity. Signatures of selection are the result of these selection events that have contributed to the adaptation of breeds to different environments and production systems. In this study, we analysed the genome variability of 19 European autochthonous pig breeds (Alentejana, Bísara, Majorcan Black, Basque, Gascon, Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda, Krškopolje pig, Black Slavonian, Turopolje, Moravka, Swallow-Bellied Mangalitsa, Schwäbisch-Hällisches Schwein, Lithuanian indigenous wattle and Lithuanian White old type) from nine countries, three European commercial breeds (Italian Large White, Italian Landrace and Italian Duroc), and European wild boars, by mining whole-genome sequencing data obtained by using a DNA-pool sequencing approach. Signatures of selection were identified by using a single-breed approach with two statistics [within-breed pooled heterozygosity (H_P) and fixation index (F_ST)] and group-based F_ST approaches, which compare groups of breeds defined according to external traits and use/specialization/type.

Results

We detected more than 22 million single nucleotide polymorphisms (SNPs) across the 23 compared populations and identified 359 chromosome regions showing signatures of selection. These regions harbour genes that are already known or new genes that are under selection and relevant for the domestication process in this species, and that affect several morphological and physiological traits (e.g. coat colours and patterns, body size, number of vertebrae and teats, ear size and conformation, reproductive traits, growth and fat deposition traits). Wild boar related signatures of selection were detected across all the genome of several autochthonous breeds, which suggests that crossbreeding (accidental or deliberate) occurred with wild boars.

Conclusions

Our findings provide a catalogue of genetic variants of many European pig populations and identify genome regions that can explain, at least in part, the phenotypic diversity of these genetic resources.

Equine Genotyping Arrays

High-quality genomic tools have been integral in understanding genomic architecture and function in the modern-day horse. The equine genetics community has a long tradition of pooling resources to develop genomic tools. Since the equine genome was sequenced in 2006, several iterations of high throughput genotyping arrays have been developed and released, enabling rapid and cost-effective genotyping. This review highlights the design considerations of each iteration, focusing on data available during development and outlining considerations in selecting the genetic variants included on each array. Additionally, we outline recent applications of equine genotyping arrays as well as future prospects and applications.

Signatures of selection analysis using whole-genome sequence data reveals novel candidate genes for pony and light horse types

Natural selection and domestication have shaped modern horse populations, resulting in a vast range of phenotypically diverse breeds. Horse breeds are classified into three types (pony, light, and draft) generally based on their body type. Understanding the genetic basis of horse type variation and selective pressures related to the evolutionary trend can be particularly important for current selection strategies. Whole-genome sequences were generated for 14 pony and 32 light horses to investigate the genetic signatures of selection of the horse type in pony and light horses. In the overlapping extremes of the fixation index and nucleotide diversity results, we found novel genomic signatures of selective sweeps near key genes previously implicated in body measurements including C4ORF33, CRB1, CPN1, FAM13A, and FGF12 that may influence variation in pony and light horse types. This study contributes to a better understanding of the genetic background of differences between pony and light horse types.

Novel Variants in the HMGA2 Gene Are Associated With Withers Height in Debao Pony

The Debao pony is a well-known dwarf horse breed in China. High-mobility group AT-hook 2 (HMGA2) gene is regarded as one of the important candidate genes regulating body height in horses. The aim of this study was to study the association between mutations in HMGA2 gene and withers height in Debao ponies. The polymorphisms in all exons and partial introns of the HMGA2 gene were screened with sequencing across 180 Debao ponies. And the association between the DNA variants and withers height was analyzed. Seven genetic variants were identified in HMGA2 gene, including six novel variants. Among them, six mutations were located in two closed linked blocks. The three novel variants (In1-1, E5-1, and E5-2) in the 1st intron and the fifth exon and a known mutation (In1-2) had significant association with withers height in Debao ponies. These results suggest that the four variants have the potential to be used as genetic markers for dwarf horse breeding activities.

Revisiting the Population Genetics of Human Height

CONTEXT

Recent advances in genetics and genomics present unique opportunities for enhancing knowledge of human physiology and disease susceptibility. An outstanding example of these new insights may be seen in the study of human height, of which it has been estimated that approximately 80% is genetically determined. Over the past decade, large-scale population analyses have led to the identification of novel variation in genes and loci individually associated with changes in adult height of as much as 2 cm.

OBJECTIVE

To assess these same variants in the genomes of 213 158 individuals compiled by the Genome Aggregation Database (GnomAD) consortium, representing different population groups from around the world.

RESULTS

The majority of these height-changing alleles are substantially less prevalent in GnomAD than found previously in other cohorts, with 4 of 5 amino acid substitution variants with the largest impact on adult height being more frequent in the European population than in other groups.

CONCLUSIONS

A larger-scale analysis of individuals from diverse backgrounds will be necessary to ensure a full and accurate understanding of the genetic underpinnings of human height throughout the world, and additional studies will be needed to discern the biochemical and molecular mechanisms governing the physiological processes that explain how these variant proteins might selectively impact the biology of the growth plate. Broader understanding of the genetics of height also should set the stage for more comprehensive investigation into the causes of prevalent polygenic human diseases.