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. Modelbased 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 (F-R0H) 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.