Multiple-trait Model Research Articles

Evaluating genotype by environment interaction for growth traits in Limousine cattle

Environmental conditions affect the growth and health of animals, making it crucial to quantify heat stress and the genetic basis of heat tolerance in animal breeding. The main objective of this study was to evaluate heat stress on beef cattle growth and investigate the genetic background of tolerance to harsh environmental conditions in the Italian Limousine beef cattle. Three growth traits were analysed: average daily gain (ADG), weaning weight (WW), and yearling weight (YW). Data were collected from animals raised between 1991 and 2022 and combined with 14 environmental covariates. Records for ADG, WW, and YW encompassed 108 205, 100 058, and 24 939 individuals, respectively, with 4 617, 4 670, and 2 048 genotyped individuals. Climatic variables were compared for inclusion in a linear mixed model using the Deviance Information Criterion. Multiple-trait models and genomic information incorporated environmental conditions with the largest impact on the studied traits Genotype by environment interaction (G × E) was detected in all the studied traits, showing substantial heterogeneity of the variance components across the different environments (Env). Heritability for WW remains constant among Env; instead, for ADG and YW decreased under uncomfortable environmental conditions. YW showed the lowest genetic correlation (0.28) between divergent conditions (Env 2 and Env 5,) for ADG and WW correlations dropped below 0.50 among Env. The values of genetic correlations indicate that growth traits are moderately to strongly affected by G × E. Eigenvalue decomposition of the additive genetic (co)variance matrix for ADG, WW, and YW indicated that three components accounted for over 0.80 of the proportion of the variance explained, suggesting different animal performances across Env. Spearman rank correlations showed potential re-ranking of genotyped sires, because ADG, WW, and YW showed correlations between Env below 0.80. The accuracy of single-step genomic EBV was higher compared to EBV for al traits. Overall, the result confirms the existence of G × E for growth traits in the Italian Limousine population. Including G × E in the model allows for more environment-aware predictions, helping breeders understand how different genetic bases respond to varying conditions. Genomic predictions incorporating G × E could accelerate genetic gains and improve response to selection for heat tolerance.

Supervised Machine Learning Techniques for Breeding Value Prediction in Horses: An Example Using Gait Visual Scores.

Gait scores are widely used in the genetic evaluation of horses. However, the nature of such measurement may limit genetic progress since there is subjectivity in phenotypic information. This study aimed to assess the application of machine learning techniques in the prediction of breeding values for five visual gait scores in Campolina horses: dissociation, comfort, style, regularity, and development. The dataset contained over 5000 phenotypic records with 107,951 horses (14 generations) in the pedigree. A fixed model was used to estimate least-square solutions for fixed effects and adjusted phenotypes. Variance components and breeding values (EBV) were obtained via a multiple-trait model (MTM). Adjusted phenotypes and fixed effects solutions were used to train machine learning models (using the EBV from MTM as target variable): artificial neural network (ANN), random forest regression (RFR) and support vector regression (SVR). To validate the models, the linear regression method was used. Accuracy was comparable across all models (but it was slightly higher for ANN). The highest bias was observed for ANN, followed by MTM. Dispersion varied according to the trait; it was higher for ANN and the lowest for MTM. Machine learning is a feasible alternative to EBV prediction; however, this method will be slightly biased and over-dispersed for young animals.

118 DHGLMF90: A software tool using double hierarchical generalized linear models for estimating the genetic heterogeneity of residual variance

Abstract In animal breeding, there is a growing interest in selecting animals that exhibit uniform responses to environmental conditions, aiming for product consistency as a desirable breeding objective. This could be achieved by modeling heteroskedastic residuals when estimating breeding values through a hierarchical single-trait model that includes a mean and dispersion part. The latter could be influenced by genetic and non-genetic features. Bayesian methods can estimate heteroskedastic residuals using the Metropolis-Hastings algorithm, which is a slow and inefficient approach. Double hierarchical generalized linear models (DHGLM) provide a faster alternative to Bayesian methods. This study aimed to develop software named DHGLMF90 to implement DHGLM using a highly efficient algorithm to accurately estimate heterogeneous residual variances and breeding values, particularly suited for handling large datasets. We improved the reweighted least squares algorithm (IRWLS) to enhance convergence properties compared with prior implementations. IRWLS iteratively computes variance components and leverages them for a bivariate model, covering both the mean and dispersion part of the original model. In DHGLMF90, variance component estimation is performed using REML through BLUPF90+. The software was tested using both simulated and real datasets, the latter from previous studies. Additionally, the simulations were replicated to ensure a comprehensive evaluation of the DHGLMF90 performance. Results indicated that, in most cases, there were no significant differences in estimation compared with results obtained from previous studies. DHGLMF90 converged in order of magnitude faster than Bayesian methods. Future developments will include analyzing datasets with genomic information and multiple-trait models.

Multiple-trait genomic prediction for swine meat quality traits using gut microbiome features as a correlated trait.

Traits such as meat quality and composition are becoming valuable in modern pork production; however, they are difficult to include in genetic evaluations because of the high phenotyping costs. Combining genomic information with multiple-trait indirect selection with cheaper indicator traits is an alternative for continued cost-effective genetic improvement. Additionally, gut microbiome information is becoming more affordable to measure using targeted rRNA sequencing, and its applications in animal breeding are becoming relevant. In this paper, we investigated the usefulness of microbial information as a correlated trait in selecting meat quality in swine. This study incorporated phenotypic data encompassing marbling, colour, tenderness, loin muscle and backfat depth, along with the characterization of gut (rectal) microbiota through 16S rRNA sequencing at three distinct time points of the animal's growth curve. Genetic progress estimation and cross-validation were employed to evaluate the utility of utilizing host genomic and gut microbiota information for selecting expensive-to-record traits in crossbred individuals. Initial steps involved variance components estimation using multiple-trait models on a training dataset, where the top 25 associated operational taxonomic units (OTU) for each meat quality trait and time point were included. The second step compared the predictive ability of multiple-trait models incorporating different numbers of OTU with single-trait models in a validation set. Results demonstrated the advantage of including genomic information for some traits, while in some instances, gut microbial information proved advantageous, namely, for marbling and pH. The study suggests further investigation into the shared genetic architecture between microbial features and traits, considering microbial data's compositional and high-dimensional nature. This research proposes a straightforward method to enhance swine breeding programs for improving costly-to-record traits like meat quality by incorporating gut microbiome information.

Equivalence of variance components between standard and recursive genetic models using LDL′ transformations

BackgroundRecursive models are a category of structural equation models that propose a causal relationship between traits. These models are more parameterized than multiple trait models, and they require imposing restrictions on the parameter space to ensure statistical identification. Nevertheless, in certain situations, the likelihood of recursive models and multiple trait models are equivalent. Consequently, the estimates of variance components derived from the multiple trait mixed model can be converted into estimates under several recursive models through LDL′ or block-LDL′ transformations.ResultsThe procedure was employed on a dataset comprising five traits (birth weight—BW, weight at 90 days—W90, weight at 210 days—W210, cold carcass weight—CCW and conformation—CON) from the Pirenaica beef cattle breed. These phenotypic records were unequally distributed among 149,029 individuals and had a high percentage of missing data. The pedigree used consisted of 343,753 individuals. A Bayesian approach involving a multiple-trait mixed model was applied using a Gibbs sampler. The variance components obtained at each iteration of the Gibbs sampler were subsequently used to estimate the variance components within three distinct recursive models.ConclusionsThe LDL′ or block-LDL′ transformations applied to the variance component estimates achieved from a multiple trait mixed model enabled inference across multiple sets of recursive models, with the sole prerequisite of being likelihood equivalent. Furthermore, the aforementioned transformations simplify the handling of missing data when conducting inference within the realm of recursive models.

Comparison between multiple-trait and random regression models for genetic evaluation of weight traits in Australian meat sheep.

Random regression (RR) models are recommended as an alternative to multiple-trait (MT) models for better capturing the variance-covariance structure over a trajectory and hence more accurate genetic evaluation of traits that are repeatedly measured and genetically change gradually over time. However, a limited number of studies have been done to empirically compare RR over a MT model to determine how much extra benefit could be achieved from one method over another. We compared the prediction accuracy of RR and MT models for growth traits of Australian meat sheep measured from 60 to 525d, using 102,579 weight records from 24,872 animals. Variance components and estimated breeding values (EBVs) estimated at specific ages were compared and validated with forward prediction. The accuracy of EBVs obtained from the MT model was 0.58, 0.51, 0.54, and 0.56 for weaning, postweaning, yearling, and hogget weight stages, respectively. RR model produced accuracy estimates of 0.56, 0.51, 0.54, and 0.54 for equivalent weight stages. Regression of adjusted phenotype on EBVs was very similar between the MT and the RR models (P > 0.05). Although the RR model did not significantly increase the accuracy of predicting future progeny performance, there are other benefits of the model such as no limit to the number of records per animal, estimation of EBVs for early and late growth, no need for age correction. Therefore, RR can be considered a more flexible method for the genetic evaluation of Australian sheep for early and late growth, and no need for age correction.

Structural equation models to infer relationships between energy-related blood metabolites and milk daily energy output in Holstein cows.

During lactation, high-yielding cows experience metabolic disturbances due to milk production. Metabolic monitoring offers valuable insights into how cows manage these challenges throughout the lactation period, making it a topic of considerable interest to breeders. In this study, we used Bayesian networks to uncover potential dependencies among various energy-related blood metabolites, i.e., glucose, urea, beta-hydroxybutyrate (BHB), non-esterified fatty acids (NEFA), cholesterol (CHOL), and daily milk energy output (dMEO) in 1,254 Holstein cows. The inferred causal structure was then incorporated into structural equation models (SEM) to estimate heritabilities and additive genetic correlations among these phenotypes using both pedigree and genotypes from a 100k chip. Dependencies among traits were determined using the Hill-Climbing algorithm, implemented with the posterior distribution of the residuals obtained from the standard multiple-trait model. These identified relationships were then used to construct the SEM, considering both direct and indirect relationships. The relevant dependencies and path coefficients obtained, expressed in units of measurement variation of 1σ, were as follows: dMEO → CHOL (0.181), dMEO → BHB (-0.149), dMEO → urea (0.038), glucose → BHB (-0.55), glucose → urea (-0.194), CHOL → urea (0.175), BHB → urea (-0.049), and NEFA → urea (-0.097). Heritabilities for traits of concern obtained with SEM ranged from 0.09 to 0.2. Genetic correlations with a minimum 95% probability (P) of the posterior mean being >0 for positive means or <0 for negative means include those between dMEO and glucose (-0.583, P = 100), dMEO and BHB (0.349, P = 99), glucose and CHOL (0.325, P = 100), glucose and NEFA (-0.388, P = 100), and NEFA and BHB (0.759, P = 100). The results of this analysis revealed the existence of recursive relationships among the energy-related blood metabolites and dMEO. Understanding these connections is paramount for establishing effective genetic selection strategies, enhancing production and animal welfare.

Bayesian inference applied to soybean grown under different shading levels using the multiple-trait model

Bayesian inference applied to soybean grown under different shading levels using the multiple-trait model

Genetic relationship of energy balance predicted from milk traits with fertility in Japanese Holsteins.

We predicted the energy balance of cows from milk traits and estimated the genetic correlations of predicted energy balance (PEB) with fertility traits for the first three lactations. Data included 9,646,606 test-day records of 576,555 Holstein cows in Japan from 2015 to 2019. Genetic parameters were estimated with a multiple-trait model in which the records among lactation stages and parities were treated as separate traits. Fertility traits were conception rate at first insemination (CR), number of inseminations (NI), and days open (DO). Heritability estimates of PEB were 0.28-0.35 (first lactation), 0.15-0.29 (second), and 0.09-0.23 (third). Estimated genetic correlations among lactation stages were 0.85-1.00 (first lactation), 0.73-1.00 (second), and 0.64-1.00 (third). Estimated genetic correlations among parities were 0.82-0.96 (between first and second), 0.97-0.99 (second and third), and 0.69-0.92 (first and third). Estimated genetic correlations of PEB in early lactation with fertility were 0.04 to 0.19 for CR, -0.03 to -0.19 for NI, and -0.01 to -0.24 for DO. Genetic improvement of PEB is possible. Lower PEB in early lactation was associated with worse fertility, suggesting that improving PEB in early lactation may improve reproductive performance.

Genetic analysis of lactation consistency in U.S. Holsteins using temporal variation in daily milk weights

The ability of a dairy cow to perform reliably over time is an interesting trait to include in dairy cattle breeding programs aimed at improving dairy cow resilience. Consistency, defined as the quality of performing as expected each day of the lactation, could be highly associated with resilience, defined as animal's ability to maintain health and performance in the presence of environmental challenges, including pathogens, heat waves, and nutritional changes. A total of 51,415,022 daily milk weights collected from 2018 to 2023 were provided for 255,191 multiparous Holstein cows milked 3 times daily in conventional parlor systems on farms in 32 states. The temporal variance (TempVar) of milk yield from 5 to 305 d postpartum was computed as the log-transformed variance of daily deviations between observed and expected individual milk weights. Lower values of TempVar imply smaller day-to-day deviations from expectations, indicating consistent performance, whereas larger values indicate inconsistent performance. Expected daily milk weights were computed using 3 nonparametric and parametric regression models: 1) LOESS regression with a 0.75 span; 2) polynomial quantile regression using the median (0.5 quantile), and 3) polynomial quantile regression using a 0.7 quantile. The univariate statistical model included age at first calving and herd-year-season as fixed effects and cow as a random effect. Heritability estimates (standard errors) of TempVar phenotypes calculated over the entire lactation ranged between 0.227 (0.011) and 0.237 (0.011), demonstrating that cows are genetically predisposed to display consistent or inconsistent performance. Estimated genetic correlations calculated using a multiple trait model between TempVar traits and between lactations were high (>0.95), indicating TempVar is repeatable across lactations and robust to the model used to compute expected daily milk yield. Higher TempVar phenotypes reflect more variation in performance, hence greater inconsistency, which is undesirable. Therefore, correlations between predicted transmitting abilities (PTAs) for TempVar and milk yield of 0.57 indicate that high-producing cows exhibit more day-to-day variation in performance. Correlations with productive life and livability were -0.38 and -0.48, respectively. Correlations between PTAs for TempVar and those of postpartum health traits were also negative, ranging from -0.41 to -0.08. Given that health traits are derived from disease resistance measurements, and higher health trait PTAs are preferred, our results indicate that more consistent cows tend to have fewer health problems and greater longevity. Overall, our findings suggest that temporal variation in daily milk weights can be used to identify consistent animals that maintain expected performance throughout the lactation, which will enable selection for greater resilience to management and environmental perturbations.

Across-country genetic and genomic analyses of foot score traits in American and Australian Angus cattle

BackgroundHoof structure and health are essential for the welfare and productivity of beef cattle. Therefore, we assessed the genetic and genomic background of foot score traits in American (US) and Australian (AU) Angus cattle and investigated the feasibility of performing genomic evaluations combining data for foot score traits recorded in US and AU Angus cattle. The traits evaluated were foot angle (FA) and claw set (CS). In total, 109,294 and ~ 1.12 million animals had phenotypic and genomic information, respectively. Four sets of analyses were performed: (1) genomic connectedness between US and AU Angus cattle populations and population structure, (2) estimation of genetic parameters, (3) single-step genomic prediction of breeding values, and (4) single-step genome-wide association studies for FA and CS.ResultsThere was no clear genetic differentiation between US and AU Angus populations. Similar heritability estimates (FA: 0.22–0.24 and CS: 0.22–0.27) and moderate-to-high genetic correlations between US and AU foot scores (FA: 0.61 and CS: 0.76) were obtained. A joint-genomic prediction using data from both populations outperformed within-country genomic evaluations. A genomic prediction model considering US and AU datasets as a single population performed similarly to the scenario accounting for genotype-by-environment interactions (i.e., multiple-trait model considering US and AU records as different traits), even though the genetic correlations between countries were lower than 0.80. Common significant genomic regions were observed between US and AU for FA and CS. Significant single nucleotide polymorphisms were identified on the Bos taurus (BTA) chromosomes BTA1, BTA5, BTA11, BTA13, BTA19, BTA20, and BTA23. The candidate genes identified were primarily from growth factor gene families, including FGF12 and GDF5, which were previously associated with bone structure and repair.ConclusionsThis study presents comprehensive population structure and genetic and genomic analyses of foot scores in US and AU Angus cattle populations, which are essential for optimizing the implementation of genomic selection for improved foot scores in Angus cattle breeding programs. We have also identified candidate genes associated with foot scores in the largest Angus cattle populations in the world and made recommendations for genomic evaluations for improved foot score traits in the US and AU.

Genomic selection of soybean (Glycine max) for genetic improvement of yield and seed composition in a breeding context.

Genomic selection has been utilized for genetic improvement in both plant and animal breeding and is a favorable technique for quantitative trait development. Within this study, genomic selection was evaluated within a breeding program, using novel validation methods in addition to plant materials and data from a commercial soybean (Glycine max) breeding program. A total of 1501 inbred lines were used to test multiple genomic selection models for multiple traits. Validation included cross-validation, inter-environment, and empirical validation. The results indicated that the extended genomic best linear unbiased prediction (EGBLUP) model was the most effective model tested for yield, protein, and oil in cross-validation with accuracies of 0.50, 0.68, and 0.64, respectively. Increasing marker number from 1000 to 3000 to 6000 single nucleotide polymorphism markers leads to statistically significant increases in accuracy. Cross-environment predictions were statistically lower than cross-validation with accuracies of 0.24, 0.54, and 0.42 for yield, protein, and oil, respectively, using the extended genomic BLUP model. Empirical validation, predicting the yield of 510 soybean lines, had a prediction accuracy of 0.34, with the inclusion of a maturity covariate leading to a notable increase in accuracy. Genomic selection identified high-performance lines in inter-environment predictions: 34% of lines within the upper quartile of yield, and 51% and 48% of the highest quartile protein and oil lines, respectively. Statistically similar results occurred comparing rankings in empirical validation and selection for advancements in yield trials. These results indicate that genomic selection is a useful tool for selection decisions.

Multiple-trait model by Bayesian inference applied to environment efficient Coffea arabica with low-nitrogen nutrient

Identifying Coffea arabica cultivars that are more efficient in the use of nitrogen is an important strategy and a necessity in the context of environmental and economic impacts attributed to excessive nitrogen fertilization. Although Coffea arabica breeding data have a multi-trait structure, they are often analyzed under a single trait structure. Thus, the objectives of this study were to use a Bayesian multitrait model, to estimate heritability in the broad sense, and to select arabica coffee cultivars with better genetic potential (desirable agronomic traits) in nitrogen-restricted cultivation. The experiment was carried out in a greenhouse with 20 arabica coffee cultivars grown in a nutrient solution with low-nitrogen content (1.5 mM). The experimental design used was in randomized blocks with three replications. Six agromorphological traits of the arabica coffee breeding program and five nutritional efficiency indices were used. The Markov Chain Monte Carlo algorithm was used to estimate genetic parameters and genetic values. The agromorphological traits were considered highly heritable, with a credibility interval (95% probability): H2 = 0.9538 – 5.89E-01. The Bayesian multitrait model presents an adequate strategy for the genetic improvement of arabica coffee grown in low-nitrogen concentrations. Coffee arabica cultivars Icatu Precoce 3282, Icatu Vermelho IAC 4045, Acaiá Cerrado MG 1474, Tupi IAC 1669-33, Catucaí 785/15, Caturra Vermelho and Obatã IAC 1669/20 demonstrated greater potential for cultivation in low-nitrogen concentration.

Genetic correlations and causal effects of fighting ability on fitness traits in cattle reveal antagonistic trade-offs

Complex genetic and phenotypic relationships are theorized to link different fitness components but revealing the correlations occurring among disparate traits requires large datasets of pedigreed populations. In particular, the association between traits beneficial to social dominance with health and fitness could be antagonistic, because of trade-offs, or positive, because of greater resource acquisition by dominant individuals. Studies investigating these relationships found some empirical evidence in support of both theories, mainly using multiple trait models (MTM). However, if a trait giving a social advantage is suspected to affect the expression of other traits, MTM could provide some bias, that structural equation models (SEM) could highlight. We used Aosta Chestnut-Black Pied cattle to investigate whether the fighting ability of cows (the capability of winning social dominance interactions) is genetically correlated with health and fitness traits. We ran both MTM and SEM using a Gibbs sampling algorithm to disentangle the possible causal effects of fighting ability from the genetic correlations that this trait shares with other traits: individual milk yield, somatic cells (representing mammary health), fertility, and longevity. We found antagonistic genetic correlations, similar under both approaches, for fighting ability vs. milk, somatic cells, and fertility, Accordingly, we found only a slight causal effects of fighting ability on these traits (–0.012 to 0.059 in standardized value). However, we found genetic correlations opposite in sign between fighting ability and longevity under MTM (0.237) and SEM (–0.183), suggesting a strong causal effect (0.386 standardized) of fighting ability on longevity. In other words, MTM found a positive correlation between longevity and fighting ability, while SEM found a negative correlation. The explanation could be that for economic reasons dominant cows are kept in this population for longer, thus attaining greater longevity: using MTM, the economic importance of competitions probably covers the true genetic correlation among traits. This artificially simulates a natural situation where an antagonistic genetic correlation between longevity and fighting ability appears positive under MTM due to a non-genetic advantage obtained by the best fighters. The use of SEM to properly assess the relationships among traits is suggested in both evolutionary studies and animal breeding.

Marker density and statistical model designs to increase accuracy of genomic selection for wool traits in Angora rabbits.

The Angora rabbit, a well-known breed for fiber production, has been undergoing traditional breeding programs relying mainly on phenotypes. Genomic selection (GS) uses genomic information and promises to accelerate genetic gain. Practically, to implement GS in Angora rabbit breeding, it is necessary to evaluate different marker densities and GS models to develop suitable strategies for an optimized breeding pipeline. Considering a lack in microarray, low-coverage sequencing combined with genotype imputation was used to boost the number of SNPs across the rabbit genome. Here, in a population of 629 Angora rabbits, a total of 18,577,154 high-quality SNPs were imputed (imputation accuracy above 98%) based on low-coverage sequencing of 3.84X genomic coverage, and wool traits and body weight were measured at 70, 140 and 210 days of age. From the original markers, 0.5K, 1K, 3K, 5K, 10K, 50K, 100K, 500K, 1M and 2M were randomly selected and evaluated, resulting in 50K markers as the baseline for the heritability estimation and genomic prediction. Comparing to the GS performance of single-trait models, the prediction accuracy of nearly all traits could be improved by multi-trait models, which might because multiple-trait models used information from genetically correlated traits. Furthermore, we observed high significant negative correlation between the increased prediction accuracy from single-trait to multiple-trait models and estimated heritability. The results indicated that low-heritability traits could borrow more information from correlated traits and hence achieve higher prediction accuracy. The research first reported heritability estimation in rabbits by using genome-wide markers, and provided 50K as an optimal marker density for further microarray design, genetic evaluation and genomic selection in Angora rabbits. We expect that the work could provide strategies for GS in early selection, and optimize breeding programs in rabbits.

Multiple-trait model through Bayesian inference applied to flood-irrigated rice (Oryza sativa L)

The objectives of this study were to use a bayesian multi-trait model, estimate genetic parameters, and select flood-irrigated rice genotypes with better genetic potentials in different evaluation environments. For this, twenty-five rice genotypes and six traits belonging to the flood-irrigated rice improvement program were evaluated. The experimental design used in all experiments was a randomized block design with three replications. The Monte Carlo Markov Chain algorithm estimated genetic parameters and genetic values. The grain thickness trait was considered highly heritable, with a credibility interval ranging from: $${h}^{2}$$ : 0.9480; 0.9440; 0.8610, in environments 1, 2, and 3, respectively. The grain yields showed a weak correlation estimate between grain thickness and 100-grain weight, in all environments, with a credibility interval ranging from ( $$\rho$$ = 0.5477; 0.5762; 0.5618 and 0.5973; 0.5247; 0.5632, grain thickness and 100-grain weight, in environments 1, 2, and 3, respectively). The Bayesian multi-trait model proved to be an adequate strategy for the genetic improvement of flood-irrigated.

Multiple-trait analyses improved the accuracy of genomic prediction and the power of genome-wide association of productivity and climate change-adaptive traits in lodgepole pine

BackgroundGenomic prediction (GP) and genome-wide association (GWA) analyses are currently being employed to accelerate breeding cycles and to identify alleles or genomic regions of complex traits in forest trees species. Here, 1490 interior lodgepole pine (Pinus contorta Dougl. ex. Loud. var. latifolia Engelm) trees from four open-pollinated progeny trials were genotyped with 25,099 SNPs, and phenotyped for 15 growth, wood quality, pest resistance, drought tolerance, and defense chemical (monoterpenes) traits. The main objectives of this study were to: (1) identify genetic markers associated with these traits and determine their genetic architecture, and to compare the marker detected by single- (ST) and multiple-trait (MT) GWA models; (2) evaluate and compare the accuracy and control of bias of the genomic predictions for these traits underlying different ST and MT parametric and non-parametric GP methods. GWA, ST and MT analyses were compared using a linear transformation of genomic breeding values from the respective genomic best linear unbiased prediction (GBLUP) model. GP, ST and MT parametric and non-parametric (Reproducing Kernel Hilbert Spaces, RKHS) models were compared in terms of prediction accuracy (PA) and control of bias.ResultsMT-GWA analyses identified more significant associations than ST. Some SNPs showed potential pleiotropic effects. Averaging across traits, PA from the studied ST-GP models did not differ significantly from each other, with generally a slight superiority of the RKHS method. MT-GP models showed significantly higher PA (and lower bias) than the ST models, being generally the PA (bias) of the RKHS approach significantly higher (lower) than the GBLUP.ConclusionsThe power of GWA and the accuracy of GP were improved when MT models were used in this lodgepole pine population. Given the number of GP and GWA models fitted and the traits assessed across four progeny trials, this work has produced the most comprehensive empirical genomic study across any lodgepole pine population to date.

Effects of an herbal diet, diet x boar line and diet x genotype interactions on skin lesions and on growth performance in post-weaning pigs using a cross-classified experiment

The aims of the present study were to assess the effect of an herbal diet (HD; basal diet plus a supplement with chicory herbs) on production traits and lesions scores in post weaning pigs, and to study possible boar x diet and genotype x diet interactions considering the two feeding groups HD and control (CON). In this regard, a cross-classified research design was implemented, aiming on an equal number of boar offspring in both groups HD and CON. Grouping of piglets (993 crossbred piglets from 14 different Piétrain boars sire to German Landrace or German Edelschwein sows) was done after weaning, and the feeding experiment was conducted in post-weaning period up to the pig age of nine weeks. Studied traits included the post-weaning weight (PWW) at the age of nine week 9, average daily gain (ADG) considering the period from weaning until week 9 and a skin lesion score (LS) at week 9. Herbal diet effects and sire x diet effects were inferred via linear models for Gaussian distributed PWW and ADG, and via generalized linear mixed models with a logit link function for categorical LS. Possible genotype x diet interactions were proved in multiple-trait models via Gibbs sampling, by defining same traits recorded either in CON and HD as different traits, i.e., post-weaning weight control (PWWCON) and post-weaning weight herbal diet (PWWHD), average daily gain control (ADGCON) and average daily gain herbal diet (ADGHD), and lesion score control (LSCON) and lesion score herbal diet (LSHD). The dietary treatment had no significant effect (P > 0.05) on PWW and ADG. The dietary treatment significantly (P < 0.05) influenced the behavior of the pigs, with a lower and favorable LS score for the pigs allocated to HD. Significant sire x diet interactions were identified for LS, PWW and ADG. Pearson correlations between estimates for sire effects from both diets CON and HD (i.e., the sire x diet effect) were 0.84 for PWW, 0.71 for ADG and 0.56 for LS. Genetic correlations smaller than 0.80 between PWWCON and PWWHD, ADGCON and ADGHD, and LSCON and LSHD indicate genotype x diet interactions, especially for the functional and low heritability LS. Hence, re-rankings of sires according to their estimated breeding values (EBV) were observed in both environments CON and LD, displaying smaller genetic and EBV variations for LS in the HD group. The HD dietary supplement contributed to fewer skin lesions, probably due to the favorable effect on pig behavior, but a more pronounced genetic differentiation seems to be possible in the more challenging CON environment. Generally, the results from the present study indicate the importance of genotype x feeding interactions for pig breeding, suggesting specific boars for different feeding environments.

Improving the accuracy of genomic predictions for disease resistance traits in fish using a multiple-trait linear-threshold model

Genomic selection is effective in enhancing the selection efficiency during breeding for improved disease resistance in fish. However, given that the collection of routine disease challenge data is generally costly and that disease resistance traits are threshold traits with typically low heritability, there has been comparatively limited progress in improving the accuracy of genomic prediction. The objective of this study was to assess the advantage of a multiple-trait linear-threshold model in genomic prediction of disease resistance traits in fish. The study was based on analyses of both simulated data and real data obtained for a rainbow trout population. We simulated a linear trait (growth trait, h2 = 0.3) and a binary threshold trait (disease resistance trait, h2 = 0.1) with different genetic correlations (0.1, 0.3 and 0.5). Single- and multiple-trait models with best linear unbiased prediction (BLUP) and genomic BLUP (GBLUP) were implemented to investigate their prediction abilities. Moreover, we also assessed the impact of missing proportions of a reference population on genomic prediction for disease resistance traits. The results revealed that methods using marker information produced more accurate predictions than the pedigree-based BLUP method. Furthermore, for between-trait genetic correlations of 0.1, 0.3 and 0.5, the multiple-trait GBLUP for disease resistance traits showed 1.0%, 1.5% and 6.4% higher accuracy than single-trait GBLUP, respectively. Meanwhile, no improvement in accuracy was observed for the assessed growth trait. Moreover, with an increase in the genetic correlation between traits, we observed reductions in the consistency of estimated SNP effects using single-trait and multiple-trait GBLUP models. Our findings indicated the superiority of multiple-trait models in genomic prediction for a disease resistance trait with missing data; this superiority is more pronounced when the proportion of the reference population missing is large. Furthermore, using real data obtained for rainbow trout, we verified the advantage of the multiple-trait linear-threshold model in genomic prediction of binary survival, with the multiple-trait GBLUP model yielding a lower mean squared (absolute) error than the single-trait GBLUP model. Overall, our findings have important implications with regard to applying multiple-trait linear-threshold model in genomic selection for disease resistance in most aquaculture species.

Removing data and using metafounders alleviates biases for all traits in Lacaune dairy sheep predictions

Bias in dairy genetic evaluations, when it exists, has to be understood and properly addressed. The origin of biases is not always clear. We analyzed 40 yr of records from the Lacaune dairy sheep breeding program to evaluate the extent of bias, assess possible corrections, and emit hypotheses on its origin. The data set included 7 traits (milk yield, fat and protein contents, somatic cell score, teat angle, udder cleft, and udder depth) with records from 600,000 to 5 million depending on the trait, ∼1,900,000 animals, and ∼5,900 genotyped elite artificial insemination rams. For the ∼8% animals with missing sire, we fit 25 unknown parent groups. We used the linear regression method to compare "partial" and "whole" predictions of young rams before and after progeny testing, with 7 cut-off points, and we obtained estimates of their bias, (over)dispersion, and accuracy in early proofs. We tried (1) several scenarios as follows: multiple or single trait, the "official" (routine) evaluation, which is a mixture of both single and multiple trait, and "deletion" of data before 1990; and (2) several models as follows: BLUP and single-step genomic (SSG)BLUP with fixed unknown parent groups or metafounders, where, for metafounders, their relationship matrix gamma was estimated using either a model for inbreeding trend, or base allele frequencies estimated by peeling. The estimate of gamma obtained by modeling the inbreeding trend resulted in an estimated increase of inbreeding, based on markers, faster than the pedigree-based one. The estimated genetic trends were similar for most models and scenarios across all traits, but were shrunken when gamma was estimated by peeling. This was due to shrinking of the estimates of metafounders in the latter case. Across scenarios, all traits showed bias, generally as an overestimate of genetic trend for milk yield and an underestimate for the other traits. As for the slope, it showed overdispersion of estimated breeding values for all traits. Using multiple-trait models slightly reduced the overestimate of genetic trend and the overdispersion, as did including genomic information (i.e., SSGBLUP) when the gamma matrix was estimated by the model for inbreeding trend. However, only deletion of historical data before 1990 resulted in elimination of both kind of biases. The SSGBLUP resulted in more accurate early proofs than BLUP for all traits. We considered that a snowball effect of small errors in each genetic evaluation, combined with selection, may have resulted in biased evaluations. Improving statistical methods reduced some bias but not all, and a simple solution for this data set was to remove historical records.