High-resolution population structure and runs of homozygosity reveal the genetic architecture of complex traits in the Lipizzan horse

Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies

BACKGROUND

The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds.

RESULTS

HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20).

CONCLUSIONS

Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Genetic Diversity of Kazakhstani Equus caballus (Linnaeus, 1758) Horse Breeds Inferred from Microsatellite Markers

Understanding the genetic diversity and structure of domesticated horse (Equus caballus) populations is critical for long-term herd management and breeding programs. This study examines 435 horses from Kazakhstan, covering seven groups in three geographic areas using 11 STR markers. Identified are 136 alleles, with the mean number of alleles per locus ranging from 9 to 19. VHL20 is the most variable locus across groups, while loci HTG4, AHT4, AHT5, HTG7, and HMS3 are variable in most populations. The locus AHT5 in the Emba population shows the highest frequency of rare alleles, while the lowest frequency, 0.005, is observed in the Kulandy population. All loci were highly informative for the Kazakhstani populations of E. caballus, with PIC values higher than 0.5. Pairwise variations in Wright's FST distances show that the examined varieties have little genetic differentiation (0.05%), indicating a high degree of admixture and a continuing lineage sorting process. Phylogenetic and population structure analyses reveal three major clusters of Kazakh horses, representing (I) the Uralsk population of the Kushum breed and the monophyly of two groups: (II) the Kozhamberdy population of the Mugalzhar breed, and (III) the Mugalzhar-Kushum breed populations. Kazakhstani horse populations, while being regionally isolated, were recently in contact with each other.

The Genetic Diversity of Stallions of Different Breeds in Russia

The specifics of breeding and selection significantly affect genetic diversity and variability within a breed. We present the data obtained from the genetic analysis of 21 thoroughbred and warmblood horse breeds. The most detailed information is described from the following breeds: Arabian, Trakehner, French Trotter, Standardbred, and Soviet Heavy Horse. The analysis of 509,617 SNP variants in 87 stallions from 21 populations made it possible to estimate the genetic diversity at the genome-wide level and distinguish the studied horse breeds from each other. In this study, we searched for heterozygous and homozygous ROH regions, evaluated inbreeding using FROH analysis, and generated a population structure using Admixture 1.3 software. Our findings indicate that the Arabian breed is an ancestor of many horse breeds. The study of the full-genome architectonics of breeds is of great practical importance for preserving the genetic characteristics of breeds and managing breeding. Studies were carried out to determine homozygous regions in individual breeds and search for candidate genes in these regions. Fifty-six candidate genes for the influence of selection pressure were identified. Our research reveals genetic diversity consistent with breeding directions and the breeds' history of origin.

Identification of runs of homozygosity in Western honey bees (Apis mellifera) using whole-genome sequencing data

Runs of homozygosity (ROH) are continuous homozygous segments that arise through the transmission of haplotypes that are identical by descent. The length and distribution of ROH segments provide insights into the genetic diversity of populations and can be associated with selection signatures. Here, we analyzed reconstructed whole-genome queen genotypes, from a pool-seq data experiment including 265 Western honeybee colonies from Apis mellifera mellifera and Apis mellifera carnica. Integrating individual ROH patterns and admixture levels in a dynamic population network visualization allowed us to ascertain major differences between the two subspecies. Within A. m. mellifera, we identified well-defined substructures according to the genetic origin of the queens. Despite the current applied conservation efforts, we pinpointed 79 admixed queens. Genomic inbreeding (F ROH) strongly varied within and between the identified subpopulations. Conserved A. m. mellifera from Switzerland had the highest mean F ROH (3.39%), while queens originating from a conservation area in France, which were also highly admixed, showed significantly lower F ROH (0.45%). The majority of A. m. carnica queens were also highly admixed, except 12 purebred queens with a mean F ROH of 2.33%. Within the breed-specific ROH islands, we identified 14 coding genes for A. m. mellifera and five for A. m. carnica, respectively. Local adaption of A. m. mellifera could be suggested by the identification of genes involved in the response to ultraviolet light (Crh-BP, Uvop) and body size (Hex70a, Hex70b), while the A. m. carnica specific genes Cpr3 and Cpr4 are most likely associated with the lighter striping pattern, a morphological phenotype expected in this subspecies. We demonstrated that queen genotypes derived from pooled workers are useful tool to unravel the population dynamics in A. mellifera and provide fundamental information to conserve native honey bees.

Unlocking Horse Y Chromosome Diversity

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data

Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.

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.

Genetic Diversity and Structure of the Main Danubian Horse Paternal Genealogical Lineages Based on Microsatellite Genotyping

Simple Summary

The Danubian horse was created on the former Klementina stud farm near Pleven to satisfy the demands of the Bulgarian Army for light draft horses and to improve the working capacity of the local horse population. The privatization of the Klementina stud farm in the late 1990s and the lack of economic activity have led to a sharp reduction in the number of mares and stallions and their sale to private owners. At present, only six of the main paternal lines which participated in the creation of the Danubian horse breed are preserved: Zdravko, NONIUS XVII-30, Torpedo, Lider, Kalifa, and Hrabar. This is the first study on the genetic and population diversity of the Danubian horse paternal lines based on microsatellite markers (STRs). The results showed that the studied Danubian horse population was characterized by a high level of genetic diversity with a medium value of 0.84. The fixation index (F_ST) was 0.08 for all studied markers, which is indicative of the low genetic differentiation of the Danubian horse population. Our analysis also confirmed the low level of inbreeding and heterozygous deficiency among the animals selected from the six paternal lineages of the Danubian horse. The present research could be helpful for the development of breeding and conservation programs for the Danubian horse, as well as for making informed decisions on the management of paternal lines.

Abstract

The Danubian horse, together with the Pleven and the Eastern Bulgarian horse breeds, is one of the modern breeds in Bulgaria. The objective of this study was to compare the genetic structure and genetic diversity of six paternal genealogical lineages of the Danubian horse breed (Zdravko, NONIUS XVII-30, Torpedo, Lider, Kalifa, and Hrabar). In total, 166 individuals from the six genealogical lines were investigated, based on 15 STR markers (short tandem repeats, also known as microsatellites). In total, 184 alleles were found in the six populations, using 15 microsatellite loci. The mean number of alleles, the effective number of alleles, and the polymorphic information content (PIC) values per locus were 12.28, 9.48, and 0.73, respectively. In a comparison of the allelic diversity among sire lineages, the highest genetic diversity (Na) was observed in Lider and Kalifa (14.60 ± 0.21), while the lowest value of this parameter was observed in the Zdravko lineage 4.20 ± 0.35. The largest genetic diversity was found in loci HMS3 and HMS7, with 13 alleles, and the smallest polymorphism was noted for the locus ASB17, with 10 alleles. The level of observed heterozygosity was in the range of 0.65 ± 0.069 for the Zdravko lineage to 0.93 ± 0.01 for the Torpedo lineage. The expected heterozygosity level range was from 0.57 ± 0.048 to 0.91 ± 0.01 for all horse lineages. Structure analysis revealed three main gene pools in the study population. The first pool included the Zdravko lineage; the second had the NONIUS XVII-30, Torpedo, Lider, and Kalifa lineages; and the third defined the Hrabar lineage, which was significantly differentiated from the other genealogical lineages.

Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern

Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (H_e) and contemporary effective population size (cNe) were reported in Siciliano (H_e = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (H_e = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.

DPF3, A Putative Candidate Gene For Melanoma Etiopathogenesis in Gray Horses

Melanoma prevalence in gray horses reaches up to 50% and more. Several studies have documented a genetic melanoma predisposition which is referred to the 4.6 kb duplication in intron 6 of STX17 and its surrounding haplotype. However, the genetic background and mechanisms responsible for differences in etiopathogenesis of equine dermal melanomatosis still remain unknown. In the current study, we performed a genome wide association analysis in 141 Lipizzan horses and subsequently identified one candidate gene on chromosome 24 putatively involved in melanoma pathogenesis in gray horses. The associated SNP was located in the intronic region of DPF3, a gene which is involved in humans in cell growth, proliferation, apoptosis and motility of cancer cells. The replication study in 1210 horses from seven breeds demonstrated, that the G/G genotype of the DPF3 associated SNP exhibits putative melanoma suppression effects. As a conclusion DPF3 represents a candidate gene, which might play an essential role for gray horses coping with high genetic melanoma related tumor load.

Equine vitiligo-like depigmentation in grey horses is related to genes involved in immune response and tumor metastasis

Background

In horses, the autoimmune disease vitiligo is characterized by the loss of melanocytes and results in patchy depigmentation of the skin around the eyes, muzzle and the perianal region. Vitiligo-like depigmentation occurs predominantly in horses displaying the grey coat colour and is observed at a prevalence level of 26.0–67.0% in grey horses compared with only 0.8–3.5% in non-grey horses. While the polygenetic background of this complex disease is well documented in humans, the underlying candidate genes for this skin disorder in horses remain unknown. In this study we aim to perform a genome-wide association study (GWAS) for identifying putative candidate loci for vitiligo-like depigmentation in horses.

Methods

In the current study, we performed a GWAS analysis using high-density 670 k single nucleotide polymorphism (SNP) data from 152 Lipizzan and 104 Noriker horses, which were phenotyped for vitiligo-like depigmentation by visual inspection. After quality control 376,219 SNPs remained for analyses, the genome-wide Bonferroni corrected significance level was p < 1.33e-7.

Results

We identified seven candidate genes on four chromosomes (ECA1, ECA13, ECA17, ECA20) putatively involved in vitiligo pathogenesis in grey horses. The highlighted genes PHF11, SETDB2, CARHSP1 and LITAFD, are associated with the innate immune system, while the genes RCBTB1, LITAFD, NUBPL, PTP4A1, play a role in tumor suppression and metastasis. The antagonistic pathogenesis of vitiligo in relation to cancer specific enhanced cell motility and/or metastasis on typical melanoma predilection sites underlines a plausible involvement of RCBTB1, LITAFD, NUBPL, and PTP4A1.

Conclusions

The proposed candidate genes for equine vitiligo-like depigmentation, indicate an antagonistic relation between vitiligo and tumor metastasis in a horse population with higher incidence of melanoma. Further replication and expression studies should lead to a better understanding of this skin disorder in horses.

Genome-Wide Scanning for Signatures of Selection Revealed the Putative Genomic Regions and Candidate Genes Controlling Milk Composition and Coat Color Traits in Sahiwal Cattle

Background

In the evolutionary time scale, selection shapes the genetic variation and alters the architecture of genome in the organisms. Selection leaves detectable signatures at the genomic coordinates that provide clues about the protein-coding regions. Sahiwal is a valuable indicine cattle adapted to tropical environments with desirable milk attributes. Insights into the genomic regions under putative selection may reveal the molecular mechanisms affecting the quantitative and other important traits. To understand this, the present investigation was undertaken to explore signatures of selection in the genome of Sahiwal cattle using a medium-density genotyping INDUS chip.

Result

De-correlated composite of multiple selection signals (DCMS), which combines five different univariate statistics, was computed in the dataset to detect the signatures of selection in the Sahiwal genome. Gene annotations, Quantitative Trait Loci (QTL) enrichment, and functional analyses were carried out for the identification of significant genomic regions. A total of 117 genes were identified, which affect a number of important economic traits. The QTL enrichment analysis highlighted 14 significant [False Discovery Rate (FDR)-corrected p-value ≤ 0.05] regions on chromosomes BTA 1, 3, 6, 11, 20, and 21. The top three enriched QTLs were found on BTA 6, 20, and 23, which are associated with exterior, health, milk production, and reproduction traits. The present study on selection signatures revealed some key genes related with coat color (PDGFRA, KIT, and KDR), facial pigmentation (LEF), milk fat percent (MAP3K1, HADH, CYP2U1, and SGMS2), sperm membrane integrity (OSTC), lactation persistency (MRPS30, NNT, CCL28, HMGCS1, NIM1K, ZNF131, and CCDC152), milk yield (GHR and ZNF469), reproduction (NKX2-1 and DENND1A), and bovine tuberculosis susceptibility (RNF144B and PAPSS1). Further analysis of candidate gene prioritization identified four hub genes, viz., KIT, KDR, MAP3K1, and LEF, which play a role in coat color, facial pigmentation, and milk fat percentage in cattle. Gene enrichment analysis revealed significant Gene ontology (GO) terms related to breed-specific coat color and milk fat percent.

Conclusion

The key candidate genes and putative genomic regions associated with economic traits were identified in Sahiwal using single nucleotide polymorphism data and the DCMS method. It revealed selection for milk production, coat color, and adaptability to tropical climate. The knowledge about signatures of selection and candidate genes affecting phenotypes have provided a background information that can be further utilized to understand the underlying mechanism involved in these traits in Sahiwal cattle.

Association of inbreeding and regional equine leucocyte antigen homozygosity with the prevalence of insect bite hypersensitivity in Old Kladruber horse

Inbreeding depression is the reduction of performance caused by mating of close relatives. In livestock populations, inbreeding depression has been traditionally estimated by regression of phenotypes on pedigree inbreeding coefficients. This estimation can be improved by utilising genomic inbreeding coefficients. Here we estimate inbreeding depression for insect bite hypersensitivity (IBH) prevalence, the most common allergic horse disease worldwide, in Old Kladruber horse. In a deep pedigree with 3214 horses (187 genotyped), we used a generalised linear mixed model with IBH phenotype from 558 horses examined between 1996 and 2009 (1368 records). In addition to the classical pedigree information, we used the single-step approach that enabled joint use of pedigree and genomic information to estimate inbreeding depression overall genome and equine leucocyte antigen (ELA) class II region. Significant inbreeding depression was observed in all models fitting overall inbreeding coefficients (odds ratio between 1.018 and 1.074, P < 0.05) with the exception of Kalinowski's new inbreeding (P = 0.0516). The increase of ELA class II inbreeding was significantly associated with increased prevalence of IBH (odds ratio 1.018; P = 0.027). However, when fitted jointly with the overall inbreeding coefficient, the effect of ELA class II inbreeding was not significant (odds ratio 1.016; P = 0.062). Overall, the higher ELA class II and/or overall inbreeding (pedigree or genomic) was associated with increased prevalence of IBH in Old Kladruber horses. The single-step approach provides an efficient use of all the available pedigree, genomic, and phenotype information for estimation of overall and regional inbreeding effects.

Fine-scale estimation of inbreeding rates, runs of homozygosity and genome-wide heterozygosity levels in the Mangalarga Marchador horse breed

With the availability of high-density SNP panels and the establishment of approaches for characterizing homozygosity and heterozygosity sites, it is possible to access fine-scale information regarding genomes, providing more than just comparisons of different inbreeding coefficients. This is the first study that seeks to access such information for the Mangalarga Marchador (MM) horse breed on a genomic scale. To this end, we aimed to assess inbreeding levels using different coefficients, as well as to characterize homozygous and heterozygous runs in the population. Using Axiom ® Equine Genotyping Array-670k SNP (Thermo Fisher), 192 horses were genotyped. Our results showed different estimates: inbreeding from genomic coefficients (F_ROH ) = 0.16; pedigree-based (F_PED ) = 0.008; and a method based on excess homozygosity (F_HOM ) = 0.010. The correlations between the inbreeding coefficients were low to moderate, and some comparisons showed negative correlations, being practically null. In total, 85,295 runs of homozygosity (ROH) and 10,016 runs of heterozygosity (ROHet) were characterized for the 31 horse autosomal chromosomes. The class with the highest percentage of ROH was 0-2 Mbps, with 92.78% of the observations. In the ROHet results, only the 0-2 class presented observations, with chromosome 11 highlighted in a region with high genetic variability. Three regions from the ROHet analyses showed genes with known functions: tripartite motif-containing 37 (TRIM37), protein phosphatase, Mg²⁺ /Mn²⁺ dependent 1E (PPM1E) and carbonic anhydrase 10 (CA10). Therefore, our findings suggest moderate inbreeding, possibly attributed to breed formation, annulling possible recent inbreeding. Furthermore, regions with high variability in the MM genome were identified (ROHet), associated with the recent selection and important events in the development and performance of MM horses over generations.

The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events

Background

Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types.

Results

The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction.

Conclusions

Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00616-3.

Runs of homozygosity provide a genome landscape picture of inbreeding and genetic history of European autochthonous and commercial pig breeds

ROHs are long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROHs and their length are indicators of the level and origin of inbreeding. Frequent common ROHs within the same population define ROH islands and indicate hotspots of selection. In this work, we investigated ROHs in a total of 1131 pigs from 20 European local pig breeds and in three cosmopolitan breeds, genotyped with the GGP Porcine HD Genomic Profiler. plink software was used to identify ROHs. Size classes and genomic inbreeding parameters were evaluated. ROH islands were defined by evaluating different thresholds of homozygous SNP frequency. A functional overview of breed-specific ROH islands was obtained via over-representation analyses of GO biological processes. Mora Romagnola and Turopolje breeds had the largest proportions of genome covered with ROH (~1003 and ~955 Mb respectively), whereas Nero Siciliano and Sarda breeds had the lowest proportions (~207 and 247 Mb respectively). The highest proportion of long ROH (>16 Mb) was in Apulo-Calabrese, Mora Romagnola and Casertana. The largest number of ROH islands was identified in the Italian Landrace (n = 32), Cinta Senese (n = 26) and Lithuanian White Old Type (n = 22) breeds. Several ROH islands were in regions encompassing genes known to affect morphological traits. Comparative ROH structure analysis among breeds indicated the similar genetic structure of local breeds across Europe. This study contributed to understanding of the genetic history of the investigated pig breeds and provided information to manage these pig genetic resources.

Genetic Variability in Local and Imported Germplasm Chicken Populations as Revealed by Analyzing Runs of Homozygosity

Simple Summary

To maintain the uniqueness of conserved chicken populations of local and imported breeds is of great importance. In this study, we genotyped small populations belonging to 14 breeds and 7 crossbreds using an Illumina Chicken 60K SNP (Single Nucleotide Polymorphisms) BeadChip and looked for appropriate methods to characterize their purity/variability. It was not straightforward to identify crossbred individuals, and the best approach was based on calculating the length and number of homozygous regions, or runs of homozygosity (ROH), in the populations studied. The latter enabled most accurate identification of crossbreds and can be served as an effective tool in testing genome-wide purity of chicken breeds.

Abstract

Preserving breed uniqueness and purity is vitally important in developing conservation/breeding programs for a germplasm collection of rare and endangered chicken breeds. The present study was aimed at analyzing SNP genetic variability of 21 small local and imported purebred and F₁ crossbred populations and identifying crossbreeding events via whole-genome evaluation of runs of homozygosity (ROH). The admixture models more efficiently reflected population structure, pinpointing crossbreeding events in the presence of ancestral populations but not in their absence. Multidimensional scaling and F_ST-based analyses did not discriminate properly between purebred populations and F₁ crossbreds, especially when comparing related breeds. When applying the ROH-based approach, more and longer ROHs were revealed in purebred individuals/populations, suggesting this as an effective implement in genome-wide analysis of germplasm breed purity.

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.

Genetics of Skin Disease in Horses

Equine skin diseases are common, causing increased costs and reduced welfare of affected horses.Genetic testing, if available, can complement early detection, disease diagnosis, and clinical treatment and offers horse breeders the possibility to rule out carrier status. The mechanisms of complex disease can be investigated by using the latest state-of-the-art genomic technologies. Genome-based strategies may also serve as an efficient and cost-effective strategy for the management of the disease severity levels, with particular interest in complex traits such as insect bite hypersensitivity, chronic progressive lymphedema, and melanoma.

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.

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.

A Genome-Wide Association Analysis in Noriker Horses Identifies a SNP Associated With Roan Coat Color

The roan coat color in horses is characterized by dispersed white hair and dark points. This phenotype segregates in a broad range of horse breeds, while the underlying genetic background is still unknown. Previous studies mapped the roan locus to the KIT gene on equine chromosome 3 (ECA3). However, this association could not be validated across different horse breeds. Performing a genome-wide association analysis (GWAS) in Noriker horses, we identified a single nucleotide polymorphism (SNP) (ECA3:g.79,543.439 A > G) in the intron 17 of the KIT gene. The G -allele of the top associated SNP was present in other roan horses, namely Quarter Horse, Murgese, Slovenian, and Belgian draught horse, while it was absent in a panel of 15 breeds, including 657 non-roan horses. In further 379 gray Lipizzan horses, eight animals exhibited a heterozygous genotype (A/G). Comparative whole-genome sequence analysis of the KIT region revealed two deletions in the downstream region (ECA3:79,533,217_79,533,224delTCGTCTTC; ECA3:79,533,282_79,533,285delTTCT) and a 3 bp deletion combined with 17 bp insertion in intron 20 of KIT (ECA3:79,588,128_79,588,130delinsTTATCTCTATAGTAGTT). Within the Noriker sample, these loci were in complete linkage disequilibrium (LD) with the identified top SNP. Based upon pedigree information and historical records, we were able to trace back the genetic origin of roan coat color to a baroque gene pool. Furthermore, our data suggest allelic heterogeneity and the existence of additional roan alleles in ponies and breeds related to the English Thoroughbred. In order to study the roan phenotype segregating in those breeds, further association and verification studies are required.

Genomic Divergence in Swedish Warmblood Horses Selected for Equestrian Disciplines

The equestrian sport horse Swedish Warmblood (SWB) originates from versatile cavalry horses. Most modern SWB breeders have specialized their breeding either towards show jumping or dressage disciplines. The aim of this study was to explore the genomic structure of SWB horses to evaluate the presence of genomic subpopulations, and to search for signatures of selection in subgroups of SWB with high or low breeding values (EBVs) for show jumping. We analyzed high density genotype information from 380 SWB horses born in the period 2010-2011, and used Principal Coordinates Analysis and Discriminant Analysis of Principal Components to detect population stratification. Fixation index and Cross Population Extended Haplotype Homozygosity scores were used to scan the genome for potential signatures of selection. In accordance with current breeding practice, this study highlights the development of two separate breed subpopulations with putative signatures of selection in eleven chromosomes. These regions involve genes with known function in, e.g., mentality, endogenous reward system, development of connective tissues and muscles, motor control, body growth and development. This study shows genetic divergence, due to specialization towards different disciplines in SWB horses. This latter evidence can be of interest for SWB and other horse studbooks encountering specialized breeding.

Genome-Wide Homozygosity Patterns and Evidence for Selection in a Set of European and Near Eastern Horse Breeds

Intensive artificial and natural selection have shaped substantial variation among European horse breeds. Whereas most equine selection signature studies employ divergent genetic population structures in order to derive specific inter-breed targets of selection, we screened a total of 1476 horses originating from 12 breeds for the loss of genetic diversity by runs of homozygosity (ROH) utilizing a 670,000 single nucleotide polymorphism (SNP) genotyping array. Overlapping homozygous regions (ROH islands) indicating signatures of selection were identified by breed and similarities/dissimilarities between populations were evaluated. In the entire dataset, 180 ROH islands were identified, whilst 100 islands were breed specific, all other overlapped in 36 genomic regions with at least one ROH island of another breed. Furthermore, two ROH hot spots were determined at horse chromosome 3 (ECA3) and ECA11. Besides the confirmation of previously documented target genes involved in selection for coat color (MC1R, STX17, ASIP), body size (LCORL/NCAPG, ZFAT, LASP1, HMGA2), racing ability (PPARGC1A), behavioral traits (GRIN2B, NTM/OPCML) and gait patterns (DMRT3), several putative target genes related to embryonic morphogenesis (HOXB), energy metabolism (IGFBP-1, IGFBP-3), hair follicle morphogenesis (KRT25, KRT27, INTU) and autophagy (RALB) were highlighted. Furthermore, genes were pinpointed which might be involved in environmental adaptation of specific habitats (UVSSA, STXBP4, COX11, HLF, MMD).

Genome-Wide Association Studies Based on Equine Joint Angle Measurements Reveal New QTL Affecting the Conformation of Horses

The evaluation of conformation traits is an important part of selection for breeding stallions and mares. Some of these judged conformation traits involve joint angles that are associated with performance, health, and longevity. To improve our understanding of the genetic background of joint angles in horses, we have objectively measured the angles of the poll, elbow, carpal, fetlock (front and hind), hip, stifle, and hock joints based on one photograph of each of the 300 Franches-Montagnes (FM) and 224 Lipizzan (LIP) horses. After quality control, genome-wide association studies (GWASs) for these traits were performed on 495 horses, using 374,070 genome-wide single nucleotide polymorphisms (SNPs) in a mixed-effect model. We identified two significant quantitative trait loci (QTL) for the poll angle on ECA28 (p = 1.36 × 10⁻⁷), 50 kb downstream of the ALX1 gene, involved in cranial morphology, and for the elbow joint on ECA29 (p = 1.69 × 10⁻⁷), 49 kb downstream of the RSU1 gene, and 75 kb upstream of the PTER gene. Both genes are associated with bone mineral density in humans. Furthermore, we identified other suggestive QTL associated with the stifle joint on ECA8 (p = 3.10 × 10⁻⁷); the poll on ECA1 (p = 6.83 × 10⁻⁷); the fetlock joint of the hind limb on ECA27 (p = 5.42 × 10⁻⁷); and the carpal joint angle on ECA3 (p = 6.24 × 10⁻⁷), ECA4 (p = 6.07 × 10⁻⁷), and ECA7 (p = 8.83 × 10⁻⁷). The application of angular measurements in genetic studies may increase our understanding of the underlying genetic effects of important traits in equine breeding.