Polymorphisms in ten candidate genes are associated with conformational and locomotive traits in Spanish Purebred horses

Selection signatures and inbreeding: exploring genetic diversity in five native horse breeds

Horses have undergone extensive natural and artificial selection, shaping the diversity of breeds observed today. Native Italian breeds present unique traits influenced by natural selection, such as adaptation to harsh climates, or hoof strength, but face challenges due to population declines and the reduction of their original breeding purpose. This study focuses on five local Italian breeds: Bardigiano, Haflinger, Maremmano, Murgese, and Italian Heavy Draught Horse, to understand how selection has shaped their populations. A total of 1620 individuals were genotyped with a medium-density SNP chip and remapped to EquCab3. After quality control, where data were filtered based on missing genotypes per SNP (> 0.10) and missing SNPs per sample (> 0.10), 1498 horses and 54,825 SNPs remained for analysis. Population structure and runs of homozygosity (ROH) were identified, and genomic inbreeding coefficients were calculated based on ROH coverage of autosomal SNPs. ROH islands shared by ≥ 70% of horses were identified as selection signatures, and candidate genes within these regions were annotated. The inbreeding coefficient (FROH) ranged from 0.15 to 0.23, with Bardigiano and Haflinger showing the highest values probably due to selective breeding, while Maremmano, Murgese, and Italian Heavy Draught Horse displayed lower FROH, reflecting a broader diversity. ROH islands were identified on 12 chromosomes, with 23 islands distributed among breeds. Cold-blooded breeds (Bardigiano, Haflinger, and Italian Heavy Draught Horse) showed the majority, particularly on Equine Chromosome 3 (ECA3). These islands overlapped with 83 quantitative trait loci (QTLs) and 76 genes associated with morphology and health. Health-related traits such as osteochondrosis and hoof health were linked to ROH patterns, particularly in Bardigiano and Haflinger, highlighting selection for disease resistance. Signature of selections were found in the proximity of MC1R and ASIP genes likely due to their role for coat color; especially in the Haflinger and Italian Heavy Draught Horse the genotype frequency of the BIEC2_816499 SNP which is in the vicinity of the causative mutation for chestnut coat color is due to linkage disequilibrium between the two. In conclusion, this study offered valuable insights that breeders could utilize to make sound decisions. This issue would ensure the maintenance of breed genetic diversity, and the preservation and improvement of the breed's distinct traits and health standards.

Osteochondrosis in horses: An overview of genetic and other factors

Osteochondrosis (OC) is a frequent manifestation of developmental orthopaedic disease, and its severe clinical presentation is known as OC dissecans (OCD). OC is defined as a disruption of the endochondral ossification process in the epiphyseal cartilage, and this disease has been reported in different mammalian species, including humans, dogs, pigs, and horses. OCD is an important cause of lameness in sport horses and is a common cause of impaired orthopaedic potential, whose clinical signs may be of minimal magnitude or manifest as severe joint effusion or clinically noticeable lameness. The aetiology of OCD is unknown, although it has traditionally been considered to be multifactorial. In addition to genetic factors, associated factors include both non-genetic elements such as rapid growth, nutrition, trauma, anatomical conformation, and biomechanics. Since the prevalence of the disease varies greatly depending on the horse breed, from 13% in Swedish Warmblood to 53% in Lusitano breed, genetic factors have a great relevance in the appearance and development of OCD in horses. Many genetic modifications have been related, and the genes involved can be grouped into five clusters, related to fundamental functions for the correct development and regeneration of cartilage, such as collagen, laminin, cell signalling, matrix turnover, and transcriptional regulation. Changes in genes such as COL3A1, COL5A1, COL5A2, COL24A1, COL27A1 (collagen cluster), LAMB1 (laminin cluster), PTH, PHT receptors, and IHH (cell signalling), and genes encoding matrix metalloproteinases have been related to the occurrence and severity of diseases in different equine breeds. This review summarises the main factors associated with OC in horses, with particular emphasis on genetic factors.

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.

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.

Accurate haplotype construction and detection of selection signatures enabled by high quality pig genome sequences

High-quality whole-genome resequencing in large-scale pig populations with pedigree structure and multiple breeds would enable accurate construction of haplotype and robust selection-signature detection. Here, we sequence 740 pigs, combine with 149 of our previously published resequencing data, retrieve 207 resequencing datasets, and form a panel of worldwide distributed wild boars, aboriginal and highly selected pigs with pedigree structures, amounting to 1096 genomes from 43 breeds. Combining with their haplotype-informative reads and pedigree structure, we accurately construct a panel of 1874 haploid genomes with 41,964,356 genetic variants. We further demonstrate its valuable applications in GWAS by identifying five novel loci for intramuscular fat content, and in genomic selection by increasing the accuracy of estimated breeding value by 36.7%. In evolutionary selection, we detect MUC13 gene under a long-term balancing selection, as well as NPR3 gene under positive selection for pig stature. Our study provides abundant genomic variations for robust selection-signature detection and accurate haplotypes for deciphering complex traits in pigs.

Assessment and Distribution of Runs of Homozygosity in Horse Breeds Representing Different Utility Types

The present study reports runs of homozygosity (ROH) distribution in the genomes of six horse breeds (571 horses in total) representing three horse types (primitive, light, and draft horses) based on the 65k Equine BeadChip assay. Of major interest was the length, quantity, and frequency of ROH characteristics, as well as differences between horse breeds and types. Noticeable differences in the number, length and distribution of ROH between breeds were observed, as well as in genomic inbreeding coefficients. We also identified regions of the genome characterized by high ROH coverage, known as ROH islands, which may be signals of recent selection events. Eight to fourteen ROH islands were identified per breed, which spanned multiple genes. Many were involved in important horse breed characteristics, including WFIKNN2, CACNA1G, STXBP4, NOG, FAM184B, QDPR, LCORL, and the zinc finger protein family. Regions of the genome with zero ROH occurrences were also of major interest in specific populations. Depending on the breed, we detected between 2 to 57 no-ROH regions and identified 27 genes in these regions that were common for five breeds. These genes were involved in, e.g., muscle contractility (CACNA1A) and muscle development (miR-23, miR-24, miR-27). To sum up, the obtained results can be furthered analyzed in the topic of identification of markers unique for specific horse breed characteristics.

Genomic signatures for drylands adaptation at gene-rich regions in African zebu cattle

Background

Indigenous Sudanese cattle are mainly indicine/zebu (humped) type. They thrive in the harshest dryland environments characterised by high temperatures, long seasonal dry periods, nutritional shortages, and vector disease challenges. Here, we sequenced 60 indigenous Sudanese cattle from six indigenous breeds and analysed the data using three genomic scan approaches to unravel cattle adaptation to the African dryland region.

Results

We identified a set of gene-rich selective sweep regions, detected mostly on chromosomes 5, 7 and 19, shared across African and Gir zebu. These include genes involved in immune response, body size and conformation, and heat stress response. We also identified selective sweep regions unique to Sudanese zebu. Of these, a 250 kb selective sweep on chromosome 16 spans seven genes, including PLCH2, PEX10, PRKCZ, and SKI, which are involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

Conclusions

Our results suggest that environmental adaptation may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control, in zebu cattle.

•

Sudanese cattle thrive in the harshest environments of the African drylands.

•

Bos indicus shared selected genes are involved in immune response, conformation, and heat stress response.

•

Sudanese zebu-specific sweep includes genes involved in alternative adaptive metabolic strategies of insulin signalling, glucose homeostasis, and fat metabolism.

•

Environmental adaptation in zebu cattle may involve recent and ancient selection at gene-rich regions, which might be under a common regulatory genetic control.

The expression of the gastrin/cholecystokinin (GAST/CCK) family and their receptors (CCKAR/CCKBR) in the chicken changes in response to quantitative restriction and reveals a functional role of CCK in the crop

Gastrin and cholecystokinin peptides bind a common G-protein coupled receptor, cholecystokinin receptor B (CCKBR) whilst cholecystokinin receptor A (CCKAR) is preferentially bound by CCK. Gastrin and cholecystokinin mediate signalling from the gastrointestinal tract to regulate appetite and digestive function. In this study, expression of the cholecystokinin/gastrin family and distribution of their receptors expression was measured to understand the target organs for the peptides and how expression responds to changes in food intake. We confirmed the restricted expression of gastrin in the antrum and the abundant expression of cholecystokinin in the hypothalamus. The expression of gastrin in the antrum was significantly elevated in broiler breeders when released from feed restriction. CCKBR was most abundant in the hypothalamus and proventriculus. CCKAR was most abundant in the pancreas and crop, more than tenfold greater than the gastrointestinal tract. Cholecystokinin expression in the pancreas increased after removal of food restriction. CCKAR in the gastrointestinal tract peaks around the distal ileum, distal to the peak of cholecystokinin expression. There was virtually no cholecystokinin expression in the caecum but CCKAR expression was high. The CCKAR expression in the crop was unexpected, supporting a role of cholecystokinin in mediating crop emptying which was supported by the observation of in-vitro contraction after cholecystokinin administration. The response to changes in food intake and the expression pattern of the cholecystokinin/gastrin family and their receptors will stimulate and inform new hypotheses on their role in growth in poultry.

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.

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.

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.

Signatures of selection in the genome of Swedish warmblood horses selected for sport performance

Background

A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits.

Results

The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (F_ST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses.

Conclusions

Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.

Morphological and genetic diversity of Pura Raza Español horse with regard to the coat colour

Gene mutations influencing melanocytes also impact on physiological and behavioural functions. In this study, we investigated their association with four different coat colours in the Pura Raza Español (PRE) horse using morphological traits and molecular datasets. Four different subpopulations were identified according to individual coat colour: grey, bay, chestnut and black. Coat colour significantly associated with morphological measurements. Observed and expected heterozygosity values were low in grey compared with the other three subpopulations, suggesting the presence of unique ancestral alleles probably arisen by genetic drift and selection mechanism effects. Nei's distance demonstrated a clear division among subpopulations, the grey being the most divergent group. Gene flow estimates were similar, showing the lowest values in grey. Divergence times among subpopulations assessed with the average square distance suggested that grey was the original PRE population which diverged from bay, chestnut and black. Our results also demonstrated a clear morphological differentiation according to coat colour. The close genetic structure of bay and chestnut PRE subpopulations and the clear differences in most morphological traits of grey and chestnut PRE mares would suggest the pleiotropic effect of genomic regions determining coat colour in horses. However, further analysis including genomic information would be necessary to elucidate the mechanisms involved.

Genome-wide association study identifies loci and candidate genes for internal organ weights in Simmental beef cattle

Cattle internal organs as accessible raw materials have a long history of being widely used in beef processing, feed and pharmaceutical industry. These traits not only are of economic interest to breeders, but they are intrinsically linked to many valuable traits, such as growth, health, and productivity. Using the Illumina Bovine HD 770K SNP array, we performed a genome-wide association study for heart weight, liver weight, spleen weight, lung weight, and kidney weight in 1,217 Simmental cattle. In our research, 38 significant single nucleotide polymorphisms (SNPs) ( P < 1.49 × 10−6) were identified for five internal organ weight traits. These SNPs are within or near 13 genes, and some of them have been reported previously, including NDUFAF4, LCORL, BT.94996, SLIT2, FAM184B, LAP3, BBS12, MECOM, CD300LF, HSD17B3, TLR4, MXI1, and MB21D2. In addition, we detected four haplotype blocks on BTA6 containing 18 significant SNPs associated with spleen weight. Our results offer worthy insights into understanding the genetic mechanisms of internal organs' development, with potential application in breeding programs of Simmental beef cattle.