Genetic heterogeneity of within-family variance of body weight in Atlantic salmon (Salmo salar)

Genomic prediction for grain yield and micro-environmental sensitivity in winter wheat

Individuals within a common environment experience variations due to unique and non-identifiable micro-environmental factors. Genetic sensitivity to micro-environmental variation (i.e. micro-environmental sensitivity) can be identified in residuals, and genotypes with lower micro-environmental sensitivity can show greater resilience towards environmental perturbations. Micro-environmental sensitivity has been studied in animals; however, research on this topic is limited in plants and lacking in wheat. In this article, we aimed to (i) quantify the influence of genetic variation on residual dispersion and the genetic correlation between genetic effects on (expressed) phenotypes and residual dispersion for wheat grain yield using a double hierarchical generalized linear model (DHGLM); and (ii) evaluate the predictive performance of the proposed DHGLM for prediction of additive genetic effects on (expressed) phenotypes and its residual dispersion. Analyses were based on 2,456 advanced breeding lines tested in replicated trials within and across different environments in Denmark and genotyped with a 15K SNP-Illumina-BeadChip. We found that micro-environmental sensitivity for grain yield is heritable, and there is potential for its reduction. The genetic correlation between additive effects on (expressed) phenotypes and dispersion was investigated, and we observed an intermediate correlation. From these results, we concluded that breeding for reduced micro-environmental sensitivity is possible and can be included within breeding objectives without compromising selection for increased yield. The predictive ability and variance inflation for predictions of the DHGLM and a linear mixed model allowing heteroscedasticity of residual variance in different environments (LMM-HET) were evaluated using leave-one-line-out cross-validation. The LMM-HET and DHGLM showed good and similar performance for predicting additive effects on (expressed) phenotypes. In addition, the accuracy of predicting genetic effects on residual dispersion was sufficient to allow genetic selection for resilience. Such findings suggests that DHGLM may be a good choice to increase grain yield and reduce its micro-environmental sensitivity.

Capturing indirect genetic effects on phenotypic variability: Competition meets canalization

Phenotypic variability of a genotype is relevant both in natural and domestic populations. In the past two decades, variability has been studied as a heritable quantitative genetic trait in its own right, often referred to as inherited variability or environmental canalization. So far, studies on inherited variability have only considered genetic effects of the focal individual, that is, direct genetic effects on inherited variability. Observations from aquaculture populations and some plants, however, suggest that an additional source of genetic variation in inherited variability may be generated through competition. Social interactions, such as competition, are often a source of Indirect Genetic Effects (IGE). An IGE is a heritable effect of an individual on the trait value of another individual. IGEs may substantially affect heritable variation underlying the trait, and the direction and magnitude of response to selection. To understand the contribution of IGEs to evolution of environmental canalization in natural populations, and to exploit such inherited variability in animal and plant breeding, we need statistical models to capture this effect. To our knowledge, it is unknown to what extent the current statistical models commonly used for IGE and inherited variability capture the effect of competition on inherited variability. Here, we investigate the potential of current statistical models for inherited variability and trait values, to capture the direct and indirect genetic effects of competition on variability. Our results show that a direct model of inherited variability almost entirely captures the genetic sensitivity of individuals to competition, whereas an indirect model of inherited variability captures the cooperative genetic effects of individuals on their partners. Models for trait levels, however, capture only a small part of the genetic effects of competition. The estimation of direct and indirect genetic effects of competition, therefore, is possible with models for inherited variability but may require a two‐step analysis.

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Genetic parameters for uniformity of harvest weight in Pacific white shrimp (Litopenaeus vannamei)

Background

Uniformity of body weight is a trait of great economic importance in the production of white shrimp (Litopenaeus vannamei). A necessary condition to improve this trait through selective breeding is the existence of genetic variability for the environmental variance of body weight. Although several studies have reported such variability in other aquaculture species, to our knowledge, no estimates are available for shrimp. Our aim in this study was to estimate the genetic variance for weight uniformity in a farmed population of shrimp to determine the potential of including this trait in the selection program. We also estimated the genetic correlation of weight uniformity between two environments (selection nucleus and commercial population).

Methods

The database contained phenotypic records for body weight on 51,346 individuals from the selection nucleus and 38,297 individuals from the commercial population. A double hierarchical generalized linear model was used to analyse weight uniformity in the two environments. Fixed effects included sex and year for the nucleus data and sex and year-pond combination for the commercial data. Environmental and additive genetic effects were included as random effects.

Results

The estimated genetic variance for weight uniformity was greater than 0 (0.06 ± 0.01) in both the nucleus and commercial populations and the genetic coefficient of variation for the residual variance was 0.25 ± 0.01. The genetic correlation between weight and weight uniformity was close to zero in both environments. The estimate of the genetic correlation of weight uniformity between the two environments (selection nucleus and commercial population) was 0.64 ± 0.06.

Conclusions

The existence of genetic variance for weight uniformity suggests that genetic improvement of this trait is possible. Selection for weight uniformity should not decrease weight, given the near zero genetic correlation between these two traits. The strong genetic correlation of weight uniformity between the two environments indicates that response to selection for uniformity in the nucleus will be at least partially transmitted to the commercial population if this trait is included in the breeding goal.

Genetic parameters for canalization analysis of morphological traits in the Pura Raza Español horse

Measurements from 13 different morphological traits of importance in the Pura Raza Español (PRE) horse were used to estimate genetic and environmental parameters following a heteroscedastic model in which data were assigned to stallions. Data sets used ranged from 20,610 (height at withers) to 48,486 measurements (length of shoulder), and the number of animals analysed in the pedigrees varied from 17,662 (height at withers) to 23,962 (dorsal-sternal diameter). Results of heritabilities of the traits varied from 0.09 (width of chest and upper neck line) to 0.30 (muscular development). Further, genetic correlations between traits and their environmental variability were estimated, obtaining values from -0.56 (muscular development) to 0.69 (height at withers). Also, predicted breeding values for the mean and for the environmental variability were obtained for all horses in the pedigrees, providing individual information about not only the expected phenotypic value of their offspring but also the expected heterogeneity among them. Results proved the possibility of improving morphological traits and reducing the heterogeneity of offspring at a time by the selection of animals and levels of systematic effects.

Genetic variance for behavioural 'predictability' of stress response

Genetic factors underpinning phenotypic variation are required if natural selection is to result in adaptive evolution. However, evolutionary and behavioural ecologists typically focus on variation among individuals in their average trait values and seek to characterize genetic contributions to this. As a result, less attention has been paid to if and how genes could contribute towards within-individual variance or trait 'predictability'. In fact, phenotypic 'predictability' can vary among individuals, and emerging evidence from livestock genetics suggests this can be due to genetic factors. Here, we test this empirically using repeated measures of a behavioural stress response trait in a pedigreed population of wild-type guppies. We ask (a) whether individuals differ in behavioural predictability and (b) whether this variation is heritable and so evolvable under selection. Using statistical methodology from the field of quantitative genetics, we find support for both hypotheses and also show evidence of a genetic correlation structure between the behavioural trait mean and individual predictability. We show that investigating sources of variability in trait predictability is statistically tractable and can yield useful biological interpretation. We conclude that, if widespread, genetic variance for 'predictability' will have major implications for the evolutionary causes and consequences of phenotypic variation.

Estimation of genetic parameters for micro-environmental sensitivities of production traits in Holstein cows using two-step method

Context The request for more uniform animal products, which is motivated chiefly by economic reasons, has enhanced the interest in decreasing variability of characters via selection. In the present dairy operation, breeding dairy cows which have strong resistance against environmental changes for main traits is very important. Aims The aim of this study was to estimate genetic parameters for heterogeneity of residual variance in milk yield and composition of Iranian Holstein cows. Methods The dataset included 305-day production records of cows which were provided by the Animal Breeding Center and Promotion of Animal Products of Iran between 1983 and 2014. In two-step method, univariate analyses were conducted to estimate variance components for 305-day production traits. Then, genetic variability of residual variances was estimated. Key results Estimates of heritability for micro-environmental sensitivities of milk, fat and protein yields in the first three lactations of Holstein cows were low and equal to 0.043, 0.028 and 0.039; 0.031, 0.019 and 0.024; 0.027, 0.016 and 0.019 respectively. Considerable genetic coefficient of variations of residual variance for above mentioned traits (0.261, 0.247 and 0.218; 0.221, 0.204 and 0.194; 0.219, 0.199 and 0.178 respectively) indicated significant additive genetic variation for micro-environmental sensitivities. Conclusions The results of this study indicate that micro-environmental sensitivities were present for milk production traits of Iranian Holsteins. High genetic coefficient of variation for micro-environmental sensitivities indicated the possibility of reducing environmental variation and increase in uniformity via selection. Implications Reduction of environmental sensitivities would increase the predicted performance of animals and decreased corresponding threats for dairy farmers.

Genetics and genomics of uniformity and resilience in livestock and aquaculture species: A review

Genetic control of residual variance offers opportunities to increase uniformity and resilience of livestock and aquaculture species. Improving uniformity and resilience of animals will improve health and welfare of animals and lead to more homogenous products. Our aims in this review were to summarize the current models and methods to study genetic control of residual variance, genetic parameters and genomic results for residual variance and discuss future research directions. Typically, the genetic coefficient of variation is high (median = 0.27; range 0-0.86) and the heritability of residual variance is low (median = 0.01; range 0-0.10). Higher heritabilities can be achieved when increasing the number of records per animal. Divergent selection experiments have supported the feasibility of selecting for high or low residual variance. Genomic studies have revealed associations in regions related to stress, including those from the heat shock protein family. Although the number of studies is growing, genetic control of residual variance is still poorly understood, but big data and genomics offer great opportunities.

Economic optimization of full-sib test group size and genotyping effort in a breeding program for Atlantic salmon

BACKGROUND

Breeding companies may want to maximize the rate of genetic gain from their breeding program within a limited budget. In salmon breeding programs, full-sibs of selection candidates are subjected to performance tests for traits that cannot be recorded on selection candidates. While marginal gains in the aggregate genotype from phenotyping and genotyping more full-sibs per candidate decrease, costs increase linearly, which suggests that there is an optimum in the allocation of the budget among these activities. Here, we studied how allocation of the fixed budget to numbers of phenotyped and genotyped test individuals in performance tests can be optimized.

METHODS

Gain in the aggregate genotype was a function of the numbers of full-sibs of selection candidates that were (1) phenotyped in a challenge test for sea lice resistance (2) phenotyped in a slaughter test (3) genotyped in the challenge test, and (4) genotyped in the slaughter test. Each of these activities was subject to budget constraints. Using a grid search, we optimized allocation of the budget among activities to maximize gain in the aggregate genotype. We performed sensitivity analyses on the maximum gain in the aggregate genotype and on the relative allocation of the budget among activities at the optimum.

RESULTS

Maximum gain in the aggregate genotype was €386/ton per generation. The response surface for gain in the aggregate genotype was rather flat around the optimum, but it curved strongly near the extremes. Maximum gain was sensitive to the size of the budget and the relative emphasis on breeding goal traits, but less sensitive to the accuracy of genomic prediction and costs of phenotyping and genotyping. The relative allocation of budget among activities at the optimum was sensitive to costs of phenotyping and genotyping and the relative emphasis on breeding goal traits, but was less sensitive to the accuracy of genomic prediction and the size of the budget.

CONCLUSIONS

There is an optimum allocation of budget to the numbers of full-sibs of selection candidates that are phenotyped and genotyped in performance tests that maximizes gain in the aggregate genotype. Although potential gains from optimizing group sizes and genotyping effort may be small, they come at no extra cost.

Estimation of genetic variation for macro- and micro-environmental sensitivities of milk yield and composition in Holstein cows using double hierarchical generalized linear models

The aim of this study was to estimate genetic parameters for environmental sensitivities in milk yield and composition of Iranian Holstein cows using the double hierarchical generalized linear model (DHGLM) method. Data set included test-day productive records of cows which were provided by the Animal Breeding Center and Promotion of Animal Products of Iran during 1983 to 2014. In the DHGLM method, a random regression model was fitted which included two parts of mean and residual variance. A random regression model (mean model) and a residual variance model were used to study the genetic variation of micro-environmental sensitivities. In order to consider macro-environmental sensitivities, DHGLM was extended using a reaction norm model, and a sire model was applied. Based on the mean model, additive genetic variances for the mean were 38.25 for milk yield, 0.23 for fat yield and 0.03 for protein yield in the first lactation, respectively. Based on the residual variance model, additive genetic variances for residual variance were 0.039 for milk yield, 0.030 for fat yield and 0.020 for protein yield in the first lactation, respectively. Estimates of genetic correlation between milk yield and macro- and micro-environmental sensitivities were 0.660 and 0.597 in the first lactation, respectively. The results of this study indicated that macro- and micro-environmental sensitivities were present for milk production traits of Iranian Holsteins. High genetic coefficient of variation for micro-environmental sensitivities indicated the possibility of reducing environmental variation and increase in uniformity via selection.

Genomic regions underlying uniformity of yearling weight in Nellore cattle evaluated under different response variables

BACKGROUND

In livestock, residual variance has been studied because of the interest to improve uniformity of production. Several studies have provided evidence that residual variance is partially under genetic control; however, few investigations have elucidated genes that control it. The aim of this study was to identify genomic regions associated with within-family residual variance of yearling weight (YW; N = 423) in Nellore bulls with high density SNP data, using different response variables. For this, solutions from double hierarchical generalized linear models (DHGLM) were used to provide the response variables, as follows: a DGHLM assuming non-null genetic correlation between mean and residual variance (rmv ≠ 0) to obtain deregressed EBV for mean (dEBVm) and residual variance (dEBVv); and a DHGLM assuming rmv = 0 to obtain two alternative response variables for residual variance, dEBVv_r0 and log-transformed variance of estimated residuals (ln_[Formula: see text]).

RESULTS

The dEBVm and dEBVv were highly correlated, resulting in common regions associated with mean and residual variance of YW. However, higher effects on variance than the mean showed that these regions had effects on the variance beyond scale effects. More independent association results between mean and residual variance were obtained when null rmv was assumed. While 13 and 4 single nucleotide polymorphisms (SNPs) showed a strong association (Bayes Factor > 20) with dEBVv and ln_[Formula: see text], respectively, only suggestive signals were found for dEBVv_r0. All overlapping 1-Mb windows among top 20 between dEBVm and dEBVv were previously associated with growth traits. The potential candidate genes for uniformity are involved in metabolism, stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation.

CONCLUSIONS

It is necessary to use a strategy like assuming null rmv to obtain genomic regions associated with uniformity that are not associated with the mean. Genes involved not only in metabolism, but also stress, inflammatory and immune responses, mineralization, neuronal activity and bone formation were the most promising biological candidates for uniformity of YW. Although no clear evidence of using a specific response variable was found, we recommend consider different response variables to study uniformity to increase evidence on candidate regions and biological mechanisms behind it.

Breeding for robustness: investigating the genotype-by-environment interaction and micro-environmental sensitivity of Genetically Improved Farmed Tilapia (Oreochromis niloticus)

Robustness has become a highly desirable breeding goal in the globalized agricultural market. Both genotype‐by‐environment interaction (G × E) and micro‐environmental sensitivity are important robustness components of aquaculture production, in which breeding stock is often disseminated to different environments. The objectives of this study were (i) to quantify the degree of G × E by assessing the growth performance of Genetically Improved Farmed Tilapia (GIFT) across three countries (Malaysia, India and China) and (ii) to quantify the genetic heterogeneity of environmental variance for body weight at harvest (BW) in GIFT as a measure of micro‐environmental sensitivity. Selection for BW was carried out for 13 generations in Malaysia. Subsets of 60 full‐sib families from Malaysia were sent to China and India after five and nine generations respectively. First, a multi‐trait animal model was used to analyse the BW in different countries as different traits. The results indicate a strong G × E. Second, a genetically structured environmental variance model, implemented using Bayesian inference, was used to analyse micro‐environmental sensitivity of BW in each country. The analysis revealed the presence of genetic heterogeneity of both BW and its environmental variance in all environments. The presence of genetic variation in residual variance of BW implies that the residual variance can be modified by selection. Incorporating both G × E and micro‐environmental sensitivity information may help in selecting robust genotypes with high performance across environments and resilience to environmental fluctuations.

Modelling the co-evolution of indirect genetic effects and inherited variability

When individuals interact, their phenotypes may be affected not only by their own genes but also by genes in their social partners. This phenomenon is known as Indirect Genetic Effects (IGEs). In aquaculture species and some plants, however, competition not only affects trait levels of individuals, but also inflates variability of trait values among individuals. In the field of quantitative genetics, the variability of trait values has been studied as a quantitative trait in itself, and is often referred to as inherited variability. Such studies, however, consider only the genetic effect of the focal individual on trait variability and do not make a connection to competition. Although the observed phenotypic relationship between competition and variability suggests an underlying genetic relationship, the current quantitative genetic models of IGE and inherited variability do not allow for such a relationship. The lack of quantitative genetic models that connect IGEs to inherited variability limits our understanding of the potential of variability to respond to selection, both in nature and agriculture. Models of trait levels, for example, show that IGEs may considerably change heritable variation in trait values. Currently, we lack the tools to investigate whether this result extends to variability of trait values. Here we present a model that integrates IGEs and inherited variability. In this model, the target phenotype, say growth rate, is a function of the genetic and environmental effects of the focal individual and of the difference in trait value between the social partner and the focal individual, multiplied by a regression coefficient. The regression coefficient is a genetic trait, which is a measure of cooperation; a negative value indicates competition, a positive value cooperation, and an increasing value due to selection indicates the evolution of cooperation. In contrast to the existing quantitative genetic models, our model allows for co-evolution of IGEs and variability, as the regression coefficient can respond to selection. Our simulations show that the model results in increased variability of body weight with increasing competition. When competition decreases, i.e., cooperation evolves, variability becomes significantly smaller. Hence, our model facilitates quantitative genetic studies on the relationship between IGEs and inherited variability. Moreover, our findings suggest that we may have been overlooking an entire level of genetic variation in variability, the one due to IGEs.

Selection for environmental variance of litter size in rabbits

Background

In recent years, there has been an increasing interest in the genetic determination of environmental variance. In the case of litter size, environmental variance can be related to the capacity of animals to adapt to new environmental conditions, which can improve animal welfare.

Results

We developed a ten-generation divergent selection experiment on environmental variance. We selected one line of rabbits for litter size homogeneity and one line for litter size heterogeneity by measuring intra-doe phenotypic variance. We proved that environmental variance of litter size is genetically determined and can be modified by selection. Response to selection was 4.5% of the original environmental variance per generation. Litter size was consistently higher in the Low line than in the High line during the entire experiment.

Conclusions

We conclude that environmental variance of litter size is genetically determined based on the results of our divergent selection experiment. This has implications for animal welfare, since animals that cope better with their environment have better welfare than more sensitive animals. We also conclude that selection for reduced environmental variance of litter size does not depress litter size.

Genetic control of residual variance of yearling weight in Nellore beef cattle

There is evidence for genetic variability in residual variance of livestock traits, which offers the potential for selection for increased uniformity of production. Different statistical approaches have been employed to study this topic; however, little is known about the concordance between them. The aim of our study was to investigate the genetic heterogeneity of residual variance on yearling weight (YW; 291.15 ± 46.67) in a Nellore beef cattle population; to compare the results of the statistical approaches, the two-step approach and the double hierarchical generalized linear model (DHGLM); and to evaluate the effectiveness of power transformation to accommodate scale differences. The comparison was based on genetic parameters, accuracy of EBV for residual variance, and cross-validation to assess predictive performance of both approaches. A total of 194,628 yearling weight records from 625 sires were used in the analysis. The results supported the hypothesis of genetic heterogeneity of residual variance on YW in Nellore beef cattle and the opportunity of selection, measured through the genetic coefficient of variation of residual variance (0.10 to 0.12 for the two-step approach and 0.17 for DHGLM, using an untransformed data set). However, low estimates of genetic variance associated with positive genetic correlations between mean and residual variance (about 0.20 for two-step and 0.76 for DHGLM for an untransformed data set) limit the genetic response to selection for uniformity of production while simultaneously increasing YW itself. Moreover, large sire families are needed to obtain accurate estimates of genetic merit for residual variance, as indicated by the low heritability estimates (<0.007). Box-Cox transformation was able to decrease the dependence of the variance on the mean and decreased the estimates of genetic parameters for residual variance. The transformation reduced but did not eliminate all the genetic heterogeneity of residual variance, highlighting its presence beyond the scale effect. The DHGLM showed higher predictive ability of EBV for residual variance and therefore should be preferred over the two-step approach.

Estimation of breeding values for uniformity of growth in Atlantic salmon (Salmo salar) using pedigree relationships or single-step genomic evaluation

Background

In farmed Atlantic salmon, heritability for uniformity of body weight is low, indicating that the accuracy of estimated breeding values (EBV) may be low. The use of genomic information could be one way to increase accuracy and, hence, obtain greater response to selection. Genomic information can be merged with pedigree information to construct a combined relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H matrix) for a single-step genomic evaluation (ssGBLUP), allowing realized relationships of the genotyped animals to be exploited, in addition to numerator pedigree relationships (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A matrix). We compared the predictive ability of EBV for uniformity of body weight in Atlantic salmon, when implementing either the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H matrix in the genetic evaluation. We used double hierarchical generalized linear models (DHGLM) based either on a sire-dam (sire-dam DHGLM) or an animal model (animal DHGLM) for both body weight and its uniformity.

Results

With the animal DHGLM, the use of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H instead of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A significantly increased the correlation between the predicted EBV and adjusted phenotypes, which is a measure of predictive ability, for both body weight and its uniformity (41.1 to 78.1%). When log-transformed body weights were used to account for a scale effect, the use of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H instead of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A produced a small and non-significant increase (1.3 to 13.9%) in predictive ability. The sire-dam DHGLM had lower predictive ability for uniformity compared to the animal DHGLM.

Conclusions

Use of the combined numerator and genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{H}}$$\end{document}H) significantly increased the predictive ability of EBV for uniformity when using the animal DHGLM for untransformed body weight. The increase was only minor when using log-transformed body weights, which may be due to the lower heritability of scaled uniformity, the lower genetic correlation of transformed body weight with its uniformity compared to the untransformed traits, and the small number of genotyped animals in the reference population. This study shows that ssGBLUP increases the accuracy of EBV for uniformity of body weight and is expected to increase response to selection in uniformity.

Electronic supplementary material

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

Derivation of economic values for production traits in aquaculture species

BACKGROUND

In breeding programs for aquaculture species, breeding goal traits are often weighted based on the desired gains but economic gain would be higher if economic values were used instead. The objectives of this study were: (1) to develop a bio-economic model to derive economic values for aquaculture species, (2) to apply the model to determine the economic importance and economic values of traits in a case-study on gilthead seabream, and (3) to validate the model by comparison with a profit equation for a simplified production system.

METHODS

A bio-economic model was developed to simulate a grow-out farm for gilthead seabream, and then used to simulate gross margin at the current levels of the traits and after one genetic standard deviation change in each trait with the other traits remaining unchanged. Economic values were derived for the traits included in the breeding goal: thermal growth coefficient (TGC), thermal feed intake coefficient (TFC), mortality rate (M), and standard deviation of harvest weight ([Formula: see text]). For a simplified production system, improvement in TGC was assumed to affect harvest weight instead of growing period. Using the bio-economic model and a profit equation, economic values were derived for harvest weight, cumulative feed intake at harvest, and overall survival.

RESULTS

Changes in gross margin showed that the order of economic importance of the traits was: TGC, TFC, M, and [Formula: see text]. Economic values in € (kg production)^-1 (trait unit)^-1 were: 0.40 for TGC, -0.45 for TFC, -7.7 for M, and -0.0011 to -0.0010 for [Formula: see text]. For the simplified production system, similar economic values were obtained with the bio-economic model and the profit equation. The advantage of the profit equation is its simplicity, while that of the bio-economic model is that it can be applied to any aquaculture species, because it can include any limiting factor and/or environmental condition that affects production.

CONCLUSIONS

We confirmed the validity of the bio-economic model. TGC is the most important trait to improve, followed by TFC and M, and the effect of [Formula: see text] on gross margin is small.

Genetic parameters for uniformity of harvest weight and body size traits in the GIFT strain of Nile tilapia

Background

Animal breeding programs have been very successful in improving the mean levels of traits through selection. However, in recent decades, reducing the variability of trait levels between individuals has become a highly desirable objective. Reaching this objective through genetic selection requires that there is genetic variation in the variability of trait levels, a phenomenon known as genetic heterogeneity of environmental (residual) variance. The aim of our study was to investigate the potential for genetic improvement of uniformity of harvest weight and body size traits (length, depth, and width) in the genetically improved farmed tilapia (GIFT) strain. In order to quantify the genetic variation in uniformity of traits and estimate the genetic correlations between level and variance of the traits, double hierarchical generalized linear models were applied to individual trait values.

Results

Our results showed substantial genetic variation in uniformity of all analyzed traits, with genetic coefficients of variation for residual variance ranging from 39 to 58 %. Genetic correlation between trait level and variance was strongly positive for harvest weight (0.60 ± 0.09), moderate and positive for body depth (0.37 ± 0.13), but not significantly different from 0 for body length and width.

Conclusions

Our results on the genetic variation in uniformity of harvest weight and body size traits show good prospects for the genetic improvement of uniformity in the GIFT strain. A high and positive genetic correlation was estimated between level and variance of harvest weight, which suggests that selection for heavier fish will also result in more variation in harvest weight. Simultaneous improvement of harvest weight and its uniformity will thus require index selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0218-9) contains supplementary material, which is available to authorized users.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population

Background

In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.

Results

In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.

Conclusions

To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

A Hierarchical Generalized Linear Model in Combination with Dispersion Modeling to Improve Sib-Pair Linkage Analysis

OBJECTIVE

We explored a hierarchical generalized linear model (HGLM) in combination with dispersion modeling to improve the sib-pair linkage analysis based on the revised Haseman-Elston regression model for a quantitative trait.

METHODS

A dispersion modeling technique was investigated for sib-pair linkage analysis using simulation studies and real data applications. We considered 4 heterogeneous dispersion settings according to a signal-to-noise ratio (SNR) in the various statistical models based on the Haseman-Elston regression model.

RESULTS

Our numerical studies demonstrated that susceptibility loci could be detected well by modeling the dispersion parameter appropriately. In particular, the HGLM had better performance than the linear regression model and the ordinary linear mixed model when the SNR is low, i.e., when substantial noise was present in the data.

CONCLUSION

The study shows that the HGLM in combination with dispersion modeling can be utilized to identify multiple markers showing linkage to familial complex traits accurately. Appropriate dispersion modeling might be more powerful to identify markers closest to the major genes which determine a quantitative trait.

Genetic (co)variance of rainbow trout (Oncorhynchus mykiss) body weight and its uniformity across production environments

Background

When rainbow trout from a single breeding program are introduced into various production environments, genotype-by-environment (GxE) interaction may occur. Although growth and its uniformity are two of the most important traits for trout producers worldwide, GxE interaction on uniformity of growth has not been studied. Our objectives were to quantify the genetic variance in body weight (BW) and its uniformity and the genetic correlation (r_g) between these traits, and to investigate the degree of GxE interaction on uniformity of BW in breeding (BE) and production (PE) environments using double hierarchical generalized linear models. Log-transformed data were also used to investigate whether the genetic variance in uniformity of BW, GxE interaction on uniformity of BW, and r_g between BW and its uniformity were influenced by a scale effect.

Results

Although heritability estimates for uniformity of BW were low and of similar magnitude in BE (0.014) and PE (0.012), the corresponding coefficients of genetic variation reached 19 and 21%, which indicated a high potential for response to selection. The genetic re-ranking for uniformity of BW (r_g = 0.56) between BE and PE was moderate but greater after log-transformation, as expressed by the low r_g (-0.08) between uniformity in BE and PE, which indicated independent genetic rankings for uniformity in the two environments when the scale effect was accounted for. The r_g between BW and its uniformity were 0.30 for BE and 0.79 for PE but with log-transformed BW, these values switched to -0.83 and -0.62, respectively.

Conclusions

Genetic variance exists for uniformity of BW in both environments but its low heritability implies that a large number of relatives are needed to reach even moderate accuracy of selection. GxE interaction on uniformity is present for both environments and sib-testing in PE is recommended when the aim is to improve uniformity across environments. Positive and negative r_g between BW and its uniformity estimated with original and log-transformed BW data, respectively, indicate that increased BW is genetically associated with increased variance in BW but with a decrease in the coefficient of variation. Thus, the scale effect substantially influences the genetic parameters of uniformity, especially the sign and magnitude of its r_g.