Analysis of inbreeding depression in the first litter size of mice in a long-term selection experiment with respect to the age of the inbreeding

Pedigree diversity and implications for genetic selection of Katahdin sheep

The Katahdin hair breed gained popularity in the United States as low input and prolific, with a propensity to exhibit parasite resistance. With the introduction of genomically enhanced estimated breeding values (GEBV) to the Katahdin genetic evaluation, defining the diversity present in the breed is pertinent. Utilizing pedigree records (n = 92,030) from 1984 to 2019 from the National Sheep Improvement Program, our objectives were to (i) estimate the completeness and quality of the pedigree, (ii) calculate diversity statistics for the whole pedigree and relevant reference subpopulations and (iii) assess the impact of current diversity on genomic selection. Reference 1 was Katahdins born from 2017 to 2019 (n = 23,494), while reference 2 was a subset with at least three generations of Katahdin ancestry (n = 9327). The completeness of the whole pedigree, and the pedigrees of reference 1 and reference 2, were above 50% through the fourth, fifth and seventh generation of ancestors, respectively. Effective population size (Ne) averaged 111 animals with a range from 42.2 to 451.0. The average generation interval was 2.9 years for the whole pedigree and reference 1, and 2.8 years for reference 2. The mean individual inbreeding and average relatedness coefficients were 1.62% and 0.91%, 1.74% and 0.90% and 2.94% and 1.46% for the whole pedigree, reference 1, and reference 2, respectively. There were over 300 effective founders in the whole pedigree and reference 1, with 169 in reference 2. Effective number of ancestors were over 150 for the whole pedigree and reference 1, while there were 67 for reference 2. Prediction accuracies increased as the reference population grew from 1k to 7.5k and plateaued at 15k animals. Given the large number of founders and ancestors contributing to the base genetic variation in the breed, the Ne is sufficient to maintain diversity while achieving progress with selection. Stable low rates of inbreeding and relatedness suggest that incorporating genetic conservation in breeding decisions is currently not of high priority. Current Ne suggests that with limited genotyping, high levels of accuracy for genomic prediction can be achieved. However, intense selection on GEBV may cause loss of genetic diversity long term.

Characterizing and Eliminating the Inbreeding Load

The authors evaluated the relevant literature related to purging, which is the interaction between selection and inbreeding in which the population may eliminate its inbreeding load at least partially. According to the relevant literature, the inbreeding load and the process of purging were evaluated via pedigree methods based on ancestral inbreeding, the inbreeding-purging model, and expressed opportunity of purging, along with genomic methods. Most ancestral inbreeding-related studies were performed in zoos, where only a small proportion of the studied populations show signs of purging. The inbreeding-purging model was developed with Drosophila, and it was used to evaluate different zoo ungulates and Pannon white rabbits. Purging was detected in both studies. The expressed opportunity of purging was applied in Jersey cattle and Pannon white rabbits. In the Jersey cattle, it had an effect of 12.6% for fitness, while in the Pannon white rabbits, the inbreeding load was between 40% and 80% of its original value. The genomic studies also signalled purging, but they also made it clear that, contrary to the detected purging, the evaluated populations still suffered from inbreeding depression. Therefore, especially for domesticated animals, it can be concluded that deliberate inbreeding with the purpose of generating purging is not advocated.

The effect of new and ancestral inbreeding on milk production traits in Iranian Holstein cattle

Inbreeding depression, the reduction of fitness and performance, is due to an increase in the mating of related individuals. Based on the purge hypothesis, inbreeding and breeding over generations reduce the effect of deleterious alleles responsible for inbreeding depression. Thus, recent inbreeding is assumed to be more harmful than ancestral inbreeding. This study aimed at evaluating the effects of new and ancestral inbreeding on milk, fat and protein production in Iranian Holstein cattle. The secondary objective was to examine the changes in predicted breeding values when the inbreeding effect was included in the model's analysis. To this end, inbreeding coefficients were calculated using the pedigree of 2,394,517 Holstein cattle to achieve these goals. In addition, 419,132 records of milk, fat and protein yields of first parity cows were collected to assess inbreeding depression and breeding values. The average inbreeding coefficients were 0.83% and 1.68% for the whole population and the inbred animals, respectively. A 1% increase in classical pedigree-based inbreeding coefficient was associated with a decrease of 11.99 kg in milk, 0.39 kg in fat and 0.29 kg in protein. The effect of ancestral inbreeding was more detrimental to performance traits than the effect of new inbreeding. This result contradicted the hypothesis of purging. By including the inbreeding coefficient in the model, the rank of animals remained unchanged, but the average predicted breeding values increased. In general, inbreeding depression was observed in Iranian Holstein cows; however, no evidence of purging was observed. The average of inbreeding coefficients was not high in this population, although accounting for inbreeding coefficients in the analytical model did significantly increase the predicted breeding values. It is recommended that the analytical model incorporate the inbreeding coefficient to improve the accuracy of genetic evaluation. In future studies, inbreeding depression should be assessed using genomic data for performance and reproduction traits.

Inbreeding depression on growth and prolificacy traits in two lines of rabbit

Inbreeding depression in closed populations impairs animal fitness, health, and productivity. However, not all inbreeding is expected to be equally damaging. Recent inbreeding is thought to be more harmful than ancient inbreeding because selection decreases the frequency of unfavourable alleles with time. Accordingly, selection efficiency is improved by inbreeding in a process called purging. This research aimed to quantify inbreeding depression on growth and prolificacy traits in two lines of rabbits selected for just one growth (Caldes line) or prolificacy (Prat line) trait, and also to find some evidence of purging of deleterious alleles by selection. Caldes line comprised 51 generations and 124,371 animals in the pedigree. Prat line comprised 34 generations and 161,039 animals in the pedigree. The effects of old, intermediate, and new inbreeding (Fold, Fint, and Fnew), as well as total cumulated classical inbreeding (F) and 3 measurements of ancestral inbreeding (AHC, Fa.K, and Fa.B) were estimated for average daily gain (ADG), slaughter weight (SW), weaning weight (WW), born alive (BA), the total number of kits (NT), and the number of weaned kits (NW). There was a clear inbreeding depression for all growth and prolificacy traits in the Caldes line (-7.19 g/d, -0.45 kg, -0.25 kg, -6 kits, -4 kits, and -4 kits per unit of increase in F for ADG, SW, WW, BA, NT, and NW, respectively) and also in Prat line (-7.48 g/d, -0.31 kg, -0.11 kg, -4 kits, -5 kits, and -4 kits per unit of increase in F for ADG, SW, WW, BA, NT, and NW, respectively). The inbreeding partition appears to be a reliable alternative for assessing inbreeding depression and purging. Thus, for example, in the Caldes line and for ADG the regression coefficients were -7.61, -5.41, and 7.76 g/d per unit of increase in Fnew, Fint, and Fold, respectively. In addition, AHC and Fa.B may provide more accurate evidence of purging than Fa.K. This study confirms the existence of inbreeding depression for growth and prolificacy traits in both lines of rabbits and shows evidence of purging of deleterious recessive alleles involved both in growth and prolificacy, independently of the selection criteria established in the line.

Breeding schemes with optimum-contribution selection or truncation selection for beef cattle destined for use on dairy females

We tested the hypothesis that the size of a beef cattle population destined for use on dairy females is smaller under optimum-contribution selection (OCS) than under truncation selection (TRS) at the same genetic gain (ΔG) and the same rate of inbreeding (ΔF). We used stochastic simulation to estimate true ΔG realized at a 0.005 ΔF in breeding schemes with OCS or TRS. The schemes for the beef cattle population also differed in the number of purebred offspring per dam and the total number of purebred offspring per generation. Dams of the next generation were exclusively selected among the one-year-old heifers. All dams were donors for embryo transfer and produced a maximum of 5 or 10 offspring. The total number of purebred offspring per generation was: 400, 800, 1,600 or 4,000 calves, and it was used as a measure of population size. Rate of inbreeding was predicted and controlled using pedigree relationships. Each OCS (TRS) scheme was simulated for 10 discrete generations and replicated 100 (200) times. The OCS scheme and the TRS scheme with a maximum of 10 offspring per dam required approximately 783 and 1,257 purebred offspring per generation to realize a true ΔG of €14 and a ΔF of 0.005 per generation. Schemes with a maximum of 5 offspring per dam required more purebred offspring per generation to realize a similar true ΔG and a similar ΔF. Our results show that OCS and multiple ovulation and embryo transfer act on selection intensity through different mechanisms to achieve fewer selection candidates and fewer selected sires and dams than under TRS at the same ΔG and a fixed ΔF. Therefore, we advocate the use of a breeding scheme with OCS and multiple ovulation and embryo transfer for beef cattle destined for use on dairy females because it is favorable both from an economic perspective and a carbon footprint perspective.

Grid search approach to discriminate between old and recent inbreeding using phenotypic, pedigree and genomic information

Background

Although inbreeding caused by the mating of animals related through a recent common ancestor is expected to have more harmful effects on phenotypes than ancient inbreeding (old inbreeding), estimating these effects requires a clear definition of recent (new) and ancient (old) inbreeding. Several methods have been proposed to classify inbreeding using pedigree and genomic data. Unfortunately, these methods are largely based on heuristic criteria such as the number of generations from a common ancestor or length of runs of homozygosity (ROH) segments. To mitigate these deficiencies, this study aimed to develop a method to classify pedigree and genomic inbreeding into recent and ancient classes based on a grid search algorithm driven by the assumption that new inbreeding tends to have a more pronounced detrimental effect on traits. The proposed method was tested using a cattle population characterized by a deep pedigree.

Results

Effects of recent and ancient inbreeding were assessed on four growth traits (birth, weaning and yearling weights and average daily gain). Thresholds to classify inbreeding into recent and ancient classes were trait-specific and varied across traits and sources of information. Using pedigree information, inbreeding generated in the last 10 to 11 generations was considered as recent. When genomic information (ROH) was used, thresholds ranged between four to seven generations, indicating, in part, the ability of ROH segments to characterize the harmful effects of inbreeding in shorter periods of time. Nevertheless, using the proposed classification method, the discrimination between new and old inbreeding was less robust when ROH segments were used compared to pedigree. Using several model comparison criteria, the proposed approach was generally better than existing methods. Recent inbreeding appeared to be more harmful across the growth traits analyzed. However, both new and old inbreeding were found to be associated with decreased yearling weight and average daily gain.

Conclusions

The proposed method provided a more objective quantitative approach for the classification of inbreeding. The proposed method detected a clear divergence in the effects of old and recent inbreeding using pedigree data and it was superior to existing methods for all analyzed traits. Using ROH data, the discrimination between old and recent inbreeding was less clear and the proposed method was superior to existing approaches for two out of the four analyzed traits. Deleterious effects of recent inbreeding were detected sooner (fewer generations) using genomic information than pedigree. Difference in the results using genomic and pedigree information could be due to the dissimilarity in the number of generations to a common ancestor. Additionally, the uncertainty associated with the identification of ROH segments and associated inbreeding could have an effect on the results. Potential biases in the estimation of inbreeding effects may occur when new and old inbreeding are discriminated based on arbitrary thresholds. To minimize the impact of inbreeding, mating designs should take the different inbreeding origins into consideration.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07872-z.

Genomic and pedigree estimation of inbreeding depression for semen traits in the Basco-Béarnaise dairy sheep breed

Inbreeding depression is associated with a decrease in performance and fitness of the animals. The goal of this study was to evaluate pedigree-based and genomic methods to estimate the level of inbreeding and inbreeding depression for 3 semen traits (volume, concentration, and motility score) in the Basco-Béarnaise sheep breed. Data comprised 16,196 (or 15,071) phenotypic records from 620 rams (of which 533 rams had genotypes of 36,464 SNPs). The pedigree included 8,266 animals, composed of the 620 rams and their ancestors. The number of equivalent complete generations for the 620 rams was 7.04. Inbreeding coefficients were estimated using genomic and pedigree-based information. Genomic inbreeding coefficients were estimated from individual SNP and using segments of homozygous SNP (runs of homozygosity, ROH). Short ROH are of old origin, whereas long ROH are due to recent inbreeding. Considering that the equivalent number of generations in Basco-Béarnaise was 6, inbreeding coefficients for ROH with a length >4 Mb refer to all (recent + old) inbreeding, those with a length >17 Mb correspond to recent inbreeding, and the difference between them indicates old inbreeding. Pedigree-based inbreeding coefficients were also estimated classically, or accounting for nonzero relationships for unknown parents, or including metafounder relationships (estimated using markers) to account for missing pedigree information. Finally, inbreeding coefficients combining genotyped and nongenotyped animal information were computed from matrix H of the single-step approach, also including metafounders. Inbreeding depression was estimated differently depending on the approach used to compute inbreeding coefficients. These 8 estimators of inbreeding coefficients were included as covariates in different animal models. No inbreeding depression was detected for sperm volume or sperm concentration. Inbreeding depression was significant for the motility of spermatozoa. The effect of old and recent inbreeding on motility was null and negative, respectively, demonstrating the existence of purging by selection of deleterious recessive alleles affecting motility. A 10% increase in inbreeding would result in a reduction in mean motility ranging between 0.09 and 0.22 points in the score (from 0 to 5). Motility is unfavorably affected by increasing recent inbreeding but the impact is very small. Runs of homozygosity and metafounders allow us to accurately estimate inbreeding depression and detect recent inbreeding.

Colony aging affects the reproductive performance of Swiss Webster females used as recipients for embryo transfer

The objective was to evaluate the influence of colony aging in a Swiss Webster (SW) outbred stock used as recipients for embryo transfer. In the first study, a retrospective analysis was performed throughout several generations during a 38-month period in 2,398 embryos transferred to 108 SW recipients. A decrease in the percentage of live pups from transferred embryos was found at the end of the period. Impairment occurred due to the incidence of maternal cannibalism that increased from 0% to 67-100% (P<0.05), while pregnancy rate (pregnant/transferred recipients) and number of pups per delivered female were not affected throughout the period (P=NS). A following study was carried out to compare the reproductive performance of SW stock vs. B6D2F1 hybrid females in a 5-year interval. The study was conducted on a total of 893 embryos transferred to 40 females (20 SW and 20 B6D2F1) in Year #1, and 514 embryos transferred to 30 females (15 SW and 15 B6D2F1) in Year #5. No cases of maternal cannibalism were found on Year #1 in any of the strains (0/10 and 0/10). However, an incidence of 44,4% (4/9) was seen on Year #5 for SW, while for B6D2F1 the incidence was 0% (0/12) (P<0.05). Further examination of the uterus showed endometrial cysts and abnormal implantation sites in SW on Year #5 but not in B6D2F1 females. In conclusion, this study reports an impairment of the reproductive performance of an early aged SW outbred stock colony mainly due to the occurrence of maternal cannibalism. This finding has important implications for embryo transfer programs conducted in mouse facilities.

Phenotypically Selective Genotyping Realizes More Genetic Gains in a Rainbow Trout Breeding Program in the Presence of Genotype-by-Environment Interactions

Selective genotyping of phenotypically superior animals may lead to bias and less accurate genomic breeding values (GEBV). Performing selective genotyping based on phenotypes measured in the breeding environment (B) is not necessarily a good strategy when the aim of a breeding program is to improve animals’ performance in the commercial environment (C). Our simulation study compared different genotyping strategies for selection candidates and for fish in C in a breeding program for rainbow trout in the presence of genotype-by-environment interactions when the program had limited genotyping resources and unregistered pedigrees of individuals. For the reference population, selective genotyping of top and bottom individuals in C based on phenotypes measured in C led to the highest genetic gains, followed by random genotyping and then selective genotyping of top individuals in C. For selection candidates, selective genotyping of top individuals in B based on phenotypes measured in B led to the highest genetic gains, followed by selective genotyping of top and bottom individuals and then random genotyping. Selective genotyping led to bias in predicting GEBV. However, in scenarios that used selective genotyping of top fish in B and random genotyping of fish in C, predictions of GEBV were unbiased, with genetic correlations of 0.2 and 0.5 between traits measured in B and C. Estimates of variance components were sensitive to genotyping strategy, with an overestimation of the variance with selective genotyping of top and bottom fish and an underestimation of the variance with selective genotyping of top fish. Unbiased estimates of variance components were obtained when fish in B and C were genotyped at random. In conclusion, we recommend phenotypic genotyping of top and bottom fish in C and top fish in B for the purpose of selecting breeding animals and random genotyping of individuals in B and C for the purpose of estimating variance components when a genomic breeding program for rainbow trout aims to improve animals’ performance in C.

Effect of recent and ancient inbreeding on production and fertility traits in Canadian Holsteins

BACKGROUND

Phenotypic performances of livestock animals decline with increasing levels of inbreeding, however, the noticeable decline known as inbreeding depression, may not be due only to the total level of inbreeding, but rather could be distinctly associated with more recent or more ancient inbreeding. Therefore, splitting inbreeding into different age classes could help in assessing detrimental effects of different ages of inbreeding. Hence, this study sought to investigate the effect of recent and ancient inbreeding on production and fertility traits in Canadian Holstein cattle with both pedigree and genomic records. Furthermore, inbreeding coefficients were estimated using traditional pedigree measure (F_PED) and genomic measures using segment based (F_ROH) and marker-by-marker (F_GRM) based approaches.

RESULTS

Inbreeding depression was found for all production and most fertility traits, for example, every 1% increase in F_PED, F_ROH and F_GRM was observed to cause a - 44.71, - 40.48 and - 48.72 kg reduction in 305-day milk yield (MY), respectively. Similarly, an extension in first service to conception (FSTC) of 0.29, 0.24 and 0.31 day in heifers was found for every 1% increase in F_PED, F_ROH and F_GRM, respectively. Fertility traits that did not show significant depression were observed to move in an unfavorable direction over time. Splitting both pedigree and genomic inbreeding into age classes resulted in recent age classes showing more detrimental inbreeding effects, while more distant age classes caused more favorable effects. For example, a - 1.56 kg loss in 305-day protein yield (PY) was observed for every 1% increase in the most recent pedigree age class, whereas a 1.33 kg gain was found per 1% increase in the most distant pedigree age class.

CONCLUSIONS

Inbreeding depression was observed for production and fertility traits. In general, recent inbreeding had unfavorable effects, while ancestral inbreeding had favorable effects. Given that more negative effects were estimated from recent inbreeding when compared to ancient inbreeding suggests that recent inbreeding should be the primary focus of selection programs. Also, further work to identify specific recent homozygous regions negatively associated with phenotypic traits could be investigated.

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.

Pedigree relationships to control inbreeding in optimum-contribution selection realise more genetic gain than genomic relationships

BACKGROUND

We tested the premise that optimum-contribution selection with pedigree relationships to control inbreeding (POCS) realises at least as much true genetic gain as optimum-contribution selection with genomic relationships (GOCS) at the same rate of true inbreeding.

METHODS

We used stochastic simulation to estimate rates of true genetic gain realised by POCS and GOCS at a 0.01 rate of true inbreeding in three breeding schemes with best linear unbiased predictions of breeding values based on pedigree (PBLUP) and genomic (GBLUP) information. The three breeding schemes differed in number of matings and litter size. Selection was for a single trait with a heritability of 0.2. The trait was controlled by 7702 biallelic quantitative-trait loci (QTL) that were distributed across a 30-M genome. The genome contained 54,218 biallelic markers that were used in GOCS and GBLUP. A total of 6012 identity-by-descent loci were placed across the genome in base populations. Unique alleles at these loci were used to calculate rates of true inbreeding. Breeding schemes were run for 10 discrete generations. Selection candidates were genotyped and phenotyped before selection.

RESULTS

POCS realised more true genetic gain than GOCS at a 0.01 rate of true inbreeding in all combinations of breeding scheme and prediction method. POCS realised 14 to 33% more true genetic gain than GOCS with PBLUP in the three breeding schemes. It realised 1.5 to 5.7% more true genetic gain than GOCS with GBLUP.

CONCLUSIONS

POCS realised more true genetic gain than GOCS because it managed expected genetic drift without restricting selection at QTL. By contrast, GOCS penalised changes in allele frequencies at markers that were generated by genetic drift and selection. Because these marker alleles were in linkage disequilibrium with QTL alleles, GOCS restricted changes in allele frequencies at QTL. This provides little incentive to use GOCS and highlights that we have more to learn before we can control inbreeding using genomic relationships in selective-breeding schemes. Until we can do so, POCS remains a worthy method of optimum-contribution selection because it realises more true genetic gain than GOCS at the same rate of true inbreeding.

Inbreeding depression in line 1 Hereford cattle population using pedigree and genomic information

This study aimed at assessing inbreeding and its effect on growth and fertility traits using the longtime closed line 1 Hereford cattle population. Inbreeding was estimated based on pedigree (FPED) and genomic information. For the latter, three estimates were derived based on the diagonal elements of the genomic relationship matrix using estimated (FGRM) or fixed (FGRM0.5) minor allele frequencies or runs of homozygosity (ROH) (FROH). A pedigree containing 10,186 animals was used to calculate FPED. Genomic inbreeding was evaluated using 785 animals genotyped for 30,810 SNP. Traits analyzed were birth weight (BWT), weaning weight (WWT), yearling weight (YWT), ADG, and age at first calving (AFC). The number of ROH per animal ranged between 6 and 119 segments with an average of 83. The shortest and longest segments were 1.36 and 64.86 Mb long, respectively, reflecting both ancient and recent inbreeding occurring in the last 30 to 40 generations. The average inbreeding was 29.2%, 16.1%, 30.2%, and 22.9% for FPED, FGRM, FGRM0.5, and FROH, respectively. FROH provided the highest correlations with FPED (r = 0.66). Across paternal half-sib families, with minimal variation in FPED, there were substantial variations in their genomic inbreeding. Inbreeding depression analyses showed that a 1% increase in an animal's FPED resulted in a decrease of 1.20 kg, 2.03 kg, and 0.004 kg/d in WWT, YWT, and ADG, respectively. Maternal inbreeding showed significantly negative effects on progeny growth performance. AFC increased by 1.4 and 0.8 d for each 1% increase in FPED of the cow and her dam, respectively. Using genomic inbreeding, similar impact on growth traits was observed although the magnitude of the effect varied between methods. Across all genomic measures, WWT, YWT, and ADG decreased by 0.21 to 0.53 kg, 0.46 to 1.13 kg, and 0.002 to 0.006 kg/d for each 1% increase in genomic inbreeding, respectively. Four chromosomes (9, 12, 17, and 27) were identified to have a significant association between their homozygosity (FROH-CHR) and growth traits. Variability in genomic inbreeding could be useful when deciding between full and half-sib selection candidates. Despite the high level of inbreeding in this study, its negative impact on growth performance was not as severe as expected, which may be attributed to the purging of the deleterious alleles due to natural or artificial selection over time.

Founder-specific inbreeding depression affects racing performance in Thoroughbred horses

The Thoroughbred horse has played an important role in both sporting and economic aspects of society since the establishment of the breed in the 1700s. The extensive pedigree and phenotypic information available for the Thoroughbred horse population provides a unique opportunity to examine the effects of 300 years of selective breeding on genetic load. By analysing the relationship between inbreeding and racing performance of 135,572 individuals, we found that selective breeding has not efficiently alleviated the Australian Thoroughbred population of its genetic load. However, we found evidence for purging in the population that might have improved racing performance over time. Over 80% of inbreeding in the contemporary population is accounted for by a small number of ancestors from the foundation of the breed. Inbreeding to these ancestors has variable effects on fitness, demonstrating that an understanding of the distribution of genetic load is important in improving the phenotypic value of a population in the future. Our findings hold value not only for Thoroughbred and other domestic breeds, but also for small and endangered populations where such comprehensive information is not available.

Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein-Friesian breeding program from 1986 to 2015

Background

In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015.

Methods

Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time.

Results

Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods.

Conclusions

Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0385-y) contains supplementary material, which is available to authorized users.

Age-based partitioning of individual genomic inbreeding levels in Belgian Blue cattle

Background

Inbreeding coefficients can be estimated either from pedigree data or from genomic data, and with genomic data, they are either global or local (when the linkage map is used). Recently, we developed a new hidden Markov model (HMM) that estimates probabilities of homozygosity-by-descent (HBD) at each marker position and automatically partitions autozygosity in multiple age-related classes (based on the length of HBD segments). Our objectives were to: (1) characterize inbreeding with our model in an intensively selected population such as the Belgian Blue Beef (BBB) cattle breed; (2) compare the properties of the model at different marker densities; and (3) compare our model with other methods.

Results

When using 600 K single nucleotide polymorphisms (SNPs), the inbreeding coefficient (probability of sampling an HBD locus in an individual) was on average 0.303 (ranging from 0.258 to 0.375). HBD-classes associated to historical ancestors (with small segments ≤ 200 kb) accounted for 21.6% of the genome length (71.4% of the total length of the genome in HBD segments), whereas classes associated to more recent ancestors accounted for only 22.6% of the total length of the genome in HBD segments. However, these recent classes presented more individual variation than more ancient classes. Although inbreeding coefficients obtained with low SNP densities (7 and 32 K) were much lower (0.060 and 0.093), they were highly correlated with those obtained at higher density (r = 0.934 and 0.975, respectively), indicating that they captured most of the individual variation. At higher SNP density, smaller HBD segments are identified and, thus, more past generations can be explored. We observed very high correlations between our estimates and those based on homozygosity (r = 0.95) or on runs-of-homozygosity (r = 0.95). As expected, pedigree-based estimates were mainly correlated with recent HBD-classes (r = 0.56).

Conclusions

Although we observed high levels of autozygosity associated with small HBD segments in BBB cattle, recent inbreeding accounted for most of the individual variation. Recent autozygosity can be captured efficiently with low-density SNP arrays and relatively simple models (e.g., two HBD classes). The HMM framework provides local HBD probabilities that are still useful at lower SNP densities.

Electronic supplementary material

The online version of this article (10.1186/s12711-017-0370-x) contains supplementary material, which is available to authorized users.

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.

A model-based approach to characterize individual inbreeding at both global and local genomic scales

Inbreeding results from the mating of related individuals and may be associated with reduced fitness because it brings together deleterious variants in one individual. In general, inbreeding is estimated with respect to an arbitrary base population consisting of ancestors that are assumed unrelated. We herein propose a model-based approach to estimate and characterize individual inbreeding at both global and local genomic scales by assuming the individual genome is a mosaic of homozygous-by-descent (HBD) and non-HBD segments. The HBD segments may originate from ancestors tracing back to different periods in the past defining distinct age-related classes. The lengths of the HBD segments are exponentially distributed with class-specific parameters reflecting that inbreeding of older origin generates on average shorter stretches of observed homozygous markers. The model is implemented in a hidden Markov model framework that uses marker allele frequencies, genetic distances, genotyping error rates and the sequences of observed genotypes. Note that genotyping errors, low-fold sequencing or genotype-by-sequencing data are easily accommodated under this framework. Based on simulations under the inference model, we show that the genomewide inbreeding coefficients and the parameters of the model are accurately estimated. In addition, when several inbreeding classes are simulated, the model captures them if their ages are sufficiently different. Complementary analyses, either on data sets simulated under more realistic models or on human, dog and sheep real data, illustrate the range of applications of the approach and how it can reveal recent demographic histories among populations (e.g., very recent bottlenecks or founder effects). The method also allows to clearly identify individuals resulting from extreme consanguineous matings.

Genomic evaluation by including dominance effects and inbreeding depression for purebred and crossbred performance with an application in pigs

BACKGROUND

Improved performance of crossbred animals is partly due to heterosis. One of the major genetic bases of heterosis is dominance, but it is seldom used in pedigree-based genetic evaluation of livestock. Recently, a trivariate genomic best linear unbiased prediction (GBLUP) model including dominance was developed, which can distinguish purebreds from crossbred animals explicitly. The objectives of this study were: (1) methodological, to show that inclusion of marker-based inbreeding accounts for directional dominance and inbreeding depression in purebred and crossbred animals, to revisit variance components of additive and dominance genetic effects using this model, and to develop marker-based estimators of genetic correlations between purebred and crossbred animals and of correlations of allele substitution effects between breeds; (2) to evaluate the impact of accounting for dominance effects and inbreeding depression on predictive ability for total number of piglets born (TNB) in a pig dataset composed of two purebred populations and their crossbreds. We also developed an equivalent model that makes the estimation of variance components tractable.

RESULTS

For TNB in Danish Landrace and Yorkshire populations and their reciprocal crosses, the estimated proportions of dominance genetic variance to additive genetic variance ranged from 5 to 11%. Genetic correlations between breeding values for purebred and crossbred performances for TNB ranged from 0.79 to 0.95 for Landrace and from 0.43 to 0.54 for Yorkshire across models. The estimated correlation of allele substitution effects between Landrace and Yorkshire was low for purebred performances, but high for crossbred performances. Predictive ability for crossbred animals was similar with or without dominance. The inbreeding depression effect increased predictive ability and the estimated inbreeding depression parameter was more negative for Landrace than for Yorkshire animals and was in between for crossbred animals.

CONCLUSIONS

Methodological developments led to closed-form estimators of inbreeding depression, variance components and correlations that can be easily interpreted in a quantitative genetics context. Our results confirm that genetic correlations of breeding values between purebred and crossbred performances within breed are positive and moderate. Inclusion of dominance in the GBLUP model does not improve predictive ability for crossbred animals, whereas inclusion of inbreeding depression does.

The use of coancestry based on shared segments for maintaining genetic diversity

We have evaluated the use of genomic coancestry coefficients based on shared segments for the maintenance of genetic diversity through optimal contributions methodology for populations of three different Austrian cattle breeds. This coancestry measure has been compared with the genomic coancestry coefficient calculated on a SNP-by-SNP basis and with pedigree-based coancestry. The regressions of the shared segments coancestry on the other two coefficients suggest that the former mainly reflect Identity By Descent but with the advantage over pedigree-based coancestry of providing the realized Identity By Descent rather than an expectation. The effective population size estimated from the rate of coancestry based on shared segments was very similar to those obtained with the other coefficients and of small magnitude (from 26.24 to 111.90). This result highlights the importance of implementing active management strategies to control the increase of inbreeding and the loss of genetic diversity in livestock breeds, even when the population size is reasonably large. One problem for the implementation of coancestry based on shared segments is the need of estimating the gametic phases of the SNPs which, given the techniques used to obtain the genotypes, are a priori unknown. This study shows, through computer simulations, that using estimates of gametic phases for computing coancestry based on shared segments does not lead to a significant loss in the diversity maintained. This has been shown to be true even when the size of the population is very small as it is usually the case in populations subjected to conservation programmes.

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.

Assessing effective population size, coancestry and inbreeding effects on litter size using the pedigree and SNP data in closed lines of the Iberian pig breed

The complete pedigree of two closed Iberian pig lines (Gamito and Torbiscal), with 798 and 4077 reproducers, has been used to measure the evolution of coancestry (f) and inbreeding (F) for autosomal and X-linked genes along 16 and 28 respective equivalent discrete generations. At the last generation, the mean values of each line were f = 0.41 and 0.22, F = 0.35 and 0.18, fX  = 0.46 and 0.22 and FX  = 0.47 and 0.19, respectively. Other calculated parameters were the effective number of founders (final values, 6.8 and 35.2) and non-founders (1.5 and 2.4), founder genome equivalents (1.2 and 2.3) and effective population size (16.0 and 57.7). Measures of Torbiscal effective size based on rates of coancestry (66.1), inbreeding (65.0) and linkage disequilibrium (71.0) were estimated from whole-genome SNP genotyping data. Values of new and old inbreeding and their respective rates by generation were computed to detect purging effects of natural selection. The analysis of 6854 Torbiscal litters showed significant negative impacts of new and fast inbreeding on litter size, as expected from the purging hypothesis: -0.20 born piglets per litter by a 10% of new inbreeding, and -0.03 and -0.02 piglets by 1% of total and new inbreeding rates, respectively. The analysis performed on 1274 litters of the Gamito line failed to show purging effects. The only significant results were reductions in -0.91 and -0.17 piglets by a 10% of old and X-linked genes inbreeding, respectively. These results may be useful for some practical issues in conservation programs of farm or captive wild animals.

Ancestral Relationships Using Metafounders: Finite Ancestral Populations and Across Population Relationships

Recent use of genomic (marker-based) relationships shows that relationships exist within and across base population (breeds or lines). However, current treatment of pedigree relationships is unable to consider relationships within or across base populations, although such relationships must exist due to finite size of the ancestral population and connections between populations. This complicates the conciliation of both approaches and, in particular, combining pedigree with genomic relationships. We present a coherent theoretical framework to consider base population in pedigree relationships. We suggest a conceptual framework that considers each ancestral population as a finite-sized pool of gametes. This generates across-individual relationships and contrasts with the classical view which each population is considered as an infinite, unrelated pool. Several ancestral populations may be connected and therefore related. Each ancestral population can be represented as a "metafounder," a pseudo-individual included as founder of the pedigree and similar to an "unknown parent group." Metafounders have self- and across relationships according to a set of parameters, which measure ancestral relationships, i.e., homozygozities within populations and relationships across populations. These parameters can be estimated from existing pedigree and marker genotypes using maximum likelihood or a method based on summary statistics, for arbitrarily complex pedigrees. Equivalences of genetic variance and variance components between the classical and this new parameterization are shown. Segregation variance on crosses of populations is modeled. Efficient algorithms for computation of relationship matrices, their inverses, and inbreeding coefficients are presented. Use of metafounders leads to compatibility of genomic and pedigree relationship matrices and to simple computing algorithms. Examples and code are given.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Estimation of ancestral inbreeding effects on stillbirth, calving ease and birthweight in German Holstein dairy cattle

In this study, the effect of different measurements of ancestral inbreeding on birthweight, calving ease and stillbirth were analysed. Three models were used to estimate the effect of ancestral inbreeding, and the estimated regression coefficient of phenotypic data on different measurements of ancestral inbreeding was used to quantify the effect of ancestral inbreeding. The first model included only one measurement of inbreeding, whereas the second model included the classical inbreeding coefficients and one alternative inbreeding coefficient. The third model included the classical inbreeding coefficients, the interaction between classical inbreeding and ancestral inbreeding, and the classical inbreeding coefficients of the dam. Phenotypic data for this study were collected from February 1998 to December 2008 on three large commercial milk farms. During this time, 36,477 calving events were recorded. All calves were weighed after birth, and 8.08% of the calves died within 48 h after calving. Calving ease was recorded on a scale between 1 and 4 (1 = easy birth, 4 = surgery), and 69.95, 20.91, 8.92 and 0.21% of the calvings were scored with 1, 2, 3 and 4, respectively. The average inbreeding coefficient of inbred animals was 0.03, and average ancestral inbreeding coefficients were 0.08 and 0.01, depending on how ancestral inbreeding was calculated. Approximately 26% of classically non-inbred animals showed ancestral inbreeding. Correlations between different inbreeding coefficients ranged between 0.46 and 0.99. No significant effect of ancestral inbreeding was found for calving ease, because the number of animals with reasonable high level of ancestral inbreeding was too low. Significant effects of ancestral inbreeding were estimated for birthweight and stillbirth. Unfavourable effects of ancestral inbreeding were observed for birthweight. However, favourable purging effects were estimated for stillbirth, indicating that purging could be partly beneficial for genetic improvement of stillbirth.

Effective population size and inbreeding depression on litter size in rabbits. A case study

The purpose of this study is to use demographic and litter size data on four Spanish maternal lines of rabbits (A, V, H and LP), as a case study, in order to: (i) estimate the effective population size of the lines, as a measure of the rate of increase of inbreeding, and (ii) study whether the inbreeding effect on litter size traits depends on the pattern of its accumulation over time. The lines are being selected for litter size at weaning and are kept closed at the same selection nucleus under the same selection and management programme. The study considered 47,794 l and a pedigree of 14,622 animals. Some practices in mating and selection management allow an increase of the inbreeding coefficient lower than 0.01 per generation in these lines of around 25 males and 125 females. Their effective population size (Ne) was around 57.3, showing that the effect of selection, increasing the inbreeding, was counterbalanced by the management practices, intended to reduce the rate of inbreeding increase. The inbreeding of each individual was broken down into three components: old, intermediate and new inbreeding. The coefficients of regression of the old, intermediate and new inbreeding on total born (TB), number born alive (NBA) and number weaned (NW) per litter showed a decreasing trend from positive to negative values. Regression coefficients significantly different from zero were those for the old inbreeding on TB (6.79 ± 2.37) and NBA (5.92 ± 2.37). The contrast between the coefficients of regression between the old and new inbreeding were significant for the three litter size traits: 7.57 ± 1.72 for TB; 6.66 ± 1.73 for NBA and 5.13 ± 1.67 for NW. These results have been interpreted as the combined action of purging unfavourable genes and artificial selection favoured by the inbreeding throughout the generations of selection.

Inbreeding depression in livestock species: review and meta-analysis

Inbreeding, by virtue of its consequence on traits of interest, is a topic of major interest for geneticists and animal breeders. Based on meta-analysis conducted on 57 studies and seven livestock species considering a wide variety of selected traits, it was estimated that inbreeding depression corresponds to on average a decrease of 0.137 percent of the mean of a trait per 1 percent of inbreeding. The decrease was larger for production traits (reduction of 0.351%) than for other trait categories. For populations raised as purebreds, inbreeding depression may impact the economic income of breeders. There is a need for studies assessing the existence of an inbreeding purge phenomenon as well as the impact of inbreeding on adaptation capacities of livestock species. Promises brought by the development of dense genotyping as well as functional genomics will increase the capacities to improve our understanding and management of the phenomenon.

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.

Assessing the genetic diversity in small farm animal populations

Genetic variation is vital for the populations to adapt to varying environments and to respond to artificial selection; therefore, any conservation and development scheme should start from assessing the state of variation in the population. There are several marker-based and pedigree-based parameters to describe genetic variation. The most suitable ones are rate of inbreeding and effective population size, because they are not dependent on the amount of pedigree records. The acceptable level for effective population size can be considered from different angles leading to a conclusion that it should be at least 50 to 100. The estimates for the effective population size can be computed from the genealogical records or from demographic and marker information when pedigree data are not available. Marker information could also be used for paternity analysis and for estimation of coancestries. The sufficient accuracy in marker-based parameters would require typing thousands of markers. Across breeds, diversity is an important source of variation to rescue problematic populations and to introgress new variants. Consideration of adaptive variation brings new aspects to the estimation of the variation between populations.

Microsatellite support for active inbreeding in a cichlid fish

In wild animal populations, the degree of inbreeding differs between species and within species between populations. Because mating with kin often results in inbreeding depression, observed inbreeding is usually regarded to be caused by limited outbreeding opportunities due to demographic factors like small population size or population substructuring. However, theory predicts inclusive benefits from mating with kin, and thus part of the observed variation in inbreeding might be due to active inbreeding preferences. Although some recent studies indeed report kin mating preferences, the evidence is still highly ambiguous. Here, we investigate inbreeding in a natural population of the West African cichlid fish Pelvicachromis taeniatus which showed clear kin mating preferences in standardized laboratory experiments but no inbreeding depression. The presented microsatellite analysis reveals that the natural population has, in comparison to two reference populations, a reduced allelic diversity (A = 3) resulting in a low heterozygosity (H_o = 0.167) pointing to a highly inbred population. Furthermore, we found a significant heterozygote deficit not only at population (F_is = 0.116) but also at subpopulation level (F_is = 0.081) suggesting that inbreeding is not only a by-product of population substructuring but possibly a consequence of behavioral kin preferences.

Within-generation mutation variance for litter size in inbred mice

The mutational input of genetic variance per generation (sigma(m)(2)) is the lower limit of the genetic variability in inbred strains of mice, although greater values could be expected due to the accumulation of new mutations in successive generations. A mixed-model analysis using Bayesian methods was applied to estimate sigma(m)(2) and the across-generation accumulated genetic variability on litter size in 46 generations of a C57BL/6J inbred strain. This allowed for a separate inference on sigma(m)(2) and on the additive genetic variance in the base population (sigma(a)(2)). The additive genetic variance in the base generation was 0.151 and quickly decreased to almost null estimates in generation 10. On the other hand, sigma(m)(2) was moderate (0.035) and the within-generation mutational variance increased up to generation 14, then oscillating between 0.102 and 0.234 in remaining generations. This pattern suggested the existence of a continuous uploading of genetic variability for litter size (h(2)=0.045). Relevant genetic drift was not detected in this population. In conclusion, our approach allowed for separate estimation of sigma(a)(2) and sigma(m)(2) within the mixed-model framework, and the heritability obtained highlighted the significant and continuous influence of new genetic variability affecting the genetic stability of inbred strains.