Whole-Genome Signatures of Selection in Sport Horses Revealed Selection Footprints Related to Musculoskeletal System Development Processes

Genome-wide analysis of population structure, effective population size and inbreeding in Iranian and exotic horses

Population structure and genetic diversity are the key parameters to study the breeding history of animals. This research aimed to provide a characterization of the population structure and to compare the effective population size (Ne), LD decay, genetic diversity, and genomic inbreeding in Iranian native Caspian (n = 38), Turkmen (n = 24) and Kurdish (n = 29) breeds and some other exotic horses consisting of Arabian (n = 24), Fell pony (n = 21) and Akhal-Teke (n = 20). A variety of statistical population analysis techniques, such as principal component analysis (PCA), discriminant analysis of principal component (DAPC) and model-based method (STRUCTURE) were employed. The results of the population analysis clearly demonstrated a distinct separation of native and exotic horse breeds and clarified the relationships between studied breeds. The effective population size (Ne) for the last six generations was estimated 54, 49, 37, 35, 27 and 26 for the Caspian, Kurdish, Arabian, Turkmen, Akhal-Teke and Fell pony breeds, respectively. The Caspian breed showed the lowest LD with an average r2 value of 0.079, while the highest was observed in Fell pony (0.148). The highest and lowest average observed heterozygosity were found in the Kurdish breeds (0.346) and Fell pony (0.290) breeds, respectively. The lowest genomic inbreeding coefficient based on run of homozygosity (FROH) and excess of homozygosity (FHOM) was in the Caspian and Kurdish breeds, respectively, while based on genomic relationship matrix) FGRM) and correlation between uniting gametes) FUNI) the lowest genomic inbreeding coefficient was found in the Kurdish breed. The estimation of genomic inbreeding rates in the six breeds revealed that FROH yielded lower estimates compared to the other three methods. Additionally, the Iranian breeds displayed lower levels of inbreeding compared to the exotic breeds. Overall, the findings of this study provide valuable insights for the development of effective breeding management strategies aimed at preserving these horse breeds.

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach

BACKGROUND

The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds.

RESULTS

The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne).

CONCLUSIONS

Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

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

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

Runs of homozygosity in Sable Island feral horses reveal the genomic consequences of inbreeding and divergence from domestic breeds

BACKGROUND

Understanding inbreeding and its impact on fitness and evolutionary potential is fundamental to species conservation and agriculture. Long stretches of homozygous genotypes, known as runs of homozygosity (ROH), result from inbreeding and their number and length can provide useful population-level information on inbreeding characteristics and locations of signatures of selection. However, the utility of ROH for conservation is limited for natural populations where baseline data and genomic tools are lacking. Comparing ROH metrics in recently feral vs. domestic populations of well understood species like the horse could provide information on the genetic health of those populations and offer insight into how such metrics compare between managed and unmanaged populations. Here we characterized ROH, inbreeding coefficients, and ROH islands in a feral horse population from Sable Island, Canada, using ~41 000 SNPs and contrasted results with those from 33 domestic breeds to assess the impacts of isolation on ROH abundance, length, distribution, and ROH islands.

RESULTS

ROH number, length, and ROH-based inbreeding coefficients (F_ROH) in Sable Island horses were generally greater than in domestic breeds. Short runs, which typically coalesce many generations prior, were more abundant than long runs in all populations, but run length distributions indicated more recent population bottlenecks in Sable Island horses. Nine ROH islands were detected in Sable Island horses, exhibiting very little overlap with those found in domestic breeds. Gene ontology (GO) enrichment analysis for Sable Island ROH islands revealed enrichment for genes associated with 3 clusters of biological pathways largely associated with metabolism and immune function.

CONCLUSIONS

This study indicates that Sable Island horses tend to be more inbred than their domestic counterparts and that most of this inbreeding is due to historical bottlenecks and founder effects rather than recent mating between close relatives. Unique ROH islands in the Sable Island population suggest adaptation to local selective pressures and/or strong genetic drift and highlight the value of this population as a reservoir of equine genetic variation. This research illustrates how ROH analyses can be applied to gain insights into the population history, genetic health, and divergence of wild or feral populations of conservation concern.

Detection of common copy number of variation underlying selection pressure in Middle Eastern horse breeds using whole-genome sequence data

Dareshouri, Arabian, and Akhal-Teke are three Middle Eastern horse breeds that have been selected for endurance and adaptation to harsh climates. Deciphering the genetic characteristics of these horses by tracing selection footprints and copy number of variations will be helpful in improving our understanding of equine breeds' development and adaptation. For this purpose, we sequenced the whole-genome of four Dareshouri horses using Illumina Hiseq panels and compared them with publicly available whole-genome sequences of Arabian (n=3) and Akhal-Teke (n=3) horses . Three tests of FLK, hapFLK, and pooled heterozygosity were applied using a sliding window (window size=100kb, step size=50kb) approach to detect putative selection signals. Copy number variation analysis was applied to investigate copy number of variants (CNVs), and the results were used to suggest selection signatures involving CNVs. Whole-genome sequencing demonstrated 8,837,950 single nucleotide polymorphisms (SNPs) in autosomal chromosomes. We suggested 58 genes and three quantitative trait loci (QTLs), including some related to horse gait, insect bite hypersensitivity, and withers height, based on selective signals detected by adjusted p-value of Mahalanobis distance based on the rank-based P-values (Md-rank-P) method. We proposed 12 genomic regions under selection pressure involving CNVs which were previously reported to be associated with metabolism energy (SLC5A8), champagne dilution in horses (SLC36A1), and synthesis of polyunsaturated fatty acids (FAT2). Only 10 Middle Eastern horses were tested in this study; therefore, the conclusions are speculative. Our findings are useful to better understanding the evolution and adaptation of Middle Eastern horse breeds.

A genome-wide scan for signatures of selection in Kurdish horse breed

The genetic diversity and genomic regions being under putative natural selection in Kurdish horse population were studied. The samples from 72 horses were genotyped by using GGP Equine 70K SNP arrays. The Ne Slope (NeS) analyses revealed that a sharp decline in Ne has probably occurred around four generations ago, and high frequency of ROH with 2-4 Mbp in length suggested that the inbreeding has probably occurred around 20 generations ago. The effective population size (Ne) was 104 horses up to three generations ago and the average inbreeding (F_ROH) was 0.047(±0.045). Using de-correlated composite of multiple selection signals (DCMS) and runs of homozygosity (ROH) analyses the genomic regions being under putative selection were detected. By using DCMS, a total of 148 significant SNP (FDR < 0.05) were identified, 40% of which were located on ECA9, where the greatest peak was observed. This genomic region harbors several known QTL which are associated with withers height (body size). Also significant genomic regions (FDR<0.05), harboring QTL associated with insect bite hypersensitivity (IBH), hair density and coat texture, alternate gaits, guttural pouch tympany and temperament were identified. By using outputs of ROH analyses, two hotspot regions (i.e. 30% of individuals was considered as threshold), were identified on ECA7 (50.11-54.36 Mbp) and ECA11 (26.10-29.07 Mbp) harboring QTL associated with withers height, alternate gait and IBH. In summary, the genomic regions being under putative natural selection which harbors known QTL associated with body size and IBH, among others, were introduced. Nevertheless, additional functional and comparative studies are necessary to corroborate their effect on the observed genetic and phenotypic diversity of the Kurdish horses.

A single-nucleotide mutation within the TBX3 enhancer increased body size in Chinese horses

Chinese ponies are endemic to the mountainous areas of southwestern China and were first reported in the archaeological record at the Royal Tomb of Zhongshan King, Mancheng, dated to approximately ∼2,100 YBP.¹ Previous work has started uncovering the genetic basis of size variation in western ponies and horses, revealing a limited number of loci, including HMGA2,²LCORL/NCAPG,³ZFAT, and LASP1.^4,5 Whether the same genetic pathways also drive the small body size of Chinese ponies, which show striking anatomical differences to Shetland ponies,⁶ remains unclear.^2,7 To test this, we combined whole-genome sequences of 187 horses across China. Statistical analyses revealed top association between genetic variation at the T-box transcription factor 3 (TBX3) and the body size. Fine-scale analysis across an extended population of 189 ponies and 574 horses narrowed down the association to one A/G SNP at an enhancer region upstream of the TBX3 (ECA8:20,644,555, p = 2.34e−39). Luciferase assays confirmed the single-nucleotide G mutation upregulating TBX3 expression, and enhancer-knockout mice exhibited shorter limbs than wild-type littermates (p < 0.01). Re-analysis of ancient DNA data showed that the G allele, which is most frequent in modern horses, first occurred some ∼2,300 years ago and rose in frequency since. This supports selection for larger size in Asia from approximately the beginning of the Chinese Empire. Overall, this study characterized the causal regulatory mutation underlying small body size in Chinese ponies and revealed size as one of the main selection targets of past Chinese breeders.

•

One single A/G SNP in TBX3 enhancer region drives size variation in Chinese horses

•

The frequency of the G variant correlates positively with size in 763 horses

•

Cellular and mice models confirm it affects TBX3 transcription and the limb length

•

The G variant first occurred ∼2,300 years ago and rose in frequency since

Shared Ancestry and Signatures of Recent Selection in Gotland Sheep

Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation F_ST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the F_ST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.

Signatures of selection reveal candidate genes involved in economic traits and cold acclimation in five Swedish cattle breeds

Background

Thousands of years of natural and artificial selection have resulted in indigenous cattle breeds that are well-adapted to the environmental challenges of their local habitat and thereby are considered as valuable genetic resources. Understanding the genetic background of such adaptation processes can help us design effective breeding objectives to preserve local breeds and improve commercial cattle. To identify regions under putative selection, GGP HD 150 K single nucleotide polymorphism (SNP) arrays were used to genotype 106 individuals representing five Swedish breeds i.e. native to different regions and covering areas with a subarctic cold climate in the north and mountainous west, to those with a continental climate in the more densely populated south regions.

Results

Five statistics were incorporated within a framework, known as de-correlated composite of multiple signals (DCMS) to detect signatures of selection. The obtained p-values were adjusted for multiple testing (FDR < 5%), and significant genomic regions were identified. Annotation of genes in these regions revealed various verified and novel candidate genes that are associated with a diverse range of traits, including e.g. high altitude adaptation and response to hypoxia (DCAF8, PPP1R12A, SLC16A3, UCP2, UCP3, TIGAR), cold acclimation (AQP3, AQP7, HSPB8), body size and stature (PLAG1, KCNA6, NDUFA9, AKAP3, C5H12orf4, RAD51AP1, FGF6, TIGAR, CCND2, CSMD3), resistance to disease and bacterial infection (CHI3L2, GBP6, PPFIBP1, REP15, CYP4F2, TIGD2, PYURF, SLC10A2, FCHSD2, ARHGEF17, RELT, PRDM2, KDM5B), reproduction (PPP1R12A, ZFP36L2, CSPP1), milk yield and components (NPC1L1, NUDCD3, ACSS1, FCHSD2), growth and feed efficiency (TMEM68, TGS1, LYN, XKR4, FOXA2, GBP2, GBP5, FGD6), and polled phenotype (URB1, EVA1C).

Conclusions

We identified genomic regions that may provide background knowledge to understand the mechanisms that are involved in economic traits and adaptation to cold climate in cattle. Incorporating p-values of different statistics in a single DCMS framework may help select and prioritize candidate genes for further analyses.

Genetic Variability in the Italian Heavy Draught Horse from Pedigree Data and Genomic Information

This study aimed to investigate the genetic diversity in the Italian Heavy Horse Breed from pedigree and genomic data. Pedigree information for 64,917 individuals were used to assess inbreeding level, effective population size (Ne), and effective numbers of founders and ancestors (fa/fe). Genotypic information from SNP markers were available for 267 individuals of both sexes, and it allowed estimating genomic inbreeding in two methods (observed versus expected homozygosity and from ROH) to study the breed genomic structure and possible selection signatures. Pedigree and genomic inbreeding were greatly correlated (0.65 on average). The inbreeding trend increased over time, apart from periods in which the base population enlarged, when Ne increased also. Recent bottlenecks did not occur in the genome, as fa/fe have shown. The observed homozygosity results were on average lower than expected, which was probably due to the use of French Breton stallions to support the breed genetic variability. High homozygous regions suggested that inbreeding increased in different periods. Two subpopulations were distinguished, which was probably due to the different inclusion of French animals by breeders. Few selection signatures were found at the population level, with possible associations to disease resistance. The almost low inbreeding rate suggested that despite the small breed size, conservation actions are not yet required.

Genetic Diversity and Signatures of Selection in a Native Italian Horse Breed Based on SNP Data

Simple Summary

The Bardigiano horse is a native Italian breed bred for living in rural areas, traditionally used in agriculture. The breed counts about 3000 horses, and it is nowadays mainly used for recreational purposes. The relatively small size and the closed status of the breed raise the issue of monitoring genetic diversity. We therefore characterized the breed’s genetic diversity based on molecular data. We showed a critical reduction of genetic variability mainly driven by past bottlenecks. We also highlighted homozygous genomic regions that might be the outcome of directional selection in recent years, in line with the conversion of Bardigiano horses from agricultural to riding purposes.

Abstract

Horses are nowadays mainly used for sport and leisure activities, and several local breeds, traditionally used in agriculture, have been exposed to a dramatic loss in population size and genetic diversity. The loss of genetic diversity negatively impacts individual fitness and reduces the potential long-term survivability of a breed. Recent advances in molecular biology and bioinformatics have allowed researchers to explore biodiversity one step further. This study aimed to evaluate the loss of genetic variability and identify genomic regions under selection pressure in the Bardigiano breed based on GGP Equine70k SNP data. The effective population size based on Linkage Disequilibrium (N_e) was equal to 39 horses, and it showed a decline over time. The average inbreeding based on runs of homozygosity (ROH) was equal to 0.17 (SD = 0.03). The majority of the ROH were relatively short (91% were ≤2 Mbp long), highlighting the occurrence of older inbreeding, rather than a more recent occurrence. A total of eight ROH islands, shared among more than 70% of the Bardigiano horses, were found. Four of them mapped to known quantitative trait loci related to morphological traits (e.g., body size and coat color) and disease susceptibility. This study provided the first genome-wide scan of genetic diversity and selection signatures in an Italian native horse breed.

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.