Impacts of genotyping strategies on long‐term genetic response in genomic selection

The present study investigated the effects of the choices of animals of reference populations on long‐term responses to genomic selection. Simulated populations comprised 300 individuals and 10 generations of selection practiced for a trait with heritability of 0.1, 0.3 or 0.5. Thirty individuals were randomly selected in the first five generations and selected by estimated breeding values from best linear unbiased prediction (BLUP) and genomic BLUP in the subsequent five generations. The reference populations comprise all animals for all generations (scenario 1), all animals for 6‐10 generations (scenario 2) and 2‐6 generations (scenario 3), and half of the animals for all generations (scenario 4). For all heritability levels, the genetic gains in generation 10 were similar in scenarios 1 and 2. Among scenarios 2 to 4, the highest genetic gains were obtained in scenario 2, with heritabilities of 0.1 and 0.3 as well as scenario 4 with heritability of 0.5. The inbreeding coefficients in scenarios 1, 2 and 4 were lower than those in BLUP, especially within cases with low heritability. These results indicate an appropriate choice of reference population can improve genetic gain and restrict inbreeding even when the reference population size is limited.     

Combined use of phenotypic and genotypic information in sampling animals for genotyping in detection of quantitative trait loci.

Conventional selective genotyping which is using the extreme phenotypes (EP) was compared with alternative criteria to find the most informative animals for genotyping with respects to mapping quantitative trait loci (QTL). Alternative sampling strategies were based on minimizing the sampling error of the estimated QTL effect (MinERR) and maximizing likelihood ratio test (MaxLRT) using both phenotypic and genotypic information. In comparison, animals were randomly genotyped either within or across families. One hundred data sets were simulated each with 30 half-sib families and 120 daughters per family. The strategies were compared in these datasets with respect to estimated effect and position of a QTL within a previously defined genomic region at genotyping 10, 20 or 30% of the animals. Combined linkage disequilibrium linkage analysis (LDLA) was applied in a variance component approach. Power to detect QTL was significantly higher for both MinERR and MaxLRT compared with EP and random genotyping methods (either across or within family), for all the proportions of genotyped animals. Power to detect significant QTL (alpha = 0.01) with 20% genotyping for MinERR and MaxLRT was 80 and 75% of that obtained with complete genotyping compared with 70 and 38% genotyping for EP within and across families respectively. With 30% genotyping, the powers were 78, 83, 78 and 58% respectively. The estimated variance components were unbiased in EP strategies (within and across family), only when at least 30% was genotyped. To decrease the number of genotyped individuals either MinERR or MaxLRT could be considered. With 20% genotyping in MinERR, the estimated QTL variance components were not significant compared with complete genotype information but all studied strategies at 20% genotyping overestimated the QTL effect. Results showed that combining the phenotypic and genotypic information in selective genotyping (e.g. MinERR and MaxLRT) is better than using only the EPs and the combined methods can be considered as alternative approaches to decrease genotyping costs, with unbiased QTL effects, decreased sampling variance of the QTL variance component and also increased the power of QTL detection.     

Divergent selection for immune responsiveness in chickens: estimation of realized heritability with an animal model.

With the aim of improving general disease resistance, chickens were divergently selected for their antibody titers 5 d after immunization with sheep red blood cells for nine generations. Selected and control lines differed significantly for primary and secondary responses after three generations. Heritability of the antibody titer was estimated by REML fitting an animal model using a derivative-free algorithm. The heritability estimate using data on all lines simultaneously was .31. Realized heritability of the antibody titer in the selected lines was estimated by using either the phenotypic cumulative response as the deviation from the control line or the mean breeding values obtained with an animal model. Values from the two methods were consistent, giving a realized heritability of .21 and .25 in the high and low lines, respectively. The genetic trend was not linear and the response to selection tended to accelerate over generations.     

Effect of selection and selective genotyping for creation of reference on bias and accuracy of genomic prediction

Reference populations for genomic selection usually involve selected individuals, which may result in biased prediction of estimated genomic breeding values (GEBV). In a simulation study, bias and accuracy of GEBV were explored for various genetic models with individuals selectively genotyped in a typical nucleus breeding program. We compared the performance of three existing methods, that is, Best Linear Unbiased Prediction of breeding values using pedigree‐based relationships (PBLUP), genomic relationships for genotyped animals only (GBLUP) and a Single‐Step approach (SSGBLUP) using both. For a scenario with no‐selection and random mating (RR), prediction was unbiased. However, lower accuracy and bias were observed for scenarios with selection and random mating (SR) or selection and positive assortative mating (SA). As expected, bias disappeared when all individuals were genotyped and used in GBLUP. SSGBLUP showed higher accuracy compared to GBLUP, and bias of prediction was negligible with SR. However, PBLUP and SSGBLUP still showed bias in SA due to high inbreeding. SSGBLUP and PBLUP were unbiased provided that inbreeding was accounted for in the relationship matrices. Selective genotyping based on extreme phenotypic contrasts increased the prediction accuracy, but prediction was biased when using GBLUP. SSGBLUP could correct the biasedness while gaining higher accuracy than GBLUP. In a typical animal breeding program, where it is too expensive to genotype all animals, it would be appropriate to genotype phenotypically contrasting selection candidates and use a Single‐Step approach to obtain accurate and unbiased prediction of GEBV.     

Long-term responses, changes in genetic variances and inbreeding depression from 122 generations of selection on increased litter size in mice

Data on mice selected for litter size over 122 generations have been analysed in order to reveal the effect of long‐term selection on responses and changes in variances over a long selection period. Originally, three lines were established from the same base population, namely an H line selected for large litter size, an L line selected for small litter size and a K line without selection. In generation 122, the mean number of pups born alive (NBA) was 22 for the H line and 11 for the K line. Phenotypic response to selection is reduced over generations, but crossing of plateaued lines increased responses and realized heritabilities. Both realized heritabilities and heritabilities from residual maximal likelihood (REML) analyses were, in general, calculated from generation (?1)–44 (period 1), 45–70 (period 2) and 71–122 (period 3) separately. Realized heritabilities were in general smaller than heritabilities estimated from mixed model analysis. An overall estimate of heritability for NBA was found to be 0.19 (±0.01) by REML analysis. Additive variance is constant over all periods in the high line and the control line, but is reduced over periods in the low line. The reduction of additive variance in the low line could probably be explained by changes in gene frequencies. In all lines, environmental variances increased over periods. Inbreeding reduced the mean litter size by 0.72 (±0.10) pups per 10% increase in inbreeding, with substantial variance between periods and lines.     

Comparison of restricted maximum likelihood and method R for estimating heritability and predicting breeding value under selection

Method R and Restricted Maximum Likelihood (REML) were compared for estimating heritability (h2) and subsequent prediction of breeding values (a) with data subject to selection. A single-trait animal model was used to generate the data and to predict breeding values. The data originated from 10 sires and 100 dams and simulation progressed for 10 overlapping generations. In simulating the data, genetic evaluation used the underlying parameter values and sires and dams were chosen by truncation selection for greatest predicted breeding values. Four alternative pedigree structures were evaluated: complete pedigree information, 50% of phenotypes with sire identities missing, 50% of phenotypes with dam identities missing, and 50% of phenotypes with sire and dams identities missing. Under selection and with complete pedigree data, Method R was a slightly less consistent estimator of h2 than REML. Estimates of h2 by both methods were biased downward when there was selection and loss of pedigree information and were unbiased when no selection was practiced. The empirical mean square error (EMSE) of Method R was several times larger than the EMSE of REML. In a subsequent analysis, different combinations of generations selected and generations sampled were simulated in an effort to disentangle the effects of both factors on Method R estimates of h2. It was observed that Method R overestimated h2 when both the sampling that is intrinsic in the method and the selection occurred in generations 6 to 10. In a final experiment, BLUP(a) were predicted with h2 estimated by either Method R or REML. Subsequently, five more generations of selection were practiced, and the mean square error of prediction (MSEP) of BLUP(a) was calculated with estimated h2 by either method, or the true value of the parameter. The MSEP of empirical BLUP(a) using Method R was greater than the MSEP of empirical BLUP(a) using REML. The latter statistic was closer to prediction error variance of BLUP(a) than the MSEP of empirical BLUP(a) using Method R, indicating that empirical BLUP(a) calculated using REML produced accurate predictions of breeding values under selection. In conclusion, the variability of h2 estimates calculated with Method R was greater than the variability of h2 estimates calculated with REML, with or without selection. Also, the MSEP of EBLUP(a) calculated using estimates of h2 by Method R was larger than MSEP of EBLUP(a) calculated with REML estimates of h2.     

Effects of nonrandom parental selection on estimation of variance components

Bayesian estimation via Gibbs sampling, REML, and Method R were compared for their empirical sampling properties in estimating genetic parameters from data subject to parental selection using an infinitesimal animal model. Models with and without contemporary groups, random or nonrandom parental selection, two levels of heritability, and none or 15% randomly missing pedigree information were considered. Nonrandom parental selection caused similar effects on estimates of variance components from all three methods. When pedigree information was complete, REML and Bayesian estimation were not biased by nonrandom parental selection for models with or without contemporary groups. Method R estimates, however, were strongly biased by nonrandom parental selection when contemporary groups were in the model. The bias was empirically shown to be a consequence of not fully accounting for gametic phase disequilibrium in the subsamples. The joint effects of nonrandom parental selection and missing pedigree information caused estimates from all methods to be highly biased. Missing pedigree information did not cause biased estimates in random mating populations. Method R estimates usually had greater mean square errors than did REML and Bayesian estimates.     

Genotyping of groups of cows to improve genomic breeding values of new traits

Genotyping females and including them into the reference set for genomic predictions in dairy cattle is considered to provide gains in reliabilities of estimated breeding values for selection candidates. This should especially be true for low heritability traits. By the use of simulation, we extended a genomic reference set for an existing trait by including a fixed number of genotyped first‐crop daughters for one or two generations of reference sires. Moreover, we calculated results for the effects of a similar strategy in a situation where for a new trait the recording of phenotypes has recently started. For this case, we compared the effect of two different genotyping strategies: first, to phenotype cows but to genotype their sires only, and second, to collect phenotypes and genotypes on the same cows. We studied the effects on validation reliabilities and unbiasedness of predicted values for selection candidates. We found that by extending the reference set with genotyped daughters it is possible to increase the validation reliability of genomic breeding values. In the case of a new trait, it is always better to collect and use genotypes and phenotypes on the same animals instead of using only sire genotypes. We found that the benefits that can be achieved are sensitive to the sampling strategy used when selecting females for genotyping.     

基因组选择及其在羊育种中的应用研究进展

数量性状是羊育种中的重要性状,受微效多基因控制、遗传力低,而传统育种方法难以提高羊的育种效率。提高动物育种效率对于选种选配工作和经济生产效益至关重要。随着育种新技术的不断革新与发展,基因组选择(genomic selection,GS)方法已成为育种技术中强大的工具,且已成功运用于个体经济价值较大的物种中,其具有缩短世代间隔、提高育种准确性、减少生产成本、提高畜禽经济效益等优势。近年来,由于基因组技术的不断成熟和各个统计模型的升级优化,以及高密度SNP芯片价格的下调,报告有关于基因组选择育种的实证和模拟研究层出不穷,且基因组选择技术已在羊育种中逐步开展,特别是在羊的重要性状中已有不少报道。由于羊的品种较多,地方性状差异化较大,个体经济价值略低,尽管基因组育种的新技术已经非常成熟,但目前仍没有在羊育种中大范围普及。为了更全面地了解该技术在羊育种中的研究现状,且基于选种选配的重要地位,作者就基因组选择在羊育种中的研究进展展开综述,主要从表型测定、基因分型、不同模型方面介绍了基因组选择在羊的重要性状中的应用和现状,讨论了其优势与挑战,并展望了基因组选择的未来发展方向。     

Implications of SNP weighting on single‐step genomic predictions for different reference population sizes

We investigated the importance of SNP weighting in populations with 2,000 to 25,000 genotyped animals. Populations were simulated with two effective sizes (20 or 100) and three numbers of QTL (10, 50 or 500). Pedigree information was available for six generations; phenotypes were recorded for the four middle generations. Animals from the last three generations were genotyped for 45,000 SNP. Single‐step genomic BLUP (ssGBLUP) and weighted ssGBLUP (WssGBLUP) were used to estimate genomic EBV using a genomic relationship matrix ( G ). The WssGBLUP performed better in small genotyped populations; however, any advantage for WssGBLUP was reduced or eliminated when more animals were genotyped. WssGBLUP had greater resolution for genome‐wide association (GWA) as did increasing the number of genotyped animals. For few QTL, accuracy was greater for WssGBLUP than ssGBLUP; however, for many QTL, accuracy was the same for both methods. The largest genotyped set was used to assess the dimensionality of genomic information (number of effective SNP). The number of effective SNP was considerably less in weighted G than in unweighted G . Once the number of independent SNP is well represented in the genotyped population, the impact of SNP weighting becomes less important.     

Multiple-trait restricted maximum likelihood for simulated measures of ovulation rate with underlying multivariate normal distributions.

A data set that was used to estimate covariance components with REML for an animal model with eight measures of ovulation rate treated as separate traits was used as a template to simulate data sets of eight multivariate normal traits that were then truncated to binomial traits. The model for simulation included eight measures on 610 animals with 1,071 animals in the numerator relationship matrix. Heritabilities were equal for the eight measures, and both genetic and phenotypic correlations among the measures were equal. Ten replications for each combination of heritability (.15, .25, and .35) and genetic correlation (.50, .66, and .90) were simulated on the normal scale. For each replicate, estimates of the eight heritabilities and 28 genetic correlations were obtained by multiple-trait REML. The usual transformation of heritability estimated on the binomial scale overestimated heritability on the normal scale. Genetic correlations on the binomial scale seriously underestimated the correlations on the normal scale. Standard errors of the estimates obtained by replication were somewhat larger than the approximate SE from REMLPK (the multi-trait REML program of K. Meyer). A final set of 10 simulated replications with heritability of .25 and genetic correlation of 1.00 resulted in average estimates of .18 for heritability and of .66 for genetic correlation that agree closely with those from the analysis of measures of ovulation at eight estrous cycles used as a template; averages for heritability of .16 and for genetic correlation of .66 were obtained.     

基于元共祖的基因组联合育种模拟研究

旨在将整合元共祖的一步法(single-step genomic best linear unbiased prediction with metafounders,MF-SSGBLUP)应用到基因组联合育种中,并与其他经典基因组选择方法进行比较分析。本研究使用QMSim软件模拟3个系谱相互独立的奶牛群体;分别使用广义最小二乘法(generalized least squares,GLS)和原始方法(naïve,NAI)估计不同群体间的祖先关系矩阵Γ;将MF-SSGBLUP、SSGBLUP和BLUP用于3个模拟群体的联合育种,评估各方法在遗传参数和育种值估计方面的差异。在不同遗传力下,GLS所得的Γ矩阵在对角线元素上略低于NAI法,在非对角线元素上没有明显差异,且基因组关系矩阵与基于元共祖构建的亲缘关系矩阵对角线元素相关系数(0.750~0.775)高于基因组关系矩阵与传统的亲缘关系矩阵相关系数(0.508~0.572)。MF-SSGBLUP遗传力估计值(0.138、0.140、0.297和0.298)与当代群体遗传力(0.107和0.296)的偏差小于其余两种方法(0.145、0.173、0.273和0.340),且MF-SSGBLUP估计育种值准确性(0.888~0.908)高于SSGBLUP法(0.863~0.876)和BLUP法(0.854~0.871)。表明,MF-SSGBLUP的遗传参数估计值无偏性更好,估计育种值准确性更高。根据上述模拟数据结果表明,在联合育种中,整合元共祖的基因组选择方法优于其他经典基因组选择方法。     

Effect of genotyped bulls with different numbers of phenotyped progenies on quantitative trait loci detection and genomic evaluation in a simulated cattle population

The objective of this study was to assess the effect of genotyped bulls with different numbers of phenotyped progenies on quantitative trait loci (QTL) detection and genomic evaluation using a simulated cattle population. Twelve generations (G1–G12) were simulated from the base generation (G0). The recent population had different effective population sizes, heritability, and number of QTL. G0–G4 were used for pedigree information. A total of 300 genotyped bulls from G5–G10 were randomly selected. Their progenies were generated in G6–G11 with different numbers of progeny per bull. Scenarios were considered according to the number of progenies and whether the genotypes were possessed by the bulls or the progenies. A genome‐wide association study and genomic evaluation were performed with a single‐step genomic best linear unbiased prediction method to calculate the power of QTL detection and the genomic estimated breeding value (GEBV). We found that genotyped bulls could be available for QTL detection depending on conditions. Additionally, using a reference population, including genotyped bulls, which had more progeny phenotypes, enabled a more accurate prediction of GEBV. However, it is desirable to have more than 4,500 individuals consisting of both genotypes and phenotypes for practical genomic evaluation.     

The empirical bias of estimates by restricted maximum likelihood, Bayesian method, and method R under selection for additive, maternal, and dominance models.

Bayesian analysis via Gibbs sampling, restricted maximum likelihood (REML), and Method R were used to estimate variance components for several models of simulated data. Four simulated data sets that included direct genetic effects and different combinations of maternal, permanent environmental, and dominance effects were used. Parents were selected randomly, on phenotype across or within contemporary groups, or on BLUP of genetic value. Estimates by Bayesian analysis and REML were always empirically unbiased in large data sets. Estimates by Method R were biased only with phenotypic selection across contemporary groups; estimates of the additive variance were biased upward, and all the other estimates were biased downward. No empirical bias was observed for Method R under selection within contemporary groups or in data without contemporary group effects. The bias of Method R estimates in small data sets was evaluated using a simple direct additive model. Method R gave biased estimates in small data sets in all types of selection except BLUP. In populations where the selection is based on BLUP of genetic value or where phenotypic selection is practiced mostly within contemporary groups, estimates by Method R are likely to be unbiased. In this case, Method R is an alternative to single-trait REML and Bayesian analysis for analyses of large data sets when the other methods are too expensive to apply.     

Genomic selection for maternal traits in pigs

The aim of this study was to compare alternative designs for implementation of genomic selection to improve maternal traits in pigs, with a conventional breeding scheme and a progeny testing scheme. The comparison was done through stochastic simulation of a pig population. It was assumed that selection was performed based on a trait that could be measured on females after the first litter, with a heritability of 0.1. Genomic selection increased genetic gain and reduced the rate of inbreeding, compared with conventional selection without progeny testing. Progeny testing could also increase genetic gain and decrease the rate of inbreeding, but because of the increased generation interval, the increase in annual genetic gain was only 7%. When genomic selection was applied, genetic gain was increased by 23 to 91%, depending on which and how many animals were genotyped. Genotyping dams in addition to the male selection candidates gave increased accuracy of the genomic breeding values, increased genetic gain, and decreased rate of inbreeding. To genotype 2 or 3 males from each litter, in order to perform within-litter selection, increased genetic gain 8 to 12%, compared with schemes with the same number of genotyped females but only 1 male candidate per litter. Comparing schemes with the same total number of genotyped animals revealed that genotyping more females caused a greater increase in genetic gain than genotyping more males because greater accuracy of selection was more advantageous than increasing the number of male selection candidates. When more than 1 male per litter was genotyped, and thereby included as selection candidates, rate of inbreeding increased because of coselection of full sibs. The conclusion is that genomic selection can increase genetic gain for traits that are measured on females, which includes several traits with economic importance in maternal pig breeds. Genotyping females is essential to obtain a high accuracy of selection.     

Detecting effective starting point of genomic selection by divergent trends from best linear unbiased prediction and single-step genomic best linear unbiased prediction in pigs,beef cattle,and broilers

Genomic selection has been adopted nationally and internationally in different livestock and plant species. However, understanding whether genomic selection has been effective or not is an essential question for both industry and academia. Once genomic evaluation started being used, estimation of breeding values with pedigree best linear unbiased prediction (BLUP) became biased because this method does not consider selection using genomic information. Hence, the effective starting point of genomic selection can be detected in two possible ways including the divergence of genetic trends and Realized Mendelian sampling (RMS) trends obtained with BLUP and single-step genomic BLUP (ssGBLUP). This study aimed to find the start date of genomic selection for a set of economically important traits in three livestock species by comparing trends obtained using BLUP and ssGBLUP. Three datasets were used for this purpose: 1) a pig dataset with 117k genotypes and 1.3M animals in pedigree, 2) an Angus cattle dataset consisted of ~842k genotypes and 11.5M animals in pedigree, and 3) a purebred broiler chicken dataset included ~154k genotypes and 1.3M birds in pedigree were used. The genetic trends for pigs diverged for the genotyped animals born in 2014 for average daily gain (ADG) and backfat (BF). In beef cattle, the trends started diverging in 2009 for weaning weight (WW) and in 2016 for postweaning gain (PWG), with little divergence for birth weight (BTW). In broiler chickens, the genetic trends estimated by ssGBLUP and BLUP diverged at breeding cycle 6 for two out of the three production traits. The RMS trends for the genotyped pigs diverged for animals born in 2014, more for ADG than for BF. In beef cattle, the RMS trends started diverging in 2009 for WW and in 2016 for PWG, with a trivial trend for BTW. In broiler chickens, the RMS trends from ssGBLUP and BLUP diverged strongly for two production traits at breeding cycle 6, with a slight divergence for another trait. Divergence of the genetic trends from ssGBLUP and BLUP indicates the onset of the genomic selection. The presence of trends for RMS indicates selective genotyping, with or without the genomic selection. The onset of genomic selection and genotyping strategies agrees with industry practices across the three species. In summary, the effective start of genomic selection can be detected by the divergence between genetic and RMS trends from BLUP and ssGBLUP.     

Comparison between genomic predictions using daughter yield deviation and conventional estimated breeding value as response variables

This study compared genomic predictions using conventional estimated breeding values (EBV) and daughter yield deviations (DYD) as response variables based on simulated data. Eight scenarios were simulated in regard to heritability (0.05 and 0.30), number of daughters per sire (30, 100, and unequal numbers with an average of 100 per sire) and numbers of genotyped sires (all or half of sires were genotyped). The simulated genome had a length of 1200 cM with 15,000 equally spaced Single-nucleotide polymorphism (SNP) markers and 500 randomly distributed Quantitative trait locus (QTL). In the simulated scenarios, the EBV approach was as effective as or slightly better than the DYD approach at predicting breeding value, dependent on simulated scenarios and statistical models. Applying a Bayesian common prior model (the same prior distribution of marker effect variance) and a linear mixed model (GBLUP), the EBV and DYD approaches provided similar genomic estimated breeding value (GEBV) reliabilities, except for scenarios with unequal numbers of daughters and half of sires without genotype, for which the EBV approach was superior to the DYD approach (by 1.2 and 2.4%). Using a Bayesian mixture prior model (mixture prior distribution of marker effect variance), the EBV approach resulted in slightly higher reliabilities of GEBV than the DYD approach (by 0.3-3.6% with an average of 1.9%), and more obvious in scenarios with low heritability, small or unequal numbers of daughters, and half of sires without genotype. Moreover, the results showed that the correlation between GEBV and conventional parent average (PA) was lower (corresponding to a relatively larger gain by including PA) when using the DYD approach than when using the EBV approach. Consequently, the two approaches led to similar reliability of an index combining GEBV and PA in most scenarios. These results indicate that EBV can be used as an alternative response variable for genomic prediction.     

Response to 13 generations of selection for increased 8-week body weight in lines of mice carrying a sheep growth hormone-based transgene

The purpose of this study was to evaluate selection in lines of transgenic mice. Two replicates of lines that either carried or did not carry the sheep metallothionein-1a sheep growth hormone transgene (oMt1a-oGH) were established. The host lines had been previously selected for rapid growth or selected randomly. Within-litter selection for increased 8-wk body weight was carried out for 13 generations. The frequency of oMt1a-oGH was monitored in all generations in the transgenic lines, but no genotypic information regarding the transgene was used as an aid to selection. The oMt1a-oGH was activated from weaning, at 3 wk, until 8 wk of age by adding ZnSO4 to the drinking water. Zinc stimulation of the transgene was not done during mating, gestation, or lactation. Data on body weights and weight gains were analyzed with a conventional mixed model and with an animal model. Genetic progress was achieved in all lines subjected to directional selection. In the control background, response to selection for 8-wk body weight was larger in the nontransgenic lines than in the transgenic lines, whereas no difference was found in the selected background. The frequency of the transgene was increased from the initial .5 to .62 in the randomly selected background but decreased to .04 in lines from a selected background. The REML estimates of variance components and genetic gain estimates varied greatly between the two methods. In general, there was better agreement between the realized heritability estimates and the heritability estimates obtained from the conventional mixed model analysis than between realized heritability estimates and results obtained using the animal model. Favorable correlated responses were obtained for 3- and 6-wk body weights and on 3- to 6- and 6- to 8-wk weight gains. Correlated responses to selection were larger in the selected than in the nonselected background but were not affected by the presence of the transgene. Results suggest that constructs similar to the oMt1a-oGH, which allow tight regulation, may be successfully incorporated into commercial livestock and should have larger effects in populations that have not been subject to selection.     

Optimum contribution selection using traditional best linear unbiased prediction and genomic breeding values in aquaculture breeding schemes

The aim of this study was to compare genetic gain for a traditional aquaculture sib breeding scheme with breeding values based on phenotypic data (TBLUP) with a breeding scheme with genome-wide (GW) breeding values. Both breeding schemes were closed nuclei with discrete generations modeled by stochastic simulation. Optimum contribution selection was applied to restrict pedigree-based inbreeding to either 0.5 or 1% per generation. There were 1,000 selection candidates and a sib test group of either 4,000 or 8,000 fish. The number of selected dams and sires to create full sib families in each generation was determined from the optimum contribution selection method. True breeding values for a trait were simulated by summing the number of each QTL allele and the true effect of each of the 1,000 simulated QTL. Breeding values in TBLUP were predicted from phenotypic and pedigree information, whereas genomic breeding values were computed from genetic markers whose effects were estimated using a genomic BLUP model. In generation 5, genetic gain was 70 and 74% greater for the GW scheme than for the TBLUP scheme for inbreeding rates of 0.5 and 1%. The reduction in genetic variance was, however, greater for the GW scheme than for the TBLUP scheme due to fixation of some QTL. As expected, accuracy of selection increased with increasing heritability (e.g., from 0.77 with a heritability of 0.2 to 0.87 with a heritability of 0.6 for GW, and from 0.53 and 0.58 for TBLUP in generation 5 with sib information only). When the trait was measured on the selection candidate compared with only on sibs and the heritability was 0.4, accuracy increased from 0.55 to 0.69 for TBLUP and from 0.83 to 0.86 for GW. The number of selected sires to get the desired rate of inbreeding was in general less in GW than in TBLUP and was 33 for GW and 83 for TBLUP (rate of inbreeding 1% and heritability 0.4). With truncation selection, genetic gain for the scheme with GW breeding values was nearly twice as large as a scheme with traditional BLUP breeding values. The results indicate that the benefits of applying GW breeding values compared with TBLUP are reduced when contributions are optimized. In conclusion, genetic gain in aquaculture breeding schemes with optimized contributions can increase by as much as 81% by applying genome-wide breeding values compared with traditional BLUP breeding values.     

Increased litter size in Rambouillet sheep: I. Estimation of genetic parameters.

The variance and covariance components needed to estimate heritabilities of and genetic correlations among litter size, ovulation rate, scrotal circumference, and BW in a flock of Rambouillet sheep were estimated using REML via an expectation-maximization type algorithm. The heritability estimates from univariate analyses were .14, .21, .25, .36, and .15 for litter size, ovulation rate, scrotal circumference, 180-d BW of females, and 180-d BW of males, respectively, and average heritability estimates from bivariate analyses were .19, .20, .20, .34, and .10 for litter size, ovulation rate, scrotal circumference, 180-d BW of females, and 180-d BW of males, respectively. The genetic correlation between litter size and ovulation rate was near unity. Body weight in ewes had a moderate genetic correlation with both litter size (.22) and ovulation rate (.20) and a low residual correlation with both litter size (.03) and ovulation rate (.09). The genetic correlation between BW in rams and scrotal circumference was 0, whereas the residual correlation was .71. The genetic correlations of scrotal circumference with litter size and ovulation rate were -.25 and +.20, respectively.