Search

Plant breeders want to develop cultivars that outperform existing genotypes. Some characteristics (here 'main traits') of these cultivars are categorical and difficult to measure directly. It is important to predict the main trait of newly developed genotypes accurately. In addition to marker data, breeding programs often have information on secondary traits (or 'phenotypes') that are easy to measure. Our goal is to improve prediction of main traits with interpretable relations by combining the two data types using variable selection techniques. However, the genomic characteristics can overwhelm the set of secondary traits, so a standard technique may fail to select any phenotypic variables. We develop a new statistical technique that ensures appropriate representation from both the secondary traits and the genotypic variables for optimal prediction. When two data types (markers and secondary traits) are available, we achieve improved prediction of a binary trait by two steps that are designed to ensure that a significant intrinsic effect of a phenotype is incorporated in the relation before accounting for extra effects of genotypes. First, we sparsely regress the secondary traits on the markers and replace the secondary traits by their residuals to obtain the effects of phenotypic variables as adjusted by the genotypic variables. Then, we develop a sparse logistic classifier using the markers and residuals so that the adjusted phenotypes may be selected first to avoid being overwhelmed by the genotypic variables due to their numerical advantage. This classifier uses forward selection aided by a penalty term and can be computed effectively by a technique called the one-pass method. It compares favorably with other classifiers on simulated and real data. [ABSTRACT FROM AUTHOR]

Objective To document the use of active management of the third stage of labour for preventing postpartum haemorrhage and to explore factors associated with such use in seven developing countries. Methods Nationally representative samples of facility-based deliveries were selected and observed to determine the use of active management of the third stage of labour and associated factors. Policies on active management were assessed through document review and interviews with relevant professionals. Findings Use of a uterotonic during the third or fourth stages of labour was nearly universal. Correct use of active management of the third stage of labour was found in only 0.5% to 32% of observed deliveries due to multiple deficiencies in practice. In every country except Indonesia, policies regarding active management were conflicting. Conclusion Developing countries have not targeted decreasing postpartum haemorrhage as an achievable goal; there is little use of active management of the third stage of labour, and policies regarding such management often conflict. Studies are needed to identify the most effective components of active management so that the most efficient package of practices can be promoted. Copyright © 2009 World Health Organization [ABSTRACT FROM AUTHOR]

Enhancing predictive modeling accuracy in wheat (Triticum aestivum) breeding through the integration of high-throughput phenotyping (HTP) data with genomic information is crucial for maximizing genetic gain. In this study, spanning four locations in the southeastern United States over 3 years, models to predict grain yield (GY) were investigated through different cross-validation approaches. The results demonstrate the superiority of multivariate comprehensive models that incorporate both genomic and HTP data, particularly in accurately predicting GY across diverse locations and years. These HTP-incorporating models achieve prediction accuracies ranging from 0.59 to 0.68, compared to 0.40-0.54 for genomic-only models when tested under different prediction scenarios both across years and locations. The comprehensive models exhibit superior generalization to new environments and achieve the highest accuracy when trained on diverse datasets. Predictive accuracy improves as models incorporate data from multiple years, highlighting the importance of considering temporal dynamics in modeling approaches. The study reveals that multivariate prediction outperformed genomic prediction methods in predicting lines across years and locations. The percentage of top 25% lines selected based on multivariate prediction was higher compared to genomic-only models, indicated by higher specificity, which is the proportion of correctly identified top-yielding lines that matched the observed top 25% performance across different sites and years. Additionally, the study addresses the prediction of untested locations based on other locations within the same year and in new years at previously tested locations. Findings show the comprehensive models effectively extrapolate to new environments, highlighting their potential for guiding breeding strategies., (© 2024 The Author(s). The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

The development of improved cultivars requires establishing multi-environment trials (METs) to evaluate their performance under a wide range of environmental conditions. However, the high phenotyping costs often limit the capacity to evaluate genotypes in all the target environments. Our main objective was to explore the potential of implementing sparse testing in cassava breeding programs to reduce the cost of phenotyping in METs. The population used in this study consisted of 435 cassava genotypes evaluated in five environments in Nigeria for dry matter (dm) and fresh root yield (fyld). Sparse testing designs were developed based on non-overlapping (NOL), completely overlapping (OL), and intermediates between NOL and OL genotypes. Three prediction models were assessed (one based on phenotypes only, while two had genomic data). All the three models had a higher predictive ability and a lower mean square error (MSE) when a large training set was used. Predictive ability increased and MSE reduced when genotype-by-environment interaction (G × E) was modeled for the same training set sizes and allocations. Predictive ability decreased while MSE increased with the increasing OL genotypes across the environments, suggesting that only a few OL genotypes may be required to set up METs for model training. Sparse testing using a model incorporating G × E could be implemented to reduce cost of phenotyping in cassava METs. If data were available, integrating crop growth models (CGMs) with genomic prediction holds the potential to improve predictive ability. The training population used for sparse testing could be optimized to determine the optimal size and distribution of genotypes to increase the predictive ability and reduce cost under a fixed budget., (© 2025 The Author(s). The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Integrating genomic, hyperspectral imaging (HSI), and environmental data enhances wheat yield predictions, with HSI providing detailed spectral insights for predicting complex grain yield (GY) traits. Incorporating HSI data with single nucleotide polymorphic markers (SNPs) resulted in a substantial improvement in predictive ability compared to the conventional genomic prediction models. Over the course of several years, the prediction ability varied due to diverse weather conditions. The most comprehensive parametric model tested, which included SNPs, HSI, and environmental covariates data, consistently achieved the best results, closely followed by machine learning (ML) approaches when considering the same omics data. For example, the most comprehensive model (M9), under the forward prediction cross-validation scheme, predicted the GY of the 2023 growing season using data from 2021 and 2022 for a correlation between predicted and observed values of 0.53. This model demonstrated superior performance compared to less complex models, emphasizing the advantage of integrating numerous data sources and their interactive effects. Furthermore, when comparing the top 25% of the predicted lines versus the corresponding observed lines with the highest GY, the M9 model returned a coincide index (CI) of 55% (i.e., in both sets, 55% of the top 25% values were common), whereas for the highest performing ML model (gradient boosting regression), the CI was of 46%. This study highlights the potential of multi-data source approaches to accelerate the selection of heat-tolerant wheat genotypes., (© 2025 The Author(s). The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Complex traits influenced by multiple genes pose challenges for marker-assisted selection (MAS) in breeding. Genomic selection (GS) is a promising strategy for achieving higher genetic gains in quantitative traits by stacking favorable alleles into elite cultivars. Resistance to Fusarium oxysporum f. sp. niveum (Fon) race 2 in watermelon is a polygenic trait with moderate heritability. This study evaluated GS as an additional approach to quantitative trait loci (QTL) analysis/marker-assisted selection (MAS) for enhancing Fon race 2 resistance in elite watermelon cultivars. Objectives were to: (1) assess the accuracy of genomic prediction (GP) models for predicting Fon race 2 resistance in a F 2:3 versus a recombinant inbred line (RIL) population, (2) rank and select families in each population based on genomic estimated breeding values (GEBVs) for developing testing populations, and (3) determined how many of the most superior families based on GEBV also have all QTL associated with Fon race 2 resistance. GBS-SNP data from genotyping-by-sequencing (GBS) for two populations were used, and parental line genome sequences were used as references. The GBLUP and random forest outperformed the other three parametric (GBLUP, Bayes B, Bayes LASSO) and three nonparametric AI (random forest, SVM linear, and SVM radial) models, with correlations of 0.48 and 0.68 in the F 2:3 and RIL population, respectively. Selection intensities (SI) of 10%, 20%, and 30% showed that superior families with highest GEBV can also comprise all QTL associated with Fon race 2 resistance, highlighting GP efficacy in improving elite watermelon cultivars with polygenic traits of disease resistance., Competing Interests: Declarations. Conflict of interest: The authors have no relevant financial interests to disclose. Dr. Amnon Levi was an editor and Sandra E. Branham serves as an associate editor for Theoretical and Applied Genetics., (© 2025. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Through enviromics, precision breeding leverages innovative geotechnologies to customize crop varieties to specific environments, potentially improving both crop yield and genetic selection gains. In Brazil's four southernmost states, data from 183 distinct geographic field trials (also accounting for 2017-2021) covered information on 164 genotypes: 79 phenotyped maize hybrid genotypes for grain yield and their 85 nonphenotyped parents. Additionally, 1342 envirotypic covariates from weather, soil, sensor-based, and satellite sources were collected to engineer 10 K synthetic enviromic markers via machine learning. Soil, radiation light, and surface temperature variations remarkably affect differential genotype yield, hinting at ecophysiological adjustments including evapotranspiration and photosynthesis. The enviromic ensemble-based random regression model showcases superior predictive performance and efficiency compared to the baseline and kernel models, matching the best genotypes to specific geographic coordinates. Clustering analysis has identified regions that minimize genotype-environment (G × E) interactions. These findings underscore the potential of enviromics in crafting specific parental combinations to breed new, higher-yielding hybrid crops. The adequate use of envirotypic information can enhance the precision and efficiency of maize breeding by providing important inputs about the environmental factors that affect the average crop performance. Generating enviromic markers associated with grain yield can enable a better selection of hybrids for specific environments., (© 2024 The Author(s). New Phytologist © 2024 New Phytologist Foundation.)

Fruit texture is a priority trait that guarantees the long-term economic sustainability of the cranberry industry through value-added products such as sweetened dried cranberries (SDCs). To develop a standard methodology to measure texture, we conducted a comparative analysis of 22 textural traits using five different methods under both harvest and postharvest conditions in 10 representative cranberry cultivars. A set of textural traits from the 10%-strain compression and puncture methods were identified that differentiate between cultivars primarily based on hardness/stiffness and elasticity properties. The complementary use of both methodologies allowed for a detailed evaluation by capturing the effect of key texture-determining factors such as structure, flesh, and skin. Furthermore, the high effectiveness of this approach in different conditions and its ability to capture high phenotypic variation in cultivars highlights its great potential for applicability in various areas of the value chain and research. Therefore, this study provides an informed reference for unifying future efforts to enhance cranberry fruit texture and quality., (Published 2024. This article is a U.S. Government work and is in the public domain in the USA. Journal of Texture Studiespublished by Wiley Periodicals LLC.)

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., Competing Interests: Authors AB and EC were employed by the company Embrapa Coffee, Brazilian Agricultural Research Corporation (Embrapa). The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Nascimento, Nascimento, Azevedo, Oliveira, Caixeta and Jarquin.)

Sweet corn breeding programs, like field corn, focus on the development of elite inbred lines to produce commercial hybrids. For this reason, genomic selection models can help the in silico prediction of hybrid crosses from the elite lines, which is hypothesized to improve the test cross scheme, leading to higher genetic gain in a breeding program. This study aimed to explore the potential of implementing genomic selection in a sweet corn breeding program through hybrid prediction in a within-site across-year and across-site framework. A total of 506 hybrids were evaluated in six environments (California, Florida, and Wisconsin, in the years 2020 and 2021). A total of 20 traits from three different groups were measured (plant-, ear-, and flavor-related traits) across the six environments. Eight statistical models were considered for prediction, as the combination of two genomic prediction models (GBLUP and RKHS) with two different kernels (additive and additive + dominance), and in a single- and multi-trait framework. Also, three different cross-validation schemes were tested (CV1, CV0, and CV00). The different models were then compared based on the correlation between the estimated breeding values/total genetic values and phenotypic measurements. Overall, heritabilities and correlations varied among the traits. The models implemented showed good accuracies for trait prediction. The GBLUP implementation outperformed RKHS in all cross-validation schemes and models. Models with additive plus dominance kernels presented a slight improvement over the models with only additive kernels for some of the models examined. In addition, models for within-site across-year and across-site performed better in the CV0 than the CV00 scheme, on average. Hence, GBLUP should be considered as a standard model for sweet corn hybrid prediction. In addition, we found that the implementation of genomic prediction in a sweet corn breeding program presented reliable results, which can improve the testcross stage by identifying the top candidates that will reach advanced field-testing stages., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2024 Peixoto, Leach, Jarquin, Flannery, Zystro, Tracy, Bhering and Resende.)

The aim of this study was to evaluate the accuracy of the ridge regression best linear unbiased prediction model across different traits, parent population sizes, and breeding strategies when estimating breeding values in common bean (Phaseolus vulgaris). Genomic selection was implemented to make selections within a breeding cycle and compared across five different breeding strategies (single seed descent, mass selection, pedigree method, modified pedigree method, and bulk breeding) following 10 breeding cycles. The model was trained on a simulated population of recombinant inbreds genotyped for 1010 single nucleotide polymorphism markers including 38 known quantitative trait loci identified in the literature. These QTL included 11 for seed yield, eight for white mold disease incidence, and 19 for days to flowering. Simulation results revealed that realized accuracies fluctuate depending on the factors investigated: trait genetic architecture, breeding strategy, and the number of initial parents used to begin the first breeding cycle. Trait architecture and breeding strategy appeared to have a larger impact on accuracy than the initial number of parents. Generally, maximum accuracies (in terms of the correlation between true and estimated breeding value) were consistently achieved under a mass selection strategy, pedigree method, and single seed descent method depending on the simulation parameters being tested. This study also investigated model updating, which involves retraining the prediction model with a new set of genotypes and phenotypes that have a closer relation to the population being tested. While it has been repeatedly shown that model updating generally improves prediction accuracy, it benefited some breeding strategies more than others. For low heritability traits (e.g., yield), conventional phenotype-based selection methods showed consistent rates of genetic gain, but genetic gain under genomic selection reached a plateau after fewer cycles. This plateauing is likely a cause of faster fixation of alleles and a diminishing of genetic variance when selections are made based on estimated breeding value as opposed to phenotype., (© 2024 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

The availability of high-dimensional genomic data and advancements in genome-based prediction models (GP) have revolutionized and contributed to accelerated genetic gains in soybean breeding programs. GP-based sparse testing is a promising concept that allows increasing the testing capacity of genotypes in environments, of genotypes or environments at a fixed cost, or a substantial reduction of costs at a fixed testing capacity. This study represents the first attempt to implement GP-based sparse testing in soybeans by evaluating different training set compositions going from non-overlapped RILs until almost the other extreme of having same set of genotypes observed across environments for different training set sizes. A total of 1,755 recombinant inbred lines (RILs) tested in nine environments were used in this study. RILs were derived from 39 bi-parental populations of the Soybean Nested Association Mapping (NAM) project. The predictive abilities of various models and training set sizes and compositions were investigated. Training compositions included a range of ratios of overlapping (O-RILs) and non-overlapping (NO-RILs) RILs across environments, as well as a methodology to maximize or minimize the genetic diversity in a fixed-size sample. Reducing the training set size compromised predictive ability in most training set compositions. Overall, maximizing the genetic diversity within the training set and the inclusion of O-RILs increased prediction accuracy given a fixed training set size; however, the most complex model was less affected by these factors. More testing environments in the early stages of the breeding pipeline can provide a more comprehensive assessment of genotype stability and adaptation which are fundamental for the precise selection of superior genotypes adapted to a wide range of environments., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision., (Copyright © 2023 Persa, Canella Vieira, Rios, Hoyos-Villegas, Messina, Runcie and Jarquin.)

The adoption of dicamba-tolerant (DT) soybean in the United States resulted in extensive off-target dicamba damage to non-DT vegetation across soybean-producing states. Although soybeans are highly sensitive to dicamba, the intensity of observed symptoms and yield losses are affected by the genetic background of genotypes. Thus, the objective of this study was to detect novel marker-trait associations and expand on previously identified genomic regions related to soybean response to off-target dicamba. A total of 551 non-DT advanced breeding lines derived from 232 unique bi-parental populations were phenotyped for off-target dicamba across nine environments for three years. Breeding lines were genotyped using the Illumina Infinium BARCSoySNP6K BeadChip. Filtered SNPs were included as predictors in Random Forest (RF) and Support Vector Machine (SVM) models in a forward stepwise selection loop to identify the combination of SNPs yielding the highest classification accuracy. Both RF and SVM models yielded high classification accuracies (0.76 and 0.79, respectively) with minor extreme misclassifications (observed tolerant predicted as susceptible, and vice-versa). Eight genomic regions associated with off-target dicamba tolerance were identified on chromosomes 6 [Linkage Group (LG) C2], 8 (LG A2), 9 (LG K), 10 (LG O), and 19 (LG L). Although the genetic architecture of tolerance is complex, high classification accuracies were obtained when including the major effect SNP identified on chromosome 6 as the sole predictor. In addition, candidate genes with annotated functions associated with phases II (conjugation of hydroxylated herbicides to endogenous sugar molecules) and III (transportation of herbicide conjugates into the vacuole) of herbicide detoxification in plants were co-localized with significant markers within each genomic region. Genomic prediction models, as reported in this study, can greatly facilitate the identification of genotypes with superior tolerance to off-target dicamba., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Canella Vieira, Zhou, Jarquin, Zhou, Diers, Riechers, Nguyen and Shannon.)

High-throughput phenotyping (HTP) has expanded the dimensionality of data in plant research; however, HTP has resulted in few novel biological discoveries to date. Field-based HTP (FHTP), using small unoccupied aerial vehicles (UAVs) equipped with imaging sensors, can be deployed routinely to monitor segregating plant population interactions with the environment under biologically meaningful conditions. Here, flowering dates and plant height, important phenological fitness traits, were collected on 520 segregating maize recombinant inbred lines (RILs) in both irrigated and drought stress trials in 2018. Using UAV phenomic, single nucleotide polymorphism (SNP) genomic, as well as combined data, flowering times were predicted using several scenarios. Untested genotypes were predicted with 0.58, 0.59, and 0.41 prediction ability for anthesis, silking, and terminal plant height, respectively, using genomic data, but prediction ability increased to 0.77, 0.76, and 0.58 when phenomic and genomic data were used together. Using the phenomic data in a genome-wide association study, a heat-related candidate gene (GRMZM2G083810; hsp18f) was discovered using temporal reflectance phenotypes belonging to flowering times (both irrigated and drought) trials where heat stress also peaked. Thus, a relationship between plants and abiotic stresses belonging to a specific time of growth was revealed only through use of temporal phenomic data. Overall, this study showed that (i) it is possible to predict complex traits using high dimensional phenomic data between different environments, and (ii) temporal phenomic data can reveal a time-dependent association between genotypes and abiotic stresses, which can help understand mechanisms to develop resilient plants., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Society for Experimental Biology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Early canopy coverage is a desirable trait that is a major determinant of yield in soybean (Glycine max). Variation in traits comprising shoot architecture can influence canopy coverage, canopy light interception, canopy-level photosynthesis, and source-sink partitioning efficiency. However, little is known about the extent of phenotypic diversity of shoot architecture traits and their genetic control in soybean. Thus, we sought to understand the contribution of shoot architecture traits to canopy coverage and to determine the genetic control of these traits. We examined the natural variation for shoot architecture traits in a set of 399 diverse maturity group I soybean (SoyMGI) accessions to identify relationships between traits, and to identify loci that are associated with canopy coverage and shoot architecture traits. Canopy coverage was correlated with branch angle, number of branches, plant height, and leaf shape. Using previously collected 50K single nucleotide polymorphism data, we identified quantitative trait locus (QTL) associated with branch angle, number of branches, branch density, leaflet shape, days to flowering, maturity, plant height, number of nodes, and stem termination. In many cases, QTL intervals overlapped with previously described genes or QTL. We also found QTL associated with branch angle and leaflet shape located on chromosomes 19 and 4, respectively, and these QTL overlapped with QTL associated with canopy coverage, suggesting the importance of branch angle and leaflet shape in determining canopy coverage. Our results highlight the role individual architecture traits play in canopy coverage and contribute information on their genetic control that could help facilitate future efforts in their genetic manipulation., (© 2023 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Soybean [Glycine max (L.) Merr.] is a significant source of protein and oil and is also widely used as animal feed. Thus, developing lines that are superior in terms of yield, protein, and oil content is important to feed the ever-growing population. As opposed to high-cost phenotyping, genotyping is both cost and time efficient for breeders because evaluating new lines in different environments (location-year combinations) can be costly. Several genomic prediction (GP) methods have been developed to use the marker and environment data effectively to predict the yield or other relevant phenotypic traits of crops. Our study compares a conventional GP method (genomic best linear unbiased predictor [GBLUP]), a kernel method (Gaussian kernel [GK]), an artificial-intelligence (AI) method (deep learning [DL]), and a hybrid method that corresponds to the emulation of a DL model using a kernel method (an arc-cosine kernel [AK]) in terms of their prediction accuracies for predicting grain yield, oil, and protein using data from the soybean nested association mapping experiment (1,379 genotypes tested in six environments, all genotypes in all environments). The relative performance of the four methods varied with the response variable and whether the model includes the genotype × environmental interaction (G×E) effects or not. The GBLUP consistently showed better performances, whereas GK and AK followed a similar pattern to GBLUP and DL performed slightly worse than the other three methods in most of the cases; however, this may also be attributed to suboptimal hyperparameters. The DL method performed particularly worse than the other three methods in presence of the G×E effects., (© 2022 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Increasing the rate of genetic gain for seed yield remains the primary breeding objective in both public and private soybean [Glycine max (L.) Merr.] breeding programs. Genomic selection (GS) has the potential to accelerate the rate of genetic gain for soybean seed yield. Limited studies to date have validated GS accuracy and directly compared GS with phenotypic selection (PS), and none have been reported in soybean. This study conducted the first empirical validation of GS for increasing seed yield using over 1,500 lines and over 7 yr (2010-2016) of replicated experiments in the University of Nebraska-Lincoln soybean breeding program. The study was designed to capture the varying genetic relatedness of the training population to three validation sets: two large biparental populations (TBP-1 and TBP-2) and a large validation set comprised of 457 preselected advanced lines derived from 45 biparental populations (TMP). We found that prediction accuracy (.54) realized in our validation experiments was comparable with what we obtained from a series of cross-validation experiments (.64). Both GS and PS were more effective for increasing the population mean performance compared with random selection (RS). We found a selection advantage of GS over PS, where higher genetic gain and identification of top-performing lines was maximized at 10-20% selected proportion. Genomic selection led to small increases in genetic similarity when compared with PS and RS presumably because of a significant shift on allelic frequencies toward the extremes, suggesting that it could erode genetic diversity more quickly. Overall, we found that GS can perform as effectively as PS but that measures should be considered to protect against loss of genetic variance when using GS., (© 2022 The Authors. The Plant Genome published by Wiley Periodicals LLC on behalf of Crop Science Society of America.)

Key Message: A reference study for breeders aiming at maximizing genetic gain in common bean. Depending on trait heritability and genetic architecture, conventional approaches may provide an advantage over other frameworks. Dry beans (Phaseolus vulgaris L.) are a nutrient dense legume that is consumed by developed and developing nations around the world. The progress to improve this crop has been quite steady. However, with the continued rise in global populations, there are demands to expedite genetic gains. Plant breeders have been at the forefront at increasing yields in the common bean. As breeding programs are both time-consuming and resource intensive, resource allocation must be carefully considered. To assist plant breeders, computer simulations can provide useful information that may then be applied to the real world. This study evaluated multiple breeding scenarios in the common bean and involved five selection strategies, three breeding frameworks, and four different parental population sizes. In addition, the breeding scenarios were implemented in three different traits: days to flowering, white mold tolerance, and seed yield. Results from the study reflect the complexity of breeding programs, with the optimal breeding scenario varying based on trait being selected. Relative genetic gains per cycle of up to 8.69% for seed yield could be obtained under the use of the optimal breeding scenario. Principal component analyses revealed similarity between strategies, where single seed descent and the modified pedigree method would often aggregate. As well, clusters in the direction of the Hamming distance eigenvector are a good indicator of poor performance in a strategy., (© 2023. The Author(s).)

The widespread adoption of genetically modified (GM) dicamba-tolerant (DT) soybean was followed by numerous reports of off-target dicamba damage and yield losses across most soybean-producing states. In this study, a subset of the USDA Soybean Germplasm Collection consisting of 382 genetically diverse soybean accessions originating from 15 countries was used to identify genomic regions associated with soybean response to off-target dicamba exposure. Accessions were genotyped with the SoySNP50K BeadChip and visually screened for damage in environments with prolonged exposure to off-target dicamba. Two models were implemented to detect significant marker-trait associations: the Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK) and a model that allows the inclusion of population structure in interaction with the environment (G×E) to account for variable patterns of genotype responses in different environments. Most accessions (84%) showed a moderate response, either moderately tolerant or moderately susceptible, with approximately 8% showing tolerance and susceptibility. No differences in off-target dicamba damage were observed across maturity groups and centers of origin. Both models identified significant associations in regions of chromosomes 10 and 19. The BLINK model identified additional significant marker-trait associations on chromosomes 11, 14, and 18, while the G×E model identified another significant marker-trait association on chromosome 15. The significant SNPs identified by both models are located within candidate genes possessing annotated functions involving different phases of herbicide detoxification in plants. These results entertain the possibility of developing non-GM soybean cultivars with improved tolerance to off-target dicamba exposure and potentially other synthetic auxin herbicides. Identification of genetic sources of tolerance and genomic regions conferring higher tolerance to off-target dicamba may sustain and improve the production of other non-DT herbicide soybean production systems, including the growing niche markets of organic and conventional soybean., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Canella Vieira, Jarquin, do Nascimento, Lee, Zhou, Smothers, Zhou, Diers, Riechers, Xu, Shannon, Chen and Nguyen.)

The likelihood of success in developing modern cultivars depend on multiple factors, including the identification of suitable parents to initiate new crosses, and characterizations of genomic regions associated with target traits. The objectives of the present study were to (a) determine the best economic weights of four major wheat diseases (leaf spot, common bunt, leaf rust, and stripe rust) and grain yield for multi-trait restrictive linear phenotypic selection index (RLPSI), (b) select the top 10% cultivars and lines (hereafter referred as genotypes) with better resistance to combinations of the four diseases and acceptable grain yield as potential parents, and (c) map genomic regions associated with resistance to each disease using genome-wide association study (GWAS). A diversity panel of 196 spring wheat genotypes was evaluated for their reaction to stripe rust at eight environments, leaf rust at four environments, leaf spot at three environments, common bunt at two environments, and grain yield at five environments. The panel was genotyped with the Wheat 90K SNP array and a few KASP SNPs of which we used 23,342 markers for statistical analyses. The RLPSI analysis performed by restricting the expected genetic gain for yield displayed significant (p < 0.05) differences among the 3125 economic weights. Using the best four economic weights, a subset of 22 of the 196 genotypes were selected as potential parents with resistance to the four diseases and acceptable grain yield. GWAS identified 37 genomic regions, which included 12 for common bunt, 13 for leaf rust, 5 for stripe rust, and 7 for leaf spot. Each genomic region explained from 6.6 to 16.9% and together accounted for 39.4% of the stripe rust, 49.1% of the leaf spot, 94.0% of the leaf rust, and 97.9% of the common bunt phenotypic variance combined across all environments. Results from this study provide valuable information for wheat breeders selecting parental combinations for new crosses to develop improved germplasm with enhanced resistance to the four diseases as well as the physical positions of genomic regions that confer resistance, which facilitates direct comparisons for independent mapping studies in the future.

Prediction models based on pedigree and/or molecular marker information are now an inextricable part of the crop breeding programs and have led to increased genetic gains in many crops. Optimization of IRRI's rice drought breeding program is crucial for better implementation of selections based on predictions. Historical datasets with precise and robust pedigree information have been a great resource to help optimize the prediction models in the breeding programs. Here, we leveraged 17 years of historical drought data along with the pedigree information to predict the new lines or environments and dissect the G × E interactions. Seven models ranging from basic to proposed higher advanced models incorporating interactions, and genotypic specific effects were used. These models were tested with three cross-validation schemes (CV1, CV2, and CV0) to assess the predictive ability of tested and untested lines in already observed environments and tested lines in novel or new environments. In general, the highest prediction abilities were obtained when the model accounting interactions between pedigrees (additive) and environment were included. The CV0 scheme (predicting unobserved or novel environments) reveals very low predictive abilities among the three schemes. CV1 and CV2 schemes that borrow information from the target and correlated environments have much higher predictive abilities. Further, predictive ability was lower when predicting lines in non-stress conditions using drought data as training set and/or vice-versa . When predicting the lines using the data sets under the same conditions (stress or non-stress data sets), much better prediction accuracy was obtained. These results provide conclusive evidence that modeling G × E interactions are important in predictions. Thus, considering G × E interactions would help to build enhanced genomic or pedigree-based prediction models in the rice breeding program. Further, it is crucial to borrow the correlated information from other environments to improve prediction accuracy., Competing Interests: The authors declare that the research review was conducted in the absence of any commercial or economic associations that could be construed as a potential conflict of interest., (Copyright © 2022 Khanna, Anumalla, Catolos, Bhosale, Jarquin and Hussain.)

The availability of high-dimensional molecular markers has allowed plant breeding programs to maximize their efficiency through the genomic prediction of a phenotype of interest. Yield is a complex quantitative trait whose expression is sensitive to environmental stimuli. In this research, we investigated the potential of incorporating soil texture information and its interaction with molecular markers via covariance structures for enhancing predictive ability across breeding scenarios. A total of 797 soybean lines derived from 367 unique bi-parental populations were genotyped using the Illumina BARCSoySNP6K and tested for yield during 5 years in Tiptonville silt loam, Sharkey clay, and Malden fine sand environments. Four statistical models were considered, including the GBLUP model (M1), the reaction norm model (M2) including the interaction between molecular markers and the environment (G×E), an extended version of M2 that also includes soil type (S), and the interaction between soil type and molecular markers (G×S) (M3), and a parsimonious version of M3 which discards the G×E term (M4). Four cross-validation scenarios simulating progeny testing and line selection of tested-untested genotypes (TG, UG) in observed-unobserved environments [OE, UE] were implemented (CV2 [TG, OE], CV1 [UG, OE], CV0 [TG, UE], and CV00 [UG, UE]). Across environments, the addition of G×S interaction in M3 decreased the amount of variability captured by the environment (-30.4%) and residual (-39.2%) terms as compared to M1. Within environments, the G×S term in M3 reduced the variability captured by the residual term by 60 and 30% when compared to M1 and M2, respectively. M3 outperformed all the other models in CV2 (0.577), CV1 (0.480), and CV0 (0.488). In addition to the Pearson correlation, other measures were considered to assess predictive ability and these showed that the addition of soil texture seems to structure/dissect the environmental term revealing its components that could enhance or hinder the predictability of a model, especially in the most complex prediction scenario (CV00). Hence, the availability of soil texture information before the growing season could be used to optimize the efficiency of a breeding program by allowing the reconsideration of field experimental design, allocation of resources, reduction of preliminary trials, and shortening of the breeding cycle., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Canella Vieira, Persa, Chen and Jarquin.)

Both the Linear Phenotypic Selection Index (LPSI) and the Restrictive Linear Phenotypic Selection Index (RLPSI) have been widely used to select parents and progenies, but the effect of economic weights on the selection parameters (the expected genetic gain, response to selection, and the correlation between the indices and genetic merits) have not been investigated in detail. Here, we (i) assessed combinations of 2304 economic weights using four traits (maturity, plant height, grain yield and grain protein content) recorded under four organically (low nitrogen) and five conventionally (high nitrogen) managed environments, (ii) compared single-trait and multi-trait selection indices (LPSI vs. RLPSI by imposing restrictions to the expected genetic gain of either yield or grain protein content), and (iii) selected a subset of about 10% spring wheat cultivars that performed very well under organic and/or conventional management systems. The multi-trait selection indices, with and without imposing restrictions, were superior to single trait selection. However, the selection parameters differed quite a lot depending on the economic weights, which suggests the need for optimizing the weights. Twenty-two of the 196 cultivars that showed superior performance under organic and/or conventional management systems were consistently selected using all five of the selected economic weights, and at least two of the selection scenarios. The selected cultivars belonged to the Canada Western Red Spring (16 cultivars), the Canada Northern Hard Red (3), and the Canada Prairie Spring Red (3), and required 83-93 days to maturity, were 72-100 cm tall, and produced from 4.0 to 6.2 t ha -1 grain yield with 14.6-17.7% GPC. The selected cultivars would be highly useful, not only as potential trait donors for breeding under an organic management system, but also for other studies, including nitrogen use efficiency.

Some previous studies have assessed the predictive ability of genome-wide selection on stripe (yellow) rust resistance in wheat, but the effect of genotype by environment interaction (GEI) in prediction accuracies has not been well studied in diverse genetic backgrounds. Here, we compared the predictive ability of a model based on phenotypic data only (M1), the main effect of phenotype and molecular markers (M2), and a model that incorporated GEI (M3) using three cross-validations (CV1, CV2, and CV0) scenarios of interest to breeders in six spring wheat populations. Each population was evaluated at three to eight field nurseries and genotyped with either the DArTseq technology or the wheat 90K single nucleotide polymorphism arrays, of which a subset of 1,058- 23,795 polymorphic markers were used for the analyses. In the CV1 scenario, the mean prediction accuracies of the M1, M2, and M3 models across the six populations varied from -0.11 to -0.07, from 0.22 to 0.49, and from 0.19 to 0.48, respectively. Mean accuracies obtained using the M3 model in the CV1 scenario were significantly greater than the M2 model in two populations, the same in three populations, and smaller in one population. In both the CV2 and CV0 scenarios, the mean prediction accuracies of the three models varied from 0.53 to 0.84 and were not significantly different in all populations, except the Attila/CDC Go in the CV2, where the M3 model gave greater accuracy than both the M1 and M2 models. Overall, the M3 model increased prediction accuracies in some populations by up to 12.4% and decreased accuracy in others by up to 17.4%, demonstrating inconsistent results among genetic backgrounds that require considering each population separately. This is the first comprehensive genome-wide prediction study that investigated details of the effect of GEI on stripe rust resistance across diverse spring wheat populations.

Some studies have investigated the potential of genomic selection (GS) on stripe rust, leaf rust, Fusarium head blight (FHB), and leaf spot in wheat, but none of them have assessed the effect of the reaction norm model that incorporated GE interactions. In addition, the prediction accuracy on common bunt has not previously been studied. Here, we investigated within-population prediction accuracies using the baseline M1 model and two reaction norm models (M2 and M3) with three random cross-validation (CV1, CV2, and CV0) schemes. Three Canadian spring wheat populations were evaluated in up to eight field environments and genotyped with 3158, 5732, and 23,795 polymorphic markers. The M3 model that incorporated GE interactions reduced residual variance by an average of 10.2% as compared with the main effect M2 model and increased prediction accuracies on average by 2-6%. In some traits, the M3 model increased prediction accuracies up to 54% as compared with the M2 model. The average prediction accuracies of the M3 model with CV1, CV2, and CV0 schemes varied from 0.02 to 0.48, from 0.25 to 0.84, and from 0.14 to 0.87, respectively. In both CV2 and CV0 schemes, stripe rust in all three populations, common bunt and leaf rust in two populations, as well as FHB severity, FHB index, and leaf spot in one population had high to very high (0.54-0.87) prediction accuracies. This is the first comprehensive genomic selection study on five major diseases in spring wheat.

Key Message: Using phenotype data of three spring wheat populations evaluated at 6-15 environments under two management systems, we found moderate to very high prediction accuracies across seven traits. The phenotype data collected under an organic management system effectively predicted the performance of lines in the conventional management and vice versa. There is growing interest in developing wheat cultivars specifically for organic agriculture, but we are not aware of the effect of organic management on the predictive ability of genomic selection (GS). Here, we evaluated within populations prediction accuracies of four GS models, four combinations of training and testing sets, three reaction norm models, and three random cross-validations (CV) schemes in three populations phenotyped under organic and conventional management systems. Our study was based on a total of 578 recombinant inbred lines and varieties from three spring wheat populations, which were evaluated for seven traits at 3-9 conventionally and 3-6 organically managed field environments and genotyped either with the wheat 90 K SNP array or DArTseq. We predicted the management systems (CV0 M ) or environments (CV0), a subset of lines that have been evaluated in either management (CV2 M ) or some environments (CV2), and the performance of newly developed lines in either management (CV1 M ) or environments (CV1). The average prediction accuracies of the model that incorporated genotype × environment interactions with CV0 and CV2 schemes varied from 0.69 to 0.97. In the CV1 and CV1M schemes, prediction accuracies ranged from - 0.12 to 0.77 depending on the reaction norm models, the traits, and populations. In most cases, grain protein showed the highest prediction accuracies. The phenotype data collected under the organic management effectively predicted the performance of lines under conventional management and vice versa. This is the first comprehensive GS study that investigated the effect of the organic management system in wheat., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

Genomic-enabled prediction models are of paramount importance for the successful implementation of genomic selection (GS) based on breeding values. As opposed to animal breeding, plant breeding includes extensive multienvironment and multiyear field trial data. Hence, genomic-enabled prediction models should include genotype × environment (G × E) interaction, which most of the time increases the prediction performance when the response of lines are different from environment to environment. In this chapter, we describe a historical timeline since 2012 related to advances of the GS models that take into account G × E interaction. We describe theoretical and practical aspects of those GS models, including the gains in prediction performance when including G × E structures for both complex continuous and categorical scale traits. Then, we detailed and explained the main G × E genomic prediction models for complex traits measured in continuous and noncontinuous (categorical) scale. Related to G × E interaction models this review also examine the analyses of the information generated with high-throughput phenotype data (phenomic) and the joint analyses of multitrait and multienvironment field trial data that is also employed in the general assessment of multitrait G × E interaction. The inclusion of nongenomic data in increasing the accuracy and biological reliability of the G × E approach is also outlined. We show the recent advances in large-scale envirotyping (enviromics), and how the use of mechanistic computational modeling can derive the crop growth and development aspects useful for predicting phenotypes and explaining G × E., (© 2022. The Author(s).)

The advent of plant phenomics, coupled with the wealth of genotypic data generated by next-generation sequencing technologies, provides exciting new resources for investigations into and improvement of complex traits. However, these new technologies also bring new challenges in quantitative genetics, namely, a need for the development of robust frameworks that can accommodate these high-dimensional data. In this chapter, we describe methods for the statistical analysis of high-throughput phenotyping (HTP) data with the goal of enhancing the prediction accuracy of genomic selection (GS). Following the Introduction in Sec. 1, Sec. 2 discusses field-based HTP, including the use of unoccupied aerial vehicles and light detection and ranging, as well as how we can achieve increased genetic gain by utilizing image data derived from HTP. Section 3 considers extending commonly used GS models to integrate HTP data as covariates associated with the principal trait response, such as yield. Particular focus is placed on single-trait, multi-trait, and genotype by environment interaction models. One unique aspect of HTP data is that phenomics platforms often produce large-scale data with high spatial and temporal resolution for capturing dynamic growth, development, and stress responses. Section 4 discusses the utility of a random regression model for performing longitudinal modeling. The chapter concludes with a discussion of some standing issues., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Genomic selection (GS) is a methodology that revolutionized the process of breeding improved genetic materials in plant and animal breeding programs. It uses predicted genomic values of the potential of untested/unobserved genotypes as surrogates of phenotypes during the selection process. Such that the predicted genomic values are obtained using exclusively the marker profiles of the untested genotypes, and these potentially can be used by breeders for screening the genotypes to be advanced in the breeding pipeline, to identify potential parents for next improvement cycles, or to find optimal crosses for targeting genotypes among others. Conceptually, GS initially requires a set of genotypes with both molecular marker information and phenotypic data for model calibration and then the performance of untested genotypes is predicted using their marker profiles only. Hence, it is expected that breeders would look at these values in order to conduct selections. Even though the concept of GS seems trivial, due to the high dimensional nature of the data delivered from modern sequencing technologies where the number of molecular markers (p) excess by far the number of data points available for model fitting (n; p ≫ n) a complete renovated set of prediction models was needed to cope with this challenge. In this chapter, we provide a conceptual framework for comparing statistical models to overcome the "large p, small n problem." Given the very large diversity of GS models only the most popular are presented here; mainly we focused on linear regression-based models and nonparametric models that predict the genetic estimated breeding values (GEBV) in a single environment considering a single trait only, mainly in the context of plant breeding., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

The effects of climate change create formidable challenges for breeders striving to produce sufficient food quantities in rapidly changing environments. It is therefore critical to investigate the ability of multi-environment genomic prediction (GP) models to predict genomic estimated breeding values (GEBVs) in extreme environments. Exploration of the impact of training set composition on the accuracy of such GEBVs is also essential. Accordingly, we examined the influence of the number of training environments and the use of environmental covariates (ECs) in GS models on four subsets of n = 500 lines of the soybean nested association mapping (SoyNAM) panel grown in nine environments in the US-North Central Region. The ensuing analyses provided insights into the influence of both of these factors for predicting grain yield in the most and the least extreme of these environments. We found that only a subset of the available environments was needed to obtain the highest observed prediction accuracies. The inclusion of ECs in the GP model did not substantially increase prediction accuracies relative to competing models, and instead more often resulted in negative prediction accuracies. Combined with the overall low prediction accuracies for grain yield in the most extreme environment, our findings highlight weaknesses in current GP approaches for prediction in extreme environments, and point to specific areas on which to focus future research efforts., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Widener, Graef, Lipka and Jarquin.)

Engineering three-dimensional (3D) sensible tissue constructs, along with the complex microarchitecture wiring of the sensory nervous system, has been an ongoing challenge in the tissue engineering field. By combining 3D bioprinting and human pluripotent stem cell (hPSC) technologies, sensible tissue constructs could be engineered in a rapid, precise, and controllable manner to replicate 3D microarchitectures and mechanosensory functionalities of the native sensory tissue (e.g. response to external stimuli). Here, we introduce a biofabrication approach to create complex 3D microarchitecture wirings. We develop an hPSC-sensory neuron (SN) laden bioink using highly purified and functional SN populations to 3D bioprint microarchitecture wirings that demonstrate responsiveness to warm/cold sense-inducing chemicals and mechanical stress. Specifically, we tailor a conventional differentiation strategy to our purification method by utilizing p75 cell surface marker and DAPT treatment along with neuronal growth factors in order to selectively differentiate neural crest cells into SNs. To create spatial resolution in 3D architectures and grow SNs in custom patterns and directions, an induced pluripotent stem cell (iPSC)-SN-laden gelatin bioink was printed on laminin-coated substrates using extrusion-based bioprinting technique. Then the printed constructs were covered with a collagen matrix that guided SNs growing in the printed micropattern. Using a sacrificial bioprinting technique, the iPSC-SNs were seeded into the hollow microchannels created by sacrificial gelatin ink printed in the gelatin methacryloyl supporting bath, thereby demonstrating controllability over axon guidance in curved lines up to several tens of centimeters in length on 2D substrates and in straight microchannels in 3D matrices. Therefore, this biofabrication approach could be amenable to incorporate sensible SN networks into the engineered skin equivalents, regenerative skin implants, and augmented somatosensory neuro-prosthetics that have the potential to regenerate sensible functions by connecting host neuron systems in injured areas., (© 2021 IOP Publishing Ltd.)

Genomic prediction provides an efficient alternative to conventional phenotypic selection for developing improved cultivars with desirable characteristics. New and improved methods to genomic prediction are continually being developed that attempt to deal with the integration of data types beyond genomic information. Modern automated weather systems offer the opportunity to capture continuous data on a range of environmental parameters at specific field locations. In principle, this information could characterize training and target environments and enhance predictive ability by incorporating weather characteristics as part of the genotype-by-environment (G×E) interaction component in prediction models. We assessed the usefulness of including weather data variables in genomic prediction models using a naïve environmental kinship model across 30 environments comprising the Genomes to Fields (G2F) initiative in 2014 and 2015. Specifically four different prediction scenarios were evaluated (i) tested genotypes in observed environments; (ii) untested genotypes in observed environments; (iii) tested genotypes in unobserved environments; and (iv) untested genotypes in unobserved environments. A set of 1,481 unique hybrids were evaluated for grain yield. Evaluations were conducted using five different models including main effect of environments; general combining ability (GCA) effects of the maternal and paternal parents modeled using the genomic relationship matrix; specific combining ability (SCA) effects between maternal and paternal parents; interactions between genetic (GCA and SCA) effects and environmental effects; and finally interactions between the genetics effects and environmental covariates. Incorporation of the genotype-by-environment interaction term improved predictive ability across all scenarios. However, predictive ability was not improved through inclusion of naive environmental covariates in G×E models. More research should be conducted to link the observed weather conditions with important physiological aspects in plant development to improve predictive ability through the inclusion of weather data., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Jarquin, de Leon, Romay, Bohn, Buckler, Ciampitti, Edwards, Ertl, Flint-Garcia, Gore, Graham, Hirsch, Holland, Hooker, Kaeppler, Knoll, Lee, Lawrence-Dill, Lynch, Moose, Murray, Nelson, Rocheford, Schnable, Schnable, Smith, Springer, Thomison, Tuinstra, Wisser, Xu, Yu and Lorenz.)

Anther extrusion (AE) is the most important male floral trait for hybrid wheat seed production. AE is a complex quantitative trait that is difficult to phenotype reliably in field experiments not only due to high genotype-by-environment effects but also due to the short expression window in the field condition. In this study, we conducted a genome-wide association scan (GWAS) and explored the possibility of applying genomic prediction (GP) for AE in the CIMMYT hybrid wheat breeding program. An elite set of male lines ( n = 603) were phenotype for anther count (AC) and anther visual score (VS) across three field experiments in 2017-2019 and genotyped with the 20K Infinitum is elect SNP array. GWAS produced five marker trait associations with small effects. For GP, the main effects of lines (L), environment (E), genomic (G) and pedigree relationships (A), and their interaction effects with environments were used to develop seven statistical models of incremental complexity. The base model used only L and E, whereas the most complex model included L, E, G, A, and G × E and A × E. These models were evaluated in three cross-validation scenarios (CV0, CV1, and CV2). In cross-validation CV0, data from two environments were used to predict an untested environment; in random cross-validation CV1, the test set was never evaluated in any environment; and in CV2, the genotypes in the test set were evaluated in only a subset of environments. The prediction accuracies ranged from -0.03 to 0.74 for AC and -0.01 to 0.54 for VS across different models and CV schemes. For both traits, the highest prediction accuracies with low variance were observed in CV2, and inclusion of the interaction effects increased prediction accuracy for AC only. In CV0, the prediction accuracy was 0.73 and 0.45 for AC and VS, respectively, indicating the high reliability of across environment prediction. Genomic prediction appears to be a very reliable tool for AE in hybrid wheat breeding. Moreover, high prediction accuracy in CV0 demonstrates the possibility of implementing genomic selection across breeding cycles in related germplasm, aiding the rapid breeding cycle., Competing Interests: YM was employed by company Syngenta France S.A.S, France. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2020 Adhikari, Basnet, Crossa, Dreisigacker, Camarillo, Bhati, Jarquin, Manes and Ibrahim.)

Key Message: Comparative assessment identified naïve interaction model, and naïve and informed interaction GS models suitable for achieving higher prediction accuracy in groundnut keeping in mind the high genotype × environment interaction for complex traits. Genomic selection (GS) can be an efficient and cost-effective breeding approach which captures both small- and large-effect genetic factors and therefore promises to achieve higher genetic gains for complex traits such as yield and oil content in groundnut. A training population was constituted with 340 elite lines followed by genotyping with 58 K 'Axiom&#95;Arachis' SNP array and phenotyping for key agronomic traits at three locations in India. Four GS models were tested using three different random cross-validation schemes (CV0, CV1 and CV2). These models are: (1) model 1 (M1 = E + L) which includes the main effects of environment (E) and line (L); (2) model 2 (M2 = E + L + G) which includes the main effects of markers (G) in addition to E and L; (3) model 3 (M3 = E + L + G + GE), a naïve interaction model; and (4) model 4 (E + L + G + LE + GE), a naïve and informed interaction model. Prediction accuracy estimated for four models indicated clear advantage of the inclusion of marker information which was reflected in better prediction accuracy achieved with models M2, M3 and M4 as compared to M1 model. High prediction accuracies (> 0.600) were observed for days to 50% flowering, days to maturity, hundred seed weight, oleic acid, rust@90 days, rust@105 days and late leaf spot@90 days, while medium prediction accuracies (0.400-0.600) were obtained for pods/plant, shelling  %, and total yield/plant. Assessment of comparative prediction accuracy for different GS models to perform selection for untested genotypes, and unobserved and unevaluated environments provided greater insights on potential application of GS breeding in groundnut.

Genomic selection (GS) has proven to be an efficient tool for predicting crop-rank performance of untested genotypes; however, when the traits have intermediate optima (phenology stages), this implementation might not be the most convenient. GS might deliver high-rank correlations but incurring in serious bias. Days to heading (DTH) is a crucial development stage in rice for regional adaptability with a significant impact on yield potential. The objective of this research consisted in develop a novel method that accurately predicts time-related traits such as DTH in unobserved environments. For this, we propose an implementation that incorporates day length information (DL) in the prediction process for two relevant scenarios: CV0, predicting tested genotypes in unobserved environments (C method); and CV00, predicting untested genotypes in unobserved environments (CB method). The use of DL has advantages over weather data since it can be determined in advance just by knowing the location and planting date. The proposed methods showed that DL information significantly helps to improve the predictive ability of DTH in unobserved environments. Under CV0, the C method returned a root-mean-square error (RMSE) of 3.9 days, a Pearson correlation (PC) of 0.98 and the differences between the predicted and observed environmental means (EMD) ranged between -4.95 and 4.67 days. For CV00, the CB method returned an RMSE of 7.3 days, a PC of 0.93 and the EMD ranged between -6.4 and 4.1 days while the conventional GS implementation produced an RMSE of 18.1 days, a PC of 0.41 and the EMD ranged between -31.5 and 28.7 days.

"Sparse testing" refers to reduced multi-environment breeding trials in which not all genotypes of interest are grown in each environment. Using genomic-enabled prediction and a model embracing genotype × environment interaction (GE), the non-observed genotype-in-environment combinations can be predicted. Consequently, the overall costs can be reduced and the testing capacities can be increased. The accuracy of predicting the unobserved data depends on different factors including (1) how many genotypes overlap between environments, (2) in how many environments each genotype is grown, and (3) which prediction method is used. In this research, we studied the predictive ability obtained when using a fixed number of plots and different sparse testing designs. The considered designs included the extreme cases of (1) no overlap of genotypes between environments, and (2) complete overlap of the genotypes between environments. In the latter case, the prediction set fully consists of genotypes that have not been tested at all. Moreover, we gradually go from one extreme to the other considering (3) intermediates between the two previous cases with varying numbers of different or non-overlapping (NO)/overlapping (O) genotypes. The empirical study is built upon two different maize hybrid data sets consisting of different genotypes crossed to two different testers (T1 and T2) and each data set was analyzed separately. For each set, phenotypic records on yield from three different environments are available. Three different prediction models were implemented, two main effects models ( M1 and M2 ), and a model ( M3) including GE. The results showed that the genome-based model including GE ( M3 ) captured more phenotypic variation than the models that did not include this component. Also, M3 provided higher prediction accuracy than models M1 and M2 for the different allocation scenarios. Reducing the size of the calibration sets decreased the prediction accuracy under all allocation designs with M3 being the less affected model; however, using the genome-enabled models ( i.e. , M2 and M3 ) the predictive ability is recovered when more genotypes are tested across environments. Our results indicate that a substantial part of the testing resources can be saved when using genome-based models including GE for optimizing sparse testing designs., (Copyright © 2020 Jarquin et al.)

Genome-wide association mapping identifies quantitative trait loci (QTL) that influence the mean differences between the marker genotypes for a given trait. While most loci influence the mean value of a trait, certain loci, known as variance heterogeneity QTL (vQTL) determine the variability of the trait instead of the mean trait value (mQTL). In the present study, we performed a variance heterogeneity genome-wide association study (vGWAS) for grain cadmium (Cd) concentration in bread wheat. We used double generalized linear model and hierarchical generalized linear model to identify vQTL associated with grain Cd. We identified novel vQTL regions on chromosomes 2A and 2B that contribute to the Cd variation and loci that affect both mean and variance heterogeneity (mvQTL) on chromosome 5A. In addition, our results demonstrated the presence of epistatic interactions between vQTL and mvQTL, which could explain variance heterogeneity. Overall, we provide novel insights into the genetic architecture of grain Cd concentration and report the first application of vGWAS in wheat. Moreover, our findings indicated that epistasis is an important mechanism underlying natural variation for grain Cd concentration., (© 2020 The Authors. The Plant Genome published by Wiley Periodicals, Inc. on behalf of Crop Science Society of America.)

Objectives: Advanced tools and resources are needed to efficiently and sustainably produce food for an increasing world population in the context of variable environmental conditions. The maize genomes to fields (G2F) initiative is a multi-institutional initiative effort that seeks to approach this challenge by developing a flexible and distributed infrastructure addressing emerging problems. G2F has generated large-scale phenotypic, genotypic, and environmental datasets using publicly available inbred lines and hybrids evaluated through a network of collaborators that are part of the G2F's genotype-by-environment (G × E) project. This report covers the public release of datasets for 2014-2017., Data Description: Datasets include inbred genotypic information; phenotypic, climatic, and soil measurements and metadata information for each testing location across years. For a subset of inbreds in 2014 and 2015, yield component phenotypes were quantified by image analysis. Data released are accompanied by README descriptions. For genotypic and phenotypic data, both raw data and a version without outliers are reported. For climatic data, a version calibrated to the nearest airport weather station and a version without outliers are reported. The 2014 and 2015 datasets are updated versions from the previously released files [1] while 2016 and 2017 datasets are newly available to the public.