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

Genomic analysis of inbreeding and coancestry in Nordic Jersey and Holstein dairy cattle populations

In recent years, genomic selection (GS) has accelerated genetic gain in dairy cattle breeds worldwide. Despite the evident genetic progress, several dairy populations have also encountered challenges such as heightened inbreeding rates and reduced effective population sizes. The challenge has been to find a balance between achieving substantial genetic gain while managing genetic diversity within the population, thereby mitigating the negative effects of inbreeding depression. This study aims to elucidate the impact of GS on pedigree and genomic rates of inbreeding (ΔF) and coancestry (ΔC) in Nordic Jersey (NJ) and Holstein (NH) cattle populations. Furthermore, key genetic metrics, including the generation interval (L), effective population size (Ne), and future effective population size (FNe) were assessed between 2 time periods, before and after GS, and across distinct animal cohorts in both breeds: females, bulls, and approved semen-producing bulls (AI-sires). Analysis of ΔF and ΔC revealed distinct trends across the studied periods and animal groups. Notably, there was a consistent increase in yearly ΔF for most animal groups in both breeds. An exception was observed in NH AI-sires, which demonstrated a slight decrease in yearly ΔF. Moreover, NJ displayed minimal changes in yearly ΔC between the periods, whereas NH exhibited elevated ΔC values across all animal groups. Particularly striking was the substantial increase in yearly ΔC within the NH female population, surging from 0.02% to 0.39% between the periods. Implementation of GS resulted in a reduction of the generation interval across all animal cohorts in both NJ and NH breeds. However, the extent of reduction was more pronounced in males compared with females. This reduction in generation interval influenced generational changes in ΔF and ΔC. Bulls and AI-sires of both breeds exhibited reduced generational ΔF between periods, in contrast to females that demonstrated an opposing pattern. Between the periods, NJ maintained a relatively stable Ne (29.4 before and 30.3 after GS), whereas NH experienced a notable decline from 54.3 to 42.8. Female groups in both breeds displayed a negative Ne trend, whereas males demonstrated either neutral or positive Ne developments. Regarding FNe, NJ exhibited positive FNe development with an increase from 40.7 to 57.2. The opposite was observed in NH, where FNe decreased from 198.8 to 42.7. In summary, it was evident that the genomic methods could detect differences between the populations and changes in ΔF and ΔC more efficiently than pedigree methods. Implementation of GS yielded positive outcomes within the NJ population regarding the rate of coancestry but the opposite was observed with NH. Moreover, analysis of ΔC data hints at the potential to decrease future ΔF through informed mating strategies. Conversely, NH faces more pressing concerns, even though ΔF remains comparatively modest in contrast to what has been observed in other Holstein populations. These findings underscore the necessity of genomic control of inbreeding and coancestry with strategic changes in the Nordic breeding schemes for dairy to ensure long-term sustainability in the forthcoming years.

Integrated analysis strategy of genome-wide functional gene mining reveals DKK2 gene underlying meat quality in Shaziling synthesized pigs

BACKGROUND

Shaziling pig is a well-known indigenous breed in China who has superior meat quality traits. However, the genetic mechanism and genomic evidence underlying meat quality characteristics of Shaziling pigs are still unclear. To explore and investigate the germplasm characteristics of Shaziling pigs, we totally analyzed 67 individual's whole genome sequencing data for the first time (20 Shaziling pigs [S], 20 Dabasha pigs [DBS], 11 Yorkshire pigs [Y], 10 Berkshire pigs [BKX], 5 Basha pigs [BS] and 1 Warthog).

RESULTS

A total of 2,538,577 SNPs with high quality were detected and 9 candidate genes which was specifically selected in S and shared in S to DBS were precisely mined and screened using an integrated analysis strategy of identity-by-descent (IBD) and selective sweep. Of them, dickkopf WNT signaling pathway inhibitor 2 (DKK2), the antagonist of Wnt signaling pathway, was the most promising candidate gene which was not only identified an association of palmitic acid and palmitoleic acid quantitative trait locus in PigQTLdb, but also specifically selected in S compared to other 48 Chinese local pigs of 12 populations and 39 foreign pigs of 4 populations. Subsequently, a mutation at 12,726-bp of DKK2 intron 1 (g.114874954 A > C) was identified associated with intramuscular fat content using method of PCR-RFLP in 21 different pig populations. We observed DKK2 specifically expressed in adipose tissues. Overexpression of DKK2 decreased the content of triglyceride, fatty acid synthase and expression of relevant genes of adipogenic and Wnt signaling pathway, while interference of DKK2 got contrary effect during adipogenesis differentiation of porcine preadipocytes and 3T3-L1 cells.

CONCLUSIONS

Our findings provide an analysis strategy for mining functional genes of important economic traits and provide fundamental data and molecular evidence for improving pig meat quality traits and molecular breeding.

Impact of kinship matrices on genetic gain and inbreeding with optimum contribution selection in a genomic dairy cattle breeding program

BACKGROUND

Genomic selection has increased genetic gain in dairy cattle, but in some cases it has resulted in higher inbreeding rates. Therefore, there is need for research on efficient management of inbreeding in genomically-selected dairy cattle populations, especially for local breeds with a small population size. Optimum contribution selection (OCS) minimizes the increase in average kinship while it maximizes genetic gain. However, there is no consensus on how to construct the kinship matrix used for OCS and whether it should be based on pedigree or genomic information. VanRaden's method 1 (VR1) is a genomic relationship matrix in which centered genotype scores are scaled with the sum of 2p(1-p) where p is the reference allele frequency at each locus, and VanRaden's method 2 (VR2) scales each locus with 2p(1-p), thereby giving greater weight to loci with a low minor allele frequency. We compared the effects of nine kinship matrices on genetic gain, kinship, inbreeding, genetic diversity, and minor allele frequency when applying OCS in a simulated small dairy cattle population. We used VR1 and VR2, each using base animals, all genotyped animals, and the current generation of animals to compute reference allele frequencies. We also set the reference allele frequencies to 0.5 for VR1 and the pedigree-based relationship matrix. We constrained OCS to select a fixed number of sires per generation for all scenarios. Efficiency of the different matrices were compared by calculating the rate of genetic gain for a given rate of increase in average kinship.

RESULTS

We found that: (i) genomic relationships were more efficient than pedigree-based relationships at managing inbreeding, (ii) reference allele frequencies computed from base animals were more efficient compared to reference allele frequencies computed from recent animals, and (iii) VR1 was slightly more efficient than VR2, but the difference was not statistically significant.

CONCLUSIONS

Using genomic relationships for OCS realizes more genetic gain for a given amount of kinship and inbreeding than using pedigree relationships when the number of sires is fixed. For a small genomic dairy cattle breeding program, we recommend that the implementation of OCS uses VR1 with reference allele frequencies estimated either from base animals or old genotyped animals.

Assessment of long-term trends in genetic mean and variance after the introduction of genomic selection in layers: a simulation study

Nucleus-based breeding programs are characterized by intense selection that results in high genetic gain, which inevitably means reduction of genetic variation in the breeding population. Therefore, genetic variation in such breeding systems is typically managed systematically, for example, by avoiding mating the closest relatives to limit progeny inbreeding. However, intense selection requires maximum effort to make such breeding programs sustainable in the long-term. The objective of this study was to use simulation to evaluate the long-term impact of genomic selection on genetic mean and variance in an intense layer chicken breeding program. We developed a large-scale stochastic simulation of an intense layer chicken breeding program to compare conventional truncation selection to genomic truncation selection optimized with either minimization of progeny inbreeding or full-scale optimal contribution selection. We compared the programs in terms of genetic mean, genic variance, conversion efficiency, rate of inbreeding, effective population size, and accuracy of selection. Our results confirmed that genomic truncation selection has immediate benefits compared to conventional truncation selection in all specified metrics. A simple minimization of progeny inbreeding after genomic truncation selection did not provide any significant improvements. Optimal contribution selection was successful in having better conversion efficiency and effective population size compared to genomic truncation selection, but it must be fine-tuned for balance between loss of genetic variance and genetic gain. In our simulation, we measured this balance using trigonometric penalty degrees between truncation selection and a balanced solution and concluded that the best results were between 45° and 65°. This balance is specific to the breeding program and depends on how much immediate genetic gain a breeding program may risk vs. save for the future. Furthermore, our results show that the persistence of accuracy is better with optimal contribution selection compared to truncation selection. In general, our results show that optimal contribution selection can ensure long-term success in intensive breeding programs using genomic selection.

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.

Effect of clonal testing on the efficiency of genomic evaluation in forest tree breeding

Through stochastic simulations, accuracies of breeding values and response to selection were assessed under traditional pedigree-(BLUP) and genomic-based evaluation methods (GBLUP) in forest tree breeding. The latter provides a methodological foundation for genomic selection. We evaluated the impact of clonal replication in progeny testing on the response to selection realized in seed orchards under variable marker density and target effective population sizes. We found that clonal replication in progeny trials boosted selection accuracy, thus providing additional genetic gains under BLUP. While a similar trend was observed for GBLUP, however, the added gains did not surpass those under BLUP. Therefore, breeding programs deploying extensive progeny testing with clonal propagation might not benefit from the deployment of genomic information. These findings could be helpful in the context of operational breeding programs.

Long-term impact of conventional and optimal contribution conservation methods on genetic diversity and genetic gain in local pig breeds

There are rich and vast genetic resources of indigenous pig breeds in the world. Currently, great attention is paid to either crossbreeding or conservation of these indigenous pig breeds, and insufficient attention is paid to the combination of conservation and breeding along with their long-term effects on genetic diversity. Therefore, the objective of this study is to compare the long-term effects of using conventional conservation and optimal contribution selection methods on genetic diversity and genetic gain. A total of 11 different methods including conventional conservation and optimal contribution selection methods were investigated using stochastic simulations. The long-term effects of using these methods were evaluated in terms of genetic diversity metrices such as expected heterozygosity (He) and the rate of genetic gain. The results indicated that the rates of true inbreeding in these conventional conservation methods were maintained at around 0.01. The optimal contribution selection methods based either on the pedigree (POCS) or genome (GOCS) information showed more genetic gain than conventional methods, and POCS achieved the largest genetic gain. Furthermore, the effect of using GOCS methods on most of the genetic diversity metrics was slightly better than the conventional conservation methods when the rate of true inbreeding was the same, but this also required more sires used in OCS methods. According to the rate of true inbreeding, there was no significant difference among these conventional methods. In conclusion, there is no significant difference in different ways of selecting sows on inbreeding when we use different conventional conservation methods. Compared with conventional methods, POCS method could achieve the most genetic gain. However, GOCS methods can not only achieve higher genetic gain, but also maintain a relatively high level of genetic diversity. Therefore, GOCS is a better choice if we want to combine conservation and breeding in actual production in the conservation farms.

Mating allocations in Nordic Red Dairy Cattle using genomic information

In this study, we compared mating allocations in Nordic Red Dairy Cattle using genomic information. We used linear programming to optimize different economic scores within each herd, considering genetic level, semen cost, the economic impact of recessive genetic defects, and genetic relationships. We selected 9,841 genotyped females born in Denmark, Finland, or Sweden in 2019 for mating allocations. We used 2 different pedigree relationship coefficients, the first tracing the pedigree 3 generations back from the parents of the potential mating and the second based on all available pedigree information. We used 3 different genomic relationship coefficients, 1 SNP-by-SNP genomic relationship and 2 based on shared genomic segments. We found high correlations (≥0.83) between the pedigree and genomic relationship measures. The mating results showed that it was possible to reduce the different genetic relationships between parents with minimal effect on genetic level. Including the cost of known recessive genetic defects eliminated expression of genetic defects. It was possible to reduce genomic relationships between parents with pedigree measures, but it was best done with genomic measures. Linear programming maximized the economic score for all herds studied within seconds, which means that it is suitable for implementation in mating software to be used by advisors and farmers.

Pre-selection against a lethal recessive allele in breeding schemes with optimum-contribution selection or truncation selection

BACKGROUND

We tested the hypothesis that breeding schemes with a pre-selection step, in which carriers of a lethal recessive allele (LRA) were culled, and with optimum-contribution selection (OCS) reduce the frequency of a LRA, control rate of inbreeding, and realise as much genetic gain as breeding schemes without a pre-selection step.

METHODS

We used stochastic simulation to estimate true genetic gain realised at a 0.01 rate of true inbreeding (ΔF_true) by breeding schemes that combined one of four pre-selection strategies with one of three selection strategies. The four pre-selection strategies were: (1) no carriers culled, (2) male carriers culled, (3) female carriers culled, and (4) all carriers culled. Carrier-status was known prior to selection. The three selection strategies were: (1) OCS in which [Formula: see text] was predicted and controlled using pedigree relationships (POCS), (2) OCS in which [Formula: see text] was predicted and controlled using genomic relationships (GOCS), and (3) truncation selection of parents. All combinations of pre-selection strategies and selection strategies were tested for three starting frequencies of the LRA (0.05, 0.10, and 0.15) and two linkage statuses with the locus that has the LRA being on a chromosome with or without loci affecting the breeding goal trait. The breeding schemes were simulated for 10 discrete generations (t = 1, …, 10). In all breeding schemes, ΔF_true was calibrated to be 0.01 per generation in generations t = 4, …, 10. Each breeding scheme was replicated 100 times.

RESULTS

We found no significant difference in true genetic gain from generations t = 4, …, 10 between breeding schemes with or without pre-selection within selection strategy. POCS and GOCS schemes realised similar true genetic gains from generations t = 4, …, 10. POCS and GOCS schemes realised 12% more true genetic gain from generations t = 4, …, 10 than truncation selection schemes.

CONCLUSIONS

We advocate for OCS schemes with pre-selection against the LRA that cause animal suffering and high costs. At LRA frequencies of 0.10 or lower, OCS schemes in which male carriers are culled reduce the frequency of LRA, control rate of inbreeding, and realise no significant reduction in true genetic gain compared to OCS schemes without pre-selection against LRA.

Why and How to Switch to Genomic Selection: Lessons From Plant and Animal Breeding Experience

The present study is a transversal analysis of the interest in genomic selection for plant and animal species. It focuses on the arguments that may convince breeders to switch to genomic selection. The arguments are classified into three different “bricks.” The first brick considers the addition of genotyping to improve the accuracy of the prediction of breeding values. The second consists of saving costs and/or shortening the breeding cycle by replacing all or a portion of the phenotyping effort with genotyping. The third concerns population management to improve the choice of parents to either optimize crossbreeding or maintain genetic diversity. We analyse the relevance of these different bricks for a wide range of animal and plant species and sought to explain the differences between species according to their biological specificities and the organization of breeding programs.

Competitive gene flow does not necessarily maximize the genetic gain of genomic breeding programs in the presence of genotype-by-environment interaction

National and international across-population selection is often recommended and fairly common in the current breeding practice of dairy cattle, with the primary aims to increase genetic gain and genetic variability. The aim of this study was to test the hypothesis that the strategy of truncation selection of sires across populations [i.e., competitive gene flow strategy (CGF)] may not necessarily maximize genetic gain in the long term in the presence of genotype-by-environment interaction (G×E). Two alternative strategies used to be compared with CGF were forced gene flow (FGF) strategies, with 10 or 50% of domestic dams forced to be mated with foreign sires (FGF10%, FGF50%). Two equal-size populations (N_dams = 1,000) that were selected for the same breeding goal trait (h² = 0.3) under G×E correlation (r_g) of either 0.9 or 0.8 were simulated to test these 3 different strategies. Each population first experienced either 5 or 20 differentiation generations (G_d), then 15 migration generations. Discrete generations were simulated for simplicity. Each population performed a within-population conventional breeding program during differentiation generations and the 3 across-population sire selection strategies based on joint genomic prediction during migration generations. The 4 G_{d_}r_g combinations defined 4 different levels of differentiation degree between the 2 populations at the start of migration. The true rate of inbreeding over the last 10 migration generations in each scenario was constrained at 0.01 to provide a fair basis for comparison of genetic gain across scenarios. Results showed that CGF maximized the genetic gain after 15 migration generations in 5_0.9 combination only, the case of the lowest differentiation degree, with a superiority of 0.4% (0.04 genetic SD units) over the suboptimal strategy. While in 5_0.8, 20_0.9, and 20_0.8 combinations, 2 FGF strategies had a superiority in genetic gain of 2.3 to 12.5% (0.21-1.07 genetic SD units) over CGF after 15 migration generations, especially FGF50%. The superiority of FGF strategies over CGF was that they alleviated inbreeding, introduced new genetic variance in the early migration period, and improved accuracy in the entire migration period. Therefore, we concluded that CGF does not necessarily maximize the genetic gain of across-population genomic breeding programs given moderate G×E. The across-population selection strategy remains to be optimized to maximize genetic gain.

Statistical model and testing designs to increase response to selection with constrained inbreeding in genomic breeding programs for pigs affected by social genetic effects

Background

Social genetic effects (SGE) are the effects of the genotype of one animal on the phenotypes of other animals within a social group. Because SGE contribute to variation in economically important traits for pigs, the inclusion of SGE in statistical models could increase responses to selection (RS) in breeding programs. In such models, increasing the relatedness of members within groups further increases RS when using pedigree-based relationships; however, this has not been demonstrated with genomic-based relationships or with a constraint on inbreeding. In this study, we compared the use of statistical models with and without SGE and compared groups composed at random versus groups composed of families in genomic selection breeding programs with a constraint on the rate of inbreeding.

Results

When SGE were of a moderate magnitude, inclusion of SGE in the statistical model substantially increased RS when SGE were considered for selection. However, when SGE were included in the model but not considered for selection, the increase in RS and in accuracy of predicted direct genetic effects (DGE) depended on the correlation between SGE and DGE. When SGE were of a low magnitude, inclusion of SGE in the model did not increase RS, probably because of the poor separation of effects and convergence issues of the algorithms. Compared to a random group composition design, groups composed of families led to higher RS. The difference in RS between the two group compositions was slightly reduced when using genomic-based compared to pedigree-based relationships.

Conclusions

The use of a statistical model that includes SGE can substantially improve response to selection at a fixed rate of inbreeding, because it allows the heritable variation from SGE to be accounted for and capitalized on. Compared to having random groups, family groups result in greater response to selection in the presence of SGE but the advantage of using family groups decreases when genomic-based relationships are used.

Management of Genetic Diversity in the Era of Genomics

Management of genetic diversity aims to (i) maintain heterozygosity, which ameliorates inbreeding depression and loss of genetic variation at loci that may become of importance in the future; and (ii) avoid genetic drift, which prevents deleterious recessives (e.g., rare disease alleles) from drifting to high frequency, and prevents random drift of (functional) traits. In the genomics era, genomics data allow for many alternative measures of inbreeding and genomic relationships. Genomic relationships/inbreeding can be classified into (i) homozygosity/heterozygosity based (e.g., molecular kinship matrix); (ii) genetic drift-based, i.e., changes of allele frequencies; or (iii) IBD-based, i.e., SNPs are used in linkage analyses to identify IBD segments. Here, alternative measures of inbreeding/relationship were used to manage genetic diversity in genomic optimal contribution (GOC) selection schemes. Contrary to classic inbreeding theory, it was found that drift and homozygosity-based inbreeding could differ substantially in GOC schemes unless diversity management was based upon IBD. When using a homozygosity-based measure of relationship, the inbreeding management resulted in allele frequency changes toward 0.5 giving a low rate of increase in homozygosity for the panel used for management, but not for unmanaged neutral loci, at the expense of a high genetic drift. When genomic relationship matrices were based on drift, following VanRaden and as in GCTA, drift was low at the expense of a high rate of increase in homozygosity. The use of IBD-based relationship matrices for inbreeding management limited both drift and the homozygosity-based rate of inbreeding to their target values. Genetic improvement per percent of inbreeding was highest when GOC used IBD-based relationships irrespective of the inbreeding measure used. Genomic relationships based on runs of homozygosity resulted in very high initial improvement per percent of inbreeding, but also in substantial discrepancies between drift and homozygosity-based rates of inbreeding, and resulted in a drift that exceeded its target value. The discrepancy between drift and homozygosity-based rates of inbreeding was caused by a covariance between initial allele frequency and the subsequent change in frequency, which becomes stronger when using data from whole genome sequence.

Controlling Coancestry and Thereby Future Inbreeding by Optimum-Contribution Selection Using Alternative Genomic-Relationship Matrices

We tested the consequences of using alternative genomic relationship matrices to predict genomic breeding values (GEBVs) and control of coancestry in optimum contribution selection, where the relationship matrix used to calculate GEBVs was not necessarily the same as that used to control coancestry. A stochastic simulation study was carried out to investigate genetic gain and true genomic inbreeding in breeding schemes that applied genomic optimum contribution selection (GOCS) with different genomic relationship matrices. Three genomic-relationship matrices were used to predict the GEBVs based on three information sources: markers (G_M), QTL (G_Q), and markers and QTL (G_A). Strictly, G_Q is not possible to implement in practice since we do not know the quantitative trait loci (QTL) positions, but more and more information is becoming available especially about the largest QTL. Two genomic-relationship matrices were used to control coancestry: G_M and G_A. Three genetic architectures were simulated: with 7702, 1000, and 500 QTLs together with 54,218 markers. Selection was for a single trait with heritability 0.2. All selection candidates were phenotyped and genotyped before selection. With 7702 QTL, there were no significant differences in rates of genetic gain at the same rate of true inbreeding using different genomic relationship matrices in GOCS. However, as the number of QTLs was reduced to 1000, prediction of GEBVs using a genomic relationship matrix constructed based on G_Q and control of coancestry using G_M realized 29.7% higher genetic gain than using G_M for both prediction and control of coancestry. Forty-three percent of this increased rate of genetic gain was due to increased accuracies of GEBVs. These findings indicate that with large numbers of QTL, it is not critical what information, i.e., markers or QTL, is used to construct genomic-relationship matrices. However, it becomes critical with small numbers of QTL. This highlights the importance of using genomic-relationship matrices that focus on QTL regions for GEBV estimation when the number of QTL is small in GOCS. Relationships used to control coancestry are preferably based on marker data.