Population structure and genomic diversity of the Einsiedler horse.

The breeding history of the Einsiedler horse is closely connected with the Benedictine cloister Einsiedeln. In the mid-nineteenth century, it was decided to use European Warmblood stallions for cross-breeding and to abandon the selection of stallions. Since that time, it has only been possible to trace back the origin of Einsiedler horses using maternal ancestry information. Here, we collected high-density genotype data for European Warmblood horses (Selle Français, Swiss Warmblood and Einsiedler) and Franches-Montagnes horses, the last native Swiss horse breed, to unravel the current population structure of the Einsiedler horse. Using commonly applied methods to ascertain fine-scale population structures, it was not possible to clearly differentiate the Einsiedler from other European Warmblood horses. However, by means of runs of homozygosity (ROH) we were able to detect breed-specific ROH islands for the Einsiedler horse, including genes involved in domestication and adaptation to high altitude. Therefore, future breeding activities should involve the screening of these breed-specific ROH segments, the revival of cryopreserved sperm and the selection of Einsiedler stallions.

Increased heterozygosity in low-pass sequencing data allows identification of blood chimeras in cattle.

In about 90% of multiple pregnancies in cattle, shared blood circulation between fetuses leads to genetic chimerism in peripheral blood and can reduce reproductive performance in heterosexual co-twins. However, the early detection of heterosexual chimeras requires specialized tests. Here, we used low-pass sequencing data with a median coverage of 0.64× generated from blood samples of 322 F1 crosses between beef and dairy cattle and identified 20 putative blood chimeras through increased levels of genome-wide heterozygosity. In contrast, for 77 samples with routine SNP microarray data generated from hair bulbs of the same F1s, we found no evidence of chimerism, simultaneously observing high levels of genotype discordance with sequencing data. Fifteen out of 18 reported twins showed signs of blood chimerism, in line with previous reports, whereas the presence of five alleged singletons with strong signs of chimerism suggests that the in-utero death rate of co-twins is at the upper limit of former estimates. Together, our results show that low-pass sequencing data allow reliable screening for blood chimeras. They further affirm that blood is not recommended as a source of DNA for the detection of germline variants.

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.

Determining Objective Parameters to Assess Gait Quality in Franches-Montagnes Horses for Ground Coverage and Over-Tracking - Part 2: At Trot.

In gait quality assessments of horses, stride length (SL) is visually associated with spectacular movements of the front limbs, and described as ground coverage, while the movement of the hind limb under the body is supposedly essential to a longer over-tracking distance (OTD). To identify movement patterns with strong associations to SL and OTD, limb and body kinematics of 24 Franches-Montagnes (FM) stallions were measured with 3D optical motion capture (OMC) on a treadmill during an incremental speed test at trot (3.3-6.5 m/s). These measurements were correlated to the scores of ground coverage and over-tracking from six breeding experts. The amount of explained variance of parameters on SL and OTD were estimated using linear mixed-effect models in two models: a full model with all parameters measurable with OMC, and a reduced model with a subset of parameters measurable with inertial measurement units (IMUs). The front limb stance duration (16%) and OTD (7%) measured with OMC, or the OMC parameters front limb stance duration (24%) and suspension duration (14%) measurable with IMUs explained most variance in SL. However, four of six breeding experts were also significantly correlated (r>|0.41|) to front limb protraction angle. OTD variance was explained with OMC parameters suspension duration (10%) and hind limb contralateral pro-retraction angles (9%) or IMU-measurable parameters suspension duration (20%) and maximal pelvis pitch (5%). Four experts' scores for over-tracking were correlated to suspension duration. These results underscore the need for precise definitions of gait quality traits.

Estimates of Genetic Parameters for Shape Space Data in Franches-Montagnes Horses.

Conformation traits such as joint angles are important selection criteria in equine breeding, but mainly consist of subjective evaluation scores given by breeding judges, showing limited variation. The horse shape space model extracts shape data from 246 landmarks (LM) and objective joint angle measurements from triplets of LM on standardized horse photographs. The heritability was estimated for 10 joint angles (seven were measured twice with different LM placements), and relative warp components of the whole shape, in 608 Franches-Montagnes (FM) horses (480 stallions, 68 mares and 60 geldings born 1940-2018, 3-25 years old). The pedigree data comprised 6986 horses. Genetic variances and covariances were estimated by restricted maximum likelihood model (REML), including the fixed effects birth year, age (linear and quadratic), height at withers (linear and quadratic), as well as postural effects (head, neck, limb position and body alignment), together with a random additive genetic animal component and the residual effect. Estimated heritability varied from 0.08 (stifle joint) to 0.37 (poll). For the shape, the type was most heritable (0.36 to 0.37) and evolved from heavy to light over time. Image-based phenotyping can improve the selection of horses for conformation traits with moderate heritability (e.g., poll, shoulder and fetlock).

Determining Objective Parameters to Assess Gait Quality in Franches-Montagnes Horses for Ground Coverage and Over-Tracking - Part 1: At Walk.

Ground coverage and over-tracking are two gait quality traits describing the forward movement of the front respectively the hind limbs in relation to stride length and over-tracking distance. To investigate the complex interplay of different movement patterns in ground coverage and over-tracking, limb and body kinematics of 24 Franches-Montagnes (FM) stallions were measured with 3D optical motion capture (OMC) on a treadmill during an incremental speed test at the walk (1.4-2.0 m/s). The significance and amount of explained variance of kinematic parameters on stride length and over-tracking distance were estimated using linear mixed-effect models, with speed and horse as random effects. Two separate models were tested: a full model with all parameters measurable by OMC, and a reduced model with a subset of parameters also measurable with inertial measurement units (IMUs). The kinematic parameters were correlated to the subjective scores from six breeding experts to interpret their external validity. The parameter for ground coverage at the walk, explaining most of the variance in stride length, were the maximal forelimb retraction angle (11%) measured with OMC, and the range of pelvis pitch (10%) if measuring with IMUs. The latter was also the most relevant for quantifying over-tracking, explaining 24% to 33% of the variance in the over-tracking distance. The scores from most breeding experts were significantly correlated (r ≥ |0.41|) with the fore- and hind limb protraction angles, which reflect the textual definition of ground coverage and over-tracking. Both gait quality traits can be objectively quantified using either OMC or IMUs.

Gait quality scoring data of Franches-Montagnes stallions at walk and trot on a treadmill by experts of the breed and their reliability.

This article presents the data obtained from the scoring of 24 stallions of the Franches-Montagnes (FM) horse breed by six experts of this breed. The experts scored six traits at walk and eight at trot from the video recordings of these stallions walking and trotting on a treadmill during an incremental speed test. The scores were given on a scale of one to nine. All experts scored the same videos twice (two scoring tests) with a time interval of two years, and without feedback from the first scoring. Video sequences were presented in a different order between first and second scoring. The inter- and intra-rater reliability of the data was assessed using intraclass correlation coefficients (ICC) to evaluate its quality.

Decreased Mite Reproduction to Select Varroa destructor (Acari: Varroidae) Resistant Honey Bees (Hymenoptera: Apidae): Limitations and Potential Methodological Improvements.

The invasive parasitic mite, Varroa destructor (Anderson and Trueman), is the major biotic threat to the survival of European honey bees, Apis mellifera L. To improve colony survival against V. destructor, the selection of resistant lineages against this parasite is considered a sustainable solution. Among selected traits, mite fertility and fecundity, often referred to as suppressed mite reproduction are increasingly used in breeding programmes. However, the current literature leaves some gaps in the assessment of the effectiveness of selecting these traits toward achieving resistance. In the population studied here, we show a low repeatability and reproducibility of mite fertility and fecundity phenotypes, as well as a low correlation of these traits with infestation rates of colonies. Phenotyping reliability could neither be improved by increasing the number of worker brood cells screened, nor by screening drone brood, which is highly attractive for the parasite and available early in the season, theoretically allowing a reduction of generation time and thus an acceleration of genetic progress in selected lineages. Our results provide an evaluation of the potential and limitations of selecting on decreased mite reproduction traits to obtain V. destructor-resistant honeybee colonies. To allow for a more precise implementation of such selection and output reporting, we propose a refined nomenclature by introducing the terms of decreased mite reproduction and reduced mite reproduction, depending on the extent of mite reproduction targeted. We also highlight the importance of ensuring accurate phenotyping ahead of initiating long-lasting selection programmes.

A Novel QTL and a Candidate Gene Are Associated with the Progressive Motility of Franches-Montagnes Stallion Spermatozoa after Thaw.

The use of frozen-thawed semen is an important reproduction tool to preserve the biodiversity of small, native horse breeds such as the Franches-Montagnes (FM). However, not all stallions produce cryotolerant semen with a progressive motility after thaw ≥ 35%. To improve our understanding of the genetic background of male fertility traits in both fresh and frozen-thawed semen, we performed genome-wide association studies (GWAS) on gel-free volume, sperm cell concentration, total sperm count, and progressive motility in fresh and frozen-thawed semen from 109 FM stallions using 335,494 genome-wide single nucleotide polymorphisms (SNPs). We identified one significant (p < 1.69 × 10-7) quantitative trait locus (QTL) on ECA6 within the SCN8A gene for progressive motility after thaw, which was previously associated with progressive motility in boars. Homozygous stallions showed a substantial drop in progressive motility after thaw. This QTL could be used to identify cryointolerant stallions, avoiding the costly cryopreservation process. Further studies are needed to confirm whether this QTL is also present in other horse breeds.

Exploring Two Honey Bee Traits for Improving Resistance Against Varroa destructor: Development and Genetic Evaluation.

For the development of novel selection traits in honey bees, applicability under field conditions is crucial. We thus evaluated two novel traits intended to provide resistance against the ectoparasitic mite Varroa destructor and to allow for their straightforward implementation in honey bee selection. These traits are new field estimates of already-described colony traits: brood recapping rate ('Recapping') and solidness ('Solidness'). 'Recapping' refers to a specific worker characteristic wherein they reseal a capped and partly opened cell containing a pupa, whilst 'Solidness' assesses the percentage of capped brood in a predefined area. According to the literature and beekeepers' experiences, a higher recapping rate and higher solidness could be related to resistance to V. destructor. During a four-year field trial in Switzerland, the two resistance traits were assessed in a total of 121 colonies of Apis mellifera mellifera. We estimated the repeatability and the heritability of the two traits and determined their phenotypic correlations with commonly applied selection traits, including other putative resistance traits. Both traits showed low repeatability between different measurements within each year. 'Recapping' had a low heritability (h2 = 0.04 to 0.05, depending on the selected model) and a negative phenotypic correlation to non-removal of pin-killed brood (r = -0.23). The heritability of 'Solidness' was moderate (h2 = 0.24 to 0.25) and did not significantly correlate with resistance traits. The two traits did not show an association with V. destructor infestation levels. Further research is needed to confirm the results, as only a small number of colonies was evaluated.

Advances and perspectives in selecting resistance traits against the parasitic mite Varroa destructor in honey bees.

BACKGROUND: In spite of the implementation of control strategies in honey bee (Apis mellifera) keeping, the invasive parasitic mite Varroa destructor remains one of the main causes of colony losses in numerous countries. Therefore, this parasite represents a serious threat to beekeeping and agro-ecosystems that benefit from the pollination services provided by honey bees. To maintain their stocks, beekeepers have to treat their colonies with acaricides every year. Selecting lineages that are resistant to infestations is deemed to be a more sustainable approach. REVIEW: Over the last three decades, numerous selection programs have been initiated to improve the host-parasite relationship and to support honey bee survival in the presence of the parasite without the need for acaricide treatments. Although resistance traits have been included in the selection strategy of honey bees, it has not been possible to globally solve the V. destructor problem. In this study, we review the literature on the reasons that have potentially limited the success of such selection programs. We compile the available information to assess the relevance of selected traits and the potential environmental effects that distort trait expression and colony survival. Limitations to the implementation of these traits in the field are also discussed. CONCLUSIONS: Improving our knowledge of the mechanisms underlying resistance to V. destructor to increase trait relevance, optimizing selection programs to reduce environmental effects, and communicating selection outcomes are all crucial to efforts aiming at establishing a balanced relationship between the invasive parasite and its new host.

Should We Agree to Disagree? An Evaluation of the Inter-Rater Reliability of Gait Quality Traits in Franches-Montagnes Stallions.

Gait quality, that is, the way horses move according to functional and aesthetic principles, englobes many traits that are scored by experts during breeding competitions. The experts can score a trait on a subjective valuating (SV) scale or on a linear profiling (LP) scale representing the biological extremes of the population. However, the reliability of the appraisal of gait quality traits has not been extensively evaluated. In this study, seven breed experts appraised the walk and trot quality of 24 Franches-Montagnes stallions presented in hand on a sand track. Inter-rater reliabilities of six traits (five SV and one LP) at the walk and eight traits (five SV and three LP) at the trot were estimated with intraclass correlation coefficients (ICCs). The inter-rater reliabilities were poor (ICC < 0.50). The scale anchoring varied between experts, and the variance of scores was low. There were no systematic differences in inter-rater reliability between LP and SV traits. Future studies should determine whether the inter-rater reliabilities may be increased by a more precise definition of the scores within each trait to improve the absolute agreement between experts, by a more uniform scale anchoring between experts, and by decreasing the number of scale items. However, considering the inherent limitations of the human eye in observing high-speed movement, the use of a field-applicable kinematic measurement system may support breeding experts in the appraisal of gait quality traits in the future.

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.

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-7), 50 kb downstream of the ALX1 gene, involved in cranial morphology, and for the elbow joint on ECA29 (p = 1.69 × 10-7), 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-7); the poll on ECA1 (p = 6.83 × 10-7); the fetlock joint of the hind limb on ECA27 (p = 5.42 × 10-7); and the carpal joint angle on ECA3 (p = 6.24 × 10-7), ECA4 (p = 6.07 × 10-7), and ECA7 (p = 8.83 × 10-7). The application of angular measurements in genetic studies may increase our understanding of the underlying genetic effects of important traits in equine breeding.

The horse Y chromosome as an informative marker for tracing sire lines.

Analysis of the Y chromosome is the best-established way to reconstruct paternal family history in humans. Here, we applied fine-scaled Y-chromosomal haplotyping in horses with biallelic markers and demonstrate the potential of our approach to address the ancestry of sire lines. We de novo assembled a draft reference of the male-specific region of the Y chromosome from Illumina short reads and then screened 5.8 million basepairs for variants in 130 specimens from intensively selected and rural breeds and nine Przewalski's horses. Among domestic horses we confirmed the predominance of a young'crown haplogroup' in Central European and North American breeds. Within the crown, we distinguished 58 haplotypes based on 211 variants, forming three major haplogroups. In addition to two previously characterised haplogroups, one observed in Arabian/Coldblooded and the other in Turkoman/Thoroughbred horses, we uncovered a third haplogroup containing Iberian lines and a North African Barb Horse. In a genealogical showcase, we distinguished the patrilines of the three English Thoroughbred founder stallions and resolved a historic controversy over the parentage of the horse 'Galopin', born in 1872. We observed two nearly instantaneous radiations in the history of Central and Northern European Y-chromosomal lineages that both occurred after domestication 5,500 years ago.

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

BACKGROUND: The sample ascertainment bias due to complex population structures remains a major challenge in genome-wide investigations of complex traits. In this study we derived the high-resolution population structure and levels of autozygosity of 377 Lipizzan horses originating from five different European stud farms utilizing the SNP genotype information of the high density 700 k Affymetrix Axiom™ Equine genotyping array. Scanning the genome for overlapping runs of homozygosity (ROH) shared by more than 50% of horses, we identified homozygous regions (ROH islands) in order to investigate the gene content of those candidate regions by gene ontology and enrichment analyses. RESULTS: The high-resolution population network approach revealed well-defined substructures according to the origin of the horses (Austria, Slovakia, Croatia and Hungary). The highest mean genome coverage of ROH (SROH) was identified in the Austrian (SROH = 342.9), followed by Croatian (SROH = 214.7), Slovakian (SROH = 205.1) and Hungarian (SROH = 171.5) subpopulations. ROH island analysis revealed five common islands on ECA11 and ECA14, hereby confirming a closer genetic relationship between the Hungarian and Croatian as well as between the Austrian and Slovakian samples. Private islands were detected for the Hungarian and the Austrian Lipizzan subpopulations. All subpopulations shared a homozygous region on ECA11, nearly identical in position and length containing among other genes the homeobox-B cluster, which was also significantly (p < 0.001) highlighted by enrichment analysis. Gene ontology terms were mostly related to biological processes involved in embryonic morphogenesis and anterior/posterior specification. Around the STX17 gene (causative for greying), we identified a ROH island harbouring the genes NR4A3, STX17, ERP44 and INVS. Within further islands on ECA14, ECA16 and ECA20 we detected the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, whereas SPRY4 and HSP90AB1 are involved in melanoma metastasis and survival rate of melanoma patients in humans. CONCLUSIONS: We demonstrated that the assessment of high-resolution population structures within one single breed supports the downstream genetic analyses (e.g. the identification of ROH islands). By means of ROH island analyses, we identified the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, which might play an important role for further studies on equine melanoma. Furthermore, our results highlighted the impact of the homeobox-A and B cluster involved in morphogenesis of Lipizzan horses.

Repeatability, reproducibility and consistency of horse shape data and its association with linearly described conformation traits in Franches-Montagnes stallions.

Linear description (LD) of conformation traits was introduced in horse breeding to minimise subjectivity in scoring. However, recent studies have shown that LD traits show essentially the same problems as traditionally scored traits, such as data converging around the mean value with very small standard deviations. To improve the assessment of conformation traits of horses, we investigated the application of the recently described horse shape space model based upon 403 digitised photographs of 243 Franches-Montagnes (FM) stallions and extracted joint angles based on specific landmark triplets. Repeatability, reproducibility and consistency of the resulting shape data and joint angles were assessed with Procrustes ANOVA (Rep) and intra-class correlation coefficients (ICC). Furthermore, we developed a subjective score to classify the posture of the horses on each photograph. We derived relative warp scores (PCs) based upon the digitised photos conducting a principal component analysis (PCA). The PCs of the shapes and joint angles were compared to the posture scores and to the linear description data using linear mixed effect models including significant posture scores as random factors. The digitisation process was highly repeatable and reproducible for the shape (Rep = 0.72-0.99, ICC = 0.99). The consistency of the shape was limited by the age and posture (p < 0.05). The angle measurements were highly repeatable within one digitiser. Between digitisers, we found a higher variability of ICC values (ICC = 0.054-0.92), indicating digitising error in specific landmarks (e.g. shoulder point). The posture scores were highly repeatable (Fleiss' kappa = 0.713-0.857). We identified significant associations (p(X2) < 0.05) with traits describing the withers height, shoulder length and incline, overall leg conformation, walk and trot step length. The horse shape data and angles provide additional information to explore the morphology of horses and therefore can be applied to improve the knowledge of the genetic architecture of LD traits.

Population Networks Associated with Runs of Homozygosity Reveal New Insights into the Breeding History of the Haflinger Horse.

Within the scope of current genetic diversity analyses, population structure and homozygosity measures are independently analyzed and interpreted. To enhance analytical power, we combined the visualization of recently described high-resolution population networks with runs of homozygosity (ROH). In this study, we demonstrate that this approach enabled us to reveal important aspects of the breeding history of the Haflinger horse. We collected high-density genotype information of 531 horses originating from 7 populations which were involved in the formation of the Haflinger, namely 32 Italian Haflingers, 78 Austrian Haflingers, 190 Noriker, 23 Bosnian Mountain Horses, 20 Gidran, 33 Shagya Arabians, and 155 Purebred Arabians. Model-based cluster analysis identified substructures within Purebred Arabian, Haflinger, and Noriker that reflected distinct genealogy (Purebred Arabian), geographic origin (Haflinger), and coat color patterns (Noriker). Analysis of ROH revealed that the 2 Arabian populations (Purebred and Shagya Arabians), Gidran and the Bosnian Mountain Horse had the highest genome proportion covered by ROH segments (306-397 Mb). The Noriker and the Austrian Haflinger showed the lowest ROH coverage (228, 282 Mb). Our combined visualization approach made it feasible to clearly identify outbred (admixture) and inbred (ROH segments) horses. Genomic inbreeding coefficients (FROH) ranged from 10.1% (Noriker) to 17.7% (Purebred Arabian). Finally it could be demonstrated, that the Austrian Haflinger sample has a lack of longer ROH segments and a deviating ROH spectrum, which is associated with past bottleneck events and the recent mating strategy favoring out-crosses within the breed.

Population structure and genomic inbreeding in nine Swiss dairy cattle populations.

BACKGROUND: Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF). RESULTS: The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26. CONCLUSIONS: The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.