Search

Orphan perennial native species are gaining importance as sustainability in agriculture becomes crucial to mitigate climate change. Nevertheless, issues related to the undomesticated status and lack of improved germplasm impede the evolution of formal agricultural initiatives. Acrocomia aculeata - a neotropical palm with potential for oil production - is an example. Breeding efforts can aid the species to reach its full potential and increase market competitiveness. Here, we present genomic information and training set optimization as alternatives to boost orphan perennial native species breeding using Acrocomia aculeata as an example. Furthermore, we compared three SNP calling methods and, for the first time, presented the prediction accuracies of three yield-related traits. We collected data for two years from 201 wild individuals. These trees were genotyped, and three references were used for SNP calling: the oil palm genome, de novo sequencing, and the A. aculeata transcriptome. The traits analyzed were fruit dry mass (FDM), pulp dry mass (PDM), and pulp oil content (OC). We compared the predictive ability of GBLUP and BayesB models in crossand real validation procedures. Afterwards, we tested several optimization criteria regarding consistency and the ability to provide the optimized training set that yielded less risk in both targeted and untargeted scenarios. Using the oil palm genome as a reference and GBLUP models had better results for the genomic prediction of FDM, OC, and PDM (prediction accuracies of 0.46, 0.45, and 0.39, respectively). Using the criteria PEV, r-score and core collection methodology provides risk-averse decisions. Training set optimization is an alternative to improve decision-making while leveraging genomic information as a cost-saving tool to accelerate plant domestication and breeding. The optimized training set can be used as a reference for the characterization of native species populations, aiding in decisions involving germplasm collection and construction of breeding populations. [ABSTRACT FROM AUTHOR]

Cultivated potato, Solanum tuberosum L., is considered an autotetraploid with 12 chromosomes with four homologous phases. However, recent evidence found that, due to frequent large phase deletions in the genome, gene ploidy is not constant across the genome. The elite cultivar "Otava" was found to have an average gene copy number of 3.2 across all loci. Breeding programs for elite potato cultivars rely increasingly on genomic prediction tools for selection breeding and elucidation of quantitative trait loci underpinning trait genetic variance. These are typically based on anonymous single nucleotide polymorphism (SNP) markers, which are usually called from, for example, SNP array or sequencing data using a tetraploid model. In this study, we analyzed the impact of using whole genome markers genotyped as either tetraploid or observed allele frequencies from genotype-by-sequencing data on single-trait additive genomic best linear unbiased prediction (GBLUP) genomic prediction (GP) models and single-marker regression genome-wide association studies of potato to evaluate the implications of capturing varying ploidy on the statistical models employed in genomic breeding. A panel of 762 offspring of a diallel cross of 18 parents of elite breeding material was used for modeling. These were genotyped by sequencing and phenotyped for five key performance traits: chipping quality, length/width ratio, senescence, dry matter content, and yield. We also estimated the read coverage required to confidently discriminate between a heterozygous triploid and tetraploid state from simulated data. It was found that using a tetraploid model neither impaired nor improved genomic predictions compared to using the observed allele frequencies that account for true marker ploidy. In genome-wide associations studies (GWAS), very minor variations of both signal amplitude and number of SNPs supporting both minor and major quantitative trait loci (QTLs) were observed between the two data sets. However, all major QTLs were reproducible using both data sets. [ABSTRACT FROM AUTHOR]

Increasing SNP density by incorporating sequence information only marginally increases prediction accuracies of breeding values in livestock. To find out why, we used statistical models and simulations to investigate the shape of distribution of estimated SNP effects (a profile) around quantitative trait nucleotides (QTNs) in populations with a small effective population size (Ne). A QTN profile created by averaging SNP effects around each QTN was similar to the shape of expected pairwise linkage disequilibrium (PLD) based on Ne and genetic distance between SNP, with a distinct peak for the QTN. Populations with smaller Ne showed lower but wider QTN profiles. However, adding more genotyped individuals with phenotypes dragged the profile closer to the QTN. The QTN profile was higher and narrower for populations with larger compared to smaller Ne. Assuming the PLD curve for the QTN profile, 80% of the additive genetic variance explained by each QTN was contained in ± 1/Ne Morgan interval around the QTN, corresponding to 2 Mb in cattle and 5 Mb in pigs and chickens. With such large intervals, identifying QTN is difficult even if all of them are in the data and the assumed genetic architecture is simplistic. Additional complexity in QTN detection arises from confounding of QTN profiles with signals due to relationships, overlapping profiles with closely spaced QTN, and spurious signals. However, small Ne allows for accurate predictions with large data even without QTN identification because QTNs are accounted for by QTN profiles if SNP density is sufficient to saturate the segments. [ABSTRACT FROM AUTHOR]

Coffee Breeding programs have traditionally relied on observing plant characteristics over years, a slow and costly process. Genomic selection (GS) offers a DNA-based alternative for faster selection of superior cultivars. Stacking Ensemble Learning (SEL) combines multiple models for potentially even more accurate selection. This study explores SEL potential in coffee breeding, aiming to improve prediction accuracy for important traits [yield (YL), total number of the fruits (NF), leaf miner infestation (LM), and cercosporiosis incidence (Cer)] in Coffea Arabica. We analyzed data from 195 individuals genotyped for 21,211 single-nucleotide polymorphism (SNP) markers. To comprehensively assess model performance, we employed a cross-validation (CV) scheme. Genomic Best Linear Unbiased Prediction (GBLUP), multivariate adaptive regression splines (MARS), Quantile Random Forest (QRF), and Random Forest (RF) served as base learners. For the meta-learner within the SEL framework, various options were explored, including Ridge Regression, RF, GBLUP, and Single Average. The SEL method was able to predict the predictive ability (PA) of important traits in Coffea Arabica. SEL presented higher PA compared with those obtained for all base learner methods. The gains in PA in relation to GBLUP were 87.44% (the ratio between the PA obtained from best Stacking model and the GBLUP), 37.83%, 199.82%, and 14.59% for YL, NF, LM and Cer, respectively. Overall, SEL presents a promising approach for GS. By combining predictions from multiple models, SEL can potentially enhance the PA of GS for complex traits. [ABSTRACT FROM AUTHOR]

Sugarcane smut and Pachymetra root rots are two serious diseases of sugarcane, with susceptible infected crops losing over 30% of yield. A heritable component to both diseases has been demonstrated, suggesting selection could improve disease resistance. Genomic selection could accelerate gains even further, enabling early selection of resistant seedlings for breeding and clonal propagation. In this study we evaluated four types of algorithms for genomic predictions of clonal performance for disease resistance. These algorithms were: Genomic best linear unbiased prediction (GBLUP), including extensions to model dominance and epistasis, Bayesian methods including BayesC and BayesR, Machine learning methods including random forest, multilayer perceptron (MLP), modified convolutional neural network (CNN) and attention networks designed to capture epistasis across the genome-wide markers. Simple hybrid methods, that first used BayesR/GWAS to identify a subset of 1000 markers with moderate to large marginal additive effects, then used attention networks to derive predictions from these effects and their interactions, were also developed and evaluated. The hypothesis for this approach was that using a subset of markers more likely to have an effect would enable better estimation of interaction effects than when there were an extremely large number of possible interactions, especially with our limited data set size. To evaluate the methods, we applied both random five-fold cross-validation and a structured PCA based cross-validation that separated 4702 sugarcane clones (that had disease phenotypes and genotyped for 26k genome wide SNP markers) by genomic relationship. The Bayesian methods (BayesR and BayesC) gave the highest accuracy of prediction, followed closely by hybrid methods with attention networks. The hybrid methods with attention networks gave the lowest variation in accuracy of prediction across validation folds (and lowest MSE), which may be a criteria worth considering in practical breeding programs. This suggests that hybrid methods incorporating the attention mechanism could be useful for genomic prediction of clonal performance, particularly where nonadditive effects may be important. [ABSTRACT FROM AUTHOR]

Orphan perennial native species are gaining importance as sustainability in agriculture becomes crucial to mitigate climate change. Nevertheless, issues related to the undomesticated status and lack of improved germplasm impede the evolution of formal agricultural initiatives. Acrocomia aculeata - a neotropical palm with potential for oil production - is an example. Breeding efforts can aid the species to reach its full potential and increase market competitiveness. Here, we present genomic information and training set optimization as alternatives to boost orphan perennial native species breeding using Acrocomia aculeata as an example. Furthermore, we compared three SNP calling methods and, for the first time, presented the prediction accuracies of three yield-related traits. We collected data for two years from 201 wild individuals. These trees were genotyped, and three references were used for SNP calling: the oil palm genome, de novo sequencing, and the A. aculeata transcriptome. The traits analyzed were fruit dry mass (FDM), pulp dry mass (PDM), and pulp oil content (OC). We compared the predictive ability of GBLUP and BayesB models in cross- and real validation procedures. Afterwards, we tested several optimization criteria regarding consistency and the ability to provide the optimized training set that yielded less risk in both targeted and untargeted scenarios. Using the oil palm genome as a reference and GBLUP models had better results for the genomic prediction of FDM, OC, and PDM (prediction accuracies of 0.46, 0.45, and 0.39, respectively). Using the criteria PEV, r-score and core collection methodology provides risk-averse decisions. Training set optimization is an alternative to improve decision-making while leveraging genomic information as a cost-saving tool to accelerate plant domestication and breeding. The optimized training set can be used as a reference for the characterization of native species populations, aiding in decisions involving germplasm collection and construction of breeding populations

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenanceprogeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positivelywith narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed. [ABSTRACT FROM AUTHOR]

Genomic prediction and genome-wide association studies are becoming widely employed in potato key performance trait QTL identifications and to support potato breeding using genomic selection. Elite cultivars are tetraploid and highly heterozygous but also share many common ancestors and generation-spanning inbreeding events, resulting from the clonal propagation of potatoes through seed potatoes. Consequentially, many SNP markers are not in a 1:1 relationship with a single allele variant but shared over several alleles that might exert varying effects on a given trait. The impact of such redundant "diluted" predictors on the statistical models underpinning genome-wide association studies (GWAS) and genomic prediction has scarcely been evaluated despite the potential impact on model accuracy and performance. We evaluated the impact of marker location, marker type, and marker density on the genomic prediction and GWAS of five key performance traits in tetraploid potato (chipping quality, dry matter content, length/width ratio, senescence, and yield). A 762-offspring panel of a diallel cross of 18 elite cultivars was genotyped by sequencing, and markers were annotated according to a reference genome. Genomic prediction models (GBLUP) were trained on four marker subsets [nonsynonymous (29,553 SNPs), synonymous (31,229), non-coding (32,388), and a combination], and robustness to marker reduction was investigated. Singlemarker regression GWAS was performed for each trait and marker subset. The best cross-validated prediction correlation coefficients of 0.54, 0.75, 0.49, 0.35, and 0.28 were obtained for chipping quality, dry matter content, length/width ratio, senescence, and yield, respectively. The trait prediction abilities were similar across all marker types, with only non-synonymous variants improving yield predictive ability by 16%. Marker reduction response did not depend on marker type but rather on trait. Traits with high predictive abilities, e.g., dry matter content, reached a plateau using fewer markers than traits with intermediate-low correlations, such as yield. The predictions were unbiased across all traits, marker types, and all marker densities >100 SNPs. Our results suggest that using non-synonymous variants does not enhance the performance of genomic prediction of most traits. The major known QTLs were identified by GWAS and were reproducible across exonic and wholegenome variant sets for dry matter content, length/width ratio, and senescence. In contrast, minor QTL detection was marker type dependent. [ABSTRACT FROM AUTHOR]

Cultivated potato, Solanum tuberosum L., is considered an autotetraploid with 12 chromosomes with four homologous phases. However, recent evidence found that, due to frequent large phase deletions in the genome, gene ploidy is not constant across the genome. The elite cultivar “Otava” was found to have an average gene copy number of 3.2 across all loci. Breeding programs for elite potato cultivars rely increasingly on genomic prediction tools for selection breeding and elucidation of quantitative trait loci underpinning trait genetic variance. These are typically based on anonymous single nucleotide polymorphism (SNP) markers, which are usually called from, for example, SNP array or sequencing data using a tetraploid model. In this study, we analyzed the impact of using whole genome markers genotyped as either tetraploid or observed allele frequencies from genotype-by-sequencing data on single-trait additive genomic best linear unbiased prediction (GBLUP) genomic prediction (GP) models and single-marker regression genome-wide association studies of potato to evaluate the implications of capturing varying ploidy on the statistical models employed in genomic breeding. A panel of 762 offspring of a diallel cross of 18 parents of elite breeding material was used for modeling. These were genotyped by sequencing and phenotyped for five key performance traits: chipping quality, length/width ratio, senescence, dry matter content, and yield. We also estimated the read coverage required to confidently discriminate between a heterozygous triploid and tetraploid state from simulated data. It was found that using a tetraploid model neither impaired nor improved genomic predictions compared to using the observed allele frequencies that account for true marker ploidy. In genome-wide associations studies (GWAS), very minor variations of both signal amplitude and number of SNPs supporting both minor and major quantitative trait loci (QTLs) were observed between the two data sets. However, all major QTLs were reproducible using both data sets.

Coffee Breeding programs have traditionally relied on observing plant characteristics over years, a slow and costly process. Genomic selection (GS) offers a DNA-based alternative for faster selection of superior cultivars. Stacking Ensemble Learning (SEL) combines multiple models for potentially even more accurate selection. This study explores SEL potential in coffee breeding, aiming to improve prediction accuracy for important traits [yield (YL), total number of the fruits (NF), leaf miner infestation (LM), and cercosporiosis incidence (Cer)] in Coffea Arabica. We analyzed data from 195 individuals genotyped for 21,211 single-nucleotide polymorphism (SNP) markers. To comprehensively assess model performance, we employed a cross-validation (CV) scheme. Genomic Best Linear Unbiased Prediction (GBLUP), multivariate adaptive regression splines (MARS), Quantile Random Forest (QRF), and Random Forest (RF) served as base learners. For the meta-learner within the SEL framework, various options were explored, including Ridge Regression, RF, GBLUP, and Single Average. The SEL method was able to predict the predictive ability (PA) of important traits in Coffea Arabica. SEL presented higher PA compared with those obtained for all base learner methods. The gains in PA in relation to GBLUP were 87.44% (the ratio between the PA obtained from best Stacking model and the GBLUP), 37.83%, 199.82%, and 14.59% for YL, NF, LM and Cer, respectively. Overall, SEL presents a promising approach for GS. By combining predictions from multiple models, SEL can potentially enhance the PA of GS for complex traits.

Abstract Genomic dissection of genetic effects on desirable traits and the subsequent use of genomic selection hold great promise for accelerating the rate of genetic improvement of forest tree species. In this study, a total of 661 offspring trees from 66 open-pollinated families of Japanese larch (Larix kaempferi (Lam.) Carrière) were sampled at a test site. The contributions of additive and non-additive effects (dominance, imprinting and epistasis) were evaluated for nine valuable traits related to growth, wood physical and chemical properties, and competitive ability using three pedigree-based and four Genomics-based Best Linear Unbiased Predictions (GBLUP) models and used to determine the genetic model. The predictive ability (PA) of two genomic prediction methods, GBLUP and Reproducing Kernel Hilbert Spaces (RKHS), was compared. The traits could be classified into two types based on different quantitative genetic architectures: for type I, including wood chemical properties and Pilodyn penetration, additive effect is the main source of variation (38.20-67.46%); for type II, including growth, competitive ability and acoustic velocity, epistasis plays a significant role (50.76-91.26%). Dominance and imprinting showed low to moderate contributions (

Genomic dissection of genetic effects on desirable traits and the subsequent use of genomic selection hold great promise for accelerating the rate of genetic improvement of forest tree species. In this study, a total of 661 offspring trees from 66 open-pollinated families of Japanese larch (Larix kaempferi (Lam.) Carrière) were sampled at a test site. The contributions of additive and non-additive effects (dominance, imprinting and epistasis) were evaluated for nine valuable traits related to growth, wood physical and chemical properties, and competitive ability using three pedigree-based and four Genomics-based Best Linear Unbiased Predictions (GBLUP) models and used to determine the genetic model. The predictive ability (PA) of two genomic prediction methods, GBLUP and Reproducing Kernel Hilbert Spaces (RKHS), was compared. The traits could be classified into two types based on different quantitative genetic architectures: for type I, including wood chemical properties and Pilodyn penetration, additive effect is the main source of variation (38.20-67.46%); for type II, including growth, competitive ability and acoustic velocity, epistasis plays a significant role (50.76-91.26%). Dominance and imprinting showed low to moderate contributions (< 36.26%). GBLUP was more suitable for traits of type I (PAs = 0.37–0.39 vs. 0.14–0.25), and RKHS was more suitable for traits of type II (PAs = 0.23–0.37 vs. 0.07–0.23). Non-additive effects make no meaningful contribution to the enhancement of PA of GBLUP method for all traits. These findings enhance our current understanding of the architecture of quantitative traits and lay the foundation for the development of genomic selection strategies in Japanese larch. [ABSTRACT FROM AUTHOR]

Simple Summary: Historically, Korea has been an agrarian society where cattle played a significant role in Korean culture as working animals. However, in recent decades, there has been an increase in meat demand and the expansion of the Korean economy. This shift has led to the predominant use of native Korean Hanwoo cattle for beef production since the 1960s. Genetic improvement programs for Hanwoo cattle traditionally focused on carcass and growth qualities due to the accessibility of trait information and the simplicity of analysis techniques, rather than reproductive traits. However, there has been a recent surge of interest in genetic analyses of the reproductive traits of Hanwoo cattle because these traits significantly impact calf productivity. To enhance beef production and the overall profitability of Hanwoo farming, it is imperative to implement genomic predictions for age at first calving (AFC), calving interval (CI), gestation length (GL), and number of artificial inseminations per conception (NAIPC) to better understand and improve their response to selection. This study aimed to estimate heritability and the accuracy of genomic estimated breeding values (GEBVs) using genomic best linear unbiased prediction (GBLUP) and Bayesian methods (BayesB, BayesLASSO, and BayesR) for traits under study. Our analysis revealed relatively lower heritability values for these traits and indicated that the accuracy of genomic prediction across all methods applied was similarly reduced, likely due to the inherent lower heritability of reproductive traits. As a result, the findings of this study provide valuable insights into the genetic breeding programs of the beef cattle industry. This study aimed to predict the accuracy of genomic estimated breeding values (GEBVs) for reproductive traits in Hanwoo cows using the GBLUP, BayesB, BayesLASSO, and BayesR methods. Accuracy estimates of GEBVs for reproductive traits were derived through fivefold cross-validation, analyzing a dataset comprising 11,348 animals and employing an Illumina Bovine 50K SNP chip. GBLUP showed an accuracy of 0.26 for AFC, while BayesB, BayesLASSO, and BayesR demonstrated values of 0.28, 0.29, and 0.29, respectively. For CI, GBLUP attained an accuracy of 0.19, whereas BayesB, BayesLASSO, and BayesR scored 0.21, 0.24, and 0.25, respectively. The accuracy for GL was uniform across GBLUP, BayesB, and BayesR at 0.31, whereas BayesLASSO showed a slightly higher accuracy of 0.33. For NAIPC, GBLUP showed an accuracy of 0.24, while BayesB, BayesLASSO, and BayesR recorded 0.22, 0.27, and 0.30, respectively. The variation in genomic prediction accuracy among methods indicated Bayesian approaches slightly outperformed GBLUP. The findings suggest that Bayesian methods, notably BayesLASSO and BayesR, offer improved predictive capabilities for reproductive traits. Future research may explore more advanced genomic approaches to enhance predictive accuracy and genetic gains in Hanwoo cattle breeding programs. [ABSTRACT FROM AUTHOR]

Copyright of Forest Research is the property of Forest Research Editorial Office and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Introduction: Genomic selection (GS) experiments in forest trees have largely reported estimates of predictive abilities from cross-validation among individuals in the same breeding generation. In such conditions, no effects of recombination, selection, drift, and environmental changes are accounted for. Here, we assessed the effectively realized predictive ability (RPA) for volume growth at harvest age by GS across generations in an operational reciprocal recurrent selection (RRS) program of hybrid Eucalyptus. Methods: Genomic best linear unbiased prediction with additive (GBLUP_G), additive plus dominance (GBLUP_G+D), and additive single-step (HBLUP) models were trained with different combinations of growth data of hybrids and pure species individuals (N = 17,462) of the G1 generation, 1,944 of which were genotyped with ~16,000 SNPs from SNP arrays. The hybrid G2 progeny trial (HPT267) was the GS target, with 1,400 selection candidates, 197 of which were genotyped still at the seedling stage, and genomically predicted for their breeding and genotypic values at the operational harvest age (6 years). Seedlings were then grown to harvest and measured, and their pedigree-based breeding and genotypic values were compared to their originally predicted genomic counterparts. Results: Genomic RPAs ≥0.80 were obtained as the genetic relatedness between G1 and G2 increased, especially when the direct parents of selection candidates were used in training. GBLUP_G+D reached RPAs ≥0.70 only when hybrid or pure species data of G1 were included in training. HBLUP was only marginally better than GBLUP. Correlations ≥0.80 were obtained between pedigree and genomic individual ranks. Rank coincidence of the top 2.5% selections was the highest for GBLUP_G (45% to 60%) compared to GBLUP_G+D. To advance the pure species RRS populations, GS models were best when trained on pure species than hybrid data, and HBLUP yielded ~20% higher predictive abilities than GBLUP, but was not better than ABLUP for ungenotyped trees. Discussion: We demonstrate that genomic data effectively enable accurate ranking of eucalypt hybrid seedlings for their yet-to-be observed volume growth at harvest age. Our results support a two-stage GS approach involving family selection by average genomic breeding value, followed by within-top-families individual GS, significantly increasing selection intensity, optimizing genotyping costs, and accelerating RRS breeding. [ABSTRACT FROM AUTHOR]

Key message: A new genomic prediction method (RHPP) was developed via combining randomized Haseman–Elston regression (RHE-reg), PCR based on genomic information of core population, and preconditioned conjugate gradient (PCG) algorithm. Computational efficiency is becoming a hot issue in the practical application of genomic prediction due to the large number of data generated by the high-throughput genotyping technology. In this study, we developed a fast genomic prediction method RHPP via combining randomized Haseman–Elston regression (RHE-reg), PCR based on genomic information of core population, and preconditioned conjugate gradient (PCG) algorithm. The simulation results demonstrated similar prediction accuracy between RHPP and GBLUP, and significantly higher computational efficiency of the former with the increase of individuals. The results of real datasets of both bread wheat and loblolly pine demonstrated that RHPP had a similar or better predictive accuracy in most cases compared with GBLUP. In the future, RHPP may be an attractive choice for analyzing large-scale and high-dimensional data. [ABSTRACT FROM AUTHOR]

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Genomic prediction and genome-wide association studies are becoming widely employed in potato key performance trait QTL identifications and to support potato breeding using genomic selection. Elite cultivars are tetraploid and highly heterozygous but also share many common ancestors and generation-spanning inbreeding events, resulting from the clonal propagation of potatoes through seed potatoes. Consequentially, many SNP markers are not in a 1:1 relationship with a single allele variant but shared over several alleles that might exert varying effects on a given trait. The impact of such redundant “diluted” predictors on the statistical models underpinning genome-wide association studies (GWAS) and genomic prediction has scarcely been evaluated despite the potential impact on model accuracy and performance. We evaluated the impact of marker location, marker type, and marker density on the genomic prediction and GWAS of five key performance traits in tetraploid potato (chipping quality, dry matter content, length/width ratio, senescence, and yield). A 762-offspring panel of a diallel cross of 18 elite cultivars was genotyped by sequencing, and markers were annotated according to a reference genome. Genomic prediction models (GBLUP) were trained on four marker subsets [non-synonymous (29,553 SNPs), synonymous (31,229), non-coding (32,388), and a combination], and robustness to marker reduction was investigated. Single-marker regression GWAS was performed for each trait and marker subset. The best cross-validated prediction correlation coefficients of 0.54, 0.75, 0.49, 0.35, and 0.28 were obtained for chipping quality, dry matter content, length/width ratio, senescence, and yield, respectively. The trait prediction abilities were similar across all marker types, with only non-synonymous variants improving yield predictive ability by 16%. Marker reduction response did not depend on marker type but rather on trait. Traits with high predictive abilities, e.g., dry matter content, reached a plateau using fewer markers than traits with intermediate-low correlations, such as yield. The predictions were unbiased across all traits, marker types, and all marker densities >100 SNPs. Our results suggest that using non-synonymous variants does not enhance the performance of genomic prediction of most traits. The major known QTLs were identified by GWAS and were reproducible across exonic and whole-genome variant sets for dry matter content, length/width ratio, and senescence. In contrast, minor QTL detection was marker type dependent.

The number of animal genotypes is rapidly increasing, and a major challenge for animal models is inverting the genomic relationship matrix (G). Matrix G has a limited dimensionality, and the algorithm for proven and young (APY) makes inverting a large G possible via the inverse of a block diagonal of G with a size equivalent to the dimensionality of G. APY divides genotyped animals into core and non-core groups, and breeding values of non-core animals are conditioned on the breeding values of core animals. Therefore, there is the possibility of opting out equations for non-core animals from the model. A methodology was presented for a reduced APY genomic BLUP (GBLUP) to equations for core animals. Using a small example dataset, the method was validated by the equality of the full and the reduced model analysis results. Absorption of fixed effect equations into random effect equations was successful in reducing the number of equations to solve and producing the same random effect solutions. Extending the method to APY single-step GBLUP (ssGBLUP) was not computationally justifiable. Other reduction techniques exist for ssGBLUP (regardless of APY or non-APY) that work by reducing the number of equations for non-genotyped animals. The number of equations can further be reduced by data pruning.

Abstract Multiple trait genomic selection incorporating correlated traits can improve the predictive ability of low‐heritability traits. In this study, we evaluated genomic prediction accuracy using multi‐trait BayesCπ method (MT‐BayesCπ), which allows for a broader range of mixture priors for important traits in beef cattle. We compared the prediction performance of MT‐BayesCπ with single‐trait genomic best linear unbiased prediction (ST‐GBLUP), multi‐trait GBLUP (MT‐GBLUP), and single‐trait BayeCπ (ST‐BayesCπ) methods. We found that ribeye area (REA) and ribeye weight (REWT) showed high heritability, while slaughter weight (SWT) and carcass weight (CWT) displayed medium heritability, and slaughter rate (SR) and feedlot average daily gain (FDG) showed low heritability. Highly positive genetic correlations were observed between CWT and SWT (0.981) and SR and REWT (0.921). Notably, the MT‐BayesCπ method showed superior predictive abilities compared to other models. Using MT‐BayesCπ method, the accuracy increased from 0.272 to 0.694 for CWT compared to ST‐GBLUP and ST‐BayesCπ. MT‐GBLUP and ST‐BayesCπ showed similar prediction accuracies, while MT‐BayesCπ showed the least biased evaluations. Additionally, our results suggested that prediction accuracy of low‐heritability traits significantly increased when they were combined with traits with high genetic correlation in a multi‐trait prediction. Our study suggests that multi‐trait genomic predictions using GBLUP and Bayesian mixture prior models is feasible for genomic selection in beef cattle. Our findings indicate that MT‐BayesCπ outperforms other models (ST‐GBLUP, MT‐GBLUP and ST‐BayesCπ), especially for low‐heritability traits.

Abstract Background Genomic prediction of breeding values (GP) has been adopted in evolutionary genomic studies to uncover microevolutionary processes of wild populations or improve captive breeding strategies. While recent evolutionary studies applied GP with individual single nucleotide polymorphism (SNP), haplotype-based GP could outperform individual SNP predictions through better capturing the linkage disequilibrium (LD) between the SNP and quantitative trait loci (QTL). This study aimed to evaluate the accuracy and bias of haplotype-based GP of immunoglobulin (Ig) A (IgA), IgE, and IgG against Teladorsagia circumcincta in lambs of an unmanaged sheep population (Soay breed) based on Genomic Best Linear Unbiased Prediction (GBLUP) and five Bayesian [BayesA, BayesB, BayesCπ, Bayesian Lasso (BayesL), and BayesR] methods. Results The accuracy and bias of GPs using SNP, haplotypic pseudo-SNP from blocks with different LD thresholds (0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.00), or the combinations of pseudo-SNPs and non-LD clustered SNPs were obtained. Across methods and marker sets, higher ranges of genomic estimated breeding values (GEBV) accuracies were observed for IgA (0.20 to 0.49), followed by IgE (0.08 to 0.20) and IgG (0.05 to 0.14). Considering the methods evaluated, up to 8% gains in GP accuracy of IgG were achieved using pseudo-SNPs compared to SNPs. Up to 3% gain in GP accuracy for IgA was also obtained using the combinations of the pseudo-SNPs with non-clustered SNPs in comparison to fitting individual SNP. No improvement in GP accuracy of IgE was observed using haplotypic pseudo-SNPs or their combination with non-clustered SNPs compared to individual SNP. Bayesian methods outperformed GBLUP for all traits. Most scenarios yielded lower accuracies for all traits with an increased LD threshold. GP models using haplotypic pseudo-SNPs predicted less-biased GEBVs mainly for IgG. For this trait, lower bias was observed with higher LD thresholds, whereas no distinct trend was observed for other traits with changes in LD. Conclusions Haplotype information improves GP performance of anti-helminthic antibody traits of IgA and IgG compared to fitting individual SNP. The observed gains in the predictive performances indicate that haplotype-based methods could benefit GP of some traits in wild animal populations.

The increasing use of genomic selection (GS) in plant and animal breeding programs has led to the development of software that fits models based on unique scenarios. Accordingly, several R packages have been developed for GS. The lmekin function from the coxme R package was one of the first functions implemented in R to fit models with random family effects using the pedigree–based relationship matrix. The function allows the user to provide the covariance structures for the random effects; thus, the GBLUP model can be fitted. This fitting process consists of replacing, in the traditional BLUP model, the additive relationship matrix derived from a pedigree by the additive relationship matrix derived from markers. Thus, the objective of this study was to employ the lmekin function in the context of genomic prediction by comparing the results of this function with those obtained using five R packages for GS: rrBLUP, BGLR, sommer, lme4qtl, and lme4GS. The comparisons were performed considering the computational times and predicted values for a wheat dataset and simulated big data. In addition, we implemented a 5-fold cross-validation scheme through considering the values predicted by the lmekin function for the wheat dataset. The results indicated that the lmekin function was effective in predicting genomic breeding values considering multiple random effects and relatively small sample sizes. The rrBLUP package processed the fastest for the scenario with only one genetic random effect, and the high temporal efficiency of the sommer package was confirmed for the scenario with more than one genetic random effect. Differences in computational times occurred because of the different algorithms implemented in the packages to estimate the variance components.

Simple Summary: The objective of our study was to evaluate the predictive ability of a multi-trait genomic prediction model to estimate heritability and genetic correlations of traits such as 305-day milk yield, milk fat percentage, milk protein percentage, milk lactose percentage, and milk dry matter percentage in the Polish Holstein population. Results showed that strong accuracies for the predictions were achieved. The genetic relations with milk yield were negative, as expected, because increasing milk yield decreases the milk components percentage. In conclusion, multi-trait genomic prediction can be more beneficial than single-trait genomic prediction. The objective of our study was to evaluate the predictive ability of a multi-trait genomic prediction model that accounts for interactions between marker effects to estimate heritability and genetic correlations of traits including 305-day milk yield, milk fat percentage, milk protein percentage, milk lactose percentage, and milk dry matter percentage in the Polish Holstein Friesian cow population. For this aim, 14,742 SNP genotype records for 586 Polish Holstein Friesian dairy cows from Poland were used. Single-Trait-ssGBLUP (ST) and Multi-Trait-ssGBLUP (MT) methods were used for estimation. We examined 305-day milk yield (MY, kg), milk fat percentage (MF, %), milk protein percentage (MP, %), milk lactose percentage (ML, %), and milk dry matter percentage (MDM, %). The results showed that the highest marker effect rank correlation was found between milk fat percentage and milk dry matter. The weakest marker effect rank correlation was found between ML and all other traits. Obtained accuracies of this study were between 0.770 and 0.882, and 0.773 and 0.876 for MT and ST, respectively, which were acceptable values. All estimated bias values were positive, which is proof of underestimation. The highest heritability value was obtained for MP (0.3029) and the lowest heritability value was calculated for ML (0.2171). Estimated heritability values were low for milk yield and milk composition as expected. The strongest genetic correlation was estimated between MDM and MF (0.4990) and the weakest genetic correlation was estimated between MY and ML (0.001). The genetic relations with milk yield were negative and can be ignored as they were not significant. In conclusion, multi-trait genomic prediction can be more beneficial than single-trait genomic prediction. [ABSTRACT FROM AUTHOR]

Advances in DNA sequencing technologies allow the sequencing of whole genomes of thousands of individuals and provide several million single nucleotide polymorphisms (SNPs) per individual. These data combined with precise and high-throughput phenotyping enable genome-wide association studies (GWAS) and the identification of SNPs underlying traits with complex genetic architectures. The identified causal SNPs and estimated allelic effects could then be used for advanced marker-assisted selection (MAS) in breeding programs. But could such MAS compete with the broadly used genomic selection (GS)? This question is of particular interest for the lengthy tree breeding strategies. Here, with our new software "SNPscan breeder," we simulated a simple tree breeding program and compared the impact of different selection criteria on genetic gain and inbreeding. Further, we assessed different genetic architectures and different levels of kinship among individuals of the breeding population. Interestingly, apart from progeny testing, GS using gBLUP performed best under almost all simulated scenarios. MAS based on GWAS results outperformed GS only if the allelic effects were estimated in large populations (ca. 10,000 individuals) of unrelated individuals. Notably, GWAS using 3,000 extreme phenotypes performed as good as the use of 10,000 phenotypes. GS increased inbreeding and thus reduced genetic diversity more strongly compared to progeny testing and GWAS-based selection. We discuss the practical implications for tree breeding programs. In conclusion, our analyses further support the potential of GS for forest tree breeding and improvement, although MAS may gain relevance with decreasing sequencing costs in the future. [ABSTRACT FROM AUTHOR]

Introduction: Genomic selection (GS) has gained global importance due to its potential to accelerate genetic progress and improve the efficiency of breeding programs. Objectives of the research: In this research we proposed a method to improve the prediction accuracy of tested lines in new (untested) environments. Method-1: The new method trained the model with a modified response variable (a difference of response variables) that decreases the lack of a nonstationary distribution between the training and testing and improved the prediction accuracy. Comparing new and conventional method: We compared the prediction accuracy of the conventional genomic best linear unbiased prediction (GBLUP) model (M1) including (or not) genotype × environment interaction (GE) (M1_GE; M1_NO_GE) versus the proposed method (M2) on several data sets. Results and discussion: The gain in prediction accuracy of M2, versus M1_GE, M1_NO_GE in terms of Pearson's correlation was of at least 4.3%, while in terms of percentage of top-yielding lines captured when was selected the 10% (Best10) and 20% (Best20) of lines was at least of 19.5%, while in terms of Normalized Root Mean Squared Error (NRMSE) was of at least of 42.29%. [ABSTRACT FROM AUTHOR]

The objective of the present study was to evaluate the breeding value and accuracy of genomic estimated breeding values (GEBVs) of carcass traits in Jeju Black cattle (JBC) using Hanwoo steers and JBC as a reference population using the single‐trait animal model. Our research included genotype and phenotype information on 19,154 Hanwoo steers with 1097 JBC acting as the reference population. Likewise, the test population consisted of 418 genotyped JBC individuals with no phenotypic records for those carcass traits. For estimating the accuracy of GEBV, we divided the entire population into three groups. Hanwoo and JBC make up the first group; Hanwoo and JBC, who has both the genotype and phenotypic records, are referred to as the reference (training) population, and JBC, who lacks phenotypic information is referred to as the test (validation) population. The second group consists of the JBC (without phenotype) as the test population and Hanwoo as a reference population with phenotype and genotypic data. The only JBCs in the third group are those who have genotypic and phenotypic data on them as a reference population but no phenotypic data on them as a test population. The single‐trait animal model was used in all three groups for statistical purposes. The reference populations estimated heritabilities for carcass weight (CWT), eye muscle area (EMA), backfat thickness (BF), and marbling score (MS) as 0.30, 0.26, 0.26, and 0.34 for the Hanwoo steer and 0.42, 0.27, 0.26, and 0.48 for JBC. The average accuracy for carcass traits in Group 1 was 0.80 for the Hanwoo and JBC reference population compared with 0.73 for the JBC test population. Although the average accuracy for carcass traits in Group 2 was 0.80, it was 0.80 for the Hanwoo reference population and only 0.56 for the JBC test population. The average accuracy for the JBC reference and test populations was 0.68 and 0.50, respectively, when they were included in the accuracy comparison without the Hanwoo reference population. Groups 1 and 2 used Hanwoo as reference population, which led to a better average accuracy; however, Group 3 only used the JBC reference and test population, which led to a lower average accuracy. This might be due to the fact that Group 3 used a smaller reference size than the group that came before it and that the genetic makeup of the Hanwoo and JBC breeds differed. The GEBV accuracy for MS was higher than that of other traits across all three analysis groups, followed by CWT, EMA, and BF, which may be partially explained by the MS traits' higher heritability. This study suggests that in order to achieve more accuracy, a large reference population particular to a breed should be established. Therefore, to increase the accuracy of GEBV prediction and the genetic benefit from genomic selection in JBC, individual reference breeds, and large populations are required. [ABSTRACT FROM AUTHOR]

Simple Summary: Some of the compounds involved in sexual steroids' metabolic pathways (i.e., androstenone, indole, and skatole) might accumulate in the adipose tissue of intact male pigs after sexual maturity, resulting in boar taint (BT). With a perspective future ban on surgical castration, in pig population where early slaughtering is not a viable option, the exploitation of genomic selection procedures might prevent this off-odor and off-flavor. The accuracy provided by the use of genomic information was equal or higher than the one obtained using pedigree information. This indicates that genomic selection could be beneficial for the traits investigated as it minimizes the need to collect individual measures of BT compounds. Several genomic regions, each with a small effect on BT compound concentrations, were identified. Genes previously associated with BT, reproduction traits, and fat metabolism are located in those genomic regions. Detection of candidate genes related to fat metabolism might be explained by the relationship between sexual steroid levels and fat deposition and be ascribed to the pig line investigated, selected for ham quality and not for lean growth. With a perspective future ban on surgical castration in Europe, selecting pigs with reduced ability to accumulate boar taint (BT) compounds (androstenone, indole, skatole) in their tissues seems a promising strategy. BT compound concentrations were quantified in the adipose tissue of 1075 boars genotyped at 29,844 SNPs. Traditional and SNP-based breeding values were estimated using pedigree-based BLUP (PBLUP) and genomic BLUP (GBLUP), respectively. Heritabilities for BT compounds were moderate (0.30–0.52). The accuracies of GBLUP and PBLUP were significantly different for androstenone (0.58 and 0.36, respectively), but comparable for indole and skatole (~0.43 and ~0.47, respectively). Several SNP windows, each explaining a small percentage of the variance of BT compound concentrations, were identified in a genome-wide association study (GWAS). A total of 18 candidate genes previously associated with BT (MX1), reproduction traits (TCF21, NME5, PTGFR, KCNQ1, UMODL1), and fat metabolism (CTSD, SYT8, TNNI2, CD81, EGR1, GIPC2, MIGA1, NEGR1, CCSER1, MTMR2, LPL, ERFE) were identified in the post-GWAS analysis. The large number of genes related to fat metabolism might be explained by the relationship between sexual steroid levels and fat deposition and be partially ascribed to the pig line investigated, which is selected for ham quality and not for lean growth. [ABSTRACT FROM AUTHOR]

IntroductionGenomic selection (GS) experiments in forest trees have largely reported estimates of predictive abilities from cross-validation among individuals in the same breeding generation. In such conditions, no effects of recombination, selection, drift, and environmental changes are accounted for. Here, we assessed the effectively realized predictive ability (RPA) for volume growth at harvest age by GS across generations in an operational reciprocal recurrent selection (RRS) program of hybrid Eucalyptus.MethodsGenomic best linear unbiased prediction with additive (GBLUP_G), additive plus dominance (GBLUP_G+D), and additive single-step (HBLUP) models were trained with different combinations of growth data of hybrids and pure species individuals (N = 17,462) of the G1 generation, 1,944 of which were genotyped with ~16,000 SNPs from SNP arrays. The hybrid G2 progeny trial (HPT267) was the GS target, with 1,400 selection candidates, 197 of which were genotyped still at the seedling stage, and genomically predicted for their breeding and genotypic values at the operational harvest age (6 years). Seedlings were then grown to harvest and measured, and their pedigree-based breeding and genotypic values were compared to their originally predicted genomic counterparts.ResultsGenomic RPAs ≥0.80 were obtained as the genetic relatedness between G1 and G2 increased, especially when the direct parents of selection candidates were used in training. GBLUP_G+D reached RPAs ≥0.70 only when hybrid or pure species data of G1 were included in training. HBLUP was only marginally better than GBLUP. Correlations ≥0.80 were obtained between pedigree and genomic individual ranks. Rank coincidence of the top 2.5% selections was the highest for GBLUP_G (45% to 60%) compared to GBLUP_G+D. To advance the pure species RRS populations, GS models were best when trained on pure species than hybrid data, and HBLUP yielded ~20% higher predictive abilities than GBLUP, but was not better than ABLUP for ungenotyped trees.DiscussionWe demonstrate that genomic data effectively enable accurate ranking of eucalypt hybrid seedlings for their yet-to-be observed volume growth at harvest age. Our results support a two-stage GS approach involving family selection by average genomic breeding value, followed by within-top-families individual GS, significantly increasing selection intensity, optimizing genotyping costs, and accelerating RRS breeding.

Background: Genomic prediction of breeding values (GP) has been adopted in evolutionary genomic studies to uncover microevolutionary processes of wild populations or improve captive breeding strategies. While recent evolutionary studies applied GP with individual single nucleotide polymorphism (SNP), haplotype-based GP could outperform individual SNP predictions through better capturing the linkage disequilibrium (LD) between the SNP and quantitative trait loci (QTL). This study aimed to evaluate the accuracy and bias of haplotype-based GP of immunoglobulin (Ig) A (IgA), IgE, and IgG against Teladorsagia circumcincta in lambs of an unmanaged sheep population (Soay breed) based on Genomic Best Linear Unbiased Prediction (GBLUP) and five Bayesian [BayesA, BayesB, BayesCπ, Bayesian Lasso (BayesL), and BayesR] methods. Results: The accuracy and bias of GPs using SNP, haplotypic pseudo-SNP from blocks with different LD thresholds (0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, and 1.00), or the combinations of pseudo-SNPs and non-LD clustered SNPs were obtained. Across methods and marker sets, higher ranges of genomic estimated breeding values (GEBV) accuracies were observed for IgA (0.20 to 0.49), followed by IgE (0.08 to 0.20) and IgG (0.05 to 0.14). Considering the methods evaluated, up to 8% gains in GP accuracy of IgG were achieved using pseudo-SNPs compared to SNPs. Up to 3% gain in GP accuracy for IgA was also obtained using the combinations of the pseudo-SNPs with non-clustered SNPs in comparison to fitting individual SNP. No improvement in GP accuracy of IgE was observed using haplotypic pseudo-SNPs or their combination with non-clustered SNPs compared to individual SNP. Bayesian methods outperformed GBLUP for all traits. Most scenarios yielded lower accuracies for all traits with an increased LD threshold. GP models using haplotypic pseudo-SNPs predicted less-biased GEBVs mainly for IgG. For this trait, lower bias was observed with higher LD thresholds, whereas no distinct trend was observed for other traits with changes in LD. Conclusions: Haplotype information improves GP performance of anti-helminthic antibody traits of IgA and IgG compared to fitting individual SNP. The observed gains in the predictive performances indicate that haplotype-based methods could benefit GP of some traits in wild animal populations. [ABSTRACT FROM AUTHOR]

Abstract Genomic prediction (GP) based on haplotype alleles can capture quantitative trait loci (QTL) effects and increase predictive ability because the haplotypes are expected to be in linkage disequilibrium (LD) with QTL. In this study, we constructed haploblocks using LD‐based and the fixed number of single nucleotide polymorphisms (fixed‐SNP) methods with Illumina BovineHD chip in beef cattle. To evaluate the performance of different haplotype block partitioning methods, we constructed haploblocks based on LD thresholds (from r2 > 0.2 to r2 > 0.8) and the number of fixed‐SNPs (5, 10, 20). The performance of predictive methods for three carcass traits including liveweight (LW), dressing percentage (DP), and longissimus dorsi muscle weight (LDMW) was evaluated using three approaches (GBLUP and BayesB model based on the SNP, GHBLUP, and BayesBH models based on the haploblock, and GHBLUP+GBLUP and BayesBH+BayesB models based on the combined haploblock and the nonblocked SNPs, which were located between blocks). In this study, we found the accuracies of LD‐based and fixed‐SNP haplotype Bayesian methods outperformed the Bayesian models (up to 8.54 ± 7.44% and 5.74 ± 2.95%, respectively). GHBLUP showed a high improvement (up to 11.29 ± 9.87%) compared with GBLUP. The Bayesian models have higher accuracies than BLUP models in most scenarios. The average computing time of the BayesBH+BayesB model can reduce by 29.3% compared with the BayesB model. The prediction accuracies using the LD‐based haplotype method showed higher improvements than the fixed‐SNP haplotype method. In addition, to avoid the influence of rare haplotypes generated from haplotype construction, we compared the performance of GP by filtering four types of minor haplotype allele frequency (MHAF) (0.01, 0.025, 0.05, and 0.1) under different conditions (LD levels were set at r2 > 0.3, and the fixed number of SNPs was 5). We found the optimal MHAF threshold for LW was 0.01, and the optimal MHAF threshold for DP and LDMW was 0.025.

Abstract Western redcedar (WRC) is an ecologically and economically important forest tree species characterized by low genetic diversity with high self‐compatibility and high heartwood durability. Using sequence capture genotyping of target genic and non‐genic regions, we genotyped 44 parent trees and 1520 offspring trees representing 26 polycross (PX) families collected from three progeny test sites using 45,378 SNPs. Trees were phenotyped for eight traits related to growth, heartwood and foliar chemistry associated with wood durability and deer browse resistance. We used the genomic realized relationship matrix for paternity assignment, maternal pedigree correction, and to estimate genetic parameters. We compared genomics‐based (GBLUP) and two pedigree‐based (ABLUP: polycross and reconstructed full‐sib [FS] pedigrees) models. Models were extended to estimate dominance genetic effects. Pedigree reconstruction revealed significant unequal male contribution and separated the 26 PX families into 438 FS families. Traditional maternal PX pedigree analysis resulted in up to 51% overestimation in genetic gain and 44% in diversity. Genomic analysis resulted in up to 22% improvement in offspring breeding value (BV) theoretical accuracy, 35% increase in expected genetic gain for forward selection, and doubled selection intensity for backward selection. Overall, all traits showed low to moderate heritability (0.09–0.28), moderate genotype by environment interaction (type‐B genetic correlation: 0.51–0.80), low to high expected genetic gain (6.01%–55%), and no significant negative genetic correlation reflecting no large trade‐offs for multi‐trait selection. Only three traits showed a significant dominance effect. GBLUP resulted in smaller but more accurate heritability estimates for five traits, but larger estimates for the wood traits. Comparison between all, genic‐coding, genic‐non‐coding and intergenic SNPs showed little difference in genetic estimates. In summary, we show that GBLUP overcomes the PX limitations, successfully captures expected historical and hidden relatedness as well as linkage disequilibrium (LD), and results in increased breeding efficiency in WRC.

Simple Summary: The genomic estimated breeding value (GEBV) using data from Brangus heifers were obtained from genomic selection (GS) methods associating the single nucleotide polymorphisms (SNP) marker genotypes with phenotypic data for economically important growth (birth, weaning, and yearling weights) and carcass (depth of rib fat, and percent intramuscular fat and longissimus muscle area) traits using the linkage disequilibrium (LD) between SNP markers and quantitative trait loci (QTL) and/or the genomic relationship between animals. The heritability estimates were found similar across genomic best linear unbiased prediction (the GBLUP), and the Bayesian (BayesA, BayesB, BayesC and Lasso) GS methods for k-means and random cluster. The Bayesian methods resulted in underestimates of heritabilities and overestimates of accuracy of GEBV. However, the GBLUP method resulted in more reasonable estimates of heritabilities and accuracies of GEBV for growth and carcass traits of heifers from a composite population. The predictive abilities and accuracies of genomic best linear unbiased prediction (GBLUP) and the Bayesian (BayesA, BayesB, BayesC and Lasso) genomic selection (GS) methods for economically important growth (birth, weaning, and yearling weights) and carcass (depth of rib fat, apercent intramuscular fat and longissimus muscle area) traits were characterized by estimating the linkage disequilibrium (LD) structure in Brangus heifers using single nucleotide polymorphisms (SNP) markers. Sharp declines in LD were observed as distance among SNP markers increased. The application of the GBLUP and the Bayesian methods to obtain the GEBV for growth and carcass traits within k-means and random clusters showed that k-means and random clustering had quite similar heritability estimates, but the Bayesian methods resulted in the lower estimates of heritability between 0.06 and 0.21 for growth and carcass traits compared with those between 0.21 and 0.35 from the GBLUP methodologies. Although the prediction ability of the GBLUP and the Bayesian methods were quite similar for growth and carcass traits, the Bayesian methods overestimated the accuracies of GEBV because of the lower estimates of heritability of growth and carcass traits. However, GBLUP resulted in accuracy of GEBV for growth and carcass traits that parallels previous reports. [ABSTRACT FROM AUTHOR]

Hybrid breeding in sorghum [Sorghum bicolor (L.) Moench] utilizes the cytoplasmic-nuclear male sterility (CMS) system for seed production and subsequently harnesses heterosis. Since the cost of developing and evaluating inbred and hybrid lines in the CMS system is costly and time-consuming, genomic prediction of parental lines and hybrids is based on genetic data genotype. We generated 602 hybrids by crossing two female (A) lines with 301 diverse and elite male (R) lines from the sorghum association panel and collected phenotypic data for agronomic traits over two years. We genotyped the inbred parents using whole genome resequencing and used 2,687,342 high quality (minor allele frequency > 2%) single nucleotide polymorphisms for genomic prediction. For grain yield, the experimental hybrids exhibited an average mid-parent heterosis of 40%. Genomic best linear unbiased prediction (GBLUP) for hybrid performance yielded an average prediction accuracy of 0.76-0.93 under the prediction scenario where both parental lines in validation sets were included in the training sets (T2). However, when only female tester was shared between training and validation sets (T1F), prediction accuracies declined by 12-90%, with plant height showing the greatest decline. Mean accuracies for predicting the general combining ability of male parents ranged from 0.33 to 0.62 for all traits. Our results showed hybrid performance for agronomic traits can be predicted with high accuracy, and optimizing genomic relationship is essential for optimal training population design for genomic selection in sorghum breeding. [ABSTRACT FROM AUTHOR]

Introduction: Genomic selection is becoming a standard technique in plant breeding and is now being introduced into forest tree breeding. Despite promising results to predict the genetic merit of superior material based on their additive breeding values, many studies and operational programs still neglect non-additive effects and their potential for enhancing genetic gains. Methods: Using two large comprehensive datasets totaling 4,066 trees from 146 full-sib families of white spruce (Picea glauca (Moench) Voss), we evaluated the effect of the inclusion of dominance on the precision of genetic parameter estimates and on the accuracy of conventional pedigree-based (ABLUP-AD) and genomic-based (GBLUP-AD) models. Results: While wood quality traits were mostly additively inherited, considerable non-additive effects and lower heritabilities were detected for growth traits. For growth, GBLUP-AD better partitioned the additive and dominance effects into roughly equal variances, while ABLUP-AD strongly overestimated dominance. The predictive abilities of breeding and total genetic value estimates were similar between ABLUP-AD and GBLUP-AD when predicting individuals from the same families as those included in the training dataset. However, GBLUP-AD outperformed ABLUP-AD when predicting for new unphenotyped families that were not represented in the training dataset, with, on average, 22% and 53% higher predictive ability of breeding and genetic values, respectively. Resampling simulations showed that GBLUP-AD required smaller sample sizes than ABLUP-AD to produce precise estimates of genetic variances and accurate predictions of genetic values. Still, regardless of the method used, large training datasets were needed to estimate additive and non-additive genetic variances precisely. Discussion: This study highlights the different quantitative genetic architectures between growth and wood traits. Furthermore, the usefulness of genomic additive-dominance models for predicting new families should allow practicing mating allocation to maximize the total genetic values for the propagation of elite material. [ABSTRACT FROM AUTHOR]

Genomic selection is a promising approach for reducing the length of the selection cycle in forest tree breeding. Its efficiency must be evaluated across generations for this purpose, but such studies have been performed for multi-generational breeding programmes in only a few forest tree species to date. We analysed a subset of the Eucalyptus globulus breeding population from the Portuguese company Altri Florestal. In total, 412 genotypes from three successive breeding generations were genotyped with 14,716 SNP markers. A comparison of pedigree-based and marker-based relationship coefficients allowed to correct several documented pedigree errors. Deregressed breeding values were estimated from phenotypic records for growth traits (height and diameter) and survival for 31 field trials distributed in one breeding zone in Portugal and used as pseudo-phenotypes for genomic prediction models. Accuracy was assessed by cross-validation according to two main scenarios: (i) a scenario based on a five random fold number, not taking generation into account, and (ii) scenarios investigating progeny validation using parental generations to predict the progenies. Accuracy was highest after pedigree correction and ranged from 0.46 to 0.60 for the first scenario, from − 0.56 to 0.72 for parent/progeny scenarios and from 0.33 to 0.79 when progenies were added to the calibration population. This genomic selection study provides promising insight for the Altri Florestal Eucalyptus breeding programme. [ABSTRACT FROM AUTHOR]

IntroductionGenomic selection (GS) has gained global importance due to its potential to accelerate genetic progress and improve the efficiency of breeding programs.Objectives of the researchIn this research we proposed a method to improve the prediction accuracy of tested lines in new (untested) environments.Method-1The new method trained the model with a modified response variable (a difference of response variables) that decreases the lack of a non-stationary distribution between the training and testing and improved the prediction accuracy.Comparing new and conventional methodWe compared the prediction accuracy of the conventional genomic best linear unbiased prediction (GBLUP) model (M1) including (or not) genotype × environment interaction (GE) (M1_GE; M1_NO_GE) versus the proposed method (M2) on several data sets.Results and discussionThe gain in prediction accuracy of M2, versus M1_GE, M1_NO_GE in terms of Pearson´s correlation was of at least 4.3%, while in terms of percentage of top-yielding lines captured when was selected the 10% (Best10) and 20% (Best20) of lines was at least of 19.5%, while in terms of Normalized Root Mean Squared Error (NRMSE) was of at least of 42.29%.

This study aimed to predict the accuracy of genomic estimated breeding values (GEBVs) for reproductive traits in Hanwoo cows using the GBLUP, BayesB, BayesLASSO, and BayesR methods. Accuracy estimates of GEBVs for reproductive traits were derived through fivefold cross-validation, analyzing a dataset comprising 11,348 animals and employing an Illumina Bovine 50K SNP chip. GBLUP showed an accuracy of 0.26 for AFC, while BayesB, BayesLASSO, and BayesR demonstrated values of 0.28, 0.29, and 0.29, respectively. For CI, GBLUP attained an accuracy of 0.19, whereas BayesB, BayesLASSO, and BayesR scored 0.21, 0.24, and 0.25, respectively. The accuracy for GL was uniform across GBLUP, BayesB, and BayesR at 0.31, whereas BayesLASSO showed a slightly higher accuracy of 0.33. For NAIPC, GBLUP showed an accuracy of 0.24, while BayesB, BayesLASSO, and BayesR recorded 0.22, 0.27, and 0.30, respectively. The variation in genomic prediction accuracy among methods indicated Bayesian approaches slightly outperformed GBLUP. The findings suggest that Bayesian methods, notably BayesLASSO and BayesR, offer improved predictive capabilities for reproductive traits. Future research may explore more advanced genomic approaches to enhance predictive accuracy and genetic gains in Hanwoo cattle breeding programs.

Genomic prediction models generally assume that errors are distributed as normal, independent, and identically distributed random variables with zero mean and equal variance. This is not always true, in addition there may be phenotypes distant from the population mean, which are usually removed when making the prediction. Quantile regression (QR) deals with statistical aspects such as high dimensionality, multicollinearity and phenotypic distribution different from the normal one. The objective of this work was to compare QR using marker and pedigree information with alternative methods such as genomic best linear unbiased prediction (GBLUP) and single-step genomic best linear unbiased prediction (ssGBLUP) to analyze the birth (BW), weaning (WW) and yearling (YW) weights of Braunvieh cattle and simulated data with different degrees of asymmetry and proportion of outliers. The predictive capacity of the models was assessed by cross-validation. The predictive performance of QR both with marker information alone and with information of markers plus pedigree, with the actual dataset, was better than the GBLUP and ssGBLUP methodologies for WW and YW. For BW, GBLUP and ssGBLUP were better, however, only quantiles 0.25, 0.50 and 0.75 were used, and the BW distribution was not asymmetric. In the simulated data experiment, correlations between "true" marker effects and estimated effects, as well as "true" and estimated signal correlations were higher when QR was used compared to GBLUP. The advantages of QR were more noticeable with asymmetric distribution of phenotypes and with a higher proportion of outliers, as was the case with the simulated dataset. [ABSTRACT FROM AUTHOR]

GBLUP, the most widely used genomic prediction (GP) method, consumes large and increasing amounts of computational resources as the training population size increases due to the inverse of the genomic relationship matrix (GRM). Therefore, in this study, we developed a new genomic prediction method (RHEPCG) that avoids the direct inverse of the GRM by combining randomized Haseman-Elston (HE) regression (RHE-reg) and a preconditioned conjugate gradient (PCG). The simulation results demonstrate that RHEPCG, in most cases, not only achieves similar predictive accuracy with GBLUP but also significantly reduces computational time. As for the real data, RHEPCG shows similar or better predictive accuracy for seven traits of the Arabidopsis thaliana F2 population and four traits of the Sorghum bicolor RIL population compared with GBLUP. This indicates that RHEPCG is a practical alternative to GBLUP and has better computational efficiency. [ABSTRACT FROM AUTHOR]

Genomic prediction (GP) based on haplotype alleles can capture quantitative trait loci (QTL) effects and increase predictive ability because the haplotypes are expected to be in linkage disequilibrium (LD) with QTL. In this study, we constructed haploblocks using LD‐based and the fixed number of single nucleotide polymorphisms (fixed‐SNP) methods with Illumina BovineHD chip in beef cattle. To evaluate the performance of different haplotype block partitioning methods, we constructed haploblocks based on LD thresholds (from r2 > 0.2 to r2 > 0.8) and the number of fixed‐SNPs (5, 10, 20). The performance of predictive methods for three carcass traits including liveweight (LW), dressing percentage (DP), and longissimus dorsi muscle weight (LDMW) was evaluated using three approaches (GBLUP and BayesB model based on the SNP, GHBLUP, and BayesBH models based on the haploblock, and GHBLUP+GBLUP and BayesBH+BayesB models based on the combined haploblock and the nonblocked SNPs, which were located between blocks). In this study, we found the accuracies of LD‐based and fixed‐SNP haplotype Bayesian methods outperformed the Bayesian models (up to 8.54 ± 7.44% and 5.74 ± 2.95%, respectively). GHBLUP showed a high improvement (up to 11.29 ± 9.87%) compared with GBLUP. The Bayesian models have higher accuracies than BLUP models in most scenarios. The average computing time of the BayesBH+BayesB model can reduce by 29.3% compared with the BayesB model. The prediction accuracies using the LD‐based haplotype method showed higher improvements than the fixed‐SNP haplotype method. In addition, to avoid the influence of rare haplotypes generated from haplotype construction, we compared the performance of GP by filtering four types of minor haplotype allele frequency (MHAF) (0.01, 0.025, 0.05, and 0.1) under different conditions (LD levels were set at r2 > 0.3, and the fixed number of SNPs was 5). We found the optimal MHAF threshold for LW was 0.01, and the optimal MHAF threshold for DP and LDMW was 0.025. [ABSTRACT FROM AUTHOR]

Using data from genome-wide molecular markers, genomic selection procedures have proved useful for estimating breeding values and phenotypic prediction. The link between an individual genotype and phenotype has been modelled using a number of parametric methods to estimate individual breeding value. It has been observed that parametric methods perform satisfactorily only when the system under study has additive genetic architecture. To capture non-additive (dominance and epistasis) effects, nonparametric approaches have also been developed; however, they typically fall short of capturing additive effects. The idea behind this study is to select the most appropriate model from each parametric and nonparametric category and build an integrated model that can incorporate the best features of both models. It was observed from the results of the current study that GBLUP performed admirably under additive architecture, while SVM's performance in non-additive architecture was found to be encouraging. A robust model for genomic prediction has been developed in light of these findings, which can handle both additive and epistatic effects simultaneously by minimizing their error variance. The developed integrated model has been assessed using standard evaluation measures like predictive ability and error variance. [ABSTRACT FROM AUTHOR]

The impact of genomic epistasis effects on the accuracy of predicting the phenotypic values of residual feed intake (RFI) in U.S. Holstein cows was evaluated using 6215 Holstein cows and 78,964 SNPs. Two SNP models and seven epistasis models were initially evaluated. Heritability estimates and the accuracy of predicting the RFI phenotypic values from 10-fold cross-validation studies identified the model with SNP additive effects and additive x additive (AxA) epistasis effects (A + AxA model) to be the best prediction model. Under the A + AxA model, additive heritability was 0.141, and AxA heritability was 0.263 that consisted of 0.260 inter-chromosome AxA heritability and 0.003 intra-chromosome AxA heritability, showing that inter-chromosome AxA effects were responsible for the accuracy increases due to AxA. Under the SNP additivemodel (A-onlymodel), the additive heritability was 0.171. In the 10 validation populations, the average accuracy for predicting the RFI phenotypic values was 0.246 (with range 0.197--0.333) under A + AxA model and was 0.231 (with range of 0.188--0.319) under the A-only model. The average increase in the accuracy of predicting the RFI phenotypic values by the A + AxA model over the A-only model was 6.49% (with range of 3.02--14.29%). Results in this study showed AxA epistasis effects had a positive impact on the accuracy of predicting the RFI phenotypic values when combined with additive effects in the prediction model. [ABSTRACT FROM AUTHOR]

Computational efficiency has become a key issue in genomic prediction (GP) owing to the massive historical datasets accumulated. We developed hereby a new super-fast GP approach (SHEAPY) combining randomized Haseman-Elston regression (RHE-reg) with a modified Algorithm for Proven and Young (APY) in an additive-effect model, using the former to estimate heritability and then the latter to invert a large genomic relationship matrix for best linear prediction. In simulation results with varied sizes of training population, GBLUP, HEAPY|A and SHEAPY showed similar predictive performance when the size of a core population was half that of a large training population and the heritability was a fixed value, and the computational speed of SHEAPY was faster than that of GBLUP and HEAPY|A. In simulation results with varied heritability, SHEAPY showed better predictive ability than GBLUP in all cases and than HEAPY|A in most cases when the size of a core population was 4/5 that of a small training population and the training population size was a fixed value. As a proof of concept, SHEAPY was applied to the analysis of two real datasets. In an Arabidopsis thaliana F2 population, the predictive performance of SHEAPY was similar to or better than that of GBLUP and HEAPY|A in most cases when the size of a core population (200) was 2/3 of that of a small training population (300). In a sorghum multiparental population, SHEAPY showed higher predictive accuracy than HEAPY|A for all of three traits, and than GBLUP for two traits. SHEAPY may become the GP method of choice for large-scale genomic data.

Genomic prediction models have proven efficient for increasing genetic gains of growth and fillet traits in aquaculture. However, measurement of these traits is expensive, time-consuming, or cannot be performed on live animals. This study investigated the potential of using correlated body traits as an alternative basis for supporting genomic prediction of average daily gain (ADG) and fillet weight (FW) in Asian seabass. A total of 954 fish were genotyped for 20,430 single nucleotide polymorphisms (SNPs). Phenotypic data at 390 days, including body weight, fillet weight, 16 size and four body-part-area measurements were used. Heritability estimates were high for ADG (0.47) and FW (0.44), low to moderate (0.12–0.39) for body sizes and body areas. Genetic correlations among pairs of traits were all high and positive (0.65–0.99). Five prediction models (1, 2, 3, 4 and 5) were evaluated using the GBLUP method and were compared using five-fold cross validation. Results indicated that model 4, which fitted the secondary traits derived from correlated body sizes and body areas as the phenotypes, was the best model, showing prediction accuracy of 0.47 with small prediction bias (0.95). The results demonstrated the advantage of using correlated traits for practical genomic prediction in Asian seabass, particularly for growth and fillet traits that are difficult to measure. • Eight size and two body area traits were used as secondary traits for genomic prediction of ADG and FW. • Genomic prediction was performed on 954 fish genotyped at 20,430 SNPs using five GBLUP models. • Model 4, was the best model, showing prediction accuracy of 0.47 with small prediction bias (0.95). [ABSTRACT FROM AUTHOR]

The objective of our study was to evaluate the predictive ability of a multi-trait genomic prediction model that accounts for interactions between marker effects to estimate heritability and genetic correlations of traits including 305-day milk yield, milk fat percentage, milk protein percentage, milk lactose percentage, and milk dry matter percentage in the Polish Holstein Friesian cow population. For this aim, 14,742 SNP genotype records for 586 Polish Holstein Friesian dairy cows from Poland were used. Single-Trait-ssGBLUP (ST) and Multi-Trait-ssGBLUP (MT) methods were used for estimation. We examined 305-day milk yield (MY, kg), milk fat percentage (MF, %), milk protein percentage (MP, %), milk lactose percentage (ML, %), and milk dry matter percentage (MDM, %). The results showed that the highest marker effect rank correlation was found between milk fat percentage and milk dry matter. The weakest marker effect rank correlation was found between ML and all other traits. Obtained accuracies of this study were between 0.770 and 0.882, and 0.773 and 0.876 for MT and ST, respectively, which were acceptable values. All estimated bias values were positive, which is proof of underestimation. The highest heritability value was obtained for MP (0.3029) and the lowest heritability value was calculated for ML (0.2171). Estimated heritability values were low for milk yield and milk composition as expected. The strongest genetic correlation was estimated between MDM and MF (0.4990) and the weakest genetic correlation was estimated between MY and ML (0.001). The genetic relations with milk yield were negative and can be ignored as they were not significant. In conclusion, multi-trait genomic prediction can be more beneficial than single-trait genomic prediction.

Los modelos de predicción genómica generalmente suponen que los errores se distribuyen como variables aleatorias normales, independientes e idénticamente distribuidas con media cero e igual varianza. Esto no siempre se cumple, además puede haber fenotipos distantes de la media poblacional, los que usualmente se eliminan al realizar la predicción. La regresión cuantil (QR) afronta aspectos estadísticos como alta dimensionalidad, multicolinealidad y distribución fenotípica diferente de la normal. El objetivo de este trabajo fue comparar QR utilizando información de marcadores y pedigrí con los métodos alternativos tales como mejor predicción lineal insesgada genómica (GBLUP) y mejor predicción lineal insesgada genómica en un solo paso (ssGBLUP) para analizar los pesos al nacimiento (PN), destete (PD) y al año (PA) de bovinos Suizo Europeo y datos simulados con diferente grado de asimetría y proporción de datos atípicos. La capacidad predictiva de los modelos se evaluó mediante validación cruzada. El desempeño predictivo de QR tanto sólo con información de marcadores como con marcadores más pedigrí, con el conjunto de datos reales, fue mejor que las metodologías GBLUP y ssGBLUP para PD y PA. Para PN GBLUP y ssGBLUP fueron mejores, sin embargo, solo se utilizaron los cuantiles 0.25, 0.50 y 0.75, y la distribución de PN no fue asimétrica. En el experimento de datos simulados, las correlaciones entre efectos de marcador “verdadero” y efectos estimados, así como las correlaciones de señales “verdaderas” y estimadas fueron más altas cuando se usó QR comparado con GBLUP. Las ventajas de QR fueron más notorias con distribución asimétrica de los fenotipos y con mayor proporción de datos atípicos, como fue el caso del conjunto de datos simulados.

Using data from genome-wide molecular markers, genomic selection procedures have proved useful for estimating breeding values and phenotypic prediction. The link between an individual genotype and phenotype has been modelled using a number of parametric methods to estimate individual breeding value. It has been observed that parametric methods perform satisfactorily only when the system under study has additive genetic architecture. To capture non-additive (dominance and epistasis) effects, nonparametric approaches have also been developed; however, they typically fall short of capturing additive effects. The idea behind this study is to select the most appropriate model from each parametric and nonparametric category and build an integrated model that can incorporate the best features of both models. It was observed from the results of the current study that GBLUP performed admirably under additive architecture, while SVM’s performance in non-additive architecture was found to be encouraging. A robust model for genomic prediction has been developed in light of these findings, which can handle both additive and epistatic effects simultaneously by minimizing their error variance. The developed integrated model has been assessed using standard evaluation measures like predictive ability and error variance.