Challenges in inbreeding estimation of large populations based on Polish Holstein-Friesian cattle pedigree

The aim of this study was to evaluate observed and future inbreeding level in Polish Holstein-Friesian cattle population. In total, over 9.8 mln animals were used in the analysis coming from the pedigree of Polish Federation of Cattle Breeders and Dairy Farmers. Inbreeding level, as an average per birth year, was estimated with the method accounting for missing parent information with the assumption of year 1950 as the base year of the population. If an animal had no ancestral records, an average inbreeding level from its birth year was assigned. Twice the average inbreeding level served as relatedness of the animal to the population, which enabled estimation of inbreeding in its offspring. The future inbreeding of potential offspring was estimated as an average of animals (bulls and cows) available for mating in a certain year. It was observed that 30–50% of animals born between 1985 and 2015 had no relevant ancestral information, which is caused by a high number of new animals and/or entire farms entering the national milk recordings. For the year 2015, the observed inbreeding level was 3.30%, which was more than twice the inbreeding with the classical approach (without missing parent information) and higher by 0.4% than the future inbreeding. The average increase of inbreeding in years 2010–2015 was 0.10%, which is similar to other countries monitored by World Holstein-Friesian Federation. However, the values might be underestimated due to low pedigree completeness. The estimates of future inbreeding suggested that observed inbreeding could be even lower and also increase slower, which indicates a constant need to monitor rate of increase in inbreeding over time. The most important aspect of presented results is the necessity to advise individual farmers to keep precise recordings of the matings on their farm in order to improve the pedigree completeness of Polish Holstein-Friesian and to use suitable mating programs to avoid too rapid growth of inbreeding.     

Exact inbreeding coefficients in populations with overlapping generations

A theory is given that allows inbreeding coefficients to be calculated exactly for populations with overlapping generations. Emphasis is placed on providing equations well suited for computer iteration. Both monoecious and dioecious populations are considered and family size is not restricted to being Poisson. One-locus and two-locus inbreeding coefficients are evaluated, although the reader may omit the two-locus sections. The exact treatment is shown to be preferable to approximate treatments in that it applies to both early and late generations for all populations sizes. Inbreeding effective numbers found by the exact treatment are compared to various approximate numbers, and the approximate values are found to be generally very good.     

A method for estimating Dekkera/Brettanomyces populations in wines

Aims:  The formation of ethylphenols in wines, a consequence of Dekkera/Brettanomyces metabolism, can affect their quality. The main aims of this work were to further our knowledge of Dekkera/Brettanomyces with respect to ethylphenol production, and to develop a methodology for detecting this spoilage yeast and for estimating its population size in wines using differential-selective media and high performance liquid chromatography (HPLC).
Methods and Results:  This work examines the reduction of p -coumaric acid and the formation of 4-vinylphenol and 4-ethylphenol (recorded by HPLC-DAD) in a prepared medium because of the activities of different yeast species and populations. A regression model was constructed for estimating the population of Dekkera/Brettanomyces at the beginning of fermentation via the conversion of hydroxycinnamic acids into ethylphenols.
Conclusions:  The proposed methodology allows the populations of Dekkera/Brettanomyces at the beginning of fermentation to be estimated in problem wines. Moreover, it avoids false positives because of yeasts resistant to the effects of the selective elements of the medium.
Significance and Impact of the Study:  This may help prevent the appearance of organoleptic anomalies in wines at the winery level.     

A mating method accounting for inbreeding and multi-trait selection in dairy cattle populations

Selection in dairy cattle populations usually takes into account both the breed profiles for many traits and their overall estimated breeding values (EBV). This can result in effective contributions of breeding animals departing substantially from contributions optimised for saving future genetic variability. In this work, we propose a mating method that considers not only inbreeding but also the detailed EBV of progeny or the EBV of sires in reference to acceptance thresholds. Penalties were defined for inbreeding and for inadequate EBV profiles. Relative reductions of penalties yielded by any mating design were expressed on a scale ranging from 0 to 1. A value of 0 represented the average performance of random matings and a value of 1 represented the maximal reduction allowed by a specialized, single-penalty, mating design. The core of the method was an adaptative simulated annealing, where the maximized function was the average of both ratios, under the constraints that both relative penalty reductions should be equal and that the within-herd concentration criterion should be equal to a predefined reasonable value. The method was tested on two French dairy cattle populations originating from the same AI organization. The optimised mating design allowed substantial reductions of penalty: 70% and 64% for the Holstein and the Normandy populations, respectively. Thus, this mating method decreased inbreeding and met various demands from breeders.     

An integrated-likelihood method for estimating genetic differentiation between populations

The aim of this article is to develop an integrated-likelihood (IL) approach to estimate the genetic differentiation between populations. The conventional maximum-likelihood (ML) and pseudolikelihood (PL) methods that use sample counts of alleles may cause severe underestimations of FST, which means overestimations of theta=4Nm, when the number of sampling localities is small. To reduce such bias in the estimation of genetic differentiation, we propose an IL method in which the mean allele frequencies over populations are regarded as nuisance parameters and are eliminated by integration. To maximize the IL function, we have developed two algorithms, a Monte Carlo EM algorithm and a Laplace approximation. Our simulation studies show that the method proposed here outperforms the conventional ML and PL methods in terms of unbiasedness and precision. The IL method was applied to real data for Pacific herring and African elephants.     

Dynamics of inbreeding depression due to deleterious mutations in small populations: mutation parameters and inbreeding rate

A multilocus stochastic model is developed to simulate the dynamics of mutational load in small populations of various sizes. Old mutations sampled from a large ancestral population at mutation-selection balance and new mutations arising each generation are considered jointly, using biologically plausible lethal and deleterious mutation parameters. The results show that inbreeding depression and the number of lethal equivalents due to partially recessive mutations can be partly purged from the population by inbreeding, and that this purging mainly involves lethals or detrimentals of large effect. However, fitness decreases continuously with inbreeding, due to increased fixation and homozygosity of mildly deleterious mutants, resulting in extinctions of very small populations with low reproductive rates. No optimum inbreeding rate or population size exists for purging with respect to fitness (viability) changes, but there is an optimum inbreeding rate at a given final level of inbreeding for reducing inbreeding depression or the number of lethal equivalents. The interaction between selection against partially recessive mutations and genetic drift in small populations also influences the rate of decay of neutral variation. Weak selection against mutants relative to genetic drift results in apparent overdominance and thus an increase in effective size (Ne) at neutral loci, and strong selection relative to drift leads to a decrease in Ne due to the increased variance in family size. The simulation results and their implications are discussed in the context of biological conservation and tests for purging.     

The use of frozen semen to minimize inbreeding in small populations

In this study, we compared the average coancestry and inbreeding levels for two genetic conservation schemes in which frozen semen from a gene bank is used to reduce the inbreeding in a live population. For a simple scheme in which only semen of generation-0 (G0) sires is used, the level of inbreeding asymptotes to 1/(2N), where N is the number of newborn sires in the base generation and rate of inbreeding goes to zero. However, when only sires of G0 are selected, all genes will eventually descend from the founder sires and all genes from the founder dams are lost. We propose an alternative scheme in which N sires from generation 1 (G1), as well as the N sires from G0, have semen conserved, and the semen of G0 and G1 sires is used for dams of odd and even generation numbers, respectively. With this scheme, the level of inbreeding asymptotes to 1/(3N) and the genes of founder dams are also conserved, because 50% of the genes of sires of G1 are derived from the founder dams. A computer simulation study shows that this is the optimum design to minimize inbreeding, even if semen from later generations is available.     

Genetic diversity in Leavenworthia populations with different inbreeding levels.

Levels of neutral genetic diversity within and between populations were compared between outcrossing (self-incompatible) and inbreeding populations in the annual plant genus Leavenworthia. Two taxonomically independent comparisons are possible, since self-incompatibility has been lost twice in the group of species studied. Within inbred populations of L.uniflora and L.crassa, no DNA sequence variants were seen among the alleles sampled, but high diversity was seen in alleles from populations of the outcrosser L. stylosa, and in self-incompatible L. crassa populations. Diversity between populations was seen in all species. Although total diversity values were lower in the sets of inbreeding populations, between-population values were as high or higher, than those in the outcrossing taxa. Possible reasons for these diversity patterns are discussed. As the effect of inbreeding appears to be a greater than twofold reduction in diversity, we argue that some process such as selection for advantageous mutations, or against deleterious mutations, or bottlenecks occurring predominantly in the inbreeders, appears necessary to account for the findings. If selection for advantageous mutations is responsible, it appears that it must be some form of local adaptive selection, rather than substitution of alleles that are advantageous throughout the species. This is consistent with the finding of high between-population diversity in the inbreeding taxa.     

Optimum contribution selection in large general tree breeding populations with an application to Scots pine

Development of selection methods that optimises selection differential subject to a constraint on the increase of inbreeding (or coancestry) in a population is an important part of breeding programmes. One such method that has received much attention in animal breeding is the optimum contribution (OC) dynamic selection method. We implemented the OC algorithm and applied it to a diallel progeny trial of Pinus sylvestris L. (Scots pine) focussing on two traits (total tree height and stem diameter). The OC method resulted in a higher increase in genetic gain (8–30%) compared to the genetic gain achieved using standard restricted selection method at the same level of coancestry constraint. Genetic merit obtained at two different levels of restriction on coancestry showed that the benefit of OC was highest when restriction was strict. At the same level of genetic merit, OC decreased coancestry with 56 and 39% for diameter and height, respectively, compared to the level of coancestry obtained using unrestricted truncation selection. Inclusion of a dominance term in the statistical model resulted in changes in contribution rank of trees with 7 and 13% for diameter and height, respectively, compared to results achieved by using a pure additive model. However, the genetic gain was higher for the pure additive model than for the model including dominance for both traits.     

A method of estimating total food consumed by fish populations

The partitioning of a volumetrically investigated fish diet into food types was employed to convert cumulated food energy (Winberg) into the wet weight of food types.     

Prediction of rates of inbreeding in selected populations

A method is presented for the prediction of rate of inbreeding for populations with discrete generations. The matrix of Wright's numerator relationships is partitioned into 'contribution' matrices which describe the contribution of the Mendelian sampling of genes of ancestors in a given generation to the relationship between individuals in later generations. These contributions stabilize with time and the value to which they stabilize is shown to be related to the asymptotic rate of inbreeding and therefore also the effective population size, Ne approximately 2N/(mu 2r + sigma 2r), where N is the number of individuals per generation and mu r and sigma 2r are the mean and variance of long-term relationships or long-term contributions. These stabilized values are then predicted using a recursive equation via the concept of selective advantage for populations with hierarchical mating structures undergoing mass selection. Account is taken of the change in genetic parameters as a consequence of selection and also the increasing 'competitiveness' of contemporaries as selection proceeds. Examples are given and predicted rates of inbreeding are compared to those calculated in simulations. For populations of 20 males and 20, 40, 100 or 200 females the rate of inbreeding was found to increase by as much as 75% over the rate of inbreeding in an unselected population depending on mating ratio, selection intensity and heritability of the selected trait. The prediction presented here estimated the rate of inbreeding usually within 5% of that calculated from simulation.     

Kin selection and coefficients of relatedness in family-structured populations with inbreeding

We consider family specific fitnesses that depend on mixed strategies of two basic phenotypes or behaviours. Pairwise interactions are assumed, but they are restricted to occur between sibs. To study the change in frequency of a rare mutant allele, we consider two different forms of weak selection, one applied through small differences in genotypic values determining individual mixed strategies, the other through small differences in viabilities according to the behaviours chosen by interacting sibs. Under these two specific forms of weak selection, we deduce conditions for initial increase in frequency of a rare mutant allele for autosomal genes in the partial selfing model as well as autosomal and sex-linked genes in the partial sib-mating model with selection before mating or selection after mating. With small differences in mixed strategies, we show that conditions for protection of a mutant allele are tantamount to conditions for initial increase in frequency obtained in additive kin selection models. With particular reference to altruism versus selfishness, we provide explicit ranges of values for the selfing or sib-mating rate based on a fixed cost-benefit ratio and the dominance scheme that allow the spreading of a rare mutant allele into the population. This study confirms that more inbreeding does not necessarily promote the evolution of altruism. Under the hypothesis of small differences in viabilities, the situation is much more intricate unless an additive model is assumed. In general however, conditions for initial increase in frequency of a mutant allele can be obtained in terms of fitness effects that depend on the genotypes of interacting individuals or their mates and generalized conditional coefficients of relatedness according to the inbreeding condition of the interacting individuals.     

Predicting rates of inbreeding in populations undergoing selection

Tractable forms of predicting rates of inbreeding (DeltaF) in selected populations with general indices, nonrandom mating, and overlapping generations were developed, with the principal results assuming a period of equilibrium in the selection process. An existing theorem concerning the relationship between squared long-term genetic contributions and rates of inbreeding was extended to nonrandom mating and to overlapping generations. DeltaF was shown to be approximately (1)/(4)(1 - omega) times the expected sum of squared lifetime contributions, where omega is the deviation from Hardy-Weinberg proportions. This relationship cannot be used for prediction since it is based upon observed quantities. Therefore, the relationship was further developed to express DeltaF in terms of expected long-term contributions that are conditional on a set of selective advantages that relate the selection processes in two consecutive generations and are predictable quantities. With random mating, if selected family sizes are assumed to be independent Poisson variables then the expected long-term contribution could be substituted for the observed, providing (1)/(4) (since omega = 0) was increased to (1)/(2). Established theory was used to provide a correction term to account for deviations from the Poisson assumptions. The equations were successfully applied, using simple linear models, to the problem of predicting DeltaF with sib indices in discrete generations since previously published solutions had proved complex.     

Quantifying inbreeding in natural populations of hermaphroditic organisms

We review molecular methods for estimating selfing rates and inbreeding in populations. Two main approaches are available: the population structure approach (PSA) and progeny-array approach (PAA). The PSA approach relies on single-generation samples and produces estimates that integrate the inbreeding history over several generations, but is based on strong assumptions (for example, inbreeding equilibrium). The PSA has classically relied on single-locus inbreeding coefficients averaged over loci. Unfortunately PSA estimates are very sensitive to technical problems such as the occurrence of null alleles at one or more of the loci. Consequently inbreeding might be substantially overestimated, especially in outbred populations. However, the robustness of the PSA has recently been greatly improved by the development of multilocus methods free of such bias. The PAA, on the other hand, is based on the comparison between offspring and mother genotypes. As a consequence, PAA estimates do not reflect long-term inbreeding history but only recent mating events of the maternal individuals studied ('here and now' selfing). In addition to selfing rates, the PAA allows estimating other mating system parameters, including biparental inbreeding and the correlation of selfing among sibs. Although PAA estimates could also be biased by technical problems, incompatibilities between the mother's genotype and her offspring allow the identification and correction of such bias. For all methods, we provide guidelines on the required number of loci and sample sizes. We conclude that the PSA and PAA are equally robust, provided multilocus information is used. Although experimental constraints may make the PAA more demanding, especially in animals, the two methods provide complementary information, and can fruitfully be conducted together.     

Minimising inbreeding in small populations by rotational mating with frozen semen

Mating plans are investigated in order to minimize inbreeding in small populations when frozen semen is available. For a single dam line it was found that specific sire rotations minimized the asymptotic level of inbreeding when semen is used repeatedly from certain generations. When semen of N foundation (G0) sires is used rotationally over generations it is shown that the inbreeding level asymptotes to 1/(2(N+1) - 2). However, if only G0 sires are used then all genes will eventually descend from the founder sires. Inbreeding can be reduced further by using sires from generation one (G1) and later as this retains genes from the founder dams in the long-term gene pool. If semen from NG0 sires and N unrelated G1 sons is used rotationally then inbreeding asymptotes to (2(N-1) + 1)/(2(2N+1) - 2). When there are more founder dams than sires, the asymptotic inbreeding can be reduced even further by using the semen of half-sib G1 sires in rotation. Optimal rotations using full-sib G1 sires or generation 2 (or later) sires will lower the asymptotic inbreeding also, but generally not by much. It was found that when unlimited frozen semen from a specified group of sires was available, the optimal mating plan was achieved by selecting each generation the sire with the least co-ancestry with the current female of the dam line.     

A simple method to approximate gene content in large pedigree populations: application to the myostatin gene in dual-purpose Belgian Blue cattle

Gene content is the number of copies of a particular allele in a genotype of an animal. Gene content can be used to study additive gene action of candidate gene. Usually genotype data are available only for a part of population and for the rest gene contents have to be calculated based on typed relatives. Methods to calculate expected gene content for animals on large complex pedigrees are relatively complex. In this paper we proposed a practical method to calculate gene content using a linear regression. The method does not estimate genotype probabilities but these can be approximated from gene content assuming Hardy-Weinberg proportions. The approach was compared with other methods on multiple simulated data sets for real bovine pedigrees of 1 082 and 907 903 animals. Different allelic frequencies (0.4 and 0.2) and proportions of the missing genotypes (90, 70, and 50%) were considered in simulation. The simulation showed that the proposed method has similar capability to predict gene content as the iterative peeling method, however it requires less time and can be more practical for large pedigrees. The method was also applied to real data on the bovine myostatin locus on a large dual-purpose Belgian Blue pedigree of 235 133 animals. It was demonstrated that the proposed method can be easily adapted for particular pedigrees.