Open Access

Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data

  • María Saura1Email author,
  • Almudena Fernández1,
  • Luis Varona2,
  • Ana I Fernández1,
  • Maria Ángeles R de Cara3,
  • Carmen Barragán1 and
  • Beatriz Villanueva1
Genetics Selection Evolution201547:1

Received: 25 July 2014

Accepted: 5 December 2014

Published: 17 January 2015



The current availability of genotypes for very large numbers of single nucleotide polymorphisms (SNPs) is leading to more accurate estimates of inbreeding coefficients and more detailed approaches for detecting inbreeding depression. In the present study, genome-wide information was used to detect inbreeding depression for two reproductive traits (total number of piglets born and number of piglets born alive) in an ancient strain of Iberian pigs (the Guadyerbas strain) that is currently under serious danger of extinction.


A total of 109 sows with phenotypic records were genotyped with the PorcineSNP60 BeadChip v1. Inbreeding depression was estimated using a bivariate animal model in which the inbreeding coefficient was included as a covariate. We used two different measures of genomic inbreeding to perform the analyses: inbreeding estimated on a SNP-by-SNP basis and inbreeding estimated from runs of homozygosity. We also performed the analyses using pedigree-based inbreeding.


Significant inbreeding depression was detected for both traits using all three measures of inbreeding. Genome-wide information allowed us to identify one region on chromosome 13 associated with inbreeding depression. This region spans from 27 to 54 Mb and overlaps with a previously detected quantitative trait locus and includes the inter-alpha-trypsin inhibitor gene cluster that is involved with embryo implantation.


Our results highlight the value of high-density SNP genotyping for providing new insights on where genes causing inbreeding depression are located in the genome. Genomic measures of inbreeding obtained on a SNP-by-SNP basis or those based on the presence/absence of runs of homozygosity represent a suitable alternative to pedigree-based measures to detect inbreeding depression, and a useful tool for mapping studies. To our knowledge, this is the first study in domesticated animals using the SNP-by-SNP inbreeding coefficient to map specific regions within chromosomes associated with inbreeding depression.


The reduction in performance due to inbreeding (i.e. inbreeding depression) has long been documented in plant and animal populations [1]. In general, inbreeding depression is most severe for traits that are closely related with fitness, but other traits can also be affected [2]. The standard approach for estimating inbreeding depression is to regress the phenotype of the trait of interest on the inbreeding coefficient (F). Typically, F, defined as the probability that both alleles at any locus within an individual are identical by descent (IBD), has been computed from pedigree information. However, the current availability of very large numbers of single nucleotide polymorphisms (SNPs) offers new opportunities to improve the accuracy of F estimates [3] and to develop more detailed approaches for detecting inbreeding depression [4-6]. Several potential advantages of using genomic F rather than pedigree-based F (F ped ) have been highlighted [4,5]. Genomic F measures homozygosity directly and thus can more accurately reflect the actual percentage of the genome that is homozygous, whereas F ped is only an expectation of that percentage. Another critical difference is that genomic F allows us to estimate inbreeding and inbreeding depression for specific genomic regions, which is not possible with F ped . In addition, genomic F can be estimated in populations where pedigree recording is difficult or impossible.

Several alternative estimates of genomic F based on SNP genotypes have been proposed. A simple estimate can be obtained on a SNP-by-SNP basis as the proportion of homozygous genotypes [3-6]. However, the drawback of this estimate is that it does not differentiate between IBD and identity by state (IBS). An alternative approach for quantifying individual homozygosity that better reflects IBD is based on runs of homozygosity (ROH). The idea is that autozygous genotypes are not evenly distributed throughout the genome but are distributed in runs that are inherited together [7,8]. This is explained by consanguineous matings causing inheritance of haplotypes that are IBD and result in homozygous stretches along the genome of the offspring [9]. It has been shown that these runs provide a good measure of individual genome-wide autozygosity (F roh ) and allow us also to distinguish between recent and ancient inbreeding [10].

The aim of this study was to detect genomic regions responsible for inbreeding depression for two reproductive traits in a highly inbred strain of Iberian pigs (the Guadyerbas strain) using different measures of genomic F.


Animals, pedigree and phenotypic data

In this study, data originated from Guadyerbas pigs maintained in a small isolated herd at the CIA ‘El Dehesón del Encinar’ (Oropesa, Toledo, Spain). Four males and 20 females founded the herd that has been maintained in isolation under a genetic conservation program [11]. Complete and very accurate genealogy is available since the foundation of the herd in 1944. It comprises about 25 generations and includes 1178 animals born in the herd between 1947 and 2011. The effective population size (N e ) estimated from the rate at which pedigree-based or SNP-based coancestry increases has been estimated to be about 10 [3].

Phenotypic data for total number of piglets born (TNB) and number of piglets born alive (NBA) in successive parities from pedigreed sows were available. Means (standard deviations) for TNB and NBA were 7.39 (2.34) and 7.06 (2.25), respectively. Farrowing facilities were improved in 2000 with a new building where piglets had ad libitum access to creep food from seven days of age onwards. No creep food was supplied before 2000. Thus, there were eight levels of management by combining season of farrowing (four seasons) and farrowing facilities (two). Boars from two strains of Iberian pigs (Guadyerbas and Torbiscal) were used in the matings. Offspring fathered by Torbiscal boars were never maintained in the herd.

Genotyping data

All animals born in the herd between 1992 and 2011 (about six generations) were genotyped. They included 86 males and 141 females. Of these females, 113 had phenotypic records for litter size. Genomic DNA extracted from blood samples was hybridized with the Illumina PorcineSNP60 BeadChip v1 and images were scanned by an external service (Universidad Autónoma de Barcelona, Spain). The SNP chip comprises 62 163 probes that are distributed across 18 autosomes and the two sex chromosomes according to the latest version of the porcine gene annotation (Sscrofa10.2). Genotype calls were obtained with the Genotyping Module of the GenomeStudio Data Analysis software (Illumina Inc.). For the purpose of increasing the power of genotype calling (i.e. to correctly determine the genotype for each individual at each SNP), we included samples from other strains of Iberian pigs. In total, 468 genotyped samples were used in this step. These comprised the 227 Guadyerbas samples (including the 113 sows with phenotypic records) and 241 samples from other strains. The extra samples were not used in any further analysis.

Quality control procedures were applied to identify problematic SNPs and samples. First, SNPs that did not satisfy the following quality control criteria were removed: Call Frequency < 0.99, GenTrainScore < 0.70, AB R Mean < 0.35 and number of inconsistencies with the genealogy > 9 (see Saura et al. [3] for further details on the filtering criteria performed). Unmapped SNPs and SNPs mapped to sex chromosomes were also excluded. A total of 51 127 SNPs remained and were used in subsequent analyses. Note that monomorphic SNPs (i.e., those with a minor allele frequency (MAF) of 0 were not removed. After filtering SNPs, the data were reanalysed and samples with a Call Rate < 0.96 and with a large number of inconsistencies with the genealogy were removed. Four samples were excluded, so the final number of genotyped Guadyerbas females available was 109. The total number of litter records from these 109 genotyped sows was 265.

Inbreeding coefficients

Different estimates of F for the sows were used for the inbreeding depression analyses:
  1. (1)

    Genealogical inbreeding coefficients (F ped ) were obtained using all pedigree information that had been recorded since the foundation of the herd.

  2. (2)

    Genomic SNP-by-SNP inbreeding coefficients (F snp ) were obtained based on the excess of SNP homozygosity, as in Keller et al. [4]. The inbreeding coefficient for individual i (F snp(i) ) was computed as F snp(i) = [(OH i  − EH)/(s − EH), where s is the number of SNPs, OH i is \( {\displaystyle {\sum}_{j=1}^s{X}_{ij},} \) where X ij is an indicator variable taking values of 1 if individual i is homozygous for SNP j, and 0 if individual i is heterozygous for SNP j, and EH is the expected homozygosity in the population. The expected homozygosity was computed as \( {\displaystyle {\sum}_{j=1}^s\left[1-2{p}_j\left(1-{p}_j\right)\right],} \) where p j is the MAF for SNP j. We also computed SNP-by-SNP inbreeding coefficients as the proportion of SNPs that are homozygous for the individual (F snp_r ). Note that F snp_r ranges from 0 to 1 but F snp can be negative (when EH is higher than OH i ).

  3. (3)

    Genomic inbreeding coefficients were also estimated based on ROH (F roh ). For a given individual i, F roh(i) was defined as the proportion of its genome that is in ROH [12]. We used our own Fortran code to detect ROH [13] that were defined using the following criteria: (i) a maximum of two missing genotypes and one heterozygous genotype were permitted in a ROH; (ii) the minimum SNP density required to define a ROH was 1 SNP per 100 kb; (iii) the maximum distance allowed between two consecutive homozygous SNPs in a ROH was 1 Mb; and (iv) the minimum number of SNPs that constitute a ROH was 30. We also performed analyses based on short and long ROH. We defined the inbreeding coefficient based on short ROH for individual i (F roh_short(i) ) as the proportion of its genome that was in ROH of lengths between 0.5 and 5 Mb and the inbreeding coefficient based on long ROH (F roh_long(i) ) as the proportion of its genome that is in ROH of lengths > 5 Mb. These thresholds were applied to assess the relative importance of distant (F roh_short ) versus recent (F roh_long ) inbreeding [13]. Long ROH are expected to be autozygous segments that originated from recent common ancestors, while short ROH are likely to have originated from more remote common ancestors [9].


Inbreeding depression analyses

Inbreeding depression was estimated by regressing the phenotype of the reproductive trait (TNB and NBA) on F. This regression was performed by including F as a covariate in a bivariate animal model analysis. The model equation for both traits was:
$$ \mathbf{y}=\mathbf{X}\beta +{\mathbf{Z}}_{\mathbf{1}}\mathbf{a}+{\mathbf{Z}}_{\mathbf{2}}\mathbf{p}+\mathbf{e}, $$
where y is the vector of observations for TNB or NBA, β is the vector of fixed effects, including the combination of season of farrowing and farrowing facilities (eight levels), parity (four levels), strain of boar (two levels) and the (linear) regression on F, a is the vector of additive genetic effects, p is the vector of permanent environmental effects associated with the sows, e is the vector of random residual effects, and X, Z 1, and Z 2 are incidence matrices relating fixed and random effects to observations. The expectation of y was assumed to be E[y] = X β, and the variances and covariances of the random effects were assumed to be V(a) = A σ 2 a, V(p) = I m σ 2 p and V(e) = I n σ 2 e, where A is the pedigree-based numerator relationship matrix of order N (number of animals in the pedigree), I m and I n are identity matrices of order m (number of sows with litter size records), and n (number of records), respectively, and σ 2 a, σ 2 p, and σ 2 e are the variances of additive genetic effects, permanent environmental effects, and residual effects, respectively. The analyses were performed using the REML/VCE 6.0 [14] and PEST [15] softwares.

Different analyses were performed by using the different inbreeding coefficients in the model (F ped , F snp , F roh , F roh_short and F roh_long ). The analysis using F ped was carried out using all available performance and pedigree records, which included 823 sows with reproductive records on 2712 litters and a total pedigree file with 1032 animals. Analyses using F snp and F roh included only records for the 109 genotyped females and used estimates of variance and covariance components obtained from the F ped analysis. Three analyses were implemented with the genomic inbreeding coefficients: (i) using F coefficients across the whole genome; (ii) using F coefficients for each autosome; and (iii) using F coefficients for specific regions within chromosomes.

Ethical statement

The current study was carried out under a Project License from the INIA Scientific Ethic Committee. Animal manipulations were performed according to the Spanish Policy for Animal Protection RD1201/05, which meets the European Union Directive 86/609 about the protection of animals used in experimentation. We hereby confirm that the INIA Scientific Ethic Committee, which is the named IACUC for the INIA, specifically approved this study.


Descriptive statistics for F ped , F snp , F snp_r and F roh are summarized in Table 1. The average F snp_r was very high but this is simply due to an effect of scale since F snp_r is not corrected for homozygosity in the base population and therefore includes IBS. The average inbreeding coefficient based on long ROH was very close to F ped , whereas the average F roh_short was about four times lower than the average F ped . Although short ROH were more abundant (about double the number) than long ROH (Figure 1), their total contribution to the autosomal genome was relatively small (11.1% for short versus 34.2% for long ROH). Chromosomes SSC1 (SSC for Sus scrofa) and SSC13 contained the longest ROH, with sizes greater than 170 Mb, while the maximum size for ROH in the other autosomes was 120 Mb (Figure 2). This was as expected, given the negative correlation between physical (Mb) chromosome size and recombination rate in pigs [16]. Fisher [17] noted that the expected length of a DNA segment that is IBD follows an exponential distribution with mean equal to 1/2 g Morgans, where g is the number of generations since the common ancestor. Given this, large ROH (>5 Mb) reflect the expected inbreeding from a common ancestor that lived less than 10 generations ago. In our data, the average F roh_long was about three times as large as the average F roh_short , thus suggesting that most inbreeding was recent.
Table 1

Mean, standard deviation (SD) and range (Min, Max) of pedigree and SNP-derived inbreeding coefficients






F ped





F snp_r





F snp





F roh





F roh_short





F roh_long





F ped = pedigree-based inbreeding; F snp_r = genomic SNP-by-SNP inbreeding estimated as the proportion of homozygous genotypes; F snp = genomic SNP-by-SNP inbreeding based on the excess of SNP homozygosity; F roh = genomic inbreeding based on all ROH; F roh_short = genomic inbreeding based on short ROH; F roh_long = genomic inbreeding based on long ROH.

Figure 1

Distribution of short and long ROH.

Figure 2

Distribution of ROH according to their size within each autosome. The number of SNPs per autosome is indicated in brackets.

Pearson correlations between the different inbreeding coefficients across individuals were positive and high (Table 2), except for correlations with F roh_short , which were smaller and negative. These negative correlations should be interpreted with caution because, although short ROH are likely to have originated from remote common ancestors, they could be covered up or included in some of the longer ROH. The high positive correlations of F ped and F snp (which are measures of overall inbreeding) with F roh_long also support most of the inbreeding to be recent. Correlations computed for each autosome followed the same pattern as those computed for the whole genome [See Additional file 1]. In magnitude, the highest positive correlations across autosomes were consistently those between F snp and F roh_long and the highest negative correlations were those between F snp and F roh_short .
Table 2

Pearson correlations (SE) between different inbreeding coefficients measured at the whole-genome level


F ped

F snp

F roh

F roh_short

F snp

0.662 (0.051)


F roh

0.631 (0.053)

0.968 (0.017)


F roh_short

−0.241 (0.066)

−0.499 (0.059)

−0.450 (0.061)


F roh_long

0.603 (0.054)

0.959 (0.019)

0.974 (0.015)

−0.641 (0.052)

F ped = pedigree-based inbreeding; F snp = genomic SNP-by-SNP inbreeding based on the excess of SNP homozygosity; F roh = genomic inbreeding based on all ROH; F roh_short = genomic inbreeding based on short ROH; F roh_long = genomic inbreeding based on long ROH.

As indicated above, estimates of genetic parameters were obtained from the pedigree-based analysis and then used in all the subsequent analyses that used genomic inbreeding coefficients. Estimates of heritability (0.052 ± 0.025 for NBA and 0.077 ± 0.029 for TNB) and permanent environmental coefficients (0.073 ± 0.023 for NBA and 0.068 ± 0.024 for TNB) were of the same order of magnitude as those previously reported for Iberian pigs [18]. These estimates led to a repeatability estimate of 0.15 for both traits. The estimated genetic correlation between traits was very high (0.964 ± 0.024) and also close to previous estimates [19,20].

A significant reduction in both NBA and TNB with increases in the inbreeding coefficient was observed when performing the pedigree-based analysis using all available data (i.e. using records from 823 sows). Estimates of inbreeding depression were −0.197 ± 0.092 for NBA and −0.211 ± 0.104 for TNB per 10% increase in F ped . Although the effect of genomic F at the whole-genome level (i.e., using F snp and F roh ) with litter size was not significant (−0.267 ± 0.186 and −0.253 ± 0.194 for NBA and TNB, respectively, when using F snp , and −0.415 ± 0.374 and −0.335 ± 0.391 for NBA and TNB, respectively, when using F roh ), significant (p < 0.05 for both traits, see Table 3) inbreeding depression was found for both traits when the genomic analyses were carried out at the autosomal level for SSC13. Estimates of inbreeding depression that were obtained from the analyses performed for each autosome using F snp are presented in Figure 3 for NBA and TNB. Only SSC13 showed a significant effect. The estimated inbreeding depression for this chromosome was −0.121 ± 0.047 and −0.117 ± 0.049 per 10% increase in F snp for both NBA and TNB, respectively.
Table 3

Estimates of inbreeding depression on SSC13 for NBA and TNB


NBA mean (SE)

TNB mean (SE)

F snp

−0.121 (0.047)*

−0.117 (0.049)*

F roh

−0.230 (0.087)**

−0.231 (0.091)*

F roh_short

0.340 (0.380)

0.562 (0.395)

F roh_long

−0.181 (0.074)*

−0.188 (0.077)*

NBA = number of piglets born alive; TNB = total number of piglets born; results are expressed as the change in phenotypic mean per 10% increase in F; F snp = genomic SNP-by-SNP inbreeding based on the excess of SNP homozygosity; F roh = genomic inbreeding based on all ROH; F roh_short = genomic inbreeding based on short ROH; F roh_long = genomic inbreeding based on long ROH; significant values are indicated with asterisk (*p < 0.05, **p < 0.01).

Figure 3

Inbreeding depression estimates expressed as the change in phenotypic mean per 10% increase in F snp and 95% confidence intervals across autosomes for number of piglets born alive (NBA) (a) and total number of piglets born (TNB) (b).

Reductions in the number of piglets per 10% increase in F snp , F roh and F roh_long for SSC13 were all significant (see Table 3) and of the same order of magnitude as those derived from the pedigree-based analyses, which were based on the entire genome. There was no significant effect of F roh_short on either NBA or TNB. This may be explained by purging of deleterious alleles in ancient generations, or by a bottleneck.

In order to detect specific genomic regions that cause inbreeding depression, all autosomes were fragmented in segments of equal size (three, five and eight segments per autosome) and three additional analyses per chromosome were carried out using F snp . When autosomes were divided into three segments, inbreeding depression was only significant for the first region of SSC13 (0.0 - 73 Mb) for both traits (Table 4). When they were divided in five segments, significance was found for the first (0–44 Mb) and the second (44–88 Mb) regions of SSC13 for both traits. Finally, when autosomes were divided into eight fragments, only the second region of SSC13 showed a significant result for both traits. This region is 32.4 Mb in size and is located between 27 and 55 Mb. No other regions showed a significant result for any of the traits.
Table 4

Estimates of inbreeding depression for different regions of SSC13 for NBA and TNB






Start (Mb)

End (Mb)

Mean (SE)

p value

Mean (SE)

p value

SSC13 divided in three segments




−0.138 (0.039)


−0.144 (0.041)





−0.035 (0.028)


−0.033 (0.029)





−0.040 (0.038)


−0.032 (0.040)


SSC13 divided in five segments




−0.131 (0.041)


−0.135 (0.043)





−0.066 (0.025)


−0.071 (0.026)





−0.039 (0.023)


−0.038 (0.023)





−0.031 (0.036)


−0.016 (0.038)





−0.019 (0.035)


−0.020 (0.036)


SSC13 divided in eight segments




−0.046 (0.046)


−0.077 (0.043)





−0.072 (0.020)


−0.064 (0.020)





−0.029 (0.020)


−0.039 (0.020)





−0.028 (0.022)


−0.024 (0.023)





−0.033 (0.023)


−0.036 (0.024)





−0.022 (0.037)


−0.005 (0.038)





−0.035 (0.023)


−0.036 (0.024)





0.003 (0.027)


0.006 (0.028)


NBA = number of piglets born alive; TNB = total number of piglets born; results are expressed as the change in phenotypic mean per 10% increase in F snp ; significant values are indicated with asterisk (**p < 0.01, ***p < 0.001).

We also performed a bootstrap test by generating 1000 bootstrap replicates of the individual phenotypes and genotypes for the region that was significant in the inbreeding depression analysis. Bootstrap replicates were created by randomly resampling the individual phenotypes and inbreeding depression analyses were repeated for each replicate. When fragmenting SSC13 in eight segments, the average of the estimates for the effect across bootstrap replicates was 0, with a standard error that was lower than that obtained from the real data (0.020) in all cases except one. This is equivalent to a bootstrap p-value of 0.001. These results constitute more evidence that the signals detected are not spurious.


We have detected inbreeding depression associated with a specific region of SSC13 for two reproductive traits in a highly inbred strain of Iberian pigs using different measures of genome-wide inbreeding coefficients. It is important to note that the signal detected was significant despite the small sample size available for the study (109 sows). This indicates that it would be extremely rare that the effects detected here are spurious.

Also, our results are consistent with those obtained by Noguera et al. [19] who performed one of the first genome-wide scans for prolificacy traits in pigs. They used data from a Guadyerbas x Meishan F2 intercross using SNPs and microsatellite markers and detected a quantitative trait locus (QTL) on SSC13 for both NBA and TNB. This QTL region extends from about 38 to 194 Mb and overlaps with the region identified here. Specifically, the region detected in the inbreeding depression analysis is shorter (it spans from 27 to 54 Mb) and overlaps with the first part of the QTL region detected by Noguera et al. [19]. We examined the gene content of this common region by using the porcine genome annotation Sscrofa10.2 in BioMart tool of Ensembl ( and the Ensembl Genes 69 database and found 271 annotated genes. Interestingly three of these genes, the inter-alpha-trypsin inhibitor heavy chains 1, 3 and 4 (ITIH-1, ITIH-3 and ITIH-4) map to 38 Mb on SSC13 and play several important roles in maintaining the uterine surface glycocalyx during placental attachment in pigs [20]. Moreover, these genes have been previously associated with NBA and TNB [21]. Specifically, using the same material as Noguera et al. [19], Balcells et al. [21] analyzed the porcine ITIH-1, −3 and −4 gene sequences in order to identify polymorphisms that could explain differences in prolificacy of sows. Their results revealed significant associations with NBA and TNB for two SNPs in ITIH-1, four SNPs in ITIH-3, and four SNPs in ITIH-4. Thus, the studies of Noguera et al. [19] and Balcells et al. [21] support our findings, since genes that affect both NBA and TNB are located in the region identified on SSC13. Another recent whole-genome association study identified QTL regions for NBA and TNB that partially overlapped with the region identified here [22].

Several studies that compared the same inbreeding coefficients as used here indicated that F roh is a better measure of IBD than F snp [4,5,23]. However, Keller et al. [4] showed that as N e decreases, the similarity of both measures of molecular inbreeding (F roh and F snp ) increases. Based on this and going down to an N e of 10 (the estimate for this herd [3]), we can expect that the measures of F snp and F roh will be highly correlated (0.97 ± 0.02). In fact, the correlation between F ped and F roh found in our study was similar to that between F ped and F snp (0.63 ± 0.05 and 0.66 ± 0.05, respectively) All this information validates the use of F snp as a measure of IBD in this population, and therefore, as a suitable coefficient to perform the inbreeding depression analyses.

Previous studies aimed at detecting inbreeding depression using SNPs have focused on the whole genome level [6,23-25]. These include human studies that investigated the association between inbreeding and particular diseases [5,11,13,26,27]. Only a few of these have attempted to identify specific genomic regions that cause depression [11,13,27]. Keller et al. [8] used a ROH mapping approach to analyze the association of F roh with schizophrenia risk. The approach consisted of dividing the autosomes in a large number of segments of equal size and performing regressions of disease status on whether or not individuals had a ROH in each segment. They found significant associations between specific genomic regions and disease status. Recently Pryce et al. [28] have followed a similar approach for detecting inbreeding depression for fertility and milk production traits in dairy cattle.

In order to compare our approach for detecting genomic regions associated to inbreeding depression using F snp with the approach of Keller et al. [8], we divided the autosomes in segments of approximately 2 Mb and recorded for presence or absence of ROH within each segment. The results from these regressions showed a significant effect in the same region on SSC13 as detected with the analysis using F snp [See Additional file 2]. Although this significant effect was lost when multitest correction was applied, the results support our previous finding. Additional file 2 shows results for NBA but the same pattern was found for TNB. It should be noted that our analysis using a continuous variable (i.e., F snp ) instead of a dichotomous variable (i.e., presence/absence of ROH) has more power to detect associations between phenotypes and inbreeding within specific regions.

We also performed an association analysis for the region involved in inbreeding depression. It was conducted using the same data (265 phenotypic records on 109 sows) and statistical model as described above for inbreeding depression analyses but also included the SNP genotype as a fixed effect. Each SNP was tested separately for association with the trait and both additive and dominant effects were estimated. We found several SNPs that showed significant dominance effects and non-significant additive effects for both traits (Additional file 3 shows the results for NBA only but similar results were found for TNB [See Additional file 3]). This result is consistent with the findings from the inbreeding depression analyses and agrees with directional dominance for the genes involved.


In summary, using genome-wide SNP information, we detected inbreeding depression in a specific region that contains genes associated with litter size in the isolated population of Guadyerbas Iberian pigs. The availability of dense SNP platforms has created opportunities to estimate homozygosity without using pedigree relationships and to obtain patterns of homozygosity along an individual’s genome. Our results highlight the importance of SNP chips for providing new insights into where genes causing inbreeding depression are located in the genome and thus, offer a complementary tool to QTL analysis for mapping studies.



This work was funded by the INIA and by Ministerio de Economía y Competitividad, Spain (grants RZ2010-00009-00-00 and CGL2012-39861-C02-02). MS was supported by a Juan de la Cierva fellowship (JCI-2011-10896) from the Ministerio de Economía y Competitividad, Spain. MARDC was funded by LabEx grant ANR-10-LABX-0003-BCDiv from the Agence Nationale de la Recherche within the program Investissements d’Avenir ANR-11-IDEX-0004-02. We thank Prof MA Toro for useful comments on the manuscript, MC Rodríguez and L Silió for facilitating samples and useful discussions and F García and R Benítez for technical support.

We most especially want to thank Dr Hélène Hayes and Prof Jack Dekers for their valuable help during the final edition of the manuscript. Two anonymous referees provided helpful comments on the manuscript.

Authors’ Affiliations

Departamento de Mejora Genética Animal, INIA
Unidad de Genética Cuantitativa y Mejora Animal, Facultad de Veterinaria, Universidad de Zaragoza
Laboratoire d’Eco-anthropologie et Ethnobiologie, UMR 7206 CNRS/MNHN/Université Paris 7, Muséum National d’Histoire Naturelle


  1. Charlesworth B, Charlesworth D: The genetic basis of inbreeding depression. Genet Res 1999, 74:329–340.PubMedView ArticleGoogle Scholar
  2. DeRose MA, Roff DA: A comparison of inbreeding depression in life-history and morphological traits in animals. Evolution 1999, 53:1288–1292.View ArticleGoogle Scholar
  3. Saura M, Fernández A, Rodríguez MC, Toro MA, Barragán C, Fernández AI, Villanueva B: Genome-wide estimates of coancestry and inbreeding in a closed herd of Iberian pigs. PLoS ONE 2013, 8:e78314.PubMed CentralPubMedView ArticleGoogle Scholar
  4. Keller MC, Visscher PM, Goddard ME: Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data. Genetics 2011, 189:237–249.PubMed CentralPubMedView ArticleGoogle Scholar
  5. McQuillan R, Eklund N, Pirastu N, Kuningas M, McEvoy B, Esko T, Corre T, Davies G, Kaakinen M, Lyytikäinen LP, Kristiansson K, Havulinna AS, Gögele M, Vitart V, Tenesa A, Aulchenko Y, Hayward C, Johansson A, Boban M, Ulivi S, Robino A, Boraska V, Igl W, Wild SH, Zgaga L, Amin N, Theodoratou E, Polašek O, Girotto G, Lopez LM: Evidence of inbreeding depression in human height. PLoS Genet 2012, 8:e1002655.PubMed CentralPubMedView ArticleGoogle Scholar
  6. Bjelland DW, Weigel KA, Vukasinovic N, Nkrumah JD: Evaluation of inbreeding depression in Holstein cattle using whole-genome SNP markers and alternative measures of genomic inbreeding. J Dairy Sci 2013, 96:1–10.View ArticleGoogle Scholar
  7. Broman KW, Weber JL: Long homozygous chromosomal segments in reference families from the centre d’Etude du polymorphisme humain. Am J Hum Genet 1999, 65:1493–1500.PubMed CentralPubMedView ArticleGoogle Scholar
  8. Keller MC, Simonson MA, Ripke S, Neale BM, Gejman PV, Howrigan DP, Lee SH, Lencz T, Levinson DF, Sullivan PF, Schizophrenia Psychiatric Genome-Wide Association Study Consortium: Runs of homozygosity implicate autozygosity as a schizophrenia risk factor. PLoS Genet 2012, 8:e1002656.PubMed CentralPubMedView ArticleGoogle Scholar
  9. Pemberton TJ, Absher D, Feldman MW, Myers RM, Rosenberg NA, Li JZ: Genomic patterns of homozygosity in worldwide human populations. Am J Hum Genet 2012, 91:275–292.PubMed CentralPubMedView ArticleGoogle Scholar
  10. McQuillan R, Leutenegger AL, Abdel-Rahman R, Franklin CS, Pericic M, Barac-Lauc L, Smolej-Narancic N, Janicijevic B, Polasek O, Tenesa A, Macleod AK, Farrington SM, Rudan P, Hayward C, Vitart V, Rudan I, Wild SH, Dunlop MG, Wright AF, Campbell H, Wilson JF: Runs of homozygosity in European populations. Am J Hum Genet 2008, 83:359–372.PubMed CentralPubMedView ArticleGoogle Scholar
  11. Toro MA, Rodrigáñez J, Silió L, Rodríguez MC: Genealogical analysis of a closed herd of black hairless Iberian Pigs. Conserv Biol 2000, 14:1843–1851.View ArticleGoogle Scholar
  12. Lencz T, Lambert C, DeRosse P, Burdick KE, Morgan TV, Kane JM, Raju Kucherlapati R, Malhotra AK: Runs of homozygosity reveal highly penetrant recessive loci in schizophrenia. Proc Natl Acad Sci USA 2007, 104:19942–19947.PubMed CentralPubMedView ArticleGoogle Scholar
  13. de Cara MAR, Villanueva B, Toro MA, Fernández J: Using genomic tools to maintain diversity and fitness in conservation programmes. Mol Ecol 2013, 22:6091–6099.PubMedView ArticleGoogle Scholar
  14. Groeneveld E, Kovac M, Mielenz N: VCE Users’ Guide and Reference Manual Version 6.0. Mariensee, Germany: Institute of Farm Animal Genetics; 2008.Google Scholar
  15. Groeneveld E, Kovac M, Wang T: PEST, a General Purpose BLUP Package for Multivariate Prediction and Estimation. In Proceedings of the 4th World Congress on Genetics Applied to Livestock Production: 23–27 July 1990. Edinburgh. Edited by W. G. Hill, R. Thompson and J. A. Wooliams; 1990:488–491.Google Scholar
  16. Muñoz M, Alves E, Ramayo-Caldas Y, Casellas J, Rodríguez C, Folch JM, Silió L, Fernández AI: Recombination rates across porcine autosomes inferred from high-density linkage maps. Anim Genet 2011, 43:620–623.PubMedView ArticleGoogle Scholar
  17. Fisher RA: A fuller theory of “junctions” in inbreeding. Heredity 1954, 8:187–197.View ArticleGoogle Scholar
  18. Fernández A, Rodrigáñez J, Zuzúarregui J, Rodríguez MC, Silió L: Genetic parameters for litter size and weight at different parities in Iberian pigs. Span J Agric Res 2008, 6:98–106.View ArticleGoogle Scholar
  19. Noguera JL, Rodríguez C, Varona L, Tomàs A, Muñoz G, Ramírez O, Barragán C, Arqué M, Bidanel JP, Amills M, Ovilo C, Sánchez A: A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL. BMC Genomics 2009, 10:636.PubMed CentralPubMedView ArticleGoogle Scholar
  20. Geisert RD, Ashworth MD, Malayer JR: Expression of inter-alpha-trypsin inhibitor heavy chains in endometrium of cyclic and pregnant gilts. Reproduction 2003, 126:621–627.PubMedView ArticleGoogle Scholar
  21. Balcells I, Castelló A, Noguera JL, Fernández-Rodríguez A, Sánchez A, Tomás A: Sequencing and gene expression of the porcine ITIH SSC13 cluster and its effect on litter size in an Iberian × Meishan F2 population. Anim Reprod Sci 2011, 128:85–92.PubMedView ArticleGoogle Scholar
  22. Onteru SK, Fan B, Du Z-Q, Garrick DJ, Stalder KJ, Rothschild MF: A whole-genome association study for pig reproductive traits. Anim Genet 2011, 43:18–26.PubMedView ArticleGoogle Scholar
  23. Ferencakovic M, Hamzic E, Gredler B, Solberg TR, Klemetsdal G, Curik I, Sölkner J: Estimates of autozygosity derived from runs of homozygosity: empirical evidence from selected cattle populations. J Anim Breed Genet 2012, 130:286–293.PubMedView ArticleGoogle Scholar
  24. Walling CA, Nussey D, Morris A, Clutton-Brock TH, Kruuk LEB, Pemberton JM: Inbreeding depression in red deer calves. BMC Evol Biol 2011, 11:318.PubMed CentralPubMedView ArticleGoogle Scholar
  25. Silió L, Rodríguez MC, Fernández A, Barragán C, Benítez R, Óvilo C, Fernández AI: Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP-derived metrics. J Anim Breed Genet 2013, 130:349–360.PubMedGoogle Scholar
  26. Ku CS, Naidoo N, Teo SM, Pawitan Y: Regions of homozygosity and their impact on complex diseases and traits. Hum Genet 2011, 129:1–15.PubMedView ArticleGoogle Scholar
  27. Simón-Sánchez J, Kilarski LL, Nalls MA, Martínez M, Schulte C, Holmans P, International Parkinson’s Disease Genomics Consortium, Wellcome Trust Case Control Consortium, Gasser T, Hardy J, Singleton AB, Wood NW, Brice A, Heutink P, Williams N, Morris HR: Cooperative genome-wide analysis shows increased homozygosity in early onset Parkinson’s disease. PLoS ONE 2012, 7:e28787.PubMed CentralPubMedView ArticleGoogle Scholar
  28. Pryce JE, Haile-Mariam M, Goddard ME, Hayes BJ: Identification of genomic regions associated with inbreeding depression in Holstein and Jersey dairy cattle. Genet Sel Evol 2014, 46:71.PubMed CentralPubMedView ArticleGoogle Scholar


© Saura et al.; licensee BioMed Central. 2015

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.