Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (F_ROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (F_HOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (F_PED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. F_HOM was highly correlated (0.671 to 0.985) with F_ROH measures in all breeds. Apulo-Calabrese and Casertana had the highest F_ROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (F_ROH1) to 16 Mbp (F_ROH16)), whereas the lowest F_ROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). F_ROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all F_ROH measures and F_PED (from 0.514 to 0.523) and between F_HOM and F_PED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between F_ROH and F_PED. This might be due to the imperfect measure of F_PED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.     

Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy

Increased inbreeding is an inevitable consequence of selection in livestock populations. The analysis of high‐density single nucleotide polymorphisms (SNPs) facilitates the identification of long and uninterrupted runs of homozygosity (ROH) that can be used to identify chromosomal regions that are identical by descent. In this work, the distribution of ROH of different lengths in five Italian cattle breeds is described. A total of 4095 bulls from five cattle breeds (2093 Italian Holstein, 749 Italian Brown, 364 Piedmontese, 410 Marchigiana and 479 Italian Simmental) were genotyped at 54K SNP loci. ROH were identified and used to estimate molecular inbreeding coefficients (F_ROH), which were compared with inbreeding coefficients estimated from pedigree information (F_PED) and using the genomic relationship matrix (F_GRM). The average number of ROH per animal ranged from 54 ± 7.2 in Piedmontese to 94.6 ± 11.6 in Italian Brown. The highest number of short ROH (related to ancient consanguinity) was found in Piedmontese, followed by Simmental. The Italian Brown and Holstein had a higher proportion of longer ROH distributed across the whole genome, revealing recent inbreeding. The F_PED were moderately correlated with F_ROH > 1 Mb (0.662, 0.700 and 0.669 in Italian Brown, Italian Holstein and Italian Simmental respectively) but poorly correlated with F_GRM (0.134, 0.128 and 0.448 for Italian Brown, Italian Holstein and Italian Simmental respectively). The inclusion of ROH > 8 Mb in the inbreeding calculation improved the correlation of F_ROH with F_PED and F_GRM. ROH are a direct measure of autozygosity at the DNA level and can overcome approximations and errors resulting from incomplete pedigree data. In populations with high linkage disequilibrium (LD) and recent inbreeding (e.g. Italian Holstein and Italian Brown), a medium‐density marker panel, such as the one used here, may provide a good estimate of inbreeding. However, in populations with low LD and ancient inbreeding, marker density would have to be increased to identify short ROH that are identical by descent more precisely.     

Runs of homozygosity reveal genome‐wide autozygosity in Italian sheep breeds

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium‐density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (F_ROH = 0.099), whereas Comisana showed the lowest (F_ROH = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

An evaluation of inbreeding measures using a whole-genome sequenced cattle pedigree

The estimation of the inbreeding coefficient (F) is essential for the study of inbreeding depression (ID) or for the management of populations under conservation. Several methods have been proposed to estimate the realized F using genetic markers, but it remains unclear which one should be used. Here we used whole-genome sequence data for 245 individuals from a Holstein cattle pedigree to empirically evaluate which estimators best capture homozygosity at variants causing ID, such as rare deleterious alleles or loci presenting heterozygote advantage and segregating at intermediate frequency. Estimators relying on the correlation between uniting gametes (F_UNI) or on the genomic relationships (F_GRM) presented the highest correlations with these variants. However, homozygosity at rare alleles remained poorly captured. A second group of estimators relying on excess homozygosity (F_HOM), homozygous-by-descent segments (F_HBD), runs-of-homozygosity (F_ROH) or on the known genealogy (F_PED) was better at capturing whole-genome homozygosity, reflecting the consequences of inbreeding on all variants, and for young alleles with low to moderate frequencies (0.10 < . < 0.25). The results indicate that F_UNI and F_GRM might present a stronger association with ID. However, the situation might be different when recessive deleterious alleles reach higher frequencies, such as in populations with a small effective population size. For locus-specific inbreeding measures or at low marker density, the ranking of the methods can also change as F_HBD makes better use of the information from neighboring markers. Finally, we confirmed that genomic measures are in general superior to pedigree-based estimates. In particular, F_PED was uncorrelated with locus-specific homozygosity.Subject terms: Conservation genomics, Animal breeding, Inbreeding     

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of F_ST were higher in Wallachian sheep (F_ST = 0.14) than in Sumava sheep (F_ST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the N_e estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (F_HOM) among individuals, from ?0.04 to 0.16 in Sumava and from ?0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.     

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome‐wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.     

Low genome‐wide homozygosity in 11 Spanish ovine breeds

The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (F_ROH) are generally small (F_ROH < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.     

The footprint of recent and strong demographic decline in the genomes of Mangalitza pigs

The Mangalitza pig breed has suffered strong population reductions due to competition with more productive cosmopolitan breeds. In the current work, we aimed to investigate the effects of this sustained demographic recession on the genomic diversity of Mangalitza pigs. By using the Porcine Single Nucleotid Polymorphism BeadChip, we have characterized the genome-wide diversity of 350 individuals including 45 Red Mangalitza (number of samples; n=20 from Hungary and n=25 from Romania), 37 Blond Mangalitza, 26 Swallow-belly Mangalitza, 48 Blond Mangalitza × Duroc crossbreds, 5 Bazna swine, 143 pigs from the Hampshire, Duroc, Landrace, Large White and Pietrain breeds and 46 wild boars from Romania (n=18) and Hungary (n=28). Performance of a multidimensional scaling plot showed that Landrace, Large White and Pietrain pigs clustered independently from Mangalitza pigs and Romanian and Hungarian wild boars. The number and total length of ROH (runs of homozygosity), as well as F_ROH coefficients (proportion of the autosomal genome covered ROH) did not show major differences between Mangalitza pigs and other wild and domestic pig populations. However, Romanian and Hungarian Red Mangalitza pigs displayed an increased frequency of very long ROH (>30 Mb) when compared with other porcine breeds. These results indicate that Red Mangalitza pigs underwent recent and strong inbreeding probably as a consequence of severe reductions in census size.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

The pattern of runs of homozygosity and genomic inbreeding in world-wide sheep populations

Genome-wide pattern of runs of homozygosity (ROH) across ovine genome can provide a useful resource for studying diversity and demography history in sheep. We analyzed 50 k SNPs chip data of 2536 animals to identify pattern, distribution and level of ROHs in 68 global sheep populations. A total of 60,301 ROHs were detected in all breeds. The majority of the detected ROHs were <16 Mb and the average total number of ROHs per individual was 23.8 ± 13.8. The ROHs greater than 1 Mb covered on average 8.2% of the sheep autosomes, 1% of which was related to the ROHs with 1–4 Mb of length. The mean sum of ROH length in two-thirds of the populations was less than 250 Mb ranging from 21.7 to near 570 Mb. The level of genomic inbreeding was relatively low. The average of the inbreeding coefficients based on ROH (F_ROH) was 0.09 ± 0.05. It was rising in a stepwise manner with distance from Southwest Asia and maximum values were detected in North European breeds. A total of 465 ROH hotspots were detected in 25 different autosomes which partially surrounding 257 Refseq genes across the genome. Most of the detected genes were related to growth, body weight, meat production and quality, wool production and pigmentation. In conclusion, our analysis showed that the sheep genome, compared with other livestock species such as cattle and pig, displays low levels of homozygosity and appropriate genetic diversity for selection response and genetic merit gain.     

Preselection statistics and Random Forest classification identify population informative single nucleotide polymorphisms in cosmopolitan and autochthonous cattle breeds

Commercial single nucleotide polymorphism (SNP) arrays have been recently developed for several species and can be used to identify informative markers to differentiate breeds or populations for several downstream applications. To identify the most discriminating genetic markers among thousands of genotyped SNPs, a few statistical approaches have been proposed. In this work, we compared several methods of SNPs preselection (Delta, F_st and principal component analyses (PCA)) in addition to Random Forest classifications to analyse SNP data from six dairy cattle breeds, including cosmopolitan (Holstein, Brown and Simmental) and autochthonous Italian breeds raised in two different regions and subjected to limited or no breeding programmes (Cinisara, Modicana, raised only in Sicily and Reggiana, raised only in Emilia Romagna). From these classifications, two panels of 96 and 48 SNPs that contain the most discriminant SNPs were created for each preselection method. These panels were evaluated in terms of the ability to discriminate as a whole and breed-by-breed, as well as linkage disequilibrium within each panel. The obtained results showed that for the 48-SNP panel, the error rate increased mainly for autochthonous breeds, probably as a consequence of their admixed origin lower selection pressure and by ascertaining bias in the construction of the SNP chip. The 96-SNP panels were generally more able to discriminate all breeds. The panel derived by PCA-chrom (obtained by a preselection chromosome by chromosome) could identify informative SNPs that were particularly useful for the assignment of minor breeds that reached the lowest value of Out Of Bag error even in the Cinisara, whose value was quite high in all other panels. Moreover, this panel contained also the lowest number of SNPs in linkage disequilibrium. Several selected SNPs are located nearby genes affecting breed-specific phenotypic traits (coat colour and stature) or associated with production traits. In general, our results demonstrated the usefulness of Random Forest in combination to other reduction techniques to identify population informative SNPs.     

SNP genotypes reveal breed substructure,selection signatures and highly inbred regions in Piétrain pigs

The Piétrain pig originates from the Belgian village Piétrain some time between 1920 and 1950. Owing to its superior conformation, the Piétrain has spread worldwide since the 1960s. As initial population sizes were limited and close inbreeding was commonplace, the breed’s genetic diversity has been questioned. Therefore, this study examines Piétrain breed substructure, diversity and selection signatures using SNP data in comparison with Duroc, Landrace and Large White populations. Principal component analysis indicated three subpopulations, and F_ST analysis showed that US Piétrains differ most from European Piétrains. Average inbreeding based on runs of homozygosity (ROH) segments larger than 4 Mb ranged between 16.7 and 20.9%. The highest chromosomal inbreeding levels were found on SSC8 (42.7%). ROH islands were found on SSC8, SSC15 and SSC18 in all Piétrain populations, but numerous population-specific ROH islands were also detected. Moreover, a large ROH island on SSC8 (34–126 Mb) appears nearly fixed in all Piétrain populations, with a unique genotype. Chromosomal ROH patterns were similar between Piétrain populations. This study shows that Piétrain populations are genetically diverging, with at least three genetically distinct populations worldwide. Increasing genetic diversity in local Piétrain populations by introgression from other Piétrain populations seems to be only limited. Moreover, a unique 90 Mb region on SSC8 appeared largely fixed in the Piétrain breed, indicating that fixation was already present before the 1960s. We believe that strong selection and inbreeding during breed formation fixed these genomic regions in Piétrains. Finally, we hypothesize that independent coat color selection may have led to large ROH pattern similarities on SSC8 between unrelated pig breeds.     

Genome-wide analyses reveal population structure and identify candidate genes associated with tail fatness in local sheep from a semi-arid area

Under a climate change perspective, the genetic make-up of local livestock breeds showing adaptive traits should be explored and preserved as a priority. We used genotype data from the ovine 50 k Illumina BeadChip for assessing breed autozygosity based on runs of homozygosity (ROH) and fine-scale genetic structure and for detecting genomic regions under selection in 63 Tunisian sheep samples. The average genomic inbreeding coefficients based on ROH were estimated at 0.017, 0.021, and 0.024 for Barbarine (BAR, n = 26), Noire de Thibar (NDT, n = 23), and Queue fine de l'Ouest (QFO, n = 14) breeds, respectively. The genomic relationships among individuals based on identity by state (IBS) distance matrix highlighted a recent introgression of QFO into the BAR and a genetic differentiation of NDT samples, possibly explained by past introgression of European gene pools. Genome-wide scan for ROH across breeds and within the BAR sample set identified an outstanding signal on chromosome 13 (46.58–49.61 Mbp). These results were confirmed using F_ST index, differentiating fat vs. thin-tailed individuals. Candidate genes under selection pressure (CDS2, PROKR1, and BMP2) were associated to lipid storage and probably preferentially selected in fat-tailed BAR animals. Our findings suggest paying more attention to preserve the genetic integrity and adaptive alleles of local sheep breeds.     

Genomic characterization of Pinzgau cattle: genetic conservation and breeding perspectives

A genome-wide scan of Slovak Pinzgau cattle was prepared for the first time in order to estimate their genetic diversity at a more detailed level compared to previously published studies. The aim of this study was to describe the genetic diversity based on the runs of homozygosity (ROH_s), linkage disequilibrium (LD) and effective population size (N_eLD) using genome-wide data. Moreover, Bayesian clustering algorithms and multivariate methods were used to detect the population structure, potential admixture level and relationship between Austrian and Slovak Pinzgau cattle with respect to a large meta-population consisting of 15 European cattle breeds. The proportion of ROH segments ranged from 0.43 to 1.91% in Slovak Pinzgau, depending on the minimum size of an ROH. The genomic inbreeding coefficients were higher than the pedigree ones possibly due to the limited number of available generations in pedigree data. The observed N_eLD was close to the limit value characterizing the endangerment status, based both on genomic and pedigree data. Population structure within analyzed breeds based on the Wright’s F_ST index, Nei’s genetic distances, and unsupervised as well as supervised analysis has been established. Overall, these analyses clearly distinguished populations based on their origin. A detailed analysis of the introgression of each breed into the Pinzgau breeds prepared using a Bayesian approach showed that the contribution of Holstein cattle in Austrian as well as Slovak Pinzgau was larger than contribution of beef breeds. A possible reason is the recent usage of Holstein sires to increase milk production. There are considerable differences between well-defined regions that clearly distinguish Austrian and Slovak Pinzgau, despite their close common history. Generally, the breeding program of Austrian Pinzgau is more focused on meat production than Slovak Pinzgau, which was clearly reflected in the obtained autozygosity islands. Considering the genetic establishment of Slovak Pinzgau population the genetic potential of the breed is insufficiently used. On a long term, more global breeding program including very close populations will be more efficient providing higher genetic progress and diversity. Established methodology how to distinguish genealogically close populations on high-throughput molecular information based of Slovak and Austrian Pinzgau can be proposed as general for analysis of differences in all highly related breeds.     

Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees

Inbreeding (mating between relatives) can dramatically reduce the fitness of offspring by causing parts of the genome to be identical by descent. Thus, measuring individual inbreeding is crucial for ecology, evolution and conservation biology. We used computer simulations to test whether the realized proportion of the genome that is identical by descent (IBD_G) is predicted better by the pedigree inbreeding coefficient (F_P) or by genomic (marker-based) measures of inbreeding. Genomic estimators of IBD_G included the increase in individual homozygosity relative to mean Hardy–Weinberg expected homozygosity (F_H), and two measures (F_ROH and F_E) that use mapped genetic markers to estimate IBD_G. IBD_G was more strongly correlated with F_H, F_E and F_ROH than with F_P across a broad range of simulated scenarios when thousands of SNPs were used. For example, IBD_G was more strongly correlated with F_ROH, F_H and F_E (estimated with ⩾10 000 SNPs) than with F_P (estimated with 20 generations of complete pedigree) in populations with a recent reduction in the effective populations size (from N_e=500 to N_e=75). F_ROH, F_H and F_E generally explained >90% of the variance in IBD_G (among individuals) when 35 K or more SNPs were used. F_P explained <80% of the variation in IBD_G on average in all simulated scenarios, even when pedigrees included 20 generations. Our results demonstrate that IBD_G can be more precisely estimated with large numbers of genetic markers than with pedigrees. We encourage researchers to adopt genomic marker-based measures of IBD_G as thousands of loci can now be genotyped in any species.     

Effective Population Size,Extended Linkage Disequilibrium and Signatures of Selection in the Rare Dog Breed Lundehund

The Lundehund is an old dog breed with remarkable anatomical features including polydactyly in all four limbs and extraordinary flexibility of the spine. We genotyped 28 Lundehund using the canine Illumina high density beadchip to estimate the effective population size (N_e) and inbreeding coefficients as well as to identify potential regions of positive selection. The decay of linkage disequilibrium was slow with r² = 0.95 in 50 kb distance. The last 7-200 generations ago, Ne was at 10-13. An increase of N_e was noted in the very recent generations with a peak value of 19 for N_e at generation 4. The FROH estimated for 50-, 65- and 358-SNP windows were 0.87, 087 and 0.81, respectively. The most likely estimates for F_ROH after removing identical-by-state segments due to linkage disequilibria were at 0.80-0.81. The extreme loss of heterozygosity has been accumulated through continued inbreeding over 200 generations within a probably closed population with a small effective population size. The mean inbreeding coefficient based on pedigree data for the last 11 generations (F_Ped = 0.10) was strongly biased downwards due to the unknown coancestry of the founders in this pedigree data. The long-range haplotype test identified regions with genes involved in processes of immunity, olfaction, woundhealing and neuronal development as potential targets of selection. The genes QSOX2, BMPR1B and PRRX2 as well as MYOM1 are candidates for selection on the Lundehund characteristics small body size, increased number of digits per paw and extraordinary mobility, respectively.     

Evaluation of artificial selection in Standard Poodles using whole-genome sequencing

Identifying regions of artificial selection within dog breeds may provide insights into genetic variation that underlies breed-specific traits or diseases—particularly if these traits or disease predispositions are fixed within a breed. In this study, we searched for runs of homozygosity (ROH) and calculated the d _i statistic (which is based upon F _ST) to identify regions of artificial selection in Standard Poodles using high-coverage, whole-genome sequencing data of 15 Standard Poodles and 49 dogs across seven other breeds. We identified consensus ROH regions ≥1 Mb in length and common to at least ten Standard Poodles covering 0.6 % of the genome, and d _i regions that most distinguish Standard Poodles from other breeds covering 3.7 % of the genome. Within these regions, we identified enriched gene pathways related to olfaction, digestion, and taste, as well as pathways related to adrenal hormone biosynthesis, T cell function, and protein ubiquitination that could contribute to the pathogenesis of some Poodle-prevalent autoimmune diseases. We also validated variants related to hair coat and skull morphology that have previously been identified as being under selective pressure in Poodles, and flagged additional polymorphisms in genes such as ITGA2B, CBX4, and TNXB that may represent strong candidates for other common Poodle disorders.