Search

Eucalyptus globulus Labill., is a recognized multipurpose tree, which stands out not only for the valuable qualities of its wood but also for the medicinal applications of the essential oil extracted from its leaves. In this study, we implemented an integrated strategy comprising genomic and phenomic approaches to predict foliar essential oil content, stem quality, and growthrelated traits within a 9-year-old breeding population of E. globulus. The strategy involved evaluating Uni/Multi-trait deep learning (DL) models by incorporating genomic data related to single nucleotide polymorphisms (SNPs) and haplotypes, as well as the phenomic data from leaf near-infrared (NIR) spectroscopy. Our results showed that essential oil content (oil yield) ranged from 0.01 to 1.69% v/fw and had no significant correlation with any growth-related traits. This suggests that selection solely based on growth-related traits did n The emphases (colored text) from revisions were removed throughout the article. Confirm that this change is fine. ot influence the essential oil content. Genomic heritability estimates ranged from 0.25 (diameter at breast height (DBH) and oil yield) to 0.71 (DBH and stem straightness (ST)), while pedigree-based heritability exhibited a broader range, from 0.05 to 0.88. Notably, oil yield was found to be moderate to highly heritable, with genomic values ranging from 0.25 to 0.60, alongside a pedigree-based estimate of 0.48. The DL prediction models consistently achieved higher prediction accuracy (PA) values with a Multi-trait approach for most traits analyzed, including oil yield (0.699), tree height (0.772), DBH (0.745), slenderness coefficient (0.616), stem volume (0.757), and ST (0.764). The Uni-trait approach achieved superior PA values solely for branching quality (0.861). NIR spectral absorbance was the best omics data for CNN or MLP models with a Multi-trait approach. These results highlight considerable genetic variation within the Eucalyptus progeny trial, particularly regarding oil production. Our results contribute significantly to understanding omics-assisted deep learning models as a breeding strategy to improve growth-related traits and optimize essential oil production in this species. [ABSTRACT FROM AUTHOR]

Purpose: More accurate prediction of phenotype traits can increase the success of genomic selection in both plant and animal breeding studies and provide more reliable disease risk prediction in humans. Traditional approaches typically use regression models based on linear assumptions between the genetic markers and the traits of interest. Non-linear models have been considered as an alternative tool for modeling genomic interactions (i.e. non-additive effects) and other subtle non-linear patterns between markers and phenotype. Deep learning has become a state-of-the-art non-linear prediction method for sound, image and language data. However, genomic data is better represented in a tabular format. The existing literature on deep learning for tabular data proposes a wide range of novel architectures and reports successful results on various datasets. Tabular deep learning applications in genome-wide prediction (GWP) are still rare. In this work, we perform an overview of the main families of recent deep learning architectures for tabular data and apply them to multi-trait regression and multi-class classification for GWP on real gene datasets. Methods: The study involves an extensive overview of recent deep learning architectures for tabular data learning: NODE, TabNet, TabR, TabTransformer, FT-Transformer, AutoInt, GANDALF, SAINT and LassoNet. These architectures are applied to multi-trait GWP. Comprehensive benchmarks of various tabular deep learning methods are conducted to identify best practices and determine their effectiveness compared to traditional methods. Results: Extensive experimental results on several genomic datasets (three for multi-trait regression and two for multi-class classification) highlight LassoNet as a standout performer, surpassing both other tabular deep learning models and the highly efficient tree based LightGBM method in terms of both best prediction accuracy and computing efficiency. Conclusion: Through series of evaluations on real-world genomic datasets, the study identifies LassoNet as a standout performer, surpassing decision tree methods like LightGBM and other tabular deep learning architectures in terms of both predictive accuracy and computing efficiency. Moreover, the inherent variable selection property of LassoNet provides a systematic way to find important genetic markers that contribute to phenotype expression. [ABSTRACT FROM AUTHOR]

Abstract Purpose More accurate prediction of phenotype traits can increase the success of genomic selection in both plant and animal breeding studies and provide more reliable disease risk prediction in humans. Traditional approaches typically use regression models based on linear assumptions between the genetic markers and the traits of interest. Non-linear models have been considered as an alternative tool for modeling genomic interactions (i.e. non-additive effects) and other subtle non-linear patterns between markers and phenotype. Deep learning has become a state-of-the-art non-linear prediction method for sound, image and language data. However, genomic data is better represented in a tabular format. The existing literature on deep learning for tabular data proposes a wide range of novel architectures and reports successful results on various datasets. Tabular deep learning applications in genome-wide prediction (GWP) are still rare. In this work, we perform an overview of the main families of recent deep learning architectures for tabular data and apply them to multi-trait regression and multi-class classification for GWP on real gene datasets. Methods The study involves an extensive overview of recent deep learning architectures for tabular data learning: NODE, TabNet, TabR, TabTransformer, FT-Transformer, AutoInt, GANDALF, SAINT and LassoNet. These architectures are applied to multi-trait GWP. Comprehensive benchmarks of various tabular deep learning methods are conducted to identify best practices and determine their effectiveness compared to traditional methods. Results Extensive experimental results on several genomic datasets (three for multi-trait regression and two for multi-class classification) highlight LassoNet as a standout performer, surpassing both other tabular deep learning models and the highly efficient tree based LightGBM method in terms of both best prediction accuracy and computing efficiency. Conclusion Through series of evaluations on real-world genomic datasets, the study identifies LassoNet as a standout performer, surpassing decision tree methods like LightGBM and other tabular deep learning architectures in terms of both predictive accuracy and computing efficiency. Moreover, the inherent variable selection property of LassoNet provides a systematic way to find important genetic markers that contribute to phenotype expression.

Abstract Background Globally, the cultivation of cotton is constrained by its tendency for extended periods of growth. Early maturity plays a potential role in rainfed-based multiple cropping system especially in the current era of climate change. In the current study, a set of 20 diverse Gossypium hirsutum genotypes were evaluated in two crop seasons with three planting densities and assessed for 11 morphological traits related to early maturity. The study aimed to identify genotype(s) that mature rapidly and accomplish well under diverse environmental conditions based on the two robust multivariate techniques called multi-trait stability index (MTSI) and multi-trait genotype-ideotype distance index (MGIDI). Results MTSI analysis revealed that out of the 20 genotypes, three genotypes, viz., NNDC-30, A-2, and S-32 accomplished well in terms of early maturity traits in two seasons. Furthermore, three genotypes were selected using MGIDI method for each planting densities with a selection intensity of 15%. The strengths and weaknesses of the genotypes selected based on MGIDI method highlighted that the breeders could focus on developing early-maturing genotypes with specific traits such as days to first flower and boll opening. The selected genotypes exhibited positive genetic gains for traits related to earliness and a successful harvest during the first and second pickings. However, there were negative gains for traits related to flowering and boll opening. Conclusion The study identified three genotypes exhibiting early maturity and accomplished well under different planting densities. The multivariate methods (MTSI and MGIDI) serve as novel approaches for selecting desired genotypes in plant breeding programs, especially across various growing environments. These methods offer exclusive benefits and can easily construe and minimize multicollinearity issues.

Background: Globally, the cultivation of cotton is constrained by its tendency for extended periods of growth. Early maturity plays a potential role in rainfed-based multiple cropping system especially in the current era of climate change. In the current study, a set of 20 diverse Gossypium hirsutum genotypes were evaluated in two crop seasons with three planting densities and assessed for 11 morphological traits related to early maturity. The study aimed to identify genotype(s) that mature rapidly and accomplish well under diverse environmental conditions based on the two robust multivariate techniques called multi-trait stability index (MTSI) and multi-trait genotype-ideotype distance index (MGIDI). Results: MTSI analysis revealed that out of the 20 genotypes, three genotypes, viz., NNDC-30, A-2, and S-32 accomplished well in terms of early maturity traits in two seasons. Furthermore, three genotypes were selected using MGIDI method for each planting densities with a selection intensity of 15%. The strengths and weaknesses of the genotypes selected based on MGIDI method highlighted that the breeders could focus on developing early-maturing genotypes with specific traits such as days to first flower and boll opening. The selected genotypes exhibited positive genetic gains for traits related to earliness and a successful harvest during the first and second pickings. However, there were negative gains for traits related to flowering and boll opening. Conclusion: The study identified three genotypes exhibiting early maturity and accomplished well under different planting densities. The multivariate methods (MTSI and MGIDI) serve as novel approaches for selecting desired genotypes in plant breeding programs, especially across various growing environments. These methods offer exclusive benefits and can easily construe and minimize multicollinearity issues. [ABSTRACT FROM AUTHOR]

The size and growth patterns of nestling birds are key determinants of their survival up to fledging and long‐term fitness. However, because traits such as feathers, skeleton and body mass can follow different developmental trajectories, our understanding of the impact of adverse weather on development requires insights into trait‐specific sensitive developmental windows. We analysed data from nestling Alpine swifts in Switzerland measured throughout growth up to the age of 50 days (i.e. fledging between 50 and 70 days), for wing length and body mass (2693 nestlings in 25 years) and sternum length (2447 nestlings in 22 years). We show that the sensitive developmental windows for wing and sternum length corresponded to the periods of trait‐specific peak growth, which span almost the whole developmental period for wings and the first half for the sternum. Adverse weather conditions during these periods slowed down growth and reduced size. Although nestling body mass at 50 days showed the greatest inter‐individual variation, this was explained by weather in the two days before measurement rather than during peak growth. Interestingly, the relationship between temperature and body mass was not linear, and the initial sharp increase in body mass associated with the increase in temperature was followed by a moderate drop on hot days, likely linked to heat stress. Nestlings experiencing adverse weather conditions during wing growth had lower survival rates up to fledging and fledged at later ages, presumably to compensate for slower wing growth. Overall, our results suggest that measures of feather growth and, to some extent, skeletal growth best capture the consequences of adverse weather conditions throughout the whole development of offspring, while body mass better reflects the short, instantaneous effects of weather conditions on their body reserves (i.e. energy depletion vs. storage in unfavourable vs. favourable conditions). [ABSTRACT FROM AUTHOR]

Eucalyptus globulus Labill., is a recognized multipurpose tree, which stands out not only for the valuable qualities of its wood but also for the medicinal applications of the essential oil extracted from its leaves. In this study, we implemented an integrated strategy comprising genomic and phenomic approaches to predict foliar essential oil content, stem quality, and growth-related traits within a 9-year-old breeding population of E. globulus. The strategy involved evaluating Uni/Multi-trait deep learning (DL) models by incorporating genomic data related to single nucleotide polymorphisms (SNPs) and haplotypes, as well as the phenomic data from leaf near-infrared (NIR) spectroscopy. Our results showed that essential oil content (oil yield) ranged from 0.01 to 1.69% v/fw and had no significant correlation with any growth-related traits. This suggests that selection solely based on growth-related traits did n The emphases (colored text) from revisions were removed throughout the article. Confirm that this change is fine. ot influence the essential oil content. Genomic heritability estimates ranged from 0.25 (diameter at breast height (DBH) and oil yield) to 0.71 (DBH and stem straightness (ST)), while pedigree-based heritability exhibited a broader range, from 0.05 to 0.88. Notably, oil yield was found to be moderate to highly heritable, with genomic values ranging from 0.25 to 0.60, alongside a pedigree-based estimate of 0.48. The DL prediction models consistently achieved higher prediction accuracy (PA) values with a Multi-trait approach for most traits analyzed, including oil yield (0.699), tree height (0.772), DBH (0.745), slenderness coefficient (0.616), stem volume (0.757), and ST (0.764). The Uni-trait approach achieved superior PA values solely for branching quality (0.861). NIR spectral absorbance was the best omics data for CNN or MLP models with a Multi-trait approach. These results highlight considerable genetic variation within the Eucalyptus progeny trial, particularly regarding oil production. Our results contribute significantly to understanding omics-assisted deep learning models as a breeding strategy to improve growth-related traits and optimize essential oil production in this species.

Large-scale phenotype data are expected to increase the accuracy of genome-wide prediction and the power of genome-wide association analyses. However, genomic analyses of high-dimensional, highly correlated traits are challenging. We developed a method for implementing high-dimensional Bayesian multivariate regression to simultaneously analyze genetic variants underlying thousands of traits. As a demonstration, we implemented the BayesC prior in the R package MegaLMM. Applied to Genomic Prediction, MegaBayesC effectively integrated hyperspectral reflectance data from 620 hyperspectral wavelengths to improve the accuracy of genetic value prediction on grain yield in a wheat dataset. Applied to Genome-Wide Association Studies, we used simulations to show that MegaBayesC can accurately estimate the effect sizes of QTL across a range of genetic architectures and causes of correlations among traits. To apply MegaBayesC to a realistic scenario involving whole-genome marker data, we developed a 2-stage procedure involving a preliminary step of candidate marker selection prior to multivariate regression. We then used MegaBayesC to identify genetic associations with flowering time in Arabidopsis thaliana, leveraging expression data from 20,843 genes. MegaBayesC selected 15 single nucleotide polymorphisms as important for flowering time, with 13 located within 100 kb of known flowering-time related genes, a higher validation rate than achieved by a single-stage analysis using only the flowering time data itself. These results demonstrate that MegaBayesC can efficiently and effectively leverage high-dimensional phenotypes in genetic analyses.

Phylogenetic signal (PS) is the propensity of closely related species to resemble each other. PS has been tested across clades of terrestrial plants; however, insight for seagrasses is still lacking. Signatures of PS and models of niche (trait) evolution can help to detect phylogenetic niche conservatism (PNC), that is, close relatives live in comparable niches.The initial goal of this study was to assess the pattern of PS for the world's seagrasses, by testing the non‐independence of phylogenetic relatedness and seagrass species traits. A phylogeny of 49 seagrasses was constructed, together with a matrix of nine traits covering morphological, life‐history and reproductive attributes. PS of traits was tested through complementary indices (Pagel's λ, Blomberg's K, Moran's I and Abouheif's Cmean). Three models of niche evolution (Brownian Motion, BM; Ornstein–Uhlenbeck, OU; Early Burst, EB) were then fitted to each trait and the multivariate trait matrix.Results supported the existence of strong PS for seagrasses, with a particularly large effect size for seagrass reproductive traits, which followed an EB evolution model. Local Indicators of Phylogenetic Association (local autocorrelation metrics that can help to identify areas of large autocorrelation) supported the presence of PS across seagrass lineages/clades, supported by the dominance of OU as the most parsimonious trait evolution model.The pattern of strong PS seems to be a consequence of long‐term PNC of seagrass traits after initial radiation.Synthesis: Our study highlights the relevance of evolution from common ancestors and shared history underpinning large seagrass phylogenetic structuring. [ABSTRACT FROM AUTHOR]

Multi-environment trials (MET) are crucial for selecting genotypes that are well-suited to different environmental conditions. Incorporating multiple traits in the analysis can provide more reliable recommendations for selecting genotypes with desirable traits, including resistance to the Mungbean Yellow Mosaic Virus (MYMV) and high yield potential. The use of a Multi-Trait Stability Index (MTSI) is a good approach for analyzing the stability of genotypes across multiple traits under MYMV stress. In the present investigation, the performance of thirteen green gram genotypes were evaluated for traits such as yield, plant height, number of branches per plant, and resistance to MYMV. The main objective of the study is to identify highly productive and stable mung bean genotypes resistant to MYMV. MTSI can be calculated by combining information on the performance of genotypes across multiple traits and environmental conditions to provide a single index that indicates the overall stability of genotypes across traits and environments. The results helped to identify two green gram genotypes (Yadadri and JNG-18) that were high-yielding with stable resistance to MYMV stress across multiple environmental conditions. This can provide useful information to breeders for the development of suitable genotypes against MYMV in the affected areas.

Feeding represents the largest economic cost in meat production; therefore, selection to improve traits related to feed efficiency is a goal in most livestock breeding programs. Residual feed intake (RFI), that is, the difference between the actual and the expected feed intake based on animal's requirements, has been used as the selection criteria to improve feed efficiency since it was proposed by Kotch in 1963. In growing pigs, it is computed as the residual of the multiple regression model of daily feed intake (DFI), on average daily gain (ADG), backfat thickness (BFT), and metabolic body weight (MW). Recently, prediction using single‐output machine learning algorithms and information from SNPs as predictor variables have been proposed for genomic selection in growing pigs, but like in other species, the prediction quality achieved for RFI has been generally poor. However, it has been suggested that it could be improved through multi‐output or stacking methods. For this purpose, four strategies were implemented to predict RFI. Two of them correspond to the computation of RFI in an indirect way using the predicted values of its components obtained from (i) individual (multiple single‐output strategy) or (ii) simultaneous predictions (multi‐output strategy). The other two correspond to the direct prediction of RFI using (iii) the individual predictions of its components as predictor variables jointly with the genotype (stacking strategy), or (iv) using only the genotypes as predictors of RFI (single‐output strategy). The single‐output strategy was considered the benchmark. This research aimed to test the former three hypotheses using data recorded from 5828 growing pigs and 45,610 SNPs. For all the strategies two different learning methods were fitted: random forest (RF) and support vector regression (SVR). A nested cross‐validation (CV) with an outer 10‐folds CV and an inner threefold CV for hyperparameter tuning was implemented to test all strategies. This scheme was repeated using as predictor variables different subsets with an increasing number (from 200 to 3000) of the most informative SNPs identified with RF. Results showed that the highest prediction performance was achieved with 1000 SNPs, although the stability of feature selection was poor (0.13 points out of 1). For all SNP subsets, the benchmark showed the best prediction performance. Using the RF as a learner and the 1000 most informative SNPs as predictors, the mean (SD) of the 10 values obtained in the test sets were: 0.23 (0.04) for the Spearman correlation, 0.83 (0.04) for the zero–one loss, and 0.33 (0.03) for the rank distance loss. We conclude that the information on predicted components of RFI (DFI, ADG, MW, and BFT) does not contribute to improve the quality of the prediction of this trait in relation to the one obtained with the single‐output strategy. [ABSTRACT FROM AUTHOR]

Genomic selection (GS) plays a pivotal role in hybrid prediction. It can enhance the selection of parental lines, accurately predict hybrid performance, and harness hybrid vigor. Likewise, it can optimize breeding strategies by reducing field trial requirements, expediting hybrid development, facilitating targeted trait improvement, and enhancing adaptability to diverse environments. Leveraging genomic information empowers breeders to make informed decisions and significantly improve the efficiency and success rate of hybrid breeding programs. In order to improve the genomic ability performance, we explored the incorporation of parental phenotypic information as covariates under a multi-trait framework. Approach 1, referred to as Pmean, directly utilized parental phenotypic information without any preprocessing. While approach 2, denoted as BV, replaced the direct use of phenotypic values of both parents with their respective breeding values. While an improvement in prediction performance was observed in both approaches, with a minimum 4.24% reduction in the normalized root mean square error (NRMSE), the direct incorporation of parental phenotypic information in the Pmean approach slightly outperformed the BV approach. We also compared these two approaches using linear and nonlinear kernels, but no relevant gain was observed. Finally, our results increase empirical evidence confirming that the integration of parental phenotypic information helps increase the prediction performance of hybrids. [ABSTRACT FROM AUTHOR]

Selection for feed efficiency is the goal for many genetic breeding programs in beef cattle. Residual feed intake has been included in genetic evaluations to reduce feed intake without compromising performance traits as liveweight, body gain or carcass traits. However, measuring feed intake is expensive, and only a small percentage of selection candidates are phenotyped. Genomic selection has become a very important tool to achieve effective genetic progress in these traits. Another effective strategy has been the implementation of multi‐trait prediction using easily recordable predictor traits on both reference animals and candidates without phenotypes, and this could be another inexpensive way to increase accuracy. The objective of this work was to analyse and compare the prediction ability of two alternative different approaches to predict GEBVs for RFI. The population of inference was Hereford bulls in Uruguay that were genotyped candidates for to selection. The first model was the conventional univariate model for RFI and the second model was a multi‐trait model which included a predictor trait (weaning weight, WW), in addition to the traits used in the first one (dry matter intake, metabolic mid test weight, average daily gain and ultrasound back fat) (DMI, MWT, ADG, UBF, respectively). GEBVs from the multi‐trait model were combined using selection index theory to derive RFI values. All analyses were performed using ssGBLUP procedure. The prediction ability of both models was tested using two validation strategies (30 different replicates of random groups of animals and validation across 9 different feed intake tests). The prediction quality was assessed by the following parameters: bias, dispersion, ratio of accuracies and the relative increase in accuracy by adding phenotypic information. All parameters showed that the univariate model outperforms the multi‐trait model, regardless of the validation strategy considered. These results indicate that including WW as a proxy trait in a multi‐trait analysis does not improve the prediction ability when all animals to be predicted are genotyped. [ABSTRACT FROM AUTHOR]

Maize (Zea mays L.), the third most widely cultivated cereal crop in the world, plays a critical role in global food security. To improve the efficiency of selecting superior genotypes in breeding programs, researchers have aimed to identify key genomic regions that impact agronomic traits. In this study, the performance of multi-trait, multi-environment deep learning models was compared to that of Bayesian models (Markov Chain Monte Carlo generalized linear mixed models (MCMCglmm), Bayesian Genomic Genotype-Environment Interaction (BGGE), and Bayesian Multi-Trait and Multi-Environment (BMTME)) in terms of the prediction accuracy of flowering-related traits (Anthesis-Silking Interval: ASI, Female Flowering: FF, and Male Flowering: MF). A tropical maize panel of 258 inbred lines from Brazil was evaluated in three sites (Cambira-2018, Sabaudia-2018, and Iguatemi-2020 and 2021) using approximately 290,000 single nucleotide polymorphisms (SNPs). The results demonstrated a 14.4% increase in prediction accuracy when employing multi-trait models compared to the use of a single trait in a single environment approach. The accuracy of predictions also improved by 6.4% when using a single trait in a multi-environment scheme compared to using multi-trait analysis. Additionally, deep learning models consistently outperformed Bayesian models in both single and multiple trait and environment approaches. A complementary genome-wide association study identified associations with 26 candidate genes related to flowering time traits, and 31 marker-trait associations were identified, accounting for 37%, 37%, and 22% of the phenotypic variation of ASI, FF and MF, respectively. In conclusion, our findings suggest that deep learning models have the potential to significantly improve the accuracy of predictions, regardless of the approach used and provide support for the efficacy of this method in genomic selection for flowering-related traits in tropical maize. [ABSTRACT FROM AUTHOR]

The main objective of the study was to evaluate and select the superior barley genotypes based on grain yield and some pheno-morphological traits using a new proposed selection index (SIIG). For this purpose, one-hundred-eight pure and four local cultivars (Norouz, Auxin, Nobahar, and WB-97-11) were evaluated as reference genotypes in four warm regions of Iran, including Ahvaz, Darab, Zabol, and Gonbad, during the 2020–2021 cropping seasons. The results of REML analysis showed that the heritability of all traits (except plant height) was higher in Gonbad than in other environments, while the lowest values were estimated in Ahvaz and Zabol environments. In addition, among the measured traits, the thousand kernel weight and grain filling period showed the highest and lowest values of heritability (0.83 and 0.01, respectively). The results showed that the seed yield of genotypes 1, 108, 3, 86, 5, 87, 19, 16, 15, 56, and 18 was higher than the four reference genotypes, and, on the other hand, the SIIG index of these genotypes was greater than or equal to 0.60. Based on the SIIG discriminator index, 4, 8, 31, and 28 genotypes with values greater than or equal to 0.60 were identified as superior for Darab, Ahvaz, Zabol, and Gonbad environments, respectively. As a conclusion, our results revealed that the SIIG index has ideal potential to identify genotypes with high yield and desirable traits. Therefore, the use of this index can be beneficial in screening better genotypes in the early stages of any breeding program for any crop. [ABSTRACT FROM AUTHOR]

While sparse testing methods have been proposed by researchers to improve the efficiency of genomic selection (GS) in breeding programs, there are several factors that can hinder this. In this research, we evaluated four methods (M1–M4) for sparse testing allocation of lines to environments under multi-environmental trails for genomic prediction of unobserved lines. The sparse testing methods described in this study are applied in a two-stage analysis to build the genomic training and testing sets in a strategy that allows each location or environment to evaluate only a subset of all genotypes rather than all of them. To ensure a valid implementation, the sparse testing methods presented here require BLUEs (or BLUPs) of the lines to be computed at the first stage using an appropriate experimental design and statistical analyses in each location (or environment). The evaluation of the four cultivar allocation methods to environments of the second stage was done with four data sets (two large and two small) under a multi-trait and uni-trait framework. We found that the multi-trait model produced better genomic prediction (GP) accuracy than the uni-trait model and that methods M3 and M4 were slightly better than methods M1 and M2 for the allocation of lines to environments. Some of the most important findings, however, were that even under a scenario where we used a training-testing relation of 15–85%, the prediction accuracy of the four methods barely decreased. This indicates that genomic sparse testing methods for data sets under these scenarios can save considerable operational and financial resources with only a small loss in precision, which can be shown in our cost-benefit analysis. [ABSTRACT FROM AUTHOR]

Large-scale phenotype data are expected to increase the accuracy of genome-wide prediction and the power of genome-wide association analyses. However, genomic analyses of high-dimensional, highly correlated traits are challenging. We developed a method for implementing high-dimensional Bayesian multivariate regression to simultaneously analyze genetic variants underlying thousands of traits. As a demonstration, we implemented the BayesC prior in the R package MegaLMM. Applied to Genomic Prediction, MegaBayesC effectively integrated hyperspectral reflectance data from 620 hyperspectral wavelengths to improve the accuracy of genetic value prediction on grain yield in a wheat dataset. Applied to Genome-Wide Association Studies, we used simulations to show that MegaBayesC can accurately estimate the effect sizes of QTL across a range of genetic architectures and causes of correlations among traits. To apply MegaBayesC to a realistic scenario involving whole-genome marker data, we developed a 2-stage procedure involving a preliminary step of candidate marker selection prior to multivariate regression. We then used MegaBayesC to identify genetic associations with flowering time in Arabidopsis thaliana, leveraging expression data from 20,843 genes. MegaBayesC selected 15 single nucleotide polymorphisms as important for flowering time, with 13 located within 100 kb of known flowering-time related genes, a higher validation rate than achieved by a single-stage analysis using only the flowering time data itself. These results demonstrate that MegaBayesC can efficiently and effectively leverage high-dimensional phenotypes in genetic analyses. [ABSTRACT FROM AUTHOR]

Incorporating the genotypic and phenotypic of the correlated traits into the multi-trait model can significantly improve the prediction accuracy of the target trait in animal and plant breeding, as well as human genetics. However, in most cases, the phenotypic information of the correlated and target trait of the individual to be evaluated was null simultaneously, particularly for the newborn. Therefore, we propose a machine learning framework, MAK, to improve the prediction accuracy of the target trait by constructing the multi-target ensemble regression chains and selecting the assistant trait automatically, which predicted the genomic estimated breeding values of the target trait using genotypic information only. The prediction ability of MAK was significantly more robust than the genomic best linear unbiased prediction, BayesB, BayesRR and the multi trait Bayesian method in the four real animal and plant datasets, and the computational efficiency of MAK was roughly 100 times faster than BayesB and BayesRR. [ABSTRACT FROM AUTHOR]

Accuracy of predicted breeding values (PBV) for low heritability traits may be increased in early generations by exploiting the information available in correlated traits. We compared the accuracy of PBV for 10 correlated traits with low to medium narrow-sense heritability ( h 2) in a genetically diverse field pea (Pisum sativum L.) population after univariate or multivariate linear mixed model (MLMM) analysis with pedigree information. In the contra-season, we crossed and selfed S1 parent plants, and in the main season we evaluated spaced plants of S0 cross progeny and S2+ (S2 or higher) self progeny of parent plants for the 10 traits. Stem strength traits included stem buckling (SB) ( h 2 = 0.05), compressed stem thickness (CST) ( h 2 = 0.12), internode length (IL) ( h 2 = 0.61) and angle of the main stem above horizontal at first flower (EAngle) ( h 2 = 0.46). Significant genetic correlations of the additive effects occurred between SB and CST (0.61), IL and EAngle (−0.90) and IL and CST (−0.36). The average accuracy of PBVs in S0 progeny increased from 0.799 to 0.841 and in S2+ progeny increased from 0.835 to 0.875 in univariate vs MLMM, respectively. An optimized mating design was constructed with optimal contribution selection based on an index of PBV for the 10 traits, and predicted genetic gain in the next cycle ranged from 1.4% (SB), 5.0% (CST), 10.5% (EAngle) and −10.5% (IL), with low achieved parental coancestry of 0.12. MLMM improved the potential genetic gain in annual cycles of early generation selection in field pea by increasing the accuracy of PBV. [ABSTRACT FROM AUTHOR]

In sheep, body weight is an economically important trait. This study sought to map genetic loci related to weaning weight and yearling weight. To this end, a single-trait and multi-trait genome-wide association study (GWAS) was performed using a high-density 600 K single nucleotide polymorphism (SNP) chip. The results showed that 43 and 56 SNPs were significantly associated with weaning weight and yearling weight, respectively. A region associated with both weaning and yearling traits (OARX: 6.74–7.04 Mb) was identified, suggesting that the same genes could play a role in regulating both these traits. This region was found to contain three genes (TBL1X, SHROOM2 and GPR143). The most significant SNP was Affx-281066395, located at 6.94 Mb (p = 1.70 × 10−17), corresponding to the SHROOM2 gene. We also identified 93 novel SNPs elated to sheep weight using multi-trait GWAS analysis. A new genomic region (OAR10: 76.04–77.23 Mb) with 22 significant SNPs were discovered. Combining transcriptomic data from multiple tissues and genomic data in sheep, we found the HINT1, ASB11 and GPR143 genes may involve in sheep body weight. So, multi-omic anlaysis is a valuable strategy identifying candidate genes related to body weight.

Maize (Zea mays L.), the third most widely cultivated cereal crop in the world, plays a critical role in global food security. To improve the efficiency of selecting superior genotypes in breeding programs, researchers have aimed to identify key genomic regions that impact agronomic traits. In this study, the performance of multi-trait, multi-environment deep learning models was compared to that of Bayesian models (Markov Chain Monte Carlo generalized linear mixed models (MCMCglmm), Bayesian Genomic Genotype-Environment Interaction (BGGE), and Bayesian Multi-Trait and Multi-Environment (BMTME)) in terms of the prediction accuracy of flowering-related traits (Anthesis-Silking Interval: ASI, Female Flowering: FF, and Male Flowering: MF). A tropical maize panel of 258 inbred lines from Brazil was evaluated in three sites (Cambira-2018, Sabaudia-2018, and Iguatemi-2020 and 2021) using approximately 290,000 single nucleotide polymorphisms (SNPs). The results demonstrated a 14.4% increase in prediction accuracy when employing multi-trait models compared to the use of a single trait in a single environment approach. The accuracy of predictions also improved by 6.4% when using a single trait in a multi-environment scheme compared to using multi-trait analysis. Additionally, deep learning models consistently outperformed Bayesian models in both single and multiple trait and environment approaches. A complementary genome-wide association study identified associations with 26 candidate genes related to flowering time traits, and 31 marker-trait associations were identified, accounting for 37%, 37%, and 22% of the phenotypic variation of ASI, FF and MF, respectively. In conclusion, our findings suggest that deep learning models have the potential to significantly improve the accuracy of predictions, regardless of the approach used and provide support for the efficacy of this method in genomic selection for flowering-related traits in tropical maize.

Incorporating measurements on correlated traits into genomic prediction models can increase prediction accuracy and selection gain. However, multi-trait genomic prediction models are complex and prone to overfitting which may result in a loss of prediction accuracy relative to single-trait genomic prediction. Cross-validation is considered the gold standard method for selecting and tuning models for genomic prediction in both plant and animal breeding. When used appropriately, cross-validation gives an accurate estimate of the prediction accuracy of a genomic prediction model, and can effectively choose among disparate models based on their expected performance in real data. However, we show that a naive cross-validation strategy applied to the multi-trait prediction problem can be severely biased and lead to sub-optimal choices between single and multi-trait models when secondary traits are used to aid in the prediction of focal traits and these secondary traits are measured on the individuals to be tested. We use simulations to demonstrate the extent of the problem and propose three partial solutions: 1) a parametric solution from selection index theory, 2) a semi-parametric method for correcting the cross-validation estimates of prediction accuracy, and 3) a fully non-parametric method which we call CV2*: validating model predictions against focal trait measurements from genetically related individuals. The current excitement over high-throughput phenotyping suggests that more comprehensive phenotype measurements will be useful for accelerating breeding programs. Using an appropriate cross-validation strategy should more reliably determine if and when combining information across multiple traits is useful.

The present study was carried out to estimate the genetic parameters for direct and maternal influences on Mecheri sheep (Ovis aries) growth traits using Bayesian multi-trait animal model. The genetic parameters were calculated using data from 2825 Mecheri lambs born between 2010 and 2020 that were kept in semi-arid tropical climate. Mecheri sheep body weight (mean ± SE) at various stages, viz. BW, WW, and BW12, were 2.6 ± 0.01, 11.1 ± 0.05, and 20.7 ± 0.13 kg, respectively. The Mecheri sheep gained 71.5 percent of their body weight at the age of 6 months. With the exception of birth weight, the weights of the animals varied considerably (P < 0.01) by the year of birth. The fixed effect of sex significantly (P < 0.01) influenced all the growth traits examined. The direct estimates of heritability (± SD) for BW, WW, BW6, and BW12 was 0.21 ± 0.041, 0.21 ± 0.041, 0.12 ± 0.052, and 0.13 ± 0.053, respectively, and the maternal heritability for BW, WW, BW6, and BW12 was 0.18 ± 0.021, 0.08 ± 0.023, 0.11 ± 0.022, and 0.13 ± 0.033, respectively. Significant variance was indicated by moderately larger direct heritability estimates for BWT and WWT, indicating that there will be more opportunities for selection response during the genetic improvement programme. For the majority of the variables examined, direct heritability values were higher than maternal heritability values. The additive genetic correlation between WW and BW6, BW9, and BW12 was 0.70 ± 0.145, 0.57 ± 0.171, and 0.50 ± 0.194, respectively. The maternal genetic correlations ranged from 0.06 ± 0.152 (BW–BW12) to 0.86 ± 0.046 (BW6–BW9), and the residual correlation varied from 0.18 ± 0.034 (BW–WW) to 0.85 ± 0.013 (BW9–BW12). The BW had a stronger genetic relationship with WW and a weaker relationship with other growth traits. The WW and BW6 showed a favourable genetic relationship, even if the tendency was decreasing in the latter stages of maturation. The BW6 and BW9 demonstrated the strongest positive genetic relationship (0.90 ± 0.052) of all the variables analysed. After examining the heritabilities and genetic correlation between WW and BW6, it was established that selection based on at WW was beneficial, since it manifests early in life and would result in moderate genetic progress through selection. [ABSTRACT FROM AUTHOR]

While genomic selection (GS) began revolutionizing plant breeding when it was proposed around 20 years ago, its practical implementation is still challenging as many factors affect its accuracy. One such factor is the choice of the statistical machine learning method. For this reason, we explore the tuning process under a multi-trait framework using the Gaussian kernel with a multi-trait Bayesian Best Linear Unbiased Predictor (GBLUP) model. We explored three methods of tuning (manual, grid search and Bayesian optimization) using 5 real datasets of breeding programs. We found that using grid search and Bayesian optimization improve between 1.9 and 6.8% the prediction accuracy regarding of using manual tuning. While the improvement in prediction accuracy in some cases can be marginal, it is very important to carry out the tuning process carefully to improve the accuracy of the GS methodology, even though this entails greater computational resources. [ABSTRACT FROM AUTHOR]

Bayesian multiple-regression methods incorporating different mixture priors for marker effects are used widely in genomic prediction. Improvement in prediction accuracies from using those methods, such as BayesB, BayesC, and BayesCπ, have been shown in single-trait analyses with both simulated and real data. These methods have been extended to multi-trait analyses, but only under the restrictive assumption that a locus simultaneously affects all the traits or none of them. This assumption is not biologically meaningful, especially in multi-trait analyses involving many traits. In this paper, we develop and implement a more general multi-trait BayesC[Formula: see text] and BayesB methods allowing a broader range of mixture priors. Our methods allow a locus to affect any combination of traits, e.g., in a 5-trait analysis, the "restrictive" model only allows two situations, whereas ours allow all 32 situations. Further, we compare our methods to single-trait methods and the "restrictive" multi-trait formulation using real and simulated data. In the real data analysis, higher prediction accuracies were observed from both our new broad-based multi-trait methods and the "restrictive" formulation. The broad-based and restrictive multi-trait methods showed similar prediction accuracies. In the simulated data analysis, higher prediction accuracies to the "restrictive" method were observed from our general multi-trait methods for intermediate training population size. The software tool JWAS offers open-source routines to perform these analyses.

Multi-trait and multi-environment analyses can improve genomic prediction by exploiting between-trait correlations and genotype-by-environment interactions. In the context of reaction norm models, genotype-byenvironment interactions can be described as functions of high-dimensional sets of markers and environmental covariates. However, comprehensive multi-trait reaction norm models accounting for marker x environmental covariates interactions are lacking. In this article, we propose to extend a reaction norm model incorporating genotype-by-environment interactions through (co)variance structures of markers and environmental covariates to a multi-trait reaction norm case. To do that, we propose a novel methodology for characterizing the environment at different growth stages based on growth degree-days (GDD). The proposed models were evaluated by variance components estimation and predictive performance for winter wheat grain yield and protein content in a set of 2,015 F6-lines. Cross-validation analyses were performed using leave-one-year-location-out (CV1) and leave-onebreeding-cycle-out (CV2) strategies. The modeling of genomic [SNPs] x environmental covariates interactions significantly improved predictive ability and reduced the variance inflation of predicted genetic values for grain yield and protein content in both cross-validation schemes. Trait-assisted genomic prediction was carried out for multi-trait models, and it significantly enhanced predictive ability and reduced variance inflation in all scenarios. The genotype by environment interaction modeling via genomic [SNPs] x environmental covariates interactions, combined with trait-assisted genomic prediction, boosted the benefits in predictive performance. The proposed multi-trait reaction norm methodology is a comprehensive approach that allows capitalizing on the benefits of multi-trait models accounting for between-trait correlations and reaction norm models exploiting high-dimensional genomic and environmental information. [ABSTRACT FROM AUTHOR]

Eucalyptus globulus Labill., is a recognized multipurpose tree, which stands out not only for the valuable qualities of its wood but also for the medicinal applications of the essential oil extracted from its leaves. In this study, we implemented an integrated strategy comprising genomic and phenomic approaches to predict foliar essential oil content, stem quality, and growth-related traits within a 9-year-old breeding population of E. globulus . The strategy involved evaluating Uni/Multi-trait deep learning (DL) models by incorporating genomic data related to single nucleotide polymorphisms (SNPs) and haplotypes, as well as the phenomic data from leaf near-infrared (NIR) spectroscopy. Our results showed that essential oil content (oil yield) ranged from 0.01 to 1.69% v/fw and had no significant correlation with any growth-related traits. This suggests that selection solely based on growth-related traits did n The emphases (colored text) from revisions were removed throughout the article. Confirm that this change is fine. ot influence the essential oil content. Genomic heritability estimates ranged from 0.25 (diameter at breast height (DBH) and oil yield) to 0.71 (DBH and stem straightness (ST)), while pedigree-based heritability exhibited a broader range, from 0.05 to 0.88. Notably, oil yield was found to be moderate to highly heritable, with genomic values ranging from 0.25 to 0.60, alongside a pedigree-based estimate of 0.48. The DL prediction models consistently achieved higher prediction accuracy (PA) values with a Multi-trait approach for most traits analyzed, including oil yield (0.699), tree height (0.772), DBH (0.745), slenderness coefficient (0.616), stem volume (0.757), and ST (0.764). The Uni-trait approach achieved superior PA values solely for branching quality (0.861). NIR spectral absorbance was the best omics data for CNN or MLP models with a Multi-trait approach. These results highlight considerable genetic variation within the Eucalyptus progeny trial, particularly regarding oil production. Our results contribute significantly to understanding omics-assisted deep learning models as a breeding strategy to improve growth-related traits and optimize essential oil production in this species., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Mieres-Castro, Maldonado and Mora-Poblete.)

Multi-trait and multi-environment analyses can improve genomic prediction by exploiting between-trait correlations and genotype-by-environment interactions. In the context of reaction norm models, genotype-by-environment interactions can be described as functions of high-dimensional sets of markers and environmental covariates. However, comprehensive multi-trait reaction norm models accounting for marker × environmental covariates interactions are lacking. In this article, we propose to extend a reaction norm model incorporating genotype-by-environment interactions through (co)variance structures of markers and environmental covariates to a multi-trait reaction norm case. To do that, we propose a novel methodology for characterizing the environment at different growth stages based on growth degree-days (GDD). The proposed models were evaluated by variance components estimation and predictive performance for winter wheat grain yield and protein content in a set of 2,015 F6-lines. Cross-validation analyses were performed using leave-one-year-location-out (CV1) and leave-one-breeding-cycle-out (CV2) strategies. The modeling of genomic [SNPs] × environmental covariates interactions significantly improved predictive ability and reduced the variance inflation of predicted genetic values for grain yield and protein content in both cross-validation schemes. Trait-assisted genomic prediction was carried out for multi-trait models, and it significantly enhanced predictive ability and reduced variance inflation in all scenarios. The genotype by environment interaction modeling via genomic [SNPs] × environmental covariates interactions, combined with trait-assisted genomic prediction, boosted the benefits in predictive performance. The proposed multi-trait reaction norm methodology is a comprehensive approach that allows capitalizing on the benefits of multi-trait models accounting for between-trait correlations and reaction norm models exploiting high-dimensional genomic and environmental information.

Genomic selection (GS) changed the way plant breeders select genotypes. GS takes advantage of phenotypic and genotypic information to training a statistical machine learning model, which is used to predict phenotypic (or breeding) values of new lines for which only genotypic information is available. Therefore, many statistical machine learning methods have been proposed for this task. Multi-trait (MT) genomic prediction models take advantage of correlated traits to improve prediction accuracy. Therefore, some multivariate statistical machine learning methods are popular for GS. In this paper, we compare the prediction performance of three MT methods: the MT genomic best linear unbiased predictor (GBLUP), the MT partial least squares (PLS) and the multi-trait random forest (RF) methods. Benchmarking was performed with six real datasets. We found that the three investigated methods produce similar results, but under predictors with genotype (G) and environment (E), that is, E + G, the MT GBLUP achieved superior performance, whereas under predictors E + G + genotype × environment (GE) and G + GE, random forest achieved the best results. We also found that the best predictions were achieved under the predictors E + G and E + G + GE. Here, we also provide the R code for the implementation of these three statistical machine learning methods in the sparse kernel method (SKM) library, which offers not only options for single-trait prediction with various statistical machine learning methods but also some options for MT predictions that can help to capture improved complex patterns in datasets that are common in genomic selection. [ABSTRACT FROM AUTHOR]

The main objective of the study was to evaluate and select the superior barley genotypes based on grain yield and some pheno-morphological traits using a new proposed selection index (SIIG). For this purpose, one-hundred-eight pure and four local cultivars (Norouz, Auxin, Nobahar, and WB-97-11) were evaluated as reference genotypes in four warm regions of Iran, including Ahvaz, Darab, Zabol, and Gonbad, during the 2020–2021 cropping seasons. The results of REML analysis showed that the heritability of all traits (except plant height) was higher in Gonbad than in other environments, while the lowest values were estimated in Ahvaz and Zabol environments. In addition, among the measured traits, the thousand kernel weight and grain filling period showed the highest and lowest values of heritability (0.83 and 0.01, respectively). The results showed that the seed yield of genotypes 1, 108, 3, 86, 5, 87, 19, 16, 15, 56, and 18 was higher than the four reference genotypes, and, on the other hand, the SIIG index of these genotypes was greater than or equal to 0.60. Based on the SIIG discriminator index, 4, 8, 31, and 28 genotypes with values greater than or equal to 0.60 were identified as superior for Darab, Ahvaz, Zabol, and Gonbad environments, respectively. As a conclusion, our results revealed that the SIIG index has ideal potential to identify genotypes with high yield and desirable traits. Therefore, the use of this index can be beneficial in screening better genotypes in the early stages of any breeding program for any crop.

Heat waves are expected to become more frequent and intense, which will impact faba bean cultivation globally. Conventional breeding methods are effective but take considerable time to achieve breeding goals, and, therefore, the identification of molecular markers associated with key genes controlling heat tolerance can facilitate and accelerate efficient variety development. We phenotyped 134 accessions in six open field experiments during summer seasons at Terbol, Lebanon, at Hudeiba, Sudan, and at Central Ferry, WA, USA from 2015 to 2018. These accessions were genotyped using genotyping by sequencing (GBS), and 10,794 high quality single nucleotide polymorphisms (SNPs) were discovered. These accessions were clustered in one diverse large group, although several discrete groups may exist surrounding it. Fifteen lines belonging to different botanical groups were identified as tolerant to heat. SNPs associated with heat tolerance using single-trait (ST) and multi-trait (MT) genome-wide association studies (GWASs) showed 9 and 11 significant associations, respectively. Through the annotation of the discovered significant SNPs, we found that SNPs from transcription factor helix–loop–helix bHLH143-like S-adenosylmethionine carrier, putative pentatricopeptide repeat-containing protein At5g08310, protein NLP8-like, and photosystem II reaction center PSB28 proteins are associated with heat tolerance. [ABSTRACT FROM AUTHOR]

[Objectives]This paper aimed to explore the predictive effect of missing phenotypes by multi-trait imputation methods under vorious conditions. It is an efficient way to predict the phenotypic missing data by statistical methods, which increase the sample size and improve the accuracy of data analysis. [Methods] In the simulation studies, multi-trait imputation methods of mean, KNN ( K-nearest neighbor), decision tree, MICE (multiple imputation by chained equations), PHENIX (phenotype imputation expediated) and softImpute were used to predict the multi-trait missing data. Meanwhile, we compared the results of imputation methods under various phenotypic missing rates, number of traits, sample sizes and trait correlation. Furthermore, the real Arabidopsis thaliana phenotypic data of days to flowering under long day, days to flowering under short day, days to flowering under long day with vernalization and days to flowering under short day with vernalization were imputed by multi-trait imputation methods, and the reliability of the imputed data were verified by genome-wide association analysis. [Results]Simulation study showed that the accuracy of imputation decreased along with increasing of phenotypic missing rate; stronger trait correlation and more traits improved the accuracy of imputation; the larger the sample sizes were, the more stable the results would be. In real data analysis, the tendency of multi-trait imputation was consistent with the simulation experiment, and the results were verified by genome-wide association study and confirmed genes. [Conclusions]Phenotypic missing rate, number of traits and trait correlation play important roles in imputation. Multi-trait imputed methods, including PHENIX,decision tree and KNN easily capture the genetic structure, and thus are more accurate and efficient in simulation study and real data analysis. [ABSTRACT FROM AUTHOR]

Abstract Background Body traits are generally controlled by several genes in vertebrates (i.e. polygenes), which in turn make them difficult to identify through association mapping. Increasing the power of association studies by combining approaches such as genotype imputation and multi-trait analysis improves the ability to detect quantitative trait loci associated with polygenic traits, such as body traits. Results A multi-trait genome-wide association study (mtGWAS) was performed to identify quantitative trait loci (QTL) and genes associated with body traits in Nile tilapia (Oreochromis niloticus) using genotypes imputed to whole-genome sequences (WGS). To increase the statistical power of mtGWAS for the detection of genetic associations, summary statistics from single-trait genome-wide association studies (stGWAS) for eight different body traits recorded in 1309 animals were used. The mtGWAS increased the statistical power from the original sample size from 13 to 44%, depending on the trait analyzed. The better resolution of the WGS data, combined with the increased power of the mtGWAS approach, allowed the detection of significant markers which were not previously found in the stGWAS. Some of the lead single nucleotide polymorphisms (SNPs) were found within important functional candidate genes previously associated with growth-related traits in other terrestrial species. For instance, we identified SNP within the α1,6-fucosyltransferase (FUT8), solute carrier family 4 member 2 (SLC4A2), A disintegrin and metalloproteinase with thrombospondin motifs 9 (ADAMTS9) and heart development protein with EGF like domains 1 (HEG1) genes, which have been associated with average daily gain in sheep, osteopetrosis in cattle, chest size in goats, and growth and meat quality in sheep, respectively. Conclusions The high-resolution mtGWAS presented here allowed the identification of significant SNPs, linked to strong functional candidate genes, associated with body traits in Nile tilapia. These results provide further insights about the genetic variants and genes underlying body trait variation in cichlid fish with high accuracy and strong statistical support.

Soft white wheat is a wheat class used in foreign and domestic markets to make various end products requiring specific quality attributes. Due to associated cost, time, and amount of seed needed, phenotyping for the end-use quality trait is delayed until later generations. Previously, we explored the potential of using genomic selection (GS) for selecting superior genotypes earlier in the breeding program. Breeders typically measure multiple traits across various locations, and it opens up the avenue for exploring multi-trait–based GS models. This study's main objective was to explore the potential of using multi-trait GS models for predicting seven different end-use quality traits using cross-validation, independent prediction, and across-location predictions in a wheat breeding program. The population used consisted of 666 soft white wheat genotypes planted for 5 years at two locations in Washington, United States. We optimized and compared the performances of four uni-trait– and multi-trait–based GS models, namely, Bayes B, genomic best linear unbiased prediction (GBLUP), multilayer perceptron (MLP), and random forests. The prediction accuracies for multi-trait GS models were 5.5 and 7.9% superior to uni-trait models for the within-environment and across-location predictions. Multi-trait machine and deep learning models performed superior to GBLUP and Bayes B for across-location predictions, but their advantages diminished when the genotype by environment component was included in the model. The highest improvement in prediction accuracy, that is, 35% was obtained for flour protein content with the multi-trait MLP model. This study showed the potential of using multi-trait–based GS models to enhance prediction accuracy by using information from previously phenotyped traits. It would assist in speeding up the breeding cycle time in a cost-friendly manner. [ABSTRACT FROM AUTHOR]

Plant breeders regularly evaluate multiple traits across multiple environments, which opens an avenue for using multiple traits in genomic prediction models. We assessed the potential of multi-trait (MT) genomic prediction model through evaluating several strategies of incorporating multiple traits (eight agronomic and malting quality traits) into the prediction models with two cross-validation schemes (CV1, predicting new lines with genotypic information only and CV2, predicting partially phenotyped lines using both genotypic and phenotypic information from correlated traits) in barley. The predictive ability was similar for single (ST-CV1) and multi-trait (MT-CV1) models to predict new lines. However, the predictive ability for agronomic traits was considerably increased when partially phenotyped lines (MT-CV2) were used. The predictive ability for grain yield using the MT-CV2 model with other agronomic traits resulted in 57% and 61% higher predictive ability than ST-CV1 and MT-CV1 models, respectively. Therefore, complex traits such as grain yield are better predicted when correlated traits are used. Similarly, a considerable increase in the predictive ability of malting quality traits was observed when correlated traits were used. The predictive ability for grain protein content using the MT-CV2 model with both agronomic and malting traits resulted in a 76% higher predictive ability than ST-CV1 and MT-CV1 models. Additionally, the higher predictive ability for new environments was obtained for all traits using the MT-CV2 model compared to the MT-CV1 model. This study showed the potential of improving the genomic prediction of complex traits by incorporating the information from multiple traits (cost-friendly and easy to measure traits) collected throughout breeding programs which could assist in speeding up breeding cycles.

Soft white wheat is a wheat class used in foreign and domestic markets to make various end products requiring specific quality attributes. Due to associated cost, time, and amount of seed needed, phenotyping for the end-use quality trait is delayed until later generations. Previously, we explored the potential of using genomic selection (GS) for selecting superior genotypes earlier in the breeding program. Breeders typically measure multiple traits across various locations, and it opens up the avenue for exploring multi-trait–based GS models. This study’s main objective was to explore the potential of using multi-trait GS models for predicting seven different end-use quality traits using cross-validation, independent prediction, and across-location predictions in a wheat breeding program. The population used consisted of 666 soft white wheat genotypes planted for 5 years at two locations in Washington, United States. We optimized and compared the performances of four uni-trait– and multi-trait–based GS models, namely, Bayes B, genomic best linear unbiased prediction (GBLUP), multilayer perceptron (MLP), and random forests. The prediction accuracies for multi-trait GS models were 5.5 and 7.9% superior to uni-trait models for the within-environment and across-location predictions. Multi-trait machine and deep learning models performed superior to GBLUP and Bayes B for across-location predictions, but their advantages diminished when the genotype by environment component was included in the model. The highest improvement in prediction accuracy, that is, 35% was obtained for flour protein content with the multi-trait MLP model. This study showed the potential of using multi-trait–based GS models to enhance prediction accuracy by using information from previously phenotyped traits. It would assist in speeding up the breeding cycle time in a cost-friendly manner.

Yield prediction is the primary goal of genomic selection (GS)-assisted crop breeding. Because yield is a complex quantitative trait, making predictions from genotypic data is challenging. Transfer learning can produce an effective model for a target task by leveraging knowledge from a different, but related, source domain and is considered a great potential method for improving yield prediction by integrating multi-trait data. However, it has not previously been applied to genotype-to-phenotype prediction owing to the lack of an efficient implementation framework. We therefore developed TrG2P, a transfer-learning-based framework. TrG2P first employs convolutional neural networks (CNN) to train models using non-yield-trait phenotypic and genotypic data, thus obtaining pre-trained models. Subsequently, the convolutional layer parameters from these pre-trained models are transferred to the yield prediction task, and the fully connected layers are retrained, thus obtaining fine-tuned models. Finally, the convolutional layer and the first fully connected layer of the fine-tuned models are fused, and the last fully connected layer is trained to enhance prediction performance. We applied TrG2P to five sets of genotypic and phenotypic data from maize (Zea mays), rice (Oryza sativa), and wheat (Triticum aestivum) and compared its model precision to that of seven other popular GS tools: ridge regression best linear unbiased prediction (rrBLUP), random forest, support vector regression, light gradient boosting machine (LightGBM), CNN, DeepGS, and deep neural network for genomic prediction (DNNGP). TrG2P improved the accuracy of yield prediction by 39.9%, 6.8%, and 1.8% in rice, maize, and wheat, respectively, compared with predictions generated by the best-performing comparison model. Our work therefore demonstrates that transfer learning is an effective strategy for improving yield prediction by integrating information from non-yield-trait data. We attribute its enhanced prediction accuracy to the valuable information available from traits associated with yield and to training dataset augmentation. The Python implementation of TrG2P is available at https://github.com/lijinlong1991/TrG2P. The web-based tool is available at http://trg2p.ebreed.cn:81., (Copyright © 2024. Published by Elsevier Inc.)

Multi-environment trials (MET) are crucial for selecting genotypes that are well-suited to different environmental conditions. Incorporating multiple traits in the analysis can provide more reliable recommendations for selecting genotypes with desirable traits, including resistance to the Mungbean Yellow Mosaic Virus (MYMV) and high yield potential. The use of a Multi-Trait Stability Index (MTSI) is a good approach for analyzing the stability of genotypes across multiple traits under MYMV stress. In the present investigation, the performance of thirteen green gram genotypes were evaluated for traits such as yield, plant height, number of branches per plant, and resistance to MYMV. The main objective of the study is to identify highly productive and stable mung bean genotypes resistant to MYMV. MTSI can be calculated by combining information on the performance of genotypes across multiple traits and environmental conditions to provide a single index that indicates the overall stability of genotypes across traits and environments. The results helped to identify two green gram genotypes (Yadadri and JNG-18) that were high-yielding with stable resistance to MYMV stress across multiple environmental conditions. This can provide useful information to breeders for the development of suitable genotypes against MYMV in the affected areas., Competing Interests: The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Hosam O. El-ansary reports financial support was provided by 10.13039/501100002383King Saud University., (© 2024 The Authors.)

Summary: The root economics space is a useful framework for plant ecology but is rarely considered for crop ecophysiology. In order to understand root trait integration in winter wheat, we combined functional phenomics with trait economic theory, utilizing genetic variation, high‐throughput phenotyping, and multivariate analyses.We phenotyped a diversity panel of 276 genotypes for root respiration and architectural traits using a novel high‐throughput method for CO2 flux and the open‐source software RhizoVisionExplorer to analyze scanned images.We uncovered substantial variation in specific root respiration (SRR) and specific root length (SRL), which were primary indicators of root metabolic and structural costs. Multiple linear regression analysis indicated that lateral root tips had the greatest SRR, and the residuals from this model were used as a new trait. Specific root respiration was negatively correlated with plant mass. Network analysis, using a Gaussian graphical model, identified root weight, SRL, diameter, and SRR as hub traits. Univariate and multivariate genetic analyses identified genetic regions associated with SRR, SRL, and root branching frequency, and proposed gene candidates.Combining functional phenomics and root economics is a promising approach to improving our understanding of crop ecophysiology. We identified root traits and genomic regions that could be harnessed to breed more efficient crops for sustainable agroecosystems. See also the Commentary on this article by Li et al., 232: 5–7. [ABSTRACT FROM AUTHOR]

Estimates of variance components are needed for implementing genetic selection. This study was conducted to genetic parameters for production and reproductive traits on Indian Karan-Fries cattle using multi-trait repeatability animal model. Data collected from ICAR-National Dairy Research Institute, Karnal, India (from 1988 to 2019) were used. Single-trait and multi-trait repeatability animal models were used for parameter estimation. The posterior mean of Heritability estimates for 305-day milk yield (305-DMY), lactation milk yield (LMY), lactation length (LL) were 0.20 ± 0.03, 0.19 ± 0.03 and 0.06 ± 0.02, respectively. For age at first calving (AFC), calving interval (CI), and days open (DO), the posterior mean of heritability estimates were 0.24 ± 0.08, 0.06 ± 0.01, and 0.07 ± 0.02, respectively. The repeatability estimates for 305-DMY, LMY, LL, CI, and DO were 0.37 ± 0.02, 0.34 ± 0.02, 0.15 ± 0.02, 0.09 ± 0.02, and 0.12 ± 0.02, respectively. Genetic correlation between milk production traits (305-DMY, LMY, and LL) was positive and strong (> 0.80). However, the genetic correlation between milk production trait and AFC ranges from − 0.31 to 0.12. Unfavorable strong genetic correlations were observed between production and reproductive traits (CI and DO) with values ranged from 0.5 to 0.7. Phenotypic correlations among 305-DMY, LMY, and LL were generally positive and high. The moderate heritability estimates and potential genetic variation for 305-DMY, TMY, and AFC suggested that genetic gain can be obtained for these traits through genetic selection. Low heritability estimates found for LL, CI and DO, indicating that the possibility of changing these traits through genetic selection is small. High genetic correlation observed between productive and fertility traits were unfavorable. The existed strong genetic and phenotypic correlation estimates between CI and DO indicates that recording only one of them would be sufficient in the herd. As the multi-trait model showed slight improvements in the h2 as well as r estimates for both productive and reproductive traits over univariate analysis, future selection with a multi-trait animal model applying Bayesian approach would be recommended. [ABSTRACT FROM AUTHOR]

• Number of dams with phenotypic records influenced genetic parameter estimates. • Direct heritability estimates remained stable across different data structures. • Covariance between direct and maternal effects was sensitive to data structure. • Maternal records influenced genetic associations, impacting multiple traits. The number of dams with their own phenotypic records is crucial for estimating the genetic correlation between direct and maternal effects in beef cattle. In this study, using multi-trait analyzes, we assessed the impact of varying the percentage of dams with own records (0, 20, 35, 50, and 70 %) on genetic parameter estimates in Nelore beef cattle. The analyzes focused on key traits in the selection index: weaning weight (WW), post-weaning weight gain (PWG), scrotal circumference, and muscling score. Direct heritability estimates remained relatively stable across data structures, while maternal heritability and the covariance between direct and maternal effects showed important variations. The genetic correlation between direct and maternal effects of WW linearly shifted from negative in D0 (−0.20) to positive in D70 (0.31), emphasizing the impact of increasing the number of dams with own records. Maternal phenotypic records influenced genetic associations between traits, notably WW and PWG (0.13 in D0 and 0.36 in D70). These findings highlight that changes in data structure can impact the estimation of genetic parameters for multiple traits. The study underscores the importance of robust data structures, maternal phenotypic records (>35 %), and complete pedigree information to obtain reliable estimates of genetic parameters in beef cattle. [ABSTRACT FROM AUTHOR]

Incorporating measurements on correlated traits into genomic prediction models can increase prediction accuracy and selection gain. However, multi-trait genomic prediction models are complex and prone to overfitting which may result in a loss of prediction accuracy relative to single-trait genomic prediction. Cross-validation is considered the gold standard method for selecting and tuning models for genomic prediction in both plant and animal breeding. When used appropriately, cross-validation gives an accurate estimate of the prediction accuracy of a genomic prediction model, and can effectively choose among disparate models based on their expected performance in real data. However, we show that a naive cross-validation strategy applied to the multi-trait prediction problem can be severely biased and lead to sub-optimal choices between single and multi-trait models when secondary traits are used to aid in the prediction of focal traits and these secondary traits are measured on the individuals to be tested. We use simulations to demonstrate the extent of the problem and propose three partial solutions: 1) a parametric solution from selection index theory, 2) a semi-parametric method for correcting the cross-validation estimates of prediction accuracy, and 3) a fully non-parametric method which we call CV2*: validating model predictions against focal trait measurements from genetically related individuals. The current excitement over high-throughput phenotyping suggests that more comprehensive phenotype measurements will be useful for accelerating breeding programs. Using an appropriate cross-validation strategy should more reliably determine if and when combining information across multiple traits is useful.

In this paper we propose a Bayesian multi-output regressor stacking (BMORS) model that is a generalization of the multi-trait regressor stacking method. The proposed BMORS model consists of two stages: in the first stage, a univariate genomic best linear unbiased prediction (GBLUP including genotype × environment interaction GE) model is implemented for each of the L traits under study; then the predictions of all traits are included as covariates in the second stage, by implementing a Ridge regression model. The main objectives of this research were to study alternative models to the existing multi-trait multi-environment (BMTME) model with respect to (1) genomic-enabled prediction accuracy, and (2) potential advantages in terms of computing resources and implementation. We compared the predictions of the BMORS model to those of the univariate GBLUP model using 7 maize and wheat datasets. We found that the proposed BMORS produced similar predictions to the univariate GBLUP model and to the BMTME model in terms of prediction accuracy; however, the best predictions were obtained under the BMTME model. In terms of computing resources, we found that the BMORS is at least 9 times faster than the BMTME method. Based on our empirical findings, the proposed BMORS model is an alternative for predicting multi-trait and multi-environment data, which are very common in genomic-enabled prediction in plant and animal breeding programs.

Abstract Background Plant breeders seek to develop cultivars with maximal agronomic value, which is often assessed using numerous, often genetically correlated traits. As intervention on one trait will affect the value of another, breeding decisions should consider the relationships among traits in the context of putative causal structures (i.e., trait networks). While multi-trait genome-wide association studies (MTM-GWAS) can infer putative genetic signals at the multivariate scale, standard MTM-GWAS does not accommodate the network structure of phenotypes, and therefore does not address how the traits are interrelated. We extended the scope of MTM-GWAS by incorporating trait network structures into GWAS using structural equation models (SEM-GWAS). Here, we illustrate the utility of SEM-GWAS using a digital metric for shoot biomass, root biomass, water use, and water use efficiency in rice. Results A salient feature of SEM-GWAS is that it can partition the total single nucleotide polymorphism (SNP) effects acting on a trait into direct and indirect effects. Using this novel approach, we show that for most QTL associated with water use, total SNP effects were driven by genetic effects acting directly on water use rather that genetic effects originating from upstream traits. Conversely, total SNP effects for water use efficiency were largely due to indirect effects originating from the upstream trait, projected shoot area. Conclusions We describe a robust framework that can be applied to multivariate phenotypes to understand the interrelationships between complex traits. This framework provides novel insights into how QTL act within a phenotypic network that would otherwise not be possible with conventional multi-trait GWAS approaches. Collectively, these results suggest that the use of SEM may enhance our understanding of complex relationships among agronomic traits.

Evidence that genomic selection (GS) is a technology that is revolutionizing plant breeding continues to grow. However, it is very well documented that its success strongly depends on statistical models, which are used by GS to perform predictions of candidate genotypes that were not phenotyped. Because there is no universally better model for prediction and models for each type of response variable are needed (continuous, binary, ordinal, count, etc.), an active area of research aims to develop statistical models for the prediction of univariate and multivariate traits in GS. However, most of the models developed so far are for univariate and continuous (Gaussian) traits. Therefore, to overcome the lack of multivariate statistical models for genome-based prediction by improving the original version of the BMTME, we propose an improved Bayesian multi-trait and multi-environment (BMTME) R package for analyzing breeding data with multiple traits and multiple environments. We also introduce Bayesian multi-output regressor stacking (BMORS) functions that are considerably efficient in terms of computational resources. The package allows parameter estimation and evaluates the prediction performance of multi-trait and multi-environment data in a reliable, efficient and user-friendly way. We illustrate the use of the BMTME with real toy datasets to show all the facilities that the software offers the user. However, for large datasets, the BME() and BMTME() functions of the BMTME R package are very intense in terms of computing time; on the other hand, less intensive computing is required with BMORS functions BMORS() and BMORS_Env() that are also included in the BMTME package.

Multi-trait and multi-environment data are common in animal and plant breeding programs. However, what is lacking are more powerful statistical models that can exploit the correlation between traits to improve prediction accuracy in the context of genomic selection (GS). Multi-trait models are more complex than univariate models and usually require more computational resources, but they are preferred because they can exploit the correlation between traits, which many times helps improve prediction accuracy. For this reason, in this paper we explore the power of multi-trait deep learning (MTDL) models in terms of prediction accuracy. The prediction performance of MTDL models was compared to the performance of the Bayesian multi-trait and multi-environment (BMTME) model proposed by Montesinos-López et al. (2016), which is a multi-trait version of the genomic best linear unbiased prediction (GBLUP) univariate model. Both models were evaluated with predictors with and without the genotype×environment interaction term. The prediction performance of both models was evaluated in terms of Pearson’s correlation using cross-validation. We found that the best predictions in two of the three data sets were found under the BMTME model, but in general the predictions of both models, BTMTE and MTDL, were similar. Among models without the genotype×environment interaction, the MTDL model was the best, while among models with genotype×environment interaction, the BMTME model was superior. These results indicate that the MTDL model is very competitive for performing predictions in the context of GS, with the important practical advantage that it requires less computational resources than the BMTME model.