Multi-trait Improvement by Predicting Genetic Correlations in Breeding Crosses

Prediction of additive genetic variances of descendants for complex families based on Mendelian sampling variances

The ability to predict the outcome of selection and mating decisions enables breeders to make strategically better selection decisions. To improve genetic progress, those individuals need to be selected whose offspring can be expected to show high genetic variance next to high breeding values. Previously published approaches enable to predict the variance of descendants of 2 future generations for up to 4 founding haplotypes, or 2 outbred individuals, based on phased genotypes, allele effects, and recombination frequencies. The purpose of this study was to develop a general approach for the analytical calculation of the genetic variance in any future generation. The core development is an equation for the prediction of the variance of double haploid lines, under the assumption of no selection and negligible drift, stemming from an arbitrary number of founder haplotypes. This double haploid variance can be decomposed into gametic Mendelian sampling variances (MSVs) of ancestors of the double haploid lines allowing usage for non-double haploid genotypes that enables application in animal breeding programs as well as in plant breeding programs. Together with the breeding values of the founders, the gametic MSV may be used in new selection criteria. We present our idea of such a criterion that describes the genetic level of selected individuals in 4 generations. Since breeding programs do select, the assumption made for predicting variances is clearly violated, which decreases the accuracy of predicted gametic MSV caused by changes in allele frequency and linkage disequilibrium. Despite violating the assumption, we found high predictive correlations of our criterion to the true genetic level that was obtained by means of simulation for the "corn" and "cattle" genome models tested in this study (0.90 and 0.97). In practice, the genotype phases, genetic map, and allele effects all need to be estimated meaning inaccuracies in their estimation will lead to inaccurate variance prediction. Investigation of variance prediction accuracy when input parameters are estimated was not part of this study.

Validation of cross-progeny variance genomic prediction using simulations and experimental data in winter elite bread wheat

KEY MESSAGE

From simulations and experimental data, the quality of cross progeny variance genomic predictions may be high, but depends on trait architecture and necessitates sufficient number of progenies. Genomic predictions are used to select genitors and crosses in plant breeding. The usefulness criterion (UC) is a cross-selection criterion that necessitates the estimation of parental mean (PM) and progeny standard deviation (SD). This study evaluates the parameters that affect the predictive ability of UC and its two components using simulations. Predictive ability increased with heritability and progeny size and decreased with QTL number, most notably for SD. Comparing scenarios where marker effects were known or estimated using prediction models, SD was strongly impacted by the quality of marker effect estimates. We proposed a new algebraic formula for SD estimation that takes into account the uncertainty of the estimation of marker effects. It improved predictions when the number of QTL was superior to 300, especially when heritability was low. We also compared estimated and observed UC using experimental data for heading date, plant height, grain protein content and yield. PM and UC estimates were significantly correlated for all traits (PM: 0.38, 0.63, 0.51 and 0.91; UC: 0.45, 0.52, 0.54 and 0.74; for yield, grain protein content, plant height and heading date, respectively), while SD was correlated only for heading date and plant height (0.64 and 0.49, respectively). According to simulations, SD estimations in the field would necessitate large progenies. This pioneering study experimentally validates genomic prediction of UC but the predictive ability depends on trait architecture and precision of marker effect estimates. We advise the breeders to adjust progeny size to realize the SD potential of a cross.

Genomic predictions of genetic variances and correlations among traits for breeding crosses in soybean

Parental selection is perhaps the most critical decision a breeder makes, establishing the foundation of the entire program for years to come. Cross selection based on predicted mean and genetic variance can be further expanded to multiple-trait improvement by predicting the genetic correlation (r G ) between pairs of traits. Our objective was to empirically assess the ability to predict the family mean, genetic variance, superior progeny mean and genetic correlation through genomic prediction in a soybean population. Data made available through the Soybean Nested Association Mapping project included phenotypic data on seven traits (days to maturity, lodging, oil, plant height, protein, seed size, and seed yield) for 39 families. Training population composition followed a leave-one-family-out cross-validation scheme, with the validation family genetic parameters predicted using the remaining families as the training set. The predictive abilities for family mean and superior progeny mean were significant for all traits while predictive ability of genetic variance was significant for four traits. We were able to validate significant predictive abilities ofr G for 18 out of 21 (86%) pairwise trait combinations (P < 0.05). The findings from this study support the use of genome-wide marker effects for predictingr G in soybean biparental crosses. If successfully implemented in breeding programs, this methodology could help to increase the rate of genetic gain for multiple correlated traits.

Development and optimization of expected cross value for mate selection problems

In this study, we address the mate selection problem in the hybridization stage of a breeding pipeline, which constitutes the multi-objective breeding goal key to the performance of a variety development program. The solution framework we formulate seeks to ensure that individuals with the most desirable genomic characteristics are selected to cross in order to maximize the likelihood of the inheritance of desirable genetic materials to the progeny. Unlike approaches that use phenotypic values for parental selection and evaluate individuals separately, we use a criterion that relies on the genetic architecture of traits and evaluates combinations of genomic information of the pairs of individuals. We introduce the expected cross value (ECV) criterion that measures the expected number of desirable alleles for gametes produced by pairs of individuals sampled from a population of potential parents. We use the ECV criterion to develop an integer linear programming formulation for the parental selection problem. The formulation is capable of controlling the inbreeding level between selected mates. We evaluate the approach or two applications: (i) improving multiple target traits simultaneously, and (ii) finding a multi-parental solution to design crossing blocks. We evaluate the performance of the ECV criterion using a simulation study. Finally, we discuss how the ECV criterion and the proposed integer linear programming techniques can be applied to improve breeding efficiency while maintaining genetic diversity in a breeding program.

Improving selection decisions with mating information by accounting for Mendelian sampling variances looking two generations ahead

BACKGROUND

Breeding programs are judged by the genetic level of animals that are used to disseminate genetic progress. These animals are typically the best ones of the population. To maximise the genetic level of very good animals in the next generation, parents that are more likely to produce top performing offspring need to be selected. The ability of individuals to produce high-performing progeny differs because of differences in their breeding values and gametic variances. Differences in gametic variances among individuals are caused by differences in heterozygosity and linkage. The use of the gametic Mendelian sampling variance has been proposed before, for use in the usefulness criterion or Index5, and in this work, we extend existing approaches by not only considering the gametic Mendelian sampling variance of individuals, but also of their potential offspring. Thus, the criteria developed in this study plan one additional generation ahead. For simplicity, we assumed that the true quantitative trait loci (QTL) effects, genetic map and the haplotypes of all animals are known.

RESULTS

In this study, we propose a new selection criterion, ExpBVSelGrOff, which describes the genetic level of selected grand-offspring that are produced by selected offspring of a particular mating. We compare our criterion with other published criteria in a stochastic simulation of an ongoing breeding program for 21 generations for proof of concept. ExpBVSelGrOff performed better than all other tested criteria, like the usefulness criterion or Index5 which have been proposed in the literature, without compromising short-term gains. After only five generations, when selection is strong (1%), selection based on ExpBVSelGrOff achieved 5.8% more commercial genetic gain and retained 25% more genetic variance without compromising inbreeding rate compared to selection based only on breeding values.

CONCLUSIONS

Our proposed selection criterion offers a new tool to accelerate genetic progress for contemporary genomic breeding programs. It retains more genetic variance than previously published criteria that plan less far ahead. Considering future gametic Mendelian sampling variances in the selection process also seems promising for maintaining more genetic variance.

Death and Dying: Grapevine Survival, Cold Hardiness, and BLUPs and Winter BLUEs in North Dakota Vineyards

A total of fourteen diverse, interspecific hybrid grapevines (Vitis spp.) were evaluated for their adaptability to North Dakota winter conditions using differential thermal analysis (DTA) of low-temperature exotherms (LTE) and bud cross-sectional assessment of survival techniques. This research was conducted in two vineyard locations in eastern North Dakota. This work demonstrates the use of DTA for monitoring and selecting cultivars capable of withstanding sub-zero temperatures. These results were assessed for quantitative genetic traits. High heritability was observed for bud LTE traits and may thus be a useful target for cold hardiness breeding programs; however, it is necessary to ensure that variance is reduced when pooling multiple sample events. After DTA sampling, grapevines were assessed for survival of primary and secondary dormant buds using cross-sectional visual evaluation of death. 'Valiant' had the greatest primary bud survival (68%), followed by 'Frontenac gris', 'Crimson Pearl', and 'King of the North'. These varieties are among those with potential for production in eastern North Dakota's environment. The newly evaluated relationships between traits and the heritability of DTA results provide valuable tools to grapevine breeders for the development of cold-tolerant genotypes for future climatic challenges.

Genotype by environment and genotype by yield*trait interactions in sugar beet: analyzing yield stability and determining key traits association

The genotype by environment interaction significantly influences plant yield, making it imperative to understand its nature for the creation of breeding programs to enhance crop production. However, this is not the only obstacle in the yield improvement process. Breeders also face the significant challenge of unfavorable and negative correlations among key traits. In this study, the stability of root yield and white sugar yield, and the association between the key traits of root yield, sugar content, nitrogen, sodium, and potassium were examined in 20 sugar beet genotypes. The study was conducted using a randomized complete block design with four replications over two consecutive years across five locations. The combined analysis of variance results revealed significant main effects of year, location, and genotype on both root yield and white sugar yield. Notably, two-way and three-way interactions between these main effects on root yield and white sugar yield resulted in a significant difference. The additive main effect and multiplicative interaction analysis revealed that the first five interaction principal components significantly impacted both the root yield and white sugar yield. The linear mixed model results for root yield and white sugar yield indicated that the genotype effect and the genotype by environment interaction were significant. The weighted average absolute scores of the best linear unbiased predictions biplot demonstrated that genotypes 20, 4, 7, 2, 16, 3, 6, 1, 14, and 15 were superior in terms of root yield. For white sugar yield, genotypes 4, 16, 3, 7, 5, 1, 10, 20, 2, and 6 stood out. These genotypes were not only stable but also had a yield value higher than the total average. All key traits, which include sugar content, sodium, potassium, and alpha amino nitrogen, demonstrated a negative correlation with root yield. Based on the genotype by yield*trait analysis results, genotypes 20, 19, and 16 demonstrated optimal performance when considering the combination of root yield with sugar content, sodium, alpha amino nitrogen, and potassium. The multi-trait stability study, genotype 13 ranked first, and genotypes 10, 8, and 9 were identified as the most ideal stable genotypes across all traits. According to the multi-trait stability index, genotype 13 emerged as the top-ranking genotype. Additionally, genotypes 10, 8, and 9 were recognized as the most stable genotypes.

Variability of Glucosinolates in Pak Choy (Brassica rapa subsp. chinensis) Germplasm

Glucosinolates are sulfur-containing phytochemicals generally abundant in cruciferous vegetables such as pak choy. Glucosinolates participate in a range of biological activities essential for promoting a healthy human body. In this study, we aimed to elucidate glucosinolate variability present in pak choy germplasm that are under conservation at the Rural Development Administration Genebank, Jeonju, Republic of Korea. The Acquity Ultra-Performance Liquid Chromatography (UHPLC) analytical system was used in profiling the glucosinolate content in leaf samples of various accessions. We identified a total of 17 glucosinolates in the germplasm. Based on principal compoment analysis performed, three separate groups of the accessions were obtained. Group 1 contained the cultivar cheongsacholong which recorded high content of glucobrassicin (an indole), glucoerucin (aliphatic), gluconasturtiin (aromatic) and glucoberteroin (aliphatic). Group 2 consisted of six accessions, BRA77/72, Lu ling gaogengbai, 9041, Wuyueman, RP-75 and DH-10, predominatly high in aliphatic compounds including glucoiberin, glucocheirolin, and sinigrin. Group 3 comprised the majority of the accessions which were characterized by high content of glucoraphanin, epiprogoitrin, progoitrin, and glucotropaeolin. These results revealed the presence of variability among the pak choy germplasm based on their glucosinolate content, providing an excellent opprtunity for future breeding for improved glucosinolate content in the crop.

Genetic Variance Estimation over Time in Broiler Breeding Programmes for Growth and Reproductive Traits

Monitoring the genetic variance of traits is a key priority to ensure the sustainability of breeding programmes in populations under directional selection, since directional selection can decrease genetic variation over time. Studies monitoring changes in genetic variation have typically used long-term data from small experimental populations selected for a handful of traits. Here, we used a large dataset from a commercial breeding line spread over a period of twenty-three years. A total of 2,059,869 records and 2,062,112 animals in the pedigree were used for the estimations of variance components for the traits: body weight (BWT; 2,059,869 records) and hen-housed egg production (HHP; 45,939 records). Data were analysed with three estimation approaches: sliding overlapping windows, under frequentist (restricted maximum likelihood (REML)) and Bayesian (Gibbs sampling) methods; expected variances using coefficients of the full relationship matrix; and a "double trait covariances" analysis by computing correlations and covariances between the same trait in two distinct consecutive windows. The genetic variance showed marginal fluctuations in its estimation over time. Whereas genetic, maternal permanent environmental, and residual variances were similar for BWT in both the REML and Gibbs methods, variance components when using the Gibbs method for HHP were smaller than the variances estimated when using REML. Large data amounts were needed to estimate variance components and detect their changes. For Gibbs (REML), the changes in genetic variance from 1999-2001 to 2020-2022 were 82.29 to 93.75 (82.84 to 93.68) for BWT and 76.68 to 95.67 (98.42 to 109.04) for HHP. Heritability presented a similar pattern as the genetic variance estimation, changing from 0.32 to 0.36 (0.32 to 0.36) for BWT and 0.16 to 0.15 (0.21 to 0.18) for HHP. On the whole, genetic parameters tended slightly to increase over time. The expected variance estimates were lower than the estimates when using overlapping windows. That indicates the low effect of the drift-selection process on the genetic variance, or likely, the presence of genetic variation sources compensating for the loss. Double trait covariance analysis confirmed the maintenance of variances over time, presenting genetic correlations >0.86 for BWT and >0.82 for HHP. Monitoring genetic variance in broiler breeding programmes is important to sustain genetic progress. Although the genetic variances of both traits fluctuated over time, in some windows, particularly between 2003 and 2020, increasing trends were observed, which warrants further research on the impact of other factors, such as novel mutations, operating on the dynamics of genetic variance.

Genomic prediction of optimal cross combinations to accelerate genetic improvement of soybean (Glycine max)

Improving yield is a primary soybean breeding goal, as yield is the main determinant of soybean's profitability. Within the breeding process, selection of cross combinations is one of most important elements. Cross prediction will assist soybean breeders in identifying the best cross combinations among parental genotypes prior to crossing, increasing genetic gain and breeding efficiency. In this study optimal cross selection methods were created and applied in soybean and validated using historical data from the University of Georgia soybean breeding program, under multiple training set compositions and marker densities utilizing multiple genomic selection models for marker evaluation. Plant materials consisted of 702 advanced breeding lines evaluated in multiple environments and genotyped using SoySNP6k BeadChips. An additional marker set, the SoySNP3k marker set, was tested in this study as well. Optimal cross selection methods were used to predict the yield of 42 previously made crosses and compared to the performance of the cross's offspring in replicated field trials. The best prediction accuracy was obtained when using Extended Genomic BLUP with the SoySNP6k marker set, consisting of 3,762 polymorphic markers, with an accuracy of 0.56 with a training set maximally related to the crosses predicted and 0.4 in a training set with minimized relatedness to predicted crosses. Prediction accuracy was most significantly impacted by training set relatedness to the predicted crosses, marker density, and the genomic model used to predict marker effects. The usefulness criterion selected had an impact on prediction accuracy within training sets with low relatedness to the crosses predicted. Optimal cross prediction provides a useful method that assists plant breeders in selecting crosses in soybean breeding.

Malting quality and preharvest sprouting traits are genetically correlated in spring malting barley

Malt for craft "all-malt" brewing can have high quality, PHS resistance, and malted in normal timeframes. Canadian style adjunct malt is associated with PHS susceptibility. Expansion of malting barley production into non-traditional growing regions and erratic weather has increased the demand for preharvest sprouting (PHS) resistant, high quality malting barley cultivars. This is hindered by the relatively unknown relationships between PHS resistance and malting quality. Here we present a three-year study of malting quality and germination at different after-ripening durations post physiological maturity. Malting quality traits alpha amylase (AA) and free amino nitrogen (FAN) and germination rate at six days post PM shared a common association with a SNP in HvMKK3 on chromosome 5H in the Seed Dormancy 2 (SD2) region responsible for PHS susceptibility. Soluble protein (SP) and soluble over total protein (S/T) both shared a common association with a marker in the SD2 region. Significant genetic correlations between PHS resistance and the malting quality traits AA, FAN, SP, S/T were detected across and within HvMKK3 allele groups. High adjunct malt quality was related to PHS susceptibility. Selection for PHS resistance led to a correlated response in malting quality traits. Results strongly suggest pleiotropy of HvMKK3 on malting quality traits and that the classic "Canadian-style" malt is caused by a PHS susceptible allele of HvMKK3. PHS susceptibility appears to benefit the production of malt intended for adjunct brewing, while PHS resistance is compatible with all-malt brewing specifications. Here we present our analysis on the effect of combining complexly inherited and correlated traits with contrasting goals to inform breeding practice in malting barley, the general principles of which can be extended to other breeding programs.

Accuracy of Selection in Early Generations of Field Pea Breeding Increases by Exploiting the Information Contained in Correlated Traits

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 (h2) 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) (h2 = 0.05), compressed stem thickness (CST) (h2 = 0.12), internode length (IL) (h2 = 0.61) and angle of the main stem above horizontal at first flower (EAngle) (h2 = 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.

Multi-trait genomic prediction improves selection accuracy for enhancing seed mineral concentrations in pea

Multi-trait genomic selection (MT-GS) has the potential to improve predictive ability by maximizing the use of information across related genotypes and genetically correlated traits. In this study, we extended the use of sparse phenotyping method into the MT-GS framework by split testing of entries to maximize borrowing of information across genotypes and predict missing phenotypes for targeted traits without additional phenotyping expenditure. Using 300 advanced breeding lines from North Dakota State University (NDSU) pulse breeding program and ∼200 USDA accessions that were evaluated for 10 nutritional traits, our results show that the proposed sparse phenotyping aided MT-GS can further improve predictive ability by >12% across traits compared with univariate (UNI) genomic selection. The proposed strategy departed from the previous reports that weak genetic correlation is a limitation to the advantage of MT-GS over UNI genomic selection, which was evident in the partially balanced phenotyping-enabled MT-GS. Our results point to heritability and genetic correlation between traits as possible metrics to optimize and further improve the estimation of model parameters, and ultimately, prediction performance. Overall, our study offers a new approach to optimize the prediction performance using the MT-GS and further highlight strategy to maximize the efficiency of GS in a plant breeding program. The sparse-testing-aided MT-GS proposed in this study can be further extended to multi-environment, multi-trait GS to improve prediction performance and further reduce the cost of phenotyping and time-consuming data collection process.

Correlation scan: identifying genomic regions that affect genetic correlations applied to fertility traits

Although the genetic correlations between complex traits have been estimated for more than a century, only recently we have started to map and understand the precise localization of the genomic region(s) that underpin these correlations. Reproductive traits are often genetically correlated. Yet, we don't fully understand the complexities, synergism, or trade-offs between male and female fertility. In this study, we used reproductive traits in two cattle populations (Brahman; BB, Tropical Composite; TC) to develop a novel framework termed correlation scan (CS). This framework was used to identify local regions associated with the genetic correlations between male and female fertility traits. Animals were genotyped with bovine high-density single nucleotide polymorphisms (SNPs) chip assay. The data used consisted of ~1000 individual records measured through frequent ovarian scanning for age at first corpus luteum (AGECL) and a laboratory assay for serum levels of insulin growth hormone (IGF1 measured in bulls, IGF1b, or cows, IGF1c). The methodology developed herein used correlations of 500-SNP effects in a 100-SNPs sliding window in each chromosome to identify local genomic regions that either drive or antagonize the genetic correlations between traits. We used Fisher's Z-statistics through a permutation method to confirm which regions of the genome harboured significant correlations. About 30% of the total genomic regions were identified as driving and antagonizing genetic correlations between male and female fertility traits in the two populations. These regions confirmed the polygenic nature of the traits being studied and pointed to genes of interest. For BB, the most important chromosome in terms of local regions is often located on bovine chromosome (BTA) 14. However, the important regions are spread across few different BTA's in TC. Quantitative trait loci (QTLs) and functional enrichment analysis revealed many significant windows co-localized with known QTLs related to milk production and fertility traits, especially puberty. In general, the enriched reproductive QTLs driving the genetic correlations between male and female fertility are the same for both cattle populations, while the antagonizing regions were population specific. Moreover, most of the antagonizing regions were mapped to chromosome X. These results suggest regions of chromosome X for further investigation into the trade-offs between male and female fertility. We compared the CS with two other recently proposed methods that map local genomic correlations. Some genomic regions were significant across methods. Yet, many significant regions identified with the CS were overlooked by other methods.

High-throughput characterization, correlation, and mapping of leaf photosynthetic and functional traits in the soybean (Glycine max) nested association mapping population

Photosynthesis is a key target to improve crop production in many species including soybean [Glycine max (L.) Merr.]. A challenge is that phenotyping photosynthetic traits by traditional approaches is slow and destructive. There is proof-of-concept for leaf hyperspectral reflectance as a rapid method to model photosynthetic traits. However, the crucial step of demonstrating that hyperspectral approaches can be used to advance understanding of the genetic architecture of photosynthetic traits is untested. To address this challenge, we used full-range (500-2,400 nm) leaf reflectance spectroscopy to build partial least squares regression models to estimate leaf traits, including the rate-limiting processes of photosynthesis, maximum Rubisco carboxylation rate, and maximum electron transport. In total, 11 models were produced from a diverse population of soybean sampled over multiple field seasons to estimate photosynthetic parameters, chlorophyll content, leaf carbon and leaf nitrogen percentage, and specific leaf area (with R2 from 0.56 to 0.96 and root mean square error approximately <10% of the range of calibration data). We explore the utility of these models by applying them to the soybean nested association mapping population, which showed variability in photosynthetic and leaf traits. Genetic mapping provided insights into the underlying genetic architecture of photosynthetic traits and potential improvement in soybean. Notably, the maximum Rubisco carboxylation rate mapped to a region of chromosome 19 containing genes encoding multiple small subunits of Rubisco. We also mapped the maximum electron transport rate to a region of chromosome 10 containing a fructose 1,6-bisphosphatase gene, encoding an important enzyme in the regeneration of ribulose 1,5-bisphosphate and the sucrose biosynthetic pathway. The estimated rate-limiting steps of photosynthesis were low or negatively correlated with yield suggesting that these traits are not influenced by the same genetic mechanisms and are not limiting yield in the soybean NAM population. Leaf carbon percentage, leaf nitrogen percentage, and specific leaf area showed strong correlations with yield and may be of interest in breeding programs as a proxy for yield. This work is among the first to use hyperspectral reflectance to model and map the genetic architecture of the rate-limiting steps of photosynthesis.

Bayesian multitrait kernel methods improve multienvironment genome-based prediction

When multitrait data are available, the preferred models are those that are able to account for correlations between phenotypic traits because when the degree of correlation is moderate or large, this increases the genomic prediction accuracy. For this reason, in this article, we explore Bayesian multitrait kernel methods for genomic prediction and we illustrate the power of these models with three-real datasets. The kernels under study were the linear, Gaussian, polynomial, and sigmoid kernels; they were compared with the conventional Ridge regression and GBLUP multitrait models. The results show that, in general, the Gaussian kernel method outperformed conventional Bayesian Ridge and GBLUP multitrait linear models by 2.2–17.45% (datasets 1–3) in terms of prediction performance based on the mean square error of prediction. This improvement in terms of prediction performance of the Bayesian multitrait kernel method can be attributed to the fact that the proposed model is able to capture nonlinear patterns more efficiently than linear multitrait models. However, not all kernels perform well in the datasets used for evaluation, which is why more than one kernel should be evaluated to be able to choose the best kernel.

Genomic mating in outbred species: predicting cross usefulness with additive and total genetic covariance matrices

Diverse crops are both outbred and clonally propagated. Breeders typically use truncation selection of parents and invest significant time, land, and money evaluating the progeny of crosses to find exceptional genotypes. We developed and tested genomic mate selection criteria suitable for organisms of arbitrary homozygosity level where the full-sibling progeny are of direct interest as future parents and/or cultivars. We extended cross variance and covariance variance prediction to include dominance effects and predicted the multivariate selection index genetic variance of crosses based on haplotypes of proposed parents, marker effects, and recombination frequencies. We combined the predicted mean and variance into usefulness criteria for parent and variety development. We present an empirical study of cassava (Manihot esculenta), a staple tropical root crop. We assessed the potential to predict the multivariate genetic distribution (means, variances, and trait covariances) of 462 cassava families in terms of additive and total value using cross-validation. Most variance (89%) and covariance (70%) prediction accuracy estimates were greater than zero. The usefulness of crosses was accurately predicted with good correspondence between the predicted and the actual mean performance of family members breeders selected for advancement as new parents and candidate varieties. We also used a directional dominance model to quantify significant inbreeding depression for most traits. We predicted 47,083 possible crosses of 306 parents and contrasted them to those previously tested to show how mate selection can reveal the new potential within the germplasm. We enable breeders to consider the potential of crosses to produce future parents (progeny with top breeding values) and varieties (progeny with top own performance).

Mitigating tradeoffs in plant breeding

Tradeoffs among plant traits help maintain relative fitness under unpredictable conditions and maximize reproductive success. However, modifying tradeoffs is a breeding challenge since many genes of minor effect are involved. The intensive crosstalk and fine-tuning between growth and defense responsive phytohormones via transcription factors optimizes growth, reproduction, and stress tolerance. There are regulating genes in grain crops that deploy diverse functions to overcome tradeoffs, e.g., miR-156-IPA1 regulates crosstalk between growth and defense to achieve high disease resistance and yield, while OsALDH2B1 loss of function causes imbalance among defense, growth, and reproduction in rice. GNI-A1 regulates seed number and weight in wheat by suppressing distal florets and altering assimilate distribution of proximal seeds in spikelets. Knocking out ABA-induced transcription repressors (AITRs) enhances abiotic stress adaptation without fitness cost in Arabidopsis. Deploying AITRs homologs in grain crops may facilitate breeding. This knowledge suggests overcoming tradeoffs through breeding may expose new ones.

Assessing the correlations and selection criteria between different traits in wheat salt-tolerant genotypes

Salinity is one of the largest stresses blocking horizontal and vertical expansion in agricultural lands. Establishing salt-tolerant genotypes is a promising method to benefit from poor water quality and salinized lands. An integrated method was developed for accomplishing reliable and effective evaluation of traits stability of salt-tolerant wheat. The study aims were to estimate the genetic relationships between explanatory traits and shoot dry matter (SDM), and determine the traits stability under three salinity levels. Morphophysiological and biochemical traits were evaluated as selection criteria for SDM improvement in wheat for salinity tolerance. Three cultivars and three high-yielding doubled haploid lines (DHLs) were used. Three salt (NaCl) levels (control (washed sand), 7 and 14 dS m-1) were applied for 45 days (at the first signs of death in the sensitive genotypes). All morphophysiological traits gradually decreased as salinity levels increased, excluding the number of roots. Decreases were more visible in sensitive genotypes than in tolerant genotypes. All biochemical traits increased as salinity levels increased. Variance inflation factors (VIFs) and condition number exhibited multicollinearity for membrane stability index and polyphenol oxidase activity. After their removal, all VIFs were <10, thereby increasing path coefficient accuracy. Total chlorophyll content (CHL) and catalase (CAT) provided significant direct effects regarding genetic and phenotypic correlations for the three salinity levels and their interactions in path analysis on SDM, indicating their stability. CHL and CAT had high heritability (>0.60%) and genetic gain (>20%) and highly significant genetic correlation, co-heritability, and selection efficiencies for SDM. CHL and CAT could be used as selection criteria for salinity tolerance in wheat-breeding programs. The tolerated line (DHL21) with the check cultivar (Sakha 93) can be also recommended as novel genetic resource for improving salinity tolerance of wheat.

Genotyping crossing parents and family bulks can facilitate cost-efficient genomic prediction strategies in small-scale line breeding programs

Genomic relationship matrices based on mid-parent and family bulk genotypes represent cost-efficient alternatives to full genomic prediction approaches with individually genotyped early generation selection candidates. The routine usage of genomic selection for improving line varieties has gained an increasing popularity in recent years. Harnessing the benefits of this approach can, however, be too costly for many small-scale breeding programs, as in most genomic breeding strategies several hundred or even thousands of lines have to be genotyped each year. The aim of this study was thus to compare a full genomic prediction strategy using individually genotyped selection candidates with genomic predictions based on genotypes obtained from pooled DNA of progeny families as well as genotypes inferred from crossing parents. A population of 722 wheat lines representing 63 families tested in more than 100 multi-environment trials during 2010-2019 was for this purpose employed to conduct an empirical study, which was supplemented by a simulation with genotypic data from further 3855 lines. A similar or higher prediction ability was achieved for grain yield, protein yield, and the protein content when using mid-parent or family bulk genotypes in comparison with pedigree selection in the empirical across family prediction scenario. The difference of these methods with a full genomic prediction strategy became furthermore marginal if pre-existing phenotypic data of the selection candidates was already available. Similar observations were made in the simulation, where the usage of individually genotyped lines or family bulks was generally preferable with smaller family sizes. The proposed methods can thus be regarded as alternatives to full genomic or pedigree selection strategies, especially when pedigree information is limited like in the exchange of germplasm between breeding programs.

Identification of Genomic Regions for Carcass Quality Traits within the American Simmental Association Carcass Merit Program

USDA quality and yield grade are primary driving forces for carcass value in the United States. Carcass improvements can be achieved by making selection decisions based on the results of genetic evaluations in the form of expected progeny differences (EPD), real-time ultrasound imaging, and physical evaluation of candidate breeding animals. In an effort to advance their ability to accurately predict the breeding value of potential sires for carcass traits, the American Simmental Association launched the Carcass Merit Program as a means to collect progeny sire group carcass information. All records were extracted from the American Simmental Association database. Progeny data were organized by sire family and progeny performance phenotypes were constructed. Sire genotypes were filtered, and a multi-locus mixed linear model was used to perform an association analysis on the genotype data, while correcting for cryptic relatedness and pedigree structure. Three chromosomes were found to have genome-wide significance and this conservative approach identified putative QTL in those regions. Three hundred ninety-three novel regions were identified across all traits, as well as 290 novel positional candidate genes. Correlations between carcass characteristics and maternal traits were less unfavorable than those previously reported.

Application of multi-trait Bayesian decision theory for parental genomic selection

In all breeding programs, the decision about which individuals to select and intermate to form the next selection cycle is crucial. The improvement of genetic stocks requires considering multiple traits simultaneously, given that economic value and net genetic merits depend on many traits; therefore, with the advance of computational and statistical tools and genomic selection (GS), researchers are focusing on multi-trait selection. Selection of the best individuals is difficult, especially in traits that are antagonistically correlated, where improvement in one trait might imply a reduction in other(s). There are approaches that facilitate multi-trait selection, and recently a Bayesian decision theory (BDT) has been proposed. Parental selection using BDT has the potential to be effective in multi-trait selection given that it summarizes all relevant quantitative genetic concepts such as heritability, response to selection and the structure of dependence between traits (correlation). In this study, we applied BDT to provide a treatment for the complexity of multi-trait parental selection using three multivariate loss functions (LF), Kullback–Leibler (KL), Energy Score, and Multivariate Asymmetric Loss (MALF), to select the best-performing parents for the next breeding cycle in two extensive real wheat data sets. Results show that the high ranking lines in genomic estimated breeding value (GEBV) for certain traits did not always have low values for the posterior expected loss (PEL). For both data sets, the KL LF gave similar importance to all traits including grain yield. In contrast, the Energy Score and MALF gave a better performance in three of four traits that were different than grain yield. The BDT approach should help breeders to decide based not only on the GEBV per se of the parent to be selected, but also on the level of uncertainty according to the Bayesian paradigm.

Genetic analysis of morphological traits of the spike and reproductivity elements of speltoid chemomutant Triticum aestivum

Long use of the methods of direct intraspecies hybridization in the selective breeding of varieties of Triticum aestivum L. has led to narrowing of their gene fund and close similarity of the genetic potential of their selective breeding traits. Using the method of experimental mutagenesis, one can induce systemic mutants with features of other hexaploid Triticum species which can freely cross-breed with other hybridization offspring, contributing to extension of genetic potential of cultivated wheat and creation of varieties with new levels of manifestation of economically beneficial traits. We studied the pattern of inheritance of morphological traits of the ear, length of the stem and the elements of productivity of speltoid chemomutant of T. aestivum. For the genetic analysis we used hybrids F1 and F2 of soft wheat obtained by cross-breeding speltoid macromutant (Smuhlianka speltoid), induced by the impact of aqueous solution of N-nitroso-N-methylurea (NMU) in the concentration of 0.025% on the seeds of Smuhlianka variety, with plants of Smuhlianka variety (Erythrospermum variety) and Podoloanka (Lutescens variety). To determine the pattern of inheritance the spike morphology, length of the stem and the elements of productivity in F1 hybrids, we calculated the extent of phenotype domination. In populations of F2, we examined plants with different combinations of phenotype manifestation of ear morphology. In F1 hybrids, the speltoid shape of the ear, absence of awns and red colour of the glumes indicated the dominant pattern of inheritance. The high level of phenotype domination of length of the stem and ear, number of spikelets in the main ear indicates the inheritance of the features according to intermediate, partly dominant and over-dominant types. Taking into account the segregation according to the features of spike morphology, awnedness and colour of glumes, the plants of F2 population were divided into phenotype classes, that is 12 and 6 in the combinations of respectively Smuhlianka speltoid × Smuhlianka and Smuhlianka speltoid × Podolianka. We determined that the obtained results are the consequences of dihybrid linkage which corresponds to the theoretical proportion of 12 : 3 : 1. Segregation into non-aristate and aristate plants corresponds to the proportion of monogene segregation of 3 : 1. Within separately distinguished phenotype classes, no independent inheritance of the shape of the ear and awnedness was observed. Dihybrid segregation of F2 plants into speltoid, squarehead and varieties Lutescence/Erythrospermum with quantitative superiority of speltoid plants suggests the control of the trait by two non-allele genes with epistatic interaction. The red colour of the glumes indicates the dominant monogenic pattern of inheritance. Absence of independent inheritance of the shape of the ear and awnedness indicates localization of genes which determine these features in one chromosome.