Inbreeding effects on the parameters of the growth function in three strains of Iberian pigs

How Depressing Is Inbreeding? A Meta-Analysis of 30 Years of Research on the Effects of Inbreeding in Livestock

Inbreeding depression has been widely documented for livestock and other animal and plant populations. Inbreeding is generally expected to have a stronger unfavorable effect on fitness traits than on other traits. Traditionally, the degree of inbreeding depression in livestock has been estimated as the slope of the linear regression of phenotypic values on pedigree-based inbreeding coefficients. With the increasing availability of SNP-data, pedigree inbreeding can now be replaced by SNP-based measures. We performed a meta-analysis of 154 studies, published from 1990 to 2020 on seven livestock species, and compared the degree of inbreeding depression (1) across different trait groups, and (2) across different pedigree-based and SNP-based measures of inbreeding. Across all studies and traits, a 1% increase in pedigree inbreeding was associated with a median decrease in phenotypic value of 0.13% of a trait’s mean, or 0.59% of a trait’s standard deviation. Inbreeding had an unfavorable effect on all sorts of traits and there was no evidence for a stronger effect on primary fitness traits (e.g., reproduction/survival traits) than on other traits (e.g., production traits or morphological traits). p-values of inbreeding depression estimates were smaller for SNP-based inbreeding measures than for pedigree inbreeding, suggesting more power for SNP-based measures. There were no consistent differences in p-values for percentage of homozygous SNPs, inbreeding based on runs of homozygosity (ROH) or inbreeding based on a genomic relationship matrix. The number of studies that directly compares these different measures, however, is limited and comparisons are furthermore complicated by differences in scale and arbitrary definitions of particularly ROH-based inbreeding. To facilitate comparisons across studies in future, we provide the dataset with inbreeding depression estimates of 154 studies and stress the importance of always reporting detailed information (on traits, inbreeding coefficients, and models used) along with inbreeding depression estimates.

Inbreeding depression load for litter size in Entrepelado and Retinto Iberian pig varieties1

Individual-specific hidden inbreeding depression load (IDL) can be accounted for in livestock populations by appropriate best linear unbiased prediction approaches. This genetic effect has a recessive pattern and reveals when inherited in terms of identity-by-descent. Nevertheless, IDL inherits as a pure additive genetic background and can be selected using standard breeding values. The main target of this research was to evaluate IDL for litter size in 2 Iberian pig varieties (Entrepelado and Retinto) from a commercial breeding-stock. Analyses were performed on the total number of piglets born (both alive and dead) and used data from 3,200 (8.02 ± 0.04 piglets/litter) Entrepelado and 4,744 Retinto litters (8.40 ± 0.03 piglets/litter). Almost 50% of Entrepelado sows were inbred (1.7% to 25.0%), whereas this percentage reduced to 37.4% in the Retinto variety (0.2% to 25.0%). The analytical model was solved by Bayesian inference and accounted for 2 systematic effects (sow age and breed/variety of the artificial insemination boar), 2 permanent environmental effects (herd-year-season and sow), and 2 genetic effects (IDL and infinitesimal additive). In terms of posterior means (PM), additive genetic and IDL variances were similar in the Entrepelado variety (PM, 0.68 vs. 0.76 piglets2, respectively) and their 95% credibility intervals (95CI) overlapped, although without including zero (0.38 to 0.94 vs. 0.15 to 1.31 piglets2, respectively). The same pattern revealed in the Retinto variety, with IDL variance (PM, 0.41 piglets2; 95CI, 0.07 to 0.88 piglets2) slightly larger than the additive genetic variance (PM, 0.37 piglets2; 95CI, 0.16 to 0.59 piglets2). The relevance of IDL was also checked by a Bayes factor and the deviance information criterion, the model including this effect being clearly favored in both cases. Although the analysis assumed null genetic covariance between IDL and infinitesimal additive effects, a moderate negative correlation (-0.31) was suggested when plotting the PM of breeding values in the Entrepelado variety; a negative genetic trend for IDL was also revealed in this Iberian pig variety (-0.25 piglets for 100% inbred offspring of individuals born in 2014), whereas no trend was detected in Retinto breeding-stock. Those were the first estimates of IDL in a commercial livestock population, they giving evidence of a relevant genetic background with potential consequences on the reproductive performance of Iberian sows.

Genomic Analysis Reveals Specific Patterns of Homozygosity and Heterozygosity in Inbred Pigs

The inbred strain of miniature pig is an ideal model for biomedical research due to its high level of homozygosity. In this study, we investigated genetic diversity, relatedness, homozygosity, and heterozygosity using the Porcine SNP60K BeadChip in both inbred and non-inbred Wuzhishan pigs (WZSPs). Our results from multidimensional scaling, admixture, and phylogenetic analyses indicated that the inbred WZSP, with its unique genetic properties, can be utilized as a novel genetic resource for pig genome studies. Inbreeding depression and run of homozygosity (ROH) analyses revealed an average of 61 and 12 ROH regions in the inbred and non-inbred genomes of WZSPs, respectively. By investigating ROH number, length, and distribution across generations, we further briefly studied the impacts of recombination and demography on ROH in these WZSPs. Finally, we explored the SNPs with higher heterozygosity across generations and their potential functional implications in the inbred WZSP. We detected 56 SNPs showing constant heterozygosity with He = 1 across six generations in inbred pigs, while only one was found in the non-inbred population. Among these SNPs, we observed nine SNPs located in swine RefSeq genes, which were found to be involved in signaling and immune processes. Together, our findings indicate that the inbred-specific pattern of homozygosity and heterozygosity in inbred pigs can offer valuable insights for elucidating the mechanisms of inbreeding in farm animals.

Direct and maternal additive effects are not the main determinants of Iberian piglet perinatal mortality

Data of 127,800 Iberian piglets were used to study genetic parameters of mortality at birth at the piglet level. These records proceed from three data sets: 4,987 litter of 2,156 sows of a dam line, 2,768 litter of 817 sows of a complete diallel cross between four Iberian strains and 7,153 litter of 2,113 sows of the Torbiscal composite line. Perinatal mortality was considered as a binary trait, and Bayesian threshold animal models were fitted to separately analyse the three data sets. The posterior means of direct heritability were 0.010, 0.004 and 0.003, and those of maternal heritability were 0.034, 0.011 and 0.014 for dam line, diallel cross and Torbiscal line, respectively. Important effects of litter size and parity order were inferred in the three data sets, of within-breed cross-breeding parameters in the diallel cross and of sex and sow handling in the Torbiscal line Therefore, the inclusion of perinatal mortality in the objective of selection is questionable in this breed and strategies for reducing piglet mortality successful in other breeds should be considered.

Genomic inbreeding coefficients based on the distribution of the length of runs of homozygosity in a closed line of Iberian pigs

BACKGROUND

The increasing availability of DNA markers provides new metrics of inbreeding based on single nucleotide polymorphisms (SNPs), i.e. molecular inbreeding or the proportion of runs of homozygosity (ROH), as alternatives to traditional pedigree-based inbreeding coefficients. However, none of these metrics incorporate the length of ROH as an indicator of recent inbreeding. Novel inbreeding coefficients that incorporate length of ROH as a random variable with an associated density are investigated.

METHODS

New inbreeding metrics based on the distribution of the length of ROH are proposed: (1) the Kolmolgorov-Smirnov test, (2) a function of the quantiles of the cumulative distribution function of an individual versus the population, and (3) fitting of an exponential distribution to ROH lengths (mean, variance, and the probability of drawing at random a ROH larger than a given threshold). The new inbreeding and pedigree-based metrics were compared using 217 sows of an Iberian line that belong to three groups: C1 (conservation), C2 (conservation derived from C1), and S (selected and derived from C1), with complete pedigrees and genotyped for 35,023 SNPs.

RESULTS

Correlations between pedigree-based and the new genomic inbreeding coefficients ranged from 0.22 to 0.72 but most ranged from 0.60 to 0.70. The correlation between quantile chromosomal inbreeding coefficients (using molecular information of just one chromosome at the time) and chromosomal length was 0.84 (SE = 0.14), supporting the hypothesis that these coefficients incorporate information on ROH length as an indication of recent inbreeding. Kolmogorov-Smirnov and exponential chromosomal inbreeding coefficients were also correlated with chromosomal length (0.57). Chromosome 1 had the largest quantile ROH inbreeding coefficient (largest ROH sizes), whereas chromosome 10 had the lowest (shortest ROH sizes). Selection for lean growth increased ROH-based inbreeding coefficients for group S when compared to unselected groups C1 and C2. At the chromosomal level, this comparison showed that the level of autozygosity and the length of ROH for most of the autosomes increased in the selection line.

CONCLUSIONS

Quantile and exponential probability inbreeding coefficients using ROH length as a random variable provide additional information about recent inbreeding compared to existing inbreeding coefficients such as molecular, pedigree-based or total ROH content inbreeding coefficients.

Estimates of direct and indirect effects for early juvenile survival in captive populations maintained for conservation purposes: the case of Cuvier's gazelle

Together with the avoidance of any negative impact of inbreeding, preservation of genetic variability for life-history traits that could undergo future selective pressure is a major issue in endangered species management programmes. However, most of these programmes ignore that, apart from the direct action of genes on such traits, parents, as contributors of offspring environment, can influence offspring performance through indirect parental effects (when parental genotype and phenotype exerts environmental influences on offspring phenotype independently of additive genetic effects). Using quantitative genetic models, we estimated the additive genetic variance for juvenile survival in a population of the endangered Cuvier's gazelle kept in captivity since 1975. The dataset analyzed included performance recording for 700 calves and a total pedigree of 740 individuals. Results indicated that in this population juvenile survival harbors significant additive genetic variance. The estimates of heritability obtained were in general moderate (0.115-0.457) and not affected by the inclusion of inbreeding in the models. Maternal genetic contribution to juvenile survival seems to be of major importance in this gazelle's population as well. Indirect genetic and indirect environmental effects assigned to mothers (i.e., maternal genetic and maternal permanent environmental effects) roughly explain a quarter of the total variance estimated for the trait analyzed. These findings have major evolutionary consequences for the species as show that offspring phenotypes can evolve strictly through changes in the environment provided by mothers. They are also relevant for the captive breeding programme of the species. To take into account, the contribution that mothers have on offspring phenotype through indirect genetic effects when designing pairing strategies might serve to identify those females with better ability to recruit, and, additionally, to predict reliable responses to selection in the captive population.

Genome-wide estimates of coancestry and inbreeding in a closed herd of ancient Iberian pigs

Maintaining genetic variation and controlling the increase in inbreeding are crucial requirements in animal conservation programs. The most widely accepted strategy for achieving these objectives is to maximize the effective population size by minimizing the global coancestry obtained from a particular pedigree. However, for most natural or captive populations genealogical information is absent. In this situation, microsatellites have been traditionally the markers of choice to characterize genetic variation, and several estimators of genealogical coefficients have been developed using marker data, with unsatisfactory results. The development of high-throughput genotyping techniques states the necessity of reviewing the paradigm that genealogical coancestry is the best parameter for measuring genetic diversity. In this study, the Illumina PorcineSNP60 BeadChip was used to obtain genome-wide estimates of rates of coancestry and inbreeding and effective population size for an ancient strain of Iberian pigs that is now in serious danger of extinction and for which very accurate genealogical information is available (the Guadyerbas strain). Genome-wide estimates were compared with those obtained from microsatellite and from pedigree data. Estimates of coancestry and inbreeding computed from the SNP chip were strongly correlated with genealogical estimates and these correlations were substantially higher than those between microsatellite and genealogical coefficients. Also, molecular coancestry computed from SNP information was a better predictor of genealogical coancestry than coancestry computed from microsatellites. Rates of change in coancestry and inbreeding and effective population size estimated from molecular data were very similar to those estimated from genealogical data. However, estimates of effective population size obtained from changes in coancestry or inbreeding differed. Our results indicate that genome-wide information represents a useful alternative to genealogical information for measuring and maintaining genetic diversity.

Genetic diversity and pedigree analysis of the Finnsheep breed

Genetic diversity in the Finnsheep breed was analyzed by quantifying the demographic trends, the depth of known pedigree, effective population size, and the amount of inbreeding, as well as identifying candidate rams within the current population for future breeding and conservation purposes. Pedigree records of 148,833 animals with a pedigree completeness coefficient > or =0.60 and born from 1989 to 2006 were used to estimate the parameters. Mean inbreeding coefficient increased by 0.10% (P < 0.001) and 0.15% (P < 0.001) per annum in all animals and breeding (i.e., reproducing) animals, respectively. Average relationship coefficients among rams, among ewes, and between rams and ewes in breeding animals increased over time and reached 1.67, 1.45, and 1.46% in the 2005 cohort, respectively. The average for breeding rams was above the other 2 averages in almost all birth years. The observed generally low average relationship coefficients between rams and ewes indicate that no extra restrictions on the use of the breeding animals are needed in the near future. Average generation interval was 2.85 yr in the studied period, and the effective population size was estimated to be 119 and 122 using different methods. Relationship coefficients of rams with other breeding rams and rams with breeding ewes are suggested to aid in situ and ex situ conservation decisions on maintaining genetic diversity of Finnsheep.

Meat quality of tenderloin from Iberian pigs as affected by breed strain and crossbreeding

In 2007, a new National Quality Standard was published in Spain to regulate the products derived from the Iberian pig carcass, including for the first time fresh meat. In the same way, four different Iberian strains were recognized as official (Lampiño, Entrepelado, Retinto and Torbiscal). A batch (n=10) of each pig strain was selected using neutral DNA markers, and another batch of the most common crossbreeding pigs (Iberian×Duroc) was included into the study as a control. The main meat quality parameters of tenderloin, the most expensive meat cut for fresh consumption, from those five pig groups were analysed. Retinto and Lampiño strains showed the closest phenotypic distances, followed by Entrepelado strain. Meat from crossed and Torbiscal pigs had lower water holding capacity, L∗ and a∗, and higher SFA than meat from the other three strains. Crossbred pigs had the lowest protein, intramuscular fat and PUFA contents.

Mitochondrial genome polymorphisms associated with longissimus muscle composition in Iberian pigs

We carried out a study to investigate the associations between mitochondrial DNA polymorphisms and meat quality traits (intramuscular fat and protein content of the longissimus) in an Iberian porcine line named Torbiscal. The studied pigs (n = 319) belong to 9 maternal lineages and were previously assigned to 6 mitochondrial haplotypes (H1 to H6), based on Cytochrome b and Dloop sequences. Statistical analyses, following a bivariate mixed model, show a greater fat content and lower protein content in H3 haplotype carriers than H1, H2, H4, H5, and H6 haplotype carriers. The magnitudes of these differences are close to 1 g of fat and -0.5 g of protein per 100 g of muscle. To identify the causative mutation of these effects on intramuscular fat and protein contents, the complete mitochondrial DNA sequence of 6 individuals was determined, each one carrying a different mitochondrial haplotype. The alignments of these 6 complete mitochondrial sequences allowed identification of 32 substitutions and 2 indels. Two polymorphic positions were exclusively detected in H3 carriers: a synonymous transition 9104C > T in the gene-coding region of Cytochrome c oxidase subunit III and a substitution 715A > G in 12S rRNA. Genotyping results of a larger number of Torbiscal samples showed the exclusive presence of 9104T and 715G alleles in H3 carriers. The detected candidate substitutions are located in essential mitochondrial genes, and although they do not change the amino acid composition, we cannot disregard a potential change in the secondary structure of their corresponding mRNA. The usefulness of these polymorphisms as markers in selection programs requires validation of the consistency of these results in other Iberian pig lines.

Characterization of OCA2 cDNA in different porcine breeds and analysis of its potential effect on skin pigmentation in a red Iberian strain

Although the function of the OCA2 gene product has not been totally clarified, variation in OCA2 has been associated with skin and hair pigmentation in human and mouse. However, its contribution to skin colour in domestic species has not been reported. In this study, cDNA and intron 9 sequences of the porcine OCA2 gene have been characterized in several pig populations. The cDNA sequence alignment of 20 animals from eight porcine populations allowed the identification of 10 single nucleotide polymorphisms (SNPs); five of the 10 SNPs were non-synonymous. The intron 9 sequence alignment of 12 animals belonging to four pig populations revealed four additional SNPs. Skin colour variation was analysed in a red strain of Iberian pigs with segregation of three SNPs forming two OCA2 intragenic haplotypes. Results from this study provide evidence of a suggestive dominant effect of haplotypes on colour intensity and indicate an important contribution of additive polygenic effects (h2 = 0.56 +/- 0.21) to the variance of this trait.