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.     

Genomic inbreeding estimation in small populations: evaluation of runs of homozygosity in three local dairy cattle breeds

In the local breeds with small population size, one of the most important problems is the increase of inbreeding coefficient (F). High levels of inbreeding lead to reduced genetic diversity and inbreeding depression. The availability of high-density single nucleotide polymorphism (SNP) arrays has facilitated the quantification of F by genomic markers in farm animals. Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes and represent an estimate of the degree of autozygosity at genome-wide level. The current study aims to quantify the genomic F derived from ROH (F_ROH) in three local dairy cattle breeds. F_ROH values were compared with F estimated from the genomic relationship matrix (F_GRM), based on the difference between observed v. expected number of homozygous genotypes (F_HOM) and the genomic homozygosity of individual i (F_{MOL i}). The molecular coancestry coefficient (f_{MOL ij}) between individuals i and j was also estimated. Individuals of Cinisara (71), Modicana (72) and Reggiana (168) were genotyped with the 50K v2 Illumina BeadChip. Genotypes from 96 animals of Italian Holstein cattle breed were also included in the analysis. We used a definition of ROH as tracts of homozygous genotypes that were >4 Mb. Among breeds, 3661 ROH were identified. Modicana showed the highest mean number of ROH per individual and the highest value of F_ROH, whereas Reggiana showed the lowest ones. Differences among breeds existed for the ROH lengths. The individuals of Italian Holstein showed high number of short ROH segments, related to ancient consanguinity. Similar results showed the Reggiana with some extreme animals with segments covering 400 Mb and more of genome. Modicana and Cinisara showed similar results between them with the total length of ROH characterized by the presence of large segments. High correlation was found between F_HOM and F_ROH ranged from 0.83 in Reggiana to 0.95 in Cinisara and Modicana. The correlations among F_ROH and other estimated F coefficients were generally lower ranged from 0.45 (F_{MOL i}−F_ROH) in Cinisara to 0.17 (F_GRM−F_ROH) in Modicana. On the basis of our results, recent inbreeding was observed in local breeds, considering that 16 Mb segments are expected to present inbreeding up to three generations ago. Our results showed the necessity of implementing conservation programs to control the rise of inbreeding and coancestry in the three Italian local dairy cattle breeds.     

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.     

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.     

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.     

Long Runs of Homozygosity Are Enriched for Deleterious Variation

Exome sequencing offers the potential to study the population-genomic variables that underlie patterns of deleterious variation. Runs of homozygosity (ROH) are long stretches of consecutive homozygous genotypes probably reflecting segments shared identically by descent as the result of processes such as consanguinity, population size reduction, and natural selection. The relationship between ROH and patterns of predicted deleterious variation can provide insight into the way in which these processes contribute to the maintenance of deleterious variants. Here, we use exome sequencing to examine ROH in relation to the distribution of deleterious variation in 27 individuals of varying levels of apparent inbreeding from 6 human populations. A significantly greater fraction of all genome-wide predicted damaging homozygotes fall in ROH than would be expected from the corresponding fraction of nondamaging homozygotes in ROH (p < 0.001). This pattern is strongest for long ROH (p < 0.05). ROH, and especially long ROH, harbor disproportionately more deleterious homozygotes than would be expected on the basis of the total ROH coverage of the genome and the genomic distribution of nondamaging homozygotes. The results accord with a hypothesis that recent inbreeding, which generates long ROH, enables rare deleterious variants to exist in homozygous form. Thus, just as inbreeding can elevate the occurrence of rare recessive diseases that represent homozygotes for strongly deleterious mutations, inbreeding magnifies the occurrence of mildly deleterious variants as well.     

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     

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.     

Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Identification of genomic regions associated with inbreeding depression in Holstein and Jersey dairy cattle

Background

Inbreeding reduces the fitness of individuals by increasing the frequency of homozygous deleterious recessive alleles. Some insight into the genetic architecture of fitness, and other complex traits, can be gained by using single nucleotide polymorphism (SNP) data to identify regions of the genome which lead to reduction in performance when identical by descent (IBD). Here, we compared the effect of genome-wide and location-specific homozygosity on fertility and milk production traits in dairy cattle.

Methods

Genotype data from more than 43 000 SNPs were available for 8853 Holstein and 4138 Jersey dairy cows that were part of a much larger dataset that had pedigree records (338 696 Holstein and 64 049 Jersey animals). Measures of inbreeding were based on: (1) pedigree data; (2) genotypes to determine the realised proportion of the genome that is IBD; (3) the proportion of the total genome that is homozygous and (4) runs of homozygosity (ROH) which are stretches of the genome that are homozygous.

Results

A 1% increase in inbreeding based either on pedigree or genomic data was associated with a decrease in milk, fat and protein yields of around 0.4 to 0.6% of the phenotypic mean, and an increase in calving interval (i.e. a deterioration in fertility) of 0.02 to 0.05% of the phenotypic mean. A genome-wide association study using ROH of more than 50 SNPs revealed genomic regions that resulted in depression of up to 12.5 d and 260 L for calving interval and milk yield, respectively, when completely homozygous.

Conclusions

Genomic measures can be used instead of pedigree-based inbreeding to estimate inbreeding depression. Both the diagonal elements of the genomic relationship matrix and the proportion of homozygous SNPs can be used to measure inbreeding. Longer ROH (>3 Mb) were found to be associated with a reduction in milk yield and captured recent inbreeding independently and in addition to overall homozygosity. Inbreeding depression can be reduced by minimizing overall inbreeding but maybe also by avoiding the production of offspring that are homozygous for deleterious alleles at specific genomic regions that are associated with inbreeding depression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0071-7) contains supplementary material, which is available to authorized users.     

Runs of homozygosity and distribution of functional variants in the cattle genome

Background

Recent developments in sequencing technology have facilitated widespread investigations of genomic variants, including continuous stretches of homozygous genomic regions. For cattle, a large proportion of these runs of homozygosity (ROH) are likely the result of inbreeding due to the accumulation of elite alleles from long-term selective breeding programs. In the present study, ROH were characterized in four cattle breeds with whole genome sequence data and the distribution of predicted functional variants was detected in ROH regions and across different ROH length classes.

Results

On average, 19.5 % of the genome was located in ROH across four cattle breeds. There were an average of 715.5 ROH per genome with an average size of ~750 kbp, ranging from 10 (minimum size considered) to 49,290 kbp. There was a significant correlation between shared short ROH regions and regions putatively under selection (p < 0.001). By investigating the relationship between ROH and the predicted deleterious and non-deleterious variants, we gained insight into the distribution of functional variation in inbred (ROH) regions. Predicted deleterious variants were more enriched in ROH regions than predicted non-deleterious variants, which is consistent with observations in the human genome. We also found that increased enrichment of deleterious variants was significantly higher in short (<100 kbp) and medium (0.1 to 3 Mbp) ROH regions compared with long (>3 Mbp) ROH regions (P < 0.001), which is different than what has been observed in the human genome.

Conclusions

This study illustrates the distribution of ROH and functional variants within ROH in cattle populations. These patterns are different from those in the human genome but consistent with the natural history of cattle populations, which is confirmed by the significant correlation between shared short ROH regions and regions putatively under selection. These findings contribute to understanding the effects of inbreeding and probably selection in shaping the distribution of functional variants in the cattle genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1715-x) contains supplementary material, which is available to authorized users.     

Genomic patterns of homozygosity in worldwide human populations

Genome-wide patterns of homozygosity runs and their variation across individuals provide a valuable and often untapped resource for studying human genetic diversity and evolutionary history. Using genotype data at 577,489 autosomal SNPs, we employed a likelihood-based approach to identify runs of homozygosity (ROH) in 1,839 individuals representing 64 worldwide populations, classifying them by length into three classes—short, intermediate, and long—with a model-based clustering algorithm. For each class, the number and total length of ROH per individual show considerable variation across individuals and populations. The total lengths of short and intermediate ROH per individual increase with the distance of a population from East Africa, in agreement with similar patterns previously observed for locus-wise homozygosity and linkage disequilibrium. By contrast, total lengths of long ROH show large interindividual variations that probably reflect recent inbreeding patterns, with higher values occurring more often in populations with known high frequencies of consanguineous unions. Across the genome, distributions of ROH are not uniform, and they have distinctive continental patterns. ROH frequencies across the genome are correlated with local genomic variables such as recombination rate, as well as with signals of recent positive selection. In addition, long ROH are more frequent in genomic regions harboring genes associated with autosomal-dominant diseases than in regions not implicated in Mendelian diseases. These results provide insight into the way in which homozygosity patterns are produced, and they generate baseline homozygosity patterns that can be used to aid homozygosity mapping of genes associated with recessive diseases.     

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.     

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.     

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.     

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.     

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.