Inbreeding trend and inbreeding depression in the Danish populations of Texel, Shropshire, and Oxford Down1

Inbreeding Depression and Purging in Fertility and Longevity Traits in Sheep Breeds from Germany

In the present study, we analysed fertility and longevity traits of 22 sheep breeds from Germany with a suitable quantity of data in the national database OviCap. The data comprised merino, meat, country and milk sheep breeds with 62,198 ewes and about 173,000 lambing records, until the fifth lambing. Across-breed means of heritabilities reached estimates of 0.13, 0.17 and 0.18 for number of lambings, average number of lambs born per lambing and number of lambs per lifetime, respectively. For age at first lambing, length of lifetime and productive life, mean heritabilities over breeds were 0.34, 0.17 and 0.32, respectively. The across-breed means of the individual rate of inbreeding were significantly negative for the average number of lambs born per lambing and number of lambs born per lifetime, and for number of lambings it was close to the significance threshold. We found declining slopes for inbreeding depression for the average number of lambs born per lambing and number of lambs born per lifetime in 16 breeds, and significantly negative slopes in five and seven breeds. For lifetime and productive life, 9/22 and 8/22 breeds showed significant inbreeding depression, while for age at first lambing, only 1/22 breeds showed significant inbreeding depression. A significant reduction in inbreeding depression due to purging effects was found for eight breeds. Fitness traits may be subject to forced directional selection. Therefore, sheep breeding programmes should give special consideration to fertility and longevity traits. Fitness related traits seem to be essential in conservation of genetic diversity within sheep breeds.

Estimation of microsatellite-based autozygosity and its correlation with pedigree inbreeding coefficient in crossbred cattle

In countries where farming is largely subsistence, no pedigree records of farm animals are maintained at farmers' herd and scientific mating plans are not observed which leads to the accumulation of inbreeding and loss of production potential. Microsatellites have been widely used as reliable molecular markers to measure inbreeding. We attempted to correlate autozygosity estimated from microsatellite data with the inbreeding coefficient (F) calculated from pedigree data in Vrindavani crossbred cattle developed in India. The inbreeding coefficient was calculated from the pedigree of ninety-six Vrindavani cattle. Animals were further grouped into three groups viz. acceptable/low (F: 0-5%), moderate (F: 5-10%) and high (F: ≥10%), based on their inbreeding coefficients. The overall mean of the inbreeding coefficient was found to be 0.070 ± 0.007. A panel of twenty-five bovine-specific loci were chosen for the study according to ISAG/FAO. The mean FIS, FST, and FIT values were 0.0548 ± 0.025, 0.012 ± 0.001 and 0.0417 ± 0.025, respectively. There was no significant correlation between the FIS values obtained and the pedigree F values. The locus-wise individual autozygosity was estimated using the method-of-moments estimator (MME) formula for locus-specific autozygosity. The autozygosities ascribing to CSSM66 and TGLA53 were found to be significantly (p < .01 and p < .05, respectively) correlated with pedigree F values.

Inbreeding Depression and Purging for Meat Performance Traits in German Sheep Breeds

This study provides estimates on genetic parameters, inbreeding depression and purging for meat performance measures from 25 German sheep breeds. All German meat, merino sheep breeds and breeds of other breeding directions with a sufficient number of pedigree and performance data were included in this study. Phenotypic traits retrieved from the national database OviCap were evaluated: daily weight gain, meatiness score and ultrasound measurements for muscle and fat thickness. We employed animal models to estimate heritability, variance and covariance components for these meat performance traits as well as inbreeding depression and purging. The heritabilities, on average, reached estimates of 0.55, 0.34, 0.53 and 0.61 for daily weight gain, meatiness score and ultrasound measurements for muscle and fat thickness, respectively. We estimated the linear regression slopes for the individual rate of inbreeding, new and ancestral inbreeding, as well as the inbreeding coefficient and its interaction with the inbreeding coefficient of Ballou, employing animal models with non-genetic effects and the additive genetic effect of the animal. Across all breeds, inbreeding was only significant for daily weight gain, whereas for all other traits, estimates were not significant. Within sheep breeds, we found significant inbreeding depression for daily weight gain in German Mutton Merino and German Blackheaded Mutton as well as for the meatiness score in German Whiteheaded Mutton. Significant effects for purging, based on ancestral inbreeding and the interaction effect of the classical inbreeding coefficient with the inbreeding coefficient of Ballou, were not obvious either across or within any sheep breed. A 1% increase in inbreeding significantly decreased the phenotypic trait median of daily weight gain across all sheep breeds by 0.50% and 0.70% of phenotypic and genetic standard deviation, respectively. Purging effects due to ancestral inbreeding were not significant in any breed or across breeds. The results of this study may indicate that inbreeding depression may be more harmful in traits under stronger selection than in traits that exert low selection pressure. The results of this study demonstrate the different effects that result in meat performance traits due to inbreeding. With increasing rates of inbreeding and critical effective population sizes, selection intensity for breeding objectives has to be critically reviewed for each sheep breed. Inbreeding depression and purging should be evaluated in order to prevent a decrease in trait means due to inbreeding and to determine whether detrimental alleles are eliminated.

Inbreeding in the Zwartbles breed population and its influence on meat production in the Czech Republic

The aim of this study was to estimate the inbreeding coefficient in the Zwartbles sheep population in the Czech Republic, as well as to investigate the effect of inbreeding on meat yield indicators. The analysis used data on the entire population since 1997, with the original database containing 13 150 animals. In this population, the average inbreeding coefficient is 3.64 %. There is a significant increase in inbreeding - from zero values to more than 4 % between 1997 and 2021. The effect of inbreeding on the weight of the animals at 100 d of age, the weight at bonitation prior to breeding, the thickness of the musculus longissimus dorsi (MLD), the thickness of the subcutaneous back fat, and the fleshiness of the animals were all evaluated. The value of the inbreeding coefficient had a statistically significant effect on weight at 100 d, MLD thickness, and back fat thickness. A negative relationship between inbreeding and all of these parameters was found when the inbreeding coefficient was increased by 1 %, resulting in a 60.2 g weight reduction at 100 d, 0.06 mm reduction in MLD thickness, and 0.013 mm reduction in back fat thickness.

The genomic structure of isolation across breed, country and strain for important South African and Australian sheep populations

Background

South Africa and Australia shares multiple important sheep breeds. For some of these breeds, genomic breeding values are provided to breeders in Australia, but not yet in South Africa. Combining genomic resources could facilitate development for across country selection, but the influence of population structures could be important to the compatability of genomic data from varying origins. The genetic structure within and across breeds, countries and strains was evaluated in this study by population genomic parameters derived from SNP-marker data. Populations were first analysed by breed and country of origin and then by subpopulations of South African and Australian Merinos.

Results

Mean estimated relatedness according to the genomic relationship matrix varied by breed (-0.11 to 0.16) and bloodline (-0.08 to 0.06) groups and depended on co-ancestry as well as recent genetic links. Measures of divergence across bloodlines (F_ST: 0.04–0.12) were sometimes more distant than across some breeds (F_ST: 0.05–0.24), but the divergence of common breeds from their across-country equivalents was weak (F_ST: 0.01–0.04). According to mean relatedness, F_ST, PCA and Admixture, the Australian Ultrafine line was better connected to the SA Cradock Fine Wool flock than with other AUS bloodlines. Levels of linkage disequilibrium (LD) between adjacent markers was generally low, but also varied across breeds (r²: 0.14–0.22) as well as bloodlines (r²: 0.15–0.19). Patterns of LD decay was also unique to breeds, but bloodlines differed only at the absolute level. Estimates of effective population size (N_e) showed genetic diversity to be high for the majority of breeds (N_e: 128–418) but also for bloodlines (N_e: 137–369).

Conclusions

This study reinforced the genetic complexity and diversity of important sheep breeds, especially the Merino breed. The results also showed that implications of isolation can be highly variable and extended beyond breed structures. However, knowledge of useful links across these population substructures allows for a fine-tuned approach in the combination of genomic resources. Isolation across country rarely proved restricting compared to other structures considered. Consequently, research into the accuracy of across-country genomic prediction is recommended.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08020-3.

Some Genetic and Environmental Effects on Equine Asthma in Polish Konik Horses

Simple Summary

Equine Asthma (EA) is a blanket term covering inflammatory diseases of the lower airways in horses. It includes mild-to-moderate equine asthma, which affects horses of any age, and severe equine asthma, which is typically seen in horses older than 7 years of. Relationships of the disease’s occurrence with sex or breed have been proven. However, some authors consider genetic background a predisposing factor, due to the fact that in some bloodlines a clinical form of asthma is more frequently observed. This indicates serious breeding and economic consequences. This study aimed to identify the factors affecting predisposition to severe equine asthma in the population of Polish Konik horses and some environmental and inbreeding effects on the disease. Generally, in the observed population, EA is negligibly affected by the factors analysed. Individual inbreeding effects on asthma were not confirmed by various statistical approaches, but significant maternal inbreeding effects were observed. These results are very important from the perspective of the currently implemented genetic resource conservation programme.

Abstract

Current knowledge of the genetic and environmental backgrounds of equine asthma seems to be insufficient, especially for primitive horse breeds. The main objectives of this study were to estimate the effects of sex, birth period, stud, parentage line and inbreeding on asthma morbidity in Polish Konik horses. Records of 274 horses (housed in two studs) were analysed. These animals were allocated to maternal and paternal lines. Individual inbreeding coefficients were extracted from the additive relationship matrix. Horses underwent diagnosis based on observation of the basic symptoms (high frequency of coughing and excessive nasal discharge). Subsequently, some horses (28 individuals) were clinically examined to confirm the earlier observations. Generally, no significant effects of parentage line on heaves morbidity were identified by the use of logistic regression, although the Pearson’s chi-squared test had shown that individuals of some maternal and paternal lines had a predisposition to severe equine asthma. It was concluded that the individual inbreeding level is not associated with the incidence of EA, but a significant effect of the maternal inbreeding coefficient may be observed. It was also found that there is some variability in the incidence of this disease between studs.

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.

Litter Size of Sheep (Ovis aries): Inbreeding Depression and Homozygous Regions

Ovine litter size (LS) is an important trait showing variability within breeds. It remains largely unknown whether inbreeding depression on LS exists based on genomic homozygous regions, and whether the homozygous regions resulted from inbreeding are significantly associated with LS in sheep. We here reanalyze a set of single nucleotide polymorphism (SNP) chip of six breeds to characterize the patterns of runs of homozygosity (ROH), to evaluate inbreeding levels and inbreeding depressions on LS, and to identify candidate homozygous regions responsible for LS. Consequently, unique ROH patterns were observed among six sheep populations. Inbreeding depression on LS was only found in Hu sheep, where a significant reduction of 0.016, 0.02, and 0.02 per 1% elevated inbreeding FROH4-8, FROH>8 and the total inbreeding measure was observed, respectively. Nine significantly homozygous regions were found for LS in Hu sheep, where some promising genes for LS possibly via regulation of the development of oocytes (NGF, AKT1, and SYCP1), fertilization (SPAG17, MORC1, TDRD9, ZFYVE21, ADGRB3, and CKB), embryo implantation (PPP1R13B, INF2, and VANGL1) and development (DPPA2, DPPA4, CDCA4, CSDE1, and ADSSL1), and reproductive health (NRG3, BAG5, CKB, and XRCC3) were identified. These results from the present study would provide insights into the genetic management and complementary understandings of LS in sheep.

Genetic and nongenetic factors influencing ewe prolificacy and lamb body weight in a closed Romanov flock

The U.S. Meat Animal Research Center was the first entity in the United States to import the Romanov breed and it has been maintained as a closed flock for over 30 yr. Incorporating this super-prolific breed into crossbred and composite populations has resulted in large improvements in ewe productivity. However, few have quantified factors contributing to genetic and nongenetic variation in ewe reproduction and lamb growth within purebred Romanov populations, which were the objectives of this study. The pedigree contained a total of 8,683 lambs born to 218 and 1,600 unique sires and dams, respectively. Number of lambs born on a per ewe exposed (NLBE) and lambing (NLBL) basis were analyzed in univariate repeatability animal models. As expected, the proportion of phenotypic variance (σP2) in litter size attributable to additive genetic (0.06 to 0.08) and permanent environmental (0.05 to 0.07) effects of the ewe was low. The service sire permanent environmental effect contributed to a small but significant amount of σP2 in NLBE (0.03) but not NLBL. However, the service sire additive genetic effect did not influence σP2 in NLBE or NLBL. Lamb body weight was recorded at birth (BWB) and upon weaning from either milk replacer (~30 d; BWW-N) or their dam (~60 d; BWW-D) and were analyzed in a three-trait model with random additive direct and maternal effects. Estimated direct heritabilities were low for all body weight (BW) traits (0.07 to 0.10). Maternal heritability was moderate for BWB (0.34) but low for weaning BW (0.11 to 0.18). This was the first to report direct and maternal genetic correlations between BW of nursery- and dam-reared lambs, and both were estimated to be moderate (0.43 to 0.47). Additionally, the direct and maternal effects of BWB were more strongly correlated with BWW-N (0.74 to 0.82) than BWW-D (0.17 to 0.33). Despite inbreeding coefficients having increased at a rate of 0.33% per birth year (1986 to 2019) in this flock, they were not consistently associated with reductions in ewe or lamb performance. Parameter estimates generally agreed with those from less-prolific breeds, and results indicate that selection can be an effective means of improving subcomponents of ewe productivity.

Assessment of between-founder heterogeneity in inbreeding depression for reproductive traits in Baluchi sheep

Magnitude of inbreeding depression could be different among descendants of various founder animals, when recessive genetic composition of a population is unevenly distributed among founder animal genomes. Reproductive records of Baluchi sheep for the litter variables size at birth (LSB), size at weaning (LSW), mean weight/lamb born (LMWLB), mean weight/lamb weaned (LMWLW), total weight at birth/ewe lambing (TLWB) and total weight at weaning/ewe lambing (TLWW) were used to examine heterogeneity in inbreeding depression between founder animals. Pedigree-based inbreeding coefficients were proportioned into components coming from founder animals and Mendellian contribution from non-founder animals. Two approaches were used to assess effects of inbreeding: overall inbreeding coefficients or partial inbreeding of the four main founder animals as possible covariates included in statistical models. Among the traits evaluated, there were effects on LMWLB, LMWLW and TLWW of inbreeding with there being a -8, -48, and -95 g decrease per 1% increase in inbreeding, respectively. Linear regressions of traits on partial inbreeding coefficients due to founder animals were of different magnitudes and ranged between -0.12 and +0.128. Heterogeneous contribution of founder animals to inbreeding depression occurred for LSB, LMWLB, and TLWW. These results indicate there was uneven distribution of recessive genetic composition among genomes of founder animals or differences in selection pressures on unfavorable alleles between different founder lines. The observed variation in founder-specific inbreeding depression indicates a small number of alleles with major effects are contributing to inbreeding depression.

Genetic diversity evolution of a sheep breed reintroduced after extinction: Tracing back Christopher Columbus' first imported sheep

New World's hair sheep breeds may genetically stem from West African introgression into established ecotypes of Spanish descent presumably extinct in the XIX Century. However, present Canary non-wooled breeds have presumably regressively resulted from the absorption of primitive individuals through Venezuelan descendants. We studied genetic diversity, structure, and evolution of the Canary hair sheep since its reintroduction in the 1950s. Demographic and genetic variability were evaluated using ENDOG (v4.8). Effective population size based on individual inbreeding rate was around one third higher than when based on individual coancestry rate. Nei's distances and equivalent subpopulations number indicated a highly-structured population. Although genetic diversity loss since the founder generations could be considered small, narrower pedigree bottlenecks could result from intraflock breeding policies and excessive contribution of few ancestors. Long generation intervals could be considered when reducing inbreeding. Wright's fixation statistics indicated slight interflock inbreeding. Pedigree completeness suggested genetic parameters were reliable, hence controlling inbreeding negative effects, could indeed, be crucial preserving these animal resources, consolidating the population in the archipelago after reintroduction.

On-farm net benefit of genotyping candidate female replacement cattle and sheep

The net benefit from investing in any technology is a function of the cost of implementation and the expected return in revenue. The objective of the present study was to quantify, using deterministic equations, the net monetary benefit from investing in genotyping of commercial females. Three case studies were presented reflecting dairy cows, beef cows and ewes based on Irish population parameters; sensitivity analyses were also performed. Parameters considered in the sensitivity analyses included the accuracy of genomic evaluations, replacement rate, proportion of female selection candidates retained as replacements, the cost of genotyping, the sire parentage error rate and the age of the female when it first gave birth. Results were presented as an annualised monetary net benefit over the lifetime of an individual, after discounting for the timing of expressions. In the base scenarios, the net benefit was greatest for dairy, followed by beef and then sheep. The net benefit improved as the reliability of the genomic evaluations improved and, in fact, a negative net benefit of genotyping was less frequent when the reliability of the genomic evaluations was high. The impact of a 10% point increase in genomic reliability was, however, greatest in sheep, followed by beef and then dairy. The net benefit of genotyping female selection candidates reduced as replacement rate increased. As genotyping costs increased, the net benefit reduced irrespective of the percentage of selection candidates kept, the replacement rate or even the population considered. Nonetheless, the association between the genotyping cost and the net benefit of genotyping differed by the percentage of selection candidates kept. Across all replacement rates evaluated, retaining 25% of the selection candidates resulted in the greatest net benefit when genotyping cost was low but the lowest net benefit when genotyping cost was high. Genotyping breakeven cost was non-linearly associated with the percentage of selection candidates retained, reaching a maximum when 50% of selection candidates were retained, irrespective of replacement rate, genomic reliability or the population. The genotyping breakeven cost was also non-linearly associated with replacement rate. The approaches outlined within provide the back-end framework for a decision support tool to quantify the net benefit of genotyping, once parameterised by the relevant population metrics.

Genetic structure and admixture in sheep from terminal breeds in the United States

Selection for performance in diverse production settings has resulted in variation across sheep breeds worldwide. Although sheep are an important species to the United States, the current genetic relationship among many terminal sire breeds is not well characterized. Suffolk, Hampshire, Shropshire and Oxford (terminal) and Rambouillet (dual purpose) sheep (n = 248) sampled from different flocks were genotyped using the Applied Biosystems Axiom Ovine Genotyping Array (50K), and additional Shropshire sheep (n = 26) using the Illumina Ovine SNP50 BeadChip. Relationships were investigated by calculating observed heterozygosity, inbreeding coefficients, eigenvalues, pairwise Wright’s F_ST estimates and an identity by state matrix. The mean observed heterozygosity for each breed ranged from 0.30 to 0.35 and was consistent with data reported in other US and Australian sheep. Suffolk from two different regions of the United States (Midwest and West) clustered separately in eigenvalue plots and the rectangular cladogram. Further, divergence was detected between Suffolk from different regions with Wright’s F_ST estimate. Shropshire animals showed the greatest divergence from other terminal breeds in this study. Admixture between breeds was examined using admixture, and based on cross‐validation estimates, the best fit number of populations (clusters) was K = 6. The greatest admixture was observed within Hampshire, Suffolk, and Shropshire breeds. When plotting eigenvalues, US terminal breeds clustered separately in comparison with sheep from other locations of the world. Understanding the genetic relationships between terminal sire breeds in sheep will inform us about the potential applicability of markers derived in one breed to other breeds based on relatedness.

Estimates of genetic parameters, genetic trends, and inbreeding in a crossbred dairy sheep research flock in the United States

For the past 2 decades, the Spooner Agriculture Research Station (ARS) of the University of Wisconsin-Madison operated the only dairy sheep research flock in North America. The objectives of the present study were to 1) obtain estimates of genetic parameters for lactation and reproductive traits in dairy ewes, 2) estimate the amount of genetic change in these traits over time, and 3) quantify the level of inbreeding in this flock over the last 20 yr. Multiple-trait repeatability models (MTRM) were used to analyze ewe traits through their first 6 parities. The first MTRM jointly analyzed milk (180-d-adjusted milk yield [180d MY]), fat (180-d-adjusted fat yield [180d FY]), and protein (180-d-adjusted protein yield [180d PY]) yields adjusted to 180 d of lactation; number of lambs born per ewe lambing (NLB); and lactation average test-day somatic cell score (LSCS). A second MTRM analyzed 180d MY, NLB, LSCS, and percentage milk fat (%F) and percentage milk protein (%P). The 3 yield traits were moderately heritable (0.26 to 0.32) and strongly genetically correlated (0.91 to 0.96). Percentage milk fat and %P were highly heritable (0.53 and 0.61, respectively) and moderately genetically correlated (0.61). Milk yield adjusted to 180 d was negatively genetically correlated with %F and %P (-0.31 and -0.34, respectively). Ewe prolificacy was not significantly ( > 0.67) genetically correlated with yield traits, %P, or LSCS but lowly negatively correlated with %F (-0.26). Lactation somatic cell score was unfavorably genetically correlated with yield traits (0.28 to 0.39) but not significantly ( > 0.09) correlated with %F, %P, and NLB. Within-trait multiple-trait models through the first 4 parities revealed that 180d MY, 180d FY, 180d PY, %F, and %P were strongly genetically correlated across parity (0.67 to 1.00). However, the genetic correlations across parity for NLB and LSCS were somewhat lower (0.51 to 0.96). Regressing predicted breeding values for 180d MY, without and with the addition of breed effects, on ewe year of birth revealed a positive genetic gain of 2.30 and 6.24 kg/yr, respectively, over the past 20 yr in this flock. Inbreeding coefficients of ewes with an extended pedigree ranged from 0.0 to 0.29, with an average of 0.07. To optimize genetic gains and avoid excessive inbreeding, the development of a national genetic improvement program should be a top priority for the growing dairy sheep industry.

Imputation of non-genotyped sheep from the genotypes of their mates and resulting progeny

Accurate genomic analyses are predicated on access to a large quantity of accurately genotyped and phenotyped animals. Because the cost of genotyping is often less than the cost of phenotyping, interest is increasing in generating genotypes for phenotyped animals. In some instances this may imply the requirement to genotype older animals with greater phenotypic information content. Biological material for these older informative animals may, however, no longer exist. The objective of the present study was to quantify the ability to impute 11 129 single nucleotide polymorphism (SNP) genotypes of non-genotyped animals (in this instance sires) from the genotypes of their progeny with or without including the genotypes of the progenys' dams (i.e. mates of the sire to be imputed). The impact on the accuracy of genotype imputation by including more progeny (and their dams') genotypes in the imputation reference population was also quantified. When genotypes of the dams were not available, genotypes of 41 sires with at least 15 genotyped progeny were used for the imputation; when genotypes of the dams were available, genotypes of 21 sires with at least 10 genotyped progeny were used for the imputation. Imputation was undertaken exploiting family and population level information. The mean and variability in the proportion of genotypes per individual that could not be imputed reduced as the number of progeny genotypes used per individual increased. Little improvement in the proportion of genotypes that could not be imputed was achieved once genotypes of seven progeny and their dams were used or genotypes of 11 progeny without their respective dam's genotypes were used. Mean imputation accuracy per individual (depicted by both concordance rates and correlation between true and imputed) increased with increasing progeny group size. Moreover, the range in mean imputation accuracy per individual reduced as more progeny genotypes were used in the imputation. If the genotype of the mate of the sire was also used, high accuracy of imputation (mean genotype concordance rate per individual of 0.988), with little additional benefit thereafter, was achieved with seven genotyped progeny. In the absence of genotypes on the dam, similar imputation accuracy could not be achieved even using genotypes on up to 15 progeny. Results therefore suggest, at least for the SNP density used in the present study, that it is possible to accurately impute the genotypes of a non-genotyped parent from the genotypes of its progeny and there is a benefit of also including the genotype of the sire's mate (i.e. dam of the progeny).

Analysis of population structure and genetic variability in Iranian buffaloes (Bubalus bubalis) using pedigree information

The objective of this study was to use pedigree analysis to evaluate the population structure, genetic variability and inbreeding in Iranian buffaloes. The analysis was based on the pedigree information of 42 285 buffaloes born from 549 sires and 6376 dams within 1697 herds. Pedigree information used in this study was collected during 1976 to 2012 by the Animal Breeding Centre of Iran. The CFC program was applied to calculate pedigree statistics and genetic structure analysis of the Iranian buffaloes. Also, the INBUPGF90 program was used for calculating regular inbreeding coefficients for individuals in the pedigree. The analysis of pedigree indicated that inbreeding coefficient ranged from 0% to 31% with an average of 3.42% and the trend of inbreeding was significantly positive over the years (P < 0.0001). Average coancestry was increased in recent years and overall generation interval was 6.62 years in Iranian buffaloes. Founder genome equivalent, founder equivalent, effective number of founders and effective number of non-founders were increased from 1976 to 2002, but their values decreased from 2002 onwards. A designed mating system to avoid inbreeding may be applied to this population of buffalo to maintain genetic diversity.

Inbreeding depression for economically important traits of Mazandaran native fowls

1. The objective was to investigate inbreeding depression for some economic traits of Mazandaran native fowls using data collected from 1992 to 2012 (21 generations) using a REML 2. The mean inbreeding coefficient (F) for the whole population and dams was 4.67% and 4.12%, respectively, and most of the inbred birds (75.79%) and inbred dams (72.58%) had F < 12.5%. 3. Individual and dam inbreeding trends were 0.55% and 0.53% per year. 4. Inbreeding depression for body weight at hatch, at 8 weeks and 12 weeks of age, age at sexual maturity, weight at sexual maturity, egg weight at 1st d of laying and average egg weight at 28, 30 and 32 weeks of laying due to a 1% increase in individual inbreeding were -0.11 g, -3.1 g, -1.3 g, 0.15 d, 0.59 g, -0.05 g and -0.03 g, respectively. 5. A 1% increase in maternal inbreeding resulted in a reduction of 0.06, 0.6 and 3.6 g in body weight at hatch, 8 weeks and 12 weeks of age.

Inbreeding effects on reproductive traits in Iranian Guilan sheep

The objective of the present study was to evaluate the effects of inbreeding on reproductive traits in Guilan sheep. Traits included were litter size at birth (LSB), litter size at weaning (LSW), litter mean weight per lamb born (LMWLB), litter mean weight per lamb weaned (LMWLW), total litter weight at birth (TLWB), and total litter weight at weaning (TLWW). Data and pedigree information used in this study were collected during 1994 to 2011 by the Agriculture Organization of Guilan Province in the north of Iran. Data were comprised of 14,534 records of lambs from 136 sires and 2021 dams. All the animals were grouped into three classes according to the inbreeding coefficients (F) obtained by their pedigree: the first class included non-inbred animals (F = 0%), and the second and third classes included inbred animals (0 < F ≤ 5% and F > 5%, respectively). The regression coefficients of LSB, LSW, LMWLB, LMWLW, TLWB, and TLWW on lamb inbreeding for a change of 1% in inbreeding were estimated to be 0.035 ± 0.0038 (P < 0.001), -0.029 ± 0.0077 (P < 0.05), -0.333 ± 0.009 (P < 0.001), -2.21 ± 0.071 (P < 0.001), -0.254 ± 0.013 (P < 0.001), and -1.95 ± 0.093 (P < 0.001), respectively. Effect of inbreeding on reproductive traits in Guilan sheep was very pronounced in the flock. The utilization of a program for planned mating system, in the present flock, suggested keeping the level of inbreeding under control and avoiding appearance of its detrimental effects.

Evaluation of inbreeding and genetic variability of five pig breeds in czech republic

The complex analysis of the pedigree records of Czech Landrace (CLA), Czech Large White-dam line (CLWd), Czech Large White-sire line (CLWs), Duroc (DC), and Pietrain (PN) was performed to determine trends of genetic diversity (GD), and to find the main sources of the GD loss. The total size of the pedigree was 132,365, 391,151, 32,913, 13,299, and 7,160 animals in CLA, CLWd, CLWs, DC, and PN, respectively. Animals born in the years 2011 through 2013 were assumed as the reference population. The average pedigree completeness index for one generation back was 95.9%, 97.4%, 91.2%, 89.8%, and 94.2% for appropriate breeds. Number of ancestors explaining 100% of gene pool was 186, 373, 125, 157, and 37 in CLA, CLWd, CLWs, DC, and PN, respectively. The relative proportion of inbred animals (58%, 58%, 54%, 47%, and 25%), the average inbreeding (2.7%, 1.4%, 2.5%, 3.6%, and 1.3%) and the average co-ancestry (3.1%, 1.6%, 3.3%, 4.2%, and 3.3%) were found over the past decade in analysed breeds. The expected inbreeding under random mating increased during the last 10 years in CLWs and PN and varied from 1.27% to 3.2%. The effective population size computed on the basis of inbreeding was 76, 74, 50, 35, and 83 in 2012 in CLA, CLWd, CLWs, DC, and PN, respectively. The shortest generation interval (1.45) was observed for CLWd in sire to son selection pathway. The longest generation interval obtained PN (1.95) in sire to daughter pathway. The average relative GD loss within last generation interval was 7.05%, 4.70%, 9.81%, 7.47%, and 10.46%, respectively. The relative proportion of GD loss due to genetic drift on total GD loss was 85.04%, 84.51%, 89.46%, 86.19%, and 83.68% in CLA, CLWd, CLWs, DC, and PN, respectively. All breeds were characterized by a high proportion of inbred animals, but the average inbreeding was low. The most vulnerable breeds to loss of GD are DC and PN. Therefore, a breeding program should be more oriented to prevent the increase of GD loss in these breeds.

Population genetic structure analysis and effect of inbreeding on body weights at different ages in Iranian Mehraban sheep

The objective of this study was to describe the population structure and inbreeding, and to quantify their effects on weights at different ages of Mehraban sheep in Iran. The analysis was based on the pedigree information of 26990 animals and 10278 body weight records from birth to yearling age. Data and pedigree information were collected during 1994 to 2011 by the breeding station of Mehraban sheep. The population structure was analyzed using the CFC program. Inbreeding of all animals was calculated by INBUPGF90 program. All animals were grouped into three classes according to their inbreeding coefficients: the first class included non-inbred animals (F=0); and the second and third classes included inbred animals (0<F<0.05 and F≥0.05, respectively). The average inbreeding in Mehraban sheep was 1.69%. Founder equivalent (f_e) values were estimated to be 4244, 3116 and 2965 during 1994-1999, 2000-2005 and 2006-2011, respectively. The effective population sizes (N_e) were 363, 5080 and 5740 during 1994-1999, 2000-2005 and 2006-2011, respectively. Generation interval was 2.15 years for this breed of sheep. Regression coefficients of birth weight, 3-month weight, 6-month weight and yearling weight on lamb inbreeding were estimated to be -6.34±0.69, -14.68±5.33, 48.00±9.43 and 98.65±15.65, respectively. Both positive and negative inbreeding effects were found in the current study. The utilization of a program for designed mating system, in the present flock, could be a suitable approach to keep the level of inbreeding under control.

Heterosis and genetic variation in the litter size of purebred and crossbred mink

Crosses between different mink lines from 3 Danish mink farms that use different breeding strategies were studied to estimate heterosis and variance components for litter size. The study was designed to analyze crosses between lines of the same color type, between different color types, and between animals originating from different farms. Effect of heterosis, color type, and variance components were estimated using Average information REML (AI-REML) algorithm implemented in the DMU package for analyzing multivariate mixed models. Females from 7 generations that gave birth to at least 1 offspring were analyzed and the effects of parity and production year were included in the analyses. Genetic trend and the proportions of the total variance explained by the effects of additive genetics (h2), common environment (due to repeated litters from the same female; c2), and dam of the female (granddame of the born litter; d2) were estimated. The results showed that mink of the Black color type potentially produced smaller litters compared to mink of the other studied color types. We found significant general maternal effect of heterosis for litter size. Analyses of specific heterosis showed a significant positive effect of crossing between lines of the same color type. Estimates of variance components revealed h2 levels for farm A, B, and C of 0.15, 0.06, and 0.09, respectively; thus litter size could be selected for in the future. The effect of common environment on litter size was also considerable, with c2 values of 0.005, 0.11, and 0.15 at farms A, B, and C, respectively. In conclusion, we recommend genetic selection as a means of increasing litter size in farmed mink.

A high throughput single nucleotide polymorphism multiplex assay for parentage assignment in New Zealand sheep

Accurate pedigree information is critical to animal breeding systems to ensure the highest rate of genetic gain and management of inbreeding. The abundance of available genomic data, together with development of high throughput genotyping platforms, means that single nucleotide polymorphisms (SNPs) are now the DNA marker of choice for genomic selection studies. Furthermore the superior qualities of SNPs compared to microsatellite markers allows for standardization between laboratories; a property that is crucial for developing an international set of markers for traceability studies. The objective of this study was to develop a high throughput SNP assay for use in the New Zealand sheep industry that gives accurate pedigree assignment and will allow a reduction in breeder input over lambing. This required two phases of development- firstly, a method of extracting quality DNA from ear-punch tissue performed in a high throughput cost efficient manner and secondly a SNP assay that has the ability to assign paternity to progeny resulting from mob mating. A likelihood based approach to infer paternity was used where sires with the highest LOD score (log of the ratio of the likelihood given parentage to likelihood given non-parentage) are assigned. An 84 “parentage SNP panel” was developed that assigned, on average, 99% of progeny to a sire in a problem where there were 3,000 progeny from 120 mob mated sires that included numerous half sib sires. In only 6% of those cases was there another sire with at least a 0.02 probability of paternity. Furthermore dam information (either recorded, or by genotyping possible dams) was absent, highlighting the SNP test’s suitability for paternity testing. Utilization of this parentage SNP assay will allow implementation of progeny testing into large commercial farms where the improved accuracy of sire assignment and genetic evaluations will increase genetic gain in the sheep industry.

A comparison of approaches to estimate the inbreeding coefficient and pairwise relatedness using genomic and pedigree data in a sheep population

Genome-wide SNP data provide a powerful tool to estimate pairwise relatedness among individuals and individual inbreeding coefficient. The aim of this study was to compare methods for estimating the two parameters in a Finnsheep population based on genome-wide SNPs and genealogies, separately. This study included ninety-nine Finnsheep in Finland that differed in coat colours (white, black, brown, grey, and black/white spotted) and were from a large pedigree comprising 319 119 animals. All the individuals were genotyped with the Illumina Ovine SNP50K BeadChip by the International Sheep Genomics Consortium. We identified three genetic subpopulations that corresponded approximately with the coat colours (grey, white, and black and brown) of the sheep. We detected a significant subdivision among the colour types (F(ST) = 5.4%, P<0.05). We applied robust algorithms for the genomic estimation of individual inbreeding (F(SNP)) and pairwise relatedness (Φ(SNP)) as implemented in the programs KING and PLINK, respectively. Estimates of the two parameters from pedigrees (F(PED) and Φ(PED)) were computed using the RelaX2 program. Values of the two parameters estimated from genomic and genealogical data were mostly consistent, in particular for the highly inbred animals (e.g. inbreeding coefficient F>0.0625) and pairs of closely related animals (e.g. the full- or half-sibs). Nevertheless, we also detected differences in the two parameters between the approaches, particularly with respect to the grey Finnsheep. This could be due to the smaller sample size and relative incompleteness of the pedigree for them.We conclude that the genome-wide genomic data will provide useful information on a per sample or pairwise-samples basis in cases of complex genealogies or in the absence of genealogical data.

Epistasis for founder-specific inbreeding depression in rabbits

Inbreeding depression is a topic of main interest in experimental and domestic species, although previous studies simplified this genetically complex effect to the linear (or quadratic) regression coefficient linked to the inbreeding coefficient of each individual or, in more recent studies, to founder-specific inbreeding coefficients. Going beyond generalizing to these traditional scenarios, our research focused on the analysis of gene-by-gene interactions leading to epistasis for inbreeding depression effects. Under a Bayesian context, inbreeding depression effects were evaluated for weaning weight (WW) in a commercial rabbit population founded from 4 bucks and 1 doe (MARET population). Founder-specific inbreeding depression effects for the 4 bucks ranged between -81.1 and 38.3 g for each 1% inbreeding. More interestingly, 2 epistatic interactions between the partial inbreeding coefficient of 2 bucks were also significant and negative, showing a -1.9 and -1.0 g reduction on WW. These results provide the first evidence of epistatic inbreeding depression phenomena in domestic species, emphasizing the complexity of the genetic architecture in mammals.

Genetic diversity and population structure of the synthetic Pannon White rabbit revealed by pedigree analyses

Demographic history, current status, and efficiency of the mating strategy were analyzed using the pedigree of Pannon White (PW) rabbits born between 1992 and 2007. Potential accumulation of detrimental effects and loss of genetic diversity were also considered. Calculations and estimates were done most often for rabbits born in 2007, whereas other reference populations (REFPOPXXXX) were specified explicitly. The pedigree contained 4,749 individuals and 580 founders, and its completeness was 82.1% up to 10 and 94.5% up to 5 generations, respectively. Generation intervals through different pathways averaged 1.2 yr. When adjusted to the pedigree completeness, the amount of inbreeding (F(i)) of rabbits was comparable (5.54%) with that of other livestock populations, whereas the 10 (30) founders contributing the most to inbreeding explained a large part of the population inbreeding [i.e., 42.24% (73.18%)]. The ancestral inbreeding coefficient of REFPOP2004 (10.67%) was one-half that of REFPOP2007 (20.66%), showing its strong dependence on pedigree length. Family variance, inbreeding, and realized effective population size were 84.18 (REFPOP2006; this variable could not be calculated for the last year examined), 37.19, and 91.08, respectively. The effective numbers of ancestors, founders, and founder genomes were 48, 26, and 7.33, respectively. Although the circular mating scheme applied was generally effective, the large accumulated reduction in genetic variability indicates the need to revise and improve the current breeding strategy.

Accounting for additive genetic mutations on litter size in Ripollesa sheep

Little is known about mutational variability in livestock, among which only a few mutations with relatively large effects have been reported. In this manuscript, mutational variability was analyzed in 1,765 litter size records from 404 Ripollesa ewes to characterize the magnitude of this genetic source of variation and check the suitability of including mutational effects in genetic evaluations of this breed. Threshold animal models accounting for additive genetic mutations were preferred to models without mutational contributions, with an average difference in the deviance information criterion of more than 5 units. Moreover, the statistical relevance of the additive genetic mutation term was checked through a Bayes factor approach, which showed that the models with mutational variability were 8.5 to 22.7 times more probable than the others. The mutational heritability (percentage of the phenotypic variance accounted for by mutational variance) was 0.6 or 0.9%, depending on whether genetic dominance effects were accounted for by the analytical model. The inclusion of mutational effects in the genetic model for evaluating litter size in Ripollesa ewes called for some minor modifications in the genetic merit order of the individuals evaluated, which suggested that the continuous uploading of new additive mutations could be taken into account to optimize the selection scheme. This study is the first attempt to estimate mutational variances in a livestock species and thereby contribute to better characterization of the genetic background of productive traits of interest.

Correlated inbreeding among relatives: occurrence, magnitude, and implications

Understanding the magnitude and causes of genetic and phenotypic resemblance among relatives is key to understanding evolutionary processes. Contrary to basic expectation, individual coefficients of inbreeding (f) were recently hypothesized to be intrinsically correlated across parents and offspring in structured populations, potentially creating an additional source of phenotypic resemblance in traits that show inbreeding depression. To test this hypothesis, we used individual-based simulations to quantify the parent-offspring correlations in f arising under random mating in populations of different size, immigration rate, and mating system. Parent-offspring correlations in f were typically positive (median r approximately 0.2-0.4) in relatively small and isolated populations. Relatively inbred parents therefore produced relatively inbred offspring on average, although the magnitude of this effect varied considerably among replicate populations. Correlations were higher given more generations of random mating, greater variance in reproductive success, polygynous rather than monogamous mating, and for midparent-offspring rather than parent-offspring relationships. Furthermore, f was also positively correlated across half-siblings, and closer relatives had more similar inbreeding coefficients across entire generations. Such intrinsic resemblance in f among relatives could provide an additional genetic benefit of mate choice and bias quantitative genetic analyses that do not account for correlated inbreeding depression.

Skew distribution of founder-specific inbreeding depression effects on the longevity of Landrace sows

Inbreeding is a biological phenomenon of special relevance in domestic species in which its influence has been typically associated with reductions in animal fitness. Nevertheless, recent research suggests substantial founder-specific variability in terms of inbreeding depression on some productive traits, although this centred on a very reduced number of founders. This research focuses on the modelling of founder-specific inbreeding depression (FSID) effects from a more general point of view, characterizing the expected distribution of FSID effects on sow longevity. Under a change-point Weibull proportional hazards model solved through Bayesian inference, a skew-normala prioridistribution for the FSID effects of 19 founders was assumed. In terms of the deviance information criterion, this model was clearly preferred to other prior distributions for FSID effects as well as to a standard analysis of the overall inbreeding depression effect, although all models were consistent with an overall negative genetic effect of inbreeding on sow longevity. The joint posterior distribution of FSID effects showed a skewed pattern with substantial right-tail overexpression, in which the mean (0·036), mode (0·034), S.D. (0·032) and asymmetry parameter (0·045) reported a higher incidence of positive estimates (85·2%) with an unfavourable effect on sow longevity. The founder with the worst inbreeding depression effect reduced sow longevity by 32 days for 1% or 167 days for 10% partial inbreeding. As a whole, our analyses highlighted substantial variability in FSID effects, with unfavourable, neutral and even favourable influences on sow longevity. This heterogeneity could be related to an uneven distribution of the recessive deleterious genetic load among founder genomes, and also with the different selection pressures applied to each founder line. The implementation of skew-normal priors also proved an appealing way to bypass the strict scenario imposed by the standard symmetric-Gaussian distribution, allowing right- and left-tail overexpression as well as non-zero modal estimates.

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.