Genome-based establishment of a high-yielding heterotic pattern for hybrid wheat breeding

Phenomic Selection for Hybrid Rapeseed Breeding

Phenomic selection is a recent approach suggested as a low-cost, high-throughput alternative to genomic selection. Instead of using genetic markers, it employs spectral data to predict complex traits using equivalent statistical models. Phenomic selection has been shown to outperform genomic selection when using spectral data that was obtained within the same generation as the traits that were predicted. However, for hybrid breeding, the key question is whether spectral data from parental genotypes can be used to effectively predict traits in the hybrid generation. Here, we aimed to evaluate the potential of phenomic selection for hybrid rapeseed breeding. We performed predictions for various traits in a structured population of 410 test hybrids, grown in multiple environments, using near-infrared spectroscopy data obtained from harvested seeds of both the hybrids and their parental lines with different linear and nonlinear models. We found that phenomic selection within the hybrid generation outperformed genomic selection for seed yield and plant height, even when spectral data was collected at single locations, while being less affected by population structure. Furthermore, we demonstrate that phenomic prediction across generations is feasible, and selecting hybrids based on spectral data obtained from parental genotypes is competitive with genomic selection. We conclude that phenomic selection is a promising approach for rapeseed breeding that can be easily implemented without any additional costs or efforts as near-infrared spectroscopy is routinely assessed in rapeseed breeding.

Deciphering the differential expression patterns of yield-related negative regulators in hexaploid wheat cultivars and hybrids at different growth stages

BACKGROUND

Hexaploid bread wheat underwent a series of polyploidization events through interspecific hybridizations that conferred adaptive plasticity and resulted in duplication and neofunctionalization of major agronomic genes. The genetic architecture of polyploid wheat not only confers adaptive plasticity but also offers huge genetic diversity. However, the contribution of different gene copies (homeologs) encoded from different subgenomes (A, B, D) at different growth stages remained unexplored.

METHODS

In this study, hybrid of elite cultivars of wheat were developed via reciprocal crosses (cytoplasm swapping) and phenotypically evaluated. We assessed differential expression profiles of yield-related negative regulators in these cultivars and their F1 hybrids and identified various cis-regulatory signatures by employing bioinformatics tools. Furthermore, the preferential expression patterns of the syntenic triads encoded from A, B, and D subgenomes were assessed to decipher their functional redundancy at six different growth stages.

RESULTS

Hybrid progenies showed better heterosis such as up to 17% increase in the average number of grains and up to 50% increase in average thousand grains weight as compared to mid-parents. Based on the expression profiling, our results indicated significant dynamic transcriptional expression patterns, portraying the different homeolog-dominance at the same stage in the different cultivars and their hybrids. Albeit belonging to same syntenic triads, a dynamic trend was observed in the regulatory signatures of these genes that might be influencing their expression profiles.

CONCLUSION

These findings can substantially contribute and provide insights for the selective introduction of better cultivars into traditional and hybrid breeding programs which can be harnessed for the improvement of future wheat.

Genomic selection in plant breeding: Key factors shaping two decades of progress

Genomic selection, the application of genomic prediction (GP) models to select candidate individuals, has significantly advanced in the past two decades, effectively accelerating genetic gains in plant breeding. This article provides a holistic overview of key factors that have influenced GP in plant breeding during this period. We delved into the pivotal roles of training population size and genetic diversity, and their relationship with the breeding population, in determining GP accuracy. Special emphasis was placed on optimizing training population size. We explored its benefits and the associated diminishing returns beyond an optimum size. This was done while considering the balance between resource allocation and maximizing prediction accuracy through current optimization algorithms. The density and distribution of single-nucleotide polymorphisms, level of linkage disequilibrium, genetic complexity, trait heritability, statistical machine-learning methods, and non-additive effects are the other vital factors. Using wheat, maize, and potato as examples, we summarize the effect of these factors on the accuracy of GP for various traits. The search for high accuracy in GP-theoretically reaching one when using the Pearson's correlation as a metric-is an active research area as yet far from optimal for various traits. We hypothesize that with ultra-high sizes of genotypic and phenotypic datasets, effective training population optimization methods and support from other omics approaches (transcriptomics, metabolomics and proteomics) coupled with deep-learning algorithms could overcome the boundaries of current limitations to achieve the highest possible prediction accuracy, making genomic selection an effective tool in plant breeding.

Dynamic Phytomeric Growth Contributes to Local Adaptation in Barley

Vascular plants have segmented body axes with iterative nodes and internodes. Appropriate node initiation and internode elongation are fundamental to plant fitness and crop yield; however, how these events are spatiotemporally coordinated remains elusive. We show that in barley (Hordeum vulgare L.), selections during domestication have extended the apical meristematic phase to promote node initiation, but constrained subsequent internode elongation. In both vegetative and reproductive phases, internode elongation displays a dynamic proximal-distal gradient, and among subpopulations of domesticated barleys worldwide, node initiation and proximal internode elongation are associated with latitudinal and longitudinal gradients, respectively. Genetic and functional analyses suggest that, in addition to their converging roles in node initiation, flowering-time genes have been repurposed to specify the timing and duration of internode elongation. Our study provides an integrated view of barley node initiation and internode elongation and suggests that plant architecture should be recognized as a collection of dynamic phytomeric units in the context of crop adaptive evolution.

Simulation-based establishment of base pools for a hybrid breeding program in winter rapeseed

KEY MESSAGE

Simulation planned pre-breeding can increase the efficiency of starting a hybrid breeding program. Starting a hybrid breeding program commonly comprises a grouping of the initial germplasm in two pools and subsequent selection on general combining ability. Investigations on pre-breeding steps before starting the selection on general combining ability are not available. Our goals were (1) to use computer simulations on the basis of DNA markers and testcross data to plan crosses that separate genetically two initial germplasm pools of rapeseed, (2) to carry out the planned crosses, and (3) to verify experimentally the pool separation as well as the increase in testcross performance. We designed a crossing program consisting of four cycles of recombination. In each cycle, the experimentally generated material was used to plan the subsequent crossing cycle with computer simulations. After finishing the crossing program, the initially overlapping pools were clearly separated in principal coordinate plots. Doubled haploid lines derived from the material of crossing cycles 1 and 2 showed an increase in relative testcross performance for yield of about 5% per cycle. We conclude that simulation-designed pre-breeding crossing schemes, that were carried out before the general combining ability-based selection of a newly started hybrid breeding program, can save time and resources, and in addition conserve more of the initial genetic variation than a direct start of a hybrid breeding program with general combining ability-based selection.

Heterosis in crop improvement

Heterosis, also known as hybrid vigor, is the phenomenon wherein a progeny exhibits superior traits relative to one or both parents. In terms of crop breeding, this usually refers to the yield advantage of F1 hybrids over both inbred parents. The development of high-yielding hybrid cultivars across a wider range of crops is key to meeting future food demands. However, conventional hybrid breeding strategies are proving to be exceptionally challenging to apply commercially in many self-pollinating crops, particularly wheat and barley. Currently in these crops, the relative performance advantage of hybrids over inbred line cultivars does not outweigh the cost of hybrid seed production. Here, we review the genetic basis of heterosis, discuss the challenges in hybrid breeding, and propose a strategy to recruit multiple heterosis-associated genes to develop lines with improved agronomic characteristics. This strategy leverages modern genetic engineering tools to synthesize supergenes by fusing multiple heterotic alleles across multiple heterosis-associated loci. We outline a plan to assess the feasibility of this approach to improve line performance using barley (Hordeum vulgare) as the model self-pollinating crop species, and a few heterosis-associated genes. The proposed method can be applied to all crops for which heterotic gene combinations can be identified.

Data-driven approaches to improve water-use efficiency and drought resistance in crop plants

With the increasing population, there lies a pressing demand for food, feed and fibre, while the changing climatic conditions pose severe challenges for agricultural production worldwide. Water is the lifeline for crop production; thus, enhancing crop water-use efficiency (WUE) and improving drought resistance in crop varieties are crucial for overcoming these challenges. Genetically-driven improvements in yield, WUE and drought tolerance traits can buffer the worst effects of climate change on crop production in dry areas. While traditional crop breeding approaches have delivered impressive results in increasing yield, the methods remain time-consuming and are often limited by the existing allelic variation present in the germplasm. Significant advances in breeding and high-throughput omics technologies in parallel with smart agriculture practices have created avenues to dramatically speed up the process of trait improvement by leveraging the vast volumes of genomic and phenotypic data. For example, individual genome and pan-genome assemblies, along with transcriptomic, metabolomic and proteomic data from germplasm collections, characterised at phenotypic levels, could be utilised to identify marker-trait associations and superior haplotypes for crop genetic improvement. In addition, these omics approaches enable the identification of genes involved in pathways leading to the expression of a trait, thereby providing an understanding of the genetic, physiological and biochemical basis of trait variation. These data-driven gene discoveries and validation approaches are essential for crop improvement pipelines, including genomic breeding, speed breeding and gene editing. Herein, we provide an overview of prospects presented using big data-driven approaches (including artificial intelligence and machine learning) to harness new genetic gains for breeding programs and develop drought-tolerant crop varieties with favourable WUE and high-yield potential traits.

Maximization through optimization? On the relationship between hybrid performance and parental genetic distance

Heterosis is the improved performance of hybrids compared with their parental components and is widely exploited in agriculture. According to quantitative genetic theory, genetic distance between parents at heterotic quantitative trait loci is required for heterosis, but how heterosis varies with genetic distance has remained elusive, despite intensive research on the topic. Experimental studies have often found a positive association between heterosis and genetic distance that, however, varied in strength. Most importantly, it has remained unclear whether heterosis increases continuously with genetic distance or whether there is an optimum genetic distance after which heterosis declines again. Here, we revisit the relationship between heterosis and genetic distance and provide perspectives on how to maximize heterosis and hybrid performance in breeding, as well as the consequences for the design of heterotic groups and the utilization of more exotic material and genetic resources.

Genomic prediction in hybrid breeding: I. Optimizing the training set design

Training sets produced by maximizing the number of parent lines, each involved in one cross, had the highest prediction accuracy for H0 hybrids, but lowest for H1 and H2 hybrids. Genomic prediction holds great promise for hybrid breeding but optimum composition of the training set (TS) as determined by the number of parents (nTS) and crosses per parent (c) has received little attention. Our objective was to examine prediction accuracy ([Formula: see text]) of GCA for lines used as parents of the TS (I1 lines) or not (I0 lines), and H0, H1 and H2 hybrids, comprising crosses of type I0 × I0, I1 × I0 and I1 × I1, respectively, as function of nTS and c. In the theory, we developed estimates for [Formula: see text] of GBLUPs for hybrids: (i)[Formula: see text] based on the expected prediction accuracy, and (ii) [Formula: see text] based on [Formula: see text] of GBLUPs of GCA and SCA effects. In the simulation part, hybrid populations were generated using molecular data from two experimental maize data sets. Additive and dominance effects of QTL borrowed from literature were used to simulate six scenarios of traits differing in the proportion (τSCA = 1%, 6%, 22%) of SCA variance in σG2 and heritability (h2 = 0.4, 0.8). Values of [Formula: see text] and [Formula: see text] closely agreed with [Formula: see text] for hybrids. For given size NTS = nTS × c of TS, [Formula: see text] of H0 hybrids and GCA of I0 lines was highest for c = 1. Conversely, for GCA of I1 lines and H1 and H2 hybrids, c = 1 yielded lowest [Formula: see text] with concordant results across all scenarios for both data sets. In view of these opposite trends, the optimum choice of c for maximizing selection response across all types of hybrids depends on the size and resources of the breeding program.

A statistical package for evaluation of hybrid performance in plant breeding via genomic selection

Hybrid breeding employs heterosis, which could potentially improve the yield and quality of a crop. Genomic selection (GS) is a promising approach for the selection of quantitative traits in plant breeding. The main objectives of this study are to (i) propose a GS-based approach to identify potential parental lines and superior hybrid combinations from a breeding population, which is composed of hybrids produced by a half diallel mating design; (ii) develop a software package for users to carry out the proposed approach. An R package, designated EHPGS, was generated to facilitate the employment of the genomic best linear unbiased model considering additive plus dominance marker effects for the hybrid performance evaluation. The R package contains a Bayesian statistical algorithm for calculating genomic estimated breeding value (GEBVs), GEBV-based specific combining ability, general combining ability, mid-parent heterosis, and better-parent heterosis. Three datasets that have been published in literature, including pumpkin (Cucurbita maxima), maize (Zea mays), and wheat (Triticum aestivum L.), were reanalyzed to illustrate the use of EHPGS.

Stacked ensembles on basis of parentage information can predict hybrid performance with an accuracy comparable to marker-based GBLUP

Testcross factorials in newly established hybrid breeding programs are often highly unbalanced, incomplete, and characterized by predominance of special combining ability (SCA) over general combining ability (GCA). This results in a low efficiency of GCA-based selection. Machine learning algorithms might improve prediction of hybrid performance in such testcross factorials, as they have been successfully applied to find complex underlying patterns in sparse data. Our objective was to compare the prediction accuracy of machine learning algorithms to that of GCA-based prediction and genomic best linear unbiased prediction (GBLUP) in six unbalanced incomplete factorials from hybrid breeding programs of rapeseed, wheat, and corn. We investigated a range of machine learning algorithms with three different types of predictor variables: (a) information on parentage of hybrids, (b) in addition hybrid performance of crosses of the parental lines with other crossing partners, and (c) genotypic marker data. In two highly incomplete and unbalanced factorials from rapeseed, in which the SCA variance contributed considerably to the genetic variance, stacked ensembles of gradient boosting machines based on parentage information outperformed GCA prediction. The stacked ensembles increased prediction accuracy from 0.39 to 0.45, and from 0.48 to 0.54 compared to GCA prediction. The prediction accuracy reached by stacked ensembles without marker data reached values comparable to those of GBLUP that requires marker data. We conclude that hybrid prediction with stacked ensembles of gradient boosting machines based on parentage information is a promising approach that is worth further investigations with other data sets in which SCA variance is high.

Transcriptomic and Physiological Analyses of Two Rice Restorer Lines under Different Nitrogen Supplies Provide Novel Insights into Hybrid Rice Breeding

Improving plant nitrogen-use efficiency (NUE) has great significance for various crops, particularly in hybrid breeding. Reducing nitrogen inputs is key to achieving sustainable rice production and mitigating environmental problems. In this study, we analyzed the transcriptomic and physiological changes in two indica restorer lines (Nanhui511 [NH511] and Minghui23 [MH23]) under high nitrogen (HN) and low nitrogen (LN) conditions. Compared to MH23, NH511 was more sensitive to different nitrogen supplies and exhibited higher nitrogen uptake and NUE under HN conditions by increasing lateral root and tiller numbers in the seedling and maturation stages, respectively. NH511 also exhibited a lower survival rate than MH23 when planted in a chlorate-containing hydroponic solution, indicating its HN uptake ability under different nitrogen-supply conditions. Transcriptomic analysis showed that NH511 has 2456 differentially expressed genes, whereas MH23 had only 266. Furthermore, these genes related to nitrogen utilization showed differential expression in NH511 under HN conditions, while the opposite was observed in MH23. Our findings revealed that NH511 could be regarded as elite rice and used for breeding high-NUE restorer lines by regulating and integrating nitrogen-utilization genes, which provides novel insights for the cultivation of high-NUE hybrid rice.

Genome-wide identification of nitrate-responsive microRNAs by small RNA sequencing in the rice restorer cultivar Nanhui 511

Rice productivity relies heavily on nitrogen fertilization, and improving nitrogen use efficiency (NUE) is important for hybrid rice breeding. Reducing nitrogen inputs is the key to achieving sustainable rice production and reducing environmental problems. Here, we analyzed the genome-wide transcriptomic changes in microRNAs (miRNAs) in the indica rice restorer cultivar Nanhui 511 (NH511) under high (HN) and low nitrogen (LN) conditions. The results showed that NH511 is sensitive to nitrogen supplies and HN conditions promoted the growth its lateral roots at the seedling stage. Furthermore, we identified 483 known miRNAs and 128 novel miRNAs by small RNA sequencing in response to nitrogen in NH511. We also detected 100 differentially expressed genes (DEGs), including 75 upregulated and 25 downregulated DEGs, under HN conditions. Among these DEGs, 43 miRNAs that exhibited a 2-fold change in their expression were identified in response to HN conditions, including 28 upregulated and 15 downregulated genes. Additionally, some differentially expressed miRNAs were further validated by qPCR analysis, which showed that miR443, miR1861b, and miR166k-3p were upregulated, whereas miR395v and miR444b.1 were downregulated under HN conditions. Moreover, the degradomes of possible target genes for miR166k-3p and miR444b.1 and expression variations were analyzed by qPCR at different time points under HN conditions. Our findings revealed comprehensive expression profiles of miRNAs responsive to HN treatments in an indica rice restorer cultivar, which advances our understanding of the regulation of nitrogen signaling mediated by miRNAs and provides novel data for high-NUE hybrid rice cultivation.

Large-scale genomic and transcriptomic profiles of rice hybrids reveal a core mechanism underlying heterosis

BACKGROUND

Heterosis is widely used in agriculture. However, its molecular mechanisms are still unclear in plants. Here, we develop, sequence, and record the phenotypes of 418 hybrids from crosses between two testers and 265 rice varieties from a mini-core collection.

RESULTS

Phenotypic analysis shows that heterosis is dependent on genetic backgrounds and environments. By genome-wide association study of 418 hybrids and their parents, we find that nonadditive QTLs are the main genetic contributors to heterosis. We show that nonadditive QTLs are more sensitive to the genetic background and environment than additive ones. Further simulations and experimental analysis support a novel mechanism, homo-insufficiency under insufficient background (HoIIB), underlying heterosis. We propose heterosis in most cases is not due to heterozygote advantage but homozygote disadvantage under the insufficient genetic background.

CONCLUSION

The HoIIB model elucidates that genetic background insufficiency is the intrinsic mechanism of background dependence, and also the core mechanism of nonadditive effects and heterosis. This model can explain most known hypotheses and phenomena about heterosis, and thus provides a novel theory for hybrid rice breeding in future.

Development of the maize 5.5K loci panel for genomic prediction through genotyping by target sequencing

Genotyping platforms are important for genetic research and molecular breeding. In this study, a low-density genotyping platform containing 5.5K SNP markers was successfully developed in maize using genotyping by target sequencing (GBTS) technology with capture-in-solution. Two maize populations (Pop1 and Pop2) were used to validate the GBTS panel for genetic and molecular breeding studies. Pop1 comprised 942 hybrids derived from 250 inbred lines and four testers, and Pop2 contained 540 hybrids which were generated from 123 new-developed inbred lines and eight testers. The genetic analyses showed that the average polymorphic information content and genetic diversity values ranged from 0.27 to 0.38 in both populations using all filtered genotyping data. The mean missing rate was 1.23% across populations. The Structure and UPGMA tree analyses revealed similar genetic divergences (76-89%) in both populations. Genomic prediction analyses showed that the prediction accuracy of reproducing kernel Hilbert space (RKHS) was slightly lower than that of genomic best linear unbiased prediction (GBLUP) and three Bayesian methods for general combining ability of grain yield per plant and three yield-related traits in both populations, whereas RKHS with additive effects showed superior advantages over the other four methods in Pop1. In Pop1, the GBLUP and three Bayesian methods with additive-dominance model improved the prediction accuracies by 4.89-134.52% for the four traits in comparison to the additive model. In Pop2, the inclusion of dominance did not improve the accuracy in most cases. In general, low accuracies (0.33-0.43) were achieved for general combing ability of the four traits in Pop1, whereas moderate-to-high accuracies (0.52-0.65) were observed in Pop2. For hybrid performance prediction, the accuracies were moderate to high (0.51-0.75) for the four traits in both populations using the additive-dominance model. This study suggests a reliable genotyping platform that can be implemented in genomic selection-assisted breeding to accelerate maize new cultivar development and improvement.

Genome-wide identification of genes enabling accurate prediction of hybrid performance from parents across environments and populations for gene-based breeding in maize

Accurate prediction of hybrid offspring complex trait phenotype from parents is paramount to enhanced plant breeding, animal breeding, and human medicine. Here we report genome-wide identification of genes enabling accurate prediction of hybrid offspring complex traits from parents using maize grain yield as the target trait. We identified 181 ZmF1GY genes enabling prediction of maize (Zea mays L.) F₁ hybrid grain yield from parents and tested their utility and efficiency for predicting F₁ hybrid grain yields from parents using their expressions, genic SNPs, and number of favorable alleles (NFAs), respectively. The ZmF1GY genes predicted hybrid grain yields from parents at an accuracy of 0.86, presented by correlation coefficient between predicted and observed phenotypes, within an environment, 0.74 across environments, and 0.64 across populations, outperforming genomic prediction by 27-406%, 23%, and 40%, respectively. Furthermore, we identified nine of the ZmF1GY genes containing SNPs or InDels in parents that increased or decreased hybrid grain yields by 14-46%. When the NFAs of these nine ZmF1GY genes were used for hybrid grain yield prediction from parents, they predicted hybrid grain yields at an accuracy of 0.79, outperforming genomic prediction by 21% that was based on up to tens of thousands of genome-wide SNPs. These results demonstrate the feasibility of developing a gene toolkit for a species enabling gene-based breeding across environments and populations that is much more powerful and efficient than current breeding, thereby helping secure the world's food production. The methodology is applicable to all crops, livestock, and humans.

Improving Yield and Yield Stability in Winter Rye by Hybrid Breeding

Rye is the only cross-pollinating small-grain cereal. The unique reproduction biology results in an exceptional complexity concerning genetic improvement of rye by breeding. Rye is a close relative of wheat and has a strong adaptation potential that refers to its mating system, making this overlooked cereal readily adjustable to a changing environment. Rye breeding addresses the emerging challenges of food security associated with climate change. The systematic identification, management, and use of its valuable natural diversity became a feasible option in outbreeding rye only following the establishment of hybrid breeding late in the 20th century. In this article, we review the most recent technological advances to improve yield and yield stability in winter rye. Based on recently released reference genome sequences, SMART breeding approaches are described to counterbalance undesired linkage drag effects of major restorer genes on grain yield. We present the development of gibberellin-sensitive semidwarf hybrids as a novel plant breeding innovation based on an approach that is different from current methods of increasing productivity in rye and wheat. Breeding of new rye cultivars with improved performance and resilience is indispensable for a renaissance of this healthy minor cereal as a homogeneous commodity with cultural relevance in Europe that allows for comparatively smooth but substantial complementation of wheat with rye-based diets, supporting the necessary restoration of the balance between human action and nature.

Epistasis Activation Contributes Substantially to Heterosis in Temperate by Tropical Maize Hybrids

Epistasis strongly affects the performance of superior maize hybrids. In this study, a multiple-hybrid population, consisting of three hybrid maize sets with varied interparental divergence, was generated by crossing 28 temperate and 23 tropical inbred lines with diverse genetic backgrounds. We obtained 1,154 tested hybrids. Among these tested hybrids, heterosis increased steadily as the heterotic genetic distance increased. Mid-parent heterosis was significantly higher in the temperate by tropical hybrids than in the temperate by temperate hybrids. Genome-wide prediction and association mapping was performed for grain weight per plant (GWPP) and days to silking (DTS) using 20K high-quality SNPs, showing that epistatic effects played a more prominent role than dominance effects in temperate by tropical maize hybrids. A total of 33 and 420 epistatic QTL were identified for GWPP and DTS, respectively, in the temperate by tropical hybrids. Protein-protein interaction network and gene-set enrichment analyses showed that epistatic genes were involved in protein interactions, which play an important role in photosynthesis, biological transcription pathways, and protein synthesis. We showed that the interaction of many minor-effect genes in the hybrids could activate the transcription activators of epistatic genes, resulting in a cascade of amplified yield heterosis. The multiple-hybrid population design enhanced our understanding of heterosis in maize, providing an insight into the acceleration of hybrid maize breeding by activating epistatic effects.

Capturing Wheat Phenotypes at the Genome Level

Recent technological advances in next-generation sequencing (NGS) technologies have dramatically reduced the cost of DNA sequencing, allowing species with large and complex genomes to be sequenced. Although bread wheat (Triticum aestivum L.) is one of the world's most important food crops, efficient exploitation of molecular marker-assisted breeding approaches has lagged behind that achieved in other crop species, due to its large polyploid genome. However, an international public-private effort spanning 9 years reported over 65% draft genome of bread wheat in 2014, and finally, after more than a decade culminated in the release of a gold-standard, fully annotated reference wheat-genome assembly in 2018. Shortly thereafter, in 2020, the genome of assemblies of additional 15 global wheat accessions was released. As a result, wheat has now entered into the pan-genomic era, where basic resources can be efficiently exploited. Wheat genotyping with a few hundred markers has been replaced by genotyping arrays, capable of characterizing hundreds of wheat lines, using thousands of markers, providing fast, relatively inexpensive, and reliable data for exploitation in wheat breeding. These advances have opened up new opportunities for marker-assisted selection (MAS) and genomic selection (GS) in wheat. Herein, we review the advances and perspectives in wheat genetics and genomics, with a focus on key traits, including grain yield, yield-related traits, end-use quality, and resistance to biotic and abiotic stresses. We also focus on reported candidate genes cloned and linked to traits of interest. Furthermore, we report on the improvement in the aforementioned quantitative traits, through the use of (i) clustered regularly interspaced short-palindromic repeats/CRISPR-associated protein 9 (CRISPR/Cas9)-mediated gene-editing and (ii) positional cloning methods, and of genomic selection. Finally, we examine the utilization of genomics for the next-generation wheat breeding, providing a practical example of using in silico bioinformatics tools that are based on the wheat reference-genome sequence.

Dissection of the practical soybean breeding pipeline by developing ZDX1, a high-throughput functional array

KEY MESSAGE

We developed the ZDX1 high-throughput functional soybean array for high accuracy evaluation and selection of both parents and progeny, which can greatly accelerate soybean breeding. Microarray technology facilitates rapid, accurate, and economical genotyping. Here, using resequencing data from 2214 representative soybean accessions, we developed the high-throughput functional array ZDX1, containing 158,959 SNPs, covering 90.92% of soybean genes and sites related to important traits. By application of the array, a total of 817 accessions were genotyped, including three subpopulations of candidate parental lines, parental lines and their progeny from practical breeding. The fixed SNPs were identified in progeny, indicating artificial selection during the breeding process. By identifying functional sites of target traits, novel soybean cyst nematode-resistant progeny and maturity-related novel sources were identified by allele combinations, demonstrating that functional sites provide an efficient method for the rapid screening of desirable traits or gene sources. Notably, we found that the breeding index (BI) was a good indicator for progeny selection. Superior progeny were derived from the combination of distantly related parents, with at least one parent having a higher BI. Furthermore, new combinations based on good performance were proposed for further breeding after excluding redundant and closely related parents. Genomic best linear unbiased prediction (GBLUP) analysis was the best analysis method and achieved the highest accuracy in predicting four traits when comparing SNPs in genic regions rather than whole genomic or intergenic SNPs. The prediction accuracy was improved by 32.1% by using progeny to expand the training population. Collectively, a versatile assay demonstrated that the functional ZDX1 array provided efficient information for the design and optimization of a breeding pipeline for accelerated soybean breeding.

Assembly of yield heterosis of an elite rice hybrid is promising by manipulating dominant quantitative trait loci

In the past, rice hybrids with strong heterosis have been obtained empirically, by developing and testing thousands of combinations. Here, we aimed to determine whether heterosis of an elite hybrid could be achieved by manipulating major quantitative trait loci. We used 202 chromosome segment substitution lines from the elite hybrid Shanyou 63 to evaluate single segment heterosis (SSH) of yield per plant and identify heterotic loci. All nine detected heterotic loci acted in a dominant fashion, and no SSH exhibited overdominance. Functional alleles of key yield-related genes Ghd7, Ghd7.1, Hd1, and GS3 were dispersed in both parents. No functional alleles of three investigated genes were expressed at higher levels in the hybrids than in the more desirable parents. A hybrid pyramiding eight heterotic loci in the female parent Zhenshan 97 background had a comparable yield to Shanyou 63 and much higher yield than Zhenshan 97. Five hybrids pyramiding eight or nine heterotic loci in the combined parental genome background showed similar yield performance to that of Shanyou 63. These results suggest that dominance underlying functional complementation is an important contributor to yield heterosis and that heterosis assembly might be successfully promised by manipulating several major dominant heterotic loci.

Improvement of genomic prediction in advanced wheat breeding lines by including additive-by-additive epistasis

Including additive and additive-by-additive epistasis in a NOIA parametrization did not yield orthogonal partitioning of genetic variances, nevertheless, it improved predictive ability in a leave-one-out cross-validation for wheat grain yield. Additive-by-additive epistasis is the principal non-additive genetic effect in inbred wheat lines and is potentially useful for developing cultivars based on total genetic merit; nevertheless, its practical benefits have been highly debated. In this article, we aimed to (i) evaluate the performance of models including additive and additive-by-additive epistatic effects for variance components (VC) estimation of grain yield in a wheat-breeding population, and (ii) to investigate whether including additive-by-additive epistasis in genomic prediction enhance wheat grain yield predictive ability (PA). In total, 2060 sixth-generation (F₆) lines from Nordic Seed A/S breeding company were phenotyped in 21 year-location combinations in Denmark, and genotyped using a 15 K-Illumina-BeadChip. Three models were used to estimate VC and heritability at plot level: (i) "I-model" (baseline), (ii) "I + G_A-model", extending I-model with an additive genomic effect, and (iii) "I + G_A + G_AA-model", extending I + G_A-model with an additive-by-additive genomic effects. The I + G_A-model and I + G_A + G_AA-model were based on the Natural and Orthogonal Interactions Approach (NOIA) parametrization. The I + G_A + G_AA-model failed to achieve orthogonal partition of genetic variances, as revealed by a change in estimated additive variance of I + G_A-model when epistasis was included in the I + G_A + G_AA-model. The PA was studied using leave-one-line-out and leave-one-breeding-cycle-out cross-validations. The I + G_A + G_AA-model increased PA significantly (16.5%) compared to the I + G_A-model in leave-one-line-out cross-validation. However, the improvement due to including epistasis was not observed in leave-one-breeding-cycle-out cross-validation. We conclude that epistatic models can be useful to enhance predictions of total genetic merit. However, even though we used the NOIA parameterization, the variance partition into orthogonal genetic effects was not possible.

The performance of phenomic selection depends on the genetic architecture of the target trait

The phenomic predictive ability depends on the genetic architecture of the target trait, being high for complex traits and low for traits with major QTL. Genomic selection is a powerful tool to assist breeding of complex traits, but a limitation is the costs required for genotyping. Recently, phenomic selection has been suggested, which uses spectral data instead of molecular markers as predictors. It was shown to be competitive with genomic prediction, as it achieved predictive abilities as high or even higher than its genomic counterpart. The objective of this study was to evaluate the performance of phenomic prediction for triticale and the dependency of the predictive ability on the genetic architecture of the target trait. We found that for traits with a complex genetic architecture, like grain yield, phenomic prediction with NIRS data as predictors achieved high predictive abilities and performed better than genomic prediction. By contrast, for mono- or oligogenic traits, for example, yellow rust, marker-based approaches achieved high predictive abilities, while those of phenomic prediction were very low. Compared with molecular markers, the predictive ability obtained using NIRS data was more robust to varying degrees of genetic relatedness between the training and prediction set. Moreover, for grain yield, smaller training sets were required to achieve a similar predictive ability for phenomic prediction than for genomic prediction. In addition, our results illustrate the potential of using field-based spectral data for phenomic prediction. Overall, our result confirmed phenomic prediction as an efficient approach to improve the selection gain for complex traits in plant breeding.

Identification of Heterotic Groups and Patterns Based on Genotypic and Phenotypic Characteristics Among Rice Accessions of Diverse Origins

Identification of the right parental combinations to maximize heterosis is the major goal of hybrid breeding, which could be achieved through identification of heterotic groups. The main objective of this study was to identify promising heterotic groups for future rice breeding programs. A collection of 359 rice genotypes of diverse origins of China and abroad, composed of inbreds, maintainers, restorers, and temperature-sensitive genic male sterile (TGMS) lines were genotyped using 10K SNP chips. The SNP data set was subjected to genomic analyses for estimation of genetic divergence and diversity. Significant variations were observed in the germplasm with the identification of six different genetic groups. These lines were assigned to the genetic groups independent of their origin. Taking an account of commercially used heterotic groups present in each cluster, three cytoplasmic male sterile (CMS) lines and 14 inbred and restorer lines with moderate to high genetic distances selected from five heterotic patterns were crossed and obtained 42 F1 hybrids. A total of 14 hybrids were found with significant maximum mid- and better-parent heterosis, namely, TaifengA × Guang122, TaifengA × Wushansimiao, and TaifengA × Minghui63 for earliness; Guang8A × Huazhan for dwarf stature; and Guang8A × Huanghuzhan-1, TaifengA × Yuexiangzhan, Guang8A × Minhui3301, TianfengA × Guang122, Guang8A × Yahui2115, TianfengA × Huanghuazhan, TianfengA × Minghui63, TianfengA × Minhui3301, TaifengA × Gui99, and Guang8A × Yuenongsimiao for yield and yield-related traits. Mid-parent and better-parent heterotic F1 hybrids were in positive correlation with the genetic distances as that manifested by commercially used heterotic groups, encouraging the use of genotypic data for identification of heterotic groups. Our study provides an informative strategy for the development of early maturing, lodging resistant and high-yielding commercial hybrids and cultivars in future heterosis breeding programs.

Genomic Prediction: Progress and Perspectives for Rice Improvement

Genomic prediction can be a powerful tool to achieve greater rates of genetic gain for quantitative traits if thoroughly integrated into a breeding strategy. In rice as in other crops, the interest in genomic prediction is very strong with a number of studies addressing multiple aspects of its use, ranging from the more conceptual to the more practical. In this chapter, we review the literature on rice (Oryza sativa) and summarize important considerations for the integration of genomic prediction in breeding programs. The irrigated breeding program at the International Rice Research Institute is used as a concrete example on which we provide data and R scripts to reproduce the analysis but also to highlight practical challenges regarding the use of predictions. The adage "To someone with a hammer, everything looks like a nail" describes a common psychological pitfall that sometimes plagues the integration and application of new technologies to a discipline. We have designed this chapter to help rice breeders avoid that pitfall and appreciate the benefits and limitations of applying genomic prediction, as it is not always the best approach nor the first step to increasing the rate of genetic gain in every context.

Multi-omics prediction of oat agronomic and seed nutritional traits across environments and in distantly related populations

Integration of multi-omics data improved prediction accuracies of oat agronomic and seed nutritional traits in multi-environment trials and distantly related populations in addition to the single-environment prediction. Multi-omics prediction has been shown to be superior to genomic prediction with genome-wide DNA-based genetic markers (G) for predicting phenotypes. However, most of the existing studies were based on historical datasets from one environment; therefore, they were unable to evaluate the efficiency of multi-omics prediction in multi-environment trials and distantly related populations. To fill those gaps, we designed a systematic experiment to collect omics data and evaluate 17 traits in two oat breeding populations planted in single and multiple environments. In the single-environment trial, transcriptomic BLUP (T), metabolomic BLUP (M), G + T, G + M, and G + T + M models showed greater prediction accuracy than GBLUP for 5, 10, 11, 17, and 17 traits, respectively, and metabolites generally performed better than transcripts when combined with SNPs. In the multi-environment trial, multi-trait models with omics data outperformed both counterpart multi-trait GBLUP models and single-environment omics models, and the highest prediction accuracy was achieved when modeling genetic covariance as an unstructured covariance model. We also demonstrated that omics data can be used to prioritize loci from one population with omics data to improve genomic prediction in a distantly related population using a two-kernel linear model that accommodated both likely casual loci with large-effect and loci that explain little or no phenotypic variance. We propose that the two-kernel linear model is superior to most genomic prediction models that assume each variant is equally likely to affect the trait and can be used to improve prediction accuracy for any trait with prior knowledge of genetic architecture.

Multi-omics prediction of oat agronomic and seed nutritional traits across environments and in distantly related populations

Integration of multi-omics data improved prediction accuracies of oat agronomic and seed nutritional traits in multi-environment trials and distantly related populations in addition to the single-environment prediction. Multi-omics prediction has been shown to be superior to genomic prediction with genome-wide DNA-based genetic markers (G) for predicting phenotypes. However, most of the existing studies were based on historical datasets from one environment; therefore, they were unable to evaluate the efficiency of multi-omics prediction in multi-environment trials and distantly related populations. To fill those gaps, we designed a systematic experiment to collect omics data and evaluate 17 traits in two oat breeding populations planted in single and multiple environments. In the single-environment trial, transcriptomic BLUP (T), metabolomic BLUP (M), G + T, G + M, and G + T + M models showed greater prediction accuracy than GBLUP for 5, 10, 11, 17, and 17 traits, respectively, and metabolites generally performed better than transcripts when combined with SNPs. In the multi-environment trial, multi-trait models with omics data outperformed both counterpart multi-trait GBLUP models and single-environment omics models, and the highest prediction accuracy was achieved when modeling genetic covariance as an unstructured covariance model. We also demonstrated that omics data can be used to prioritize loci from one population with omics data to improve genomic prediction in a distantly related population using a two-kernel linear model that accommodated both likely casual loci with large-effect and loci that explain little or no phenotypic variance. We propose that the two-kernel linear model is superior to most genomic prediction models that assume each variant is equally likely to affect the trait and can be used to improve prediction accuracy for any trait with prior knowledge of genetic architecture.

Characterization of heterosis and genomic prediction-based establishment of heterotic patterns for developing better hybrids in pigeonpea

Whole-genome resequencing (WGRS) of 396 lines, consisting of 104 hybrid parental lines and 292 germplasm lines, were used to study the molecular basis of mid-parent heterosis (MPH) and to identify complementary heterotic patterns in pigeonpea [Cajanus cajan (L.) Millsp.] hybrids. The lines and hybrids were assessed for yield and yield-related traits in multiple environments. Our analysis showed positive MPH values in 78.6% of hybrids, confirming the potential of hybrid breeding in pigeonpea. By using genome-wide prediction and association mapping approaches, we identified 129 single nucleotide polymorphisms and 52 copy number variations with significant heterotic effects and also established a high-yielding heterotic pattern in pigeonpea. In summary, our study highlights the role of WGRS data in the study and use of heterosis in crops where hybrid breeding is expected to boost selection gain in order to ensure global food security.

Optimizing the setup of multienvironmental hybrid wheat yield trials for boosting the selection capability

The accuracy of genomic prediction increases with increasing heritability, and thus the challenge of optimizing the design of multienvironment yield trials under a limited budget arises. With this in mind, we aimed to find the best of several options to sparsely distribute a fixed number of plots across different environments to increase the accuracy of hybrid performance prediction. We used a comprehensive published genomic and phenotypic data set of 1,604 winter wheat (Triticum aestivum L.) hybrids and compared several commonly used biometric models for phenotypic data analysis in a resampling study to identify the one that most accurately estimated the hybrid performance in different imbalanced trials. Our results showed that when using information about genotypic relationships, genotypic values were more strongly associated with the reference values than when this information was ignored. In addition, a balanced environmental sampling resulted in an adequate characterization of each environment and increased the accuracy for estimating the hybrid performance. One promising design involved dividing the genotypes into equally sized subgroups that were tested in a subset of environments, with the constraint that the subgroups overlapped with respect to the environments. This scenario appears to be particularly appropriate, as it provided both high accuracies in the estimates of genotypic values and had low variability resulting from the data sample used. Thus, we were able to clearly demonstrate the utility for optimizing the design of multienvironment hybrid wheat yield trials in times of genomic selection.

Overview of Metabolomic Analysis and the Integration with Multi-Omics for Economic Traits in Cattle

Metabolomics has been applied to measure the dynamic metabolic responses, to understand the systematic biological networks, to reveal the potential genetic architecture, etc., for human diseases and livestock traits. For example, the current published results include the detected relevant candidate metabolites, identified metabolic pathways, potential systematic networks, etc., for different cattle traits that can be applied for further metabolomic and integrated omics studies. Therefore, summarizing the applications of metabolomics for economic traits is required in cattle. We here provide a comprehensive review about metabolomic analysis and its integration with other omics in five aspects: (1) characterization of the metabolomic profile of cattle; (2) metabolomic applications in cattle; (3) integrated metabolomic analysis with other omics; (4) methods and tools in metabolomic analysis; and (5) further potentialities. The review aims to investigate the existing metabolomic studies by highlighting the results in cattle, integrated with other omics studies, to understand the metabolic mechanisms underlying the economic traits and to provide useful information for further research and practical breeding programs in cattle.

The metabolomic landscape of rice heterosis highlights pathway biomarkers for predicting complex phenotypes

Understanding the molecular mechanisms underlying complex phenotypes requires systematic analyses of complicated metabolic networks and contributes to improvements in the breeding efficiency of staple cereal crops and diagnostic accuracy for human diseases. Here, we selected rice (Oryza sativa) heterosis as a complex phenotype and investigated the mechanisms of both vegetative and reproductive traits using an untargeted metabolomics strategy. Heterosis-associated analytes were identified, and the overlapping analytes were shown to underlie the association patterns for six agronomic traits. The heterosis-associated analytes of four yield components and plant height collectively contributed to yield heterosis, and the degree of contribution differed among the five traits. We performed dysregulated network analyses of the high- and low-better parent heterosis hybrids and found multiple types of metabolic pathways involved in heterosis. The metabolite levels of the significantly enriched pathways (especially those from amino acid and carbohydrate metabolism) were predictive of yield heterosis (area under the curve = 0.907 with 10 features), and the predictability of these pathway biomarkers was validated with hybrids across environments and populations. Our findings elucidate the metabolomic landscape of rice heterosis and highlight the potential application of pathway biomarkers in achieving accurate predictions of complex phenotypes.

Reciprocal Recurrent Genomic Selection Is Impacted by Genotype-by-Environment Interactions

Reciprocal recurrent genomic selection is a breeding strategy aimed at improving the hybrid performance of two base populations. It promises to significantly advance hybrid breeding in wheat. Against this backdrop, the main objective of this study was to empirically investigate the potential and limitations of reciprocal recurrent genomic selection. Genome-wide predictive equations were developed using genomic and phenotypic data from a comprehensive population of 1,604 single crosses between 120 female and 15 male wheat lines. Twenty superior female lines were selected for initiation of the reciprocal recurrent genomic selection program. Focusing on the female pool, one cycle was performed with genomic selection steps at the F₂ (60 out of 629 plants) and the F₅ stage (49 out of 382 plants). Selection gain for grain yield was evaluated at six locations. Analyses of the phenotypic data showed pronounced genotype-by-environment interactions with two environments that formed an outgroup compared to the environments used for the genome-wide prediction equations. Removing these two environments for further analysis resulted in a selection gain of 1.0 dt ha^-1 compared to the hybrids of the original 20 parental lines. This underscores the potential of reciprocal recurrent genomic selection to promote hybrid wheat breeding, but also highlights the need to develop robust genome-wide predictive equations.

Dissecting Bread Wheat Heterosis through the Integration of Agronomic and Physiological Traits

Simple Summary

To meet the challenge of feeding almost 10 billion people by 2050, wheat yield has to double by 2050. However, over the past 20 years, yield increase has slowed down and even stagnated in the main producing countries. Similar to what has been observed in maize, hybrids have been suggested as a solution to overcome yield stagnation in wheat. However, wheat heterosis, i.e., the fact that a progeny surpasses the performances of its parents, is still limited and poorly understood. To better characterize this phenomenon, we developed and phenotyped for physiological and agronomic traits 91 hybrids and their nineteen female and sixteen male parents. We showed that hybrids had a longer grain filling phase that led to bigger grains and an increased thousand kernel weight. This resulted in a better yield for 86% of hybrids compared to the average yield of their parents. In addition, hybrids appeared to be less affected by the negative correlation between protein content and yield compared to pure lines. These results shed light on the physiological bases underlying yield heterosis in wheat, paving new ways to breed for better wheat hybrids that can help to meet agriculture’s challenges.

Abstract

To meet the challenge of feeding almost 10 billion people by 2050, wheat yield has to double by 2050. However, over the past 20 years, yield increase has slowed down and even stagnated in the main producing countries. Following the example of maize, hybrids have been suggested as a solution to overcome yield stagnation in wheat. However, wheat heterosis is still limited and poorly understood. Gaining a better understanding of hybrid vigor holds the key to breed for better varieties. To this aim, we have developed and phenotyped for physiological and agronomic traits an incomplete factorial design consisting of 91 hybrids and their nineteen female and sixteen male parents. Monitoring the plant development with normalized difference vegetation index revealed that 89% of the hybrids including the five higher yielding hybrids had a longer grain filling phase with a delayed senescence that results in larger grain size. This average increase of 7.7% in thousand kernel weight translated to a positive mid-parent heterosis for grain yield for 86% of hybrids. In addition, hybrids displayed a positive grain protein deviation leading to a +4.7% heterosis in protein yield. These results shed light on the physiological bases underlying yield heterosis in wheat, paving new ways to breed for better wheat hybrids.

Exploring the genic resources underlying metabolites through mGWAS and mQTL in wheat: From large-scale gene identification and pathway elucidation to crop improvement

Common wheat (Triticum aestivum L.) is a leading cereal crop, but has lagged behind with respect to the interpretation of the molecular mechanisms of phenotypes compared with other major cereal crops such as rice and maize. The recently available genome sequence of wheat affords the pre-requisite information for efficiently exploiting the potential molecular resources for decoding the genetic architecture of complex traits and identifying valuable breeding targets. Meanwhile, the successful application of metabolomics as an emergent large-scale profiling methodology in several species has demonstrated this approach to be accessible for reaching the above goals. One such productive avenue is combining metabolomics approaches with genetic designs. However, this trial is not as widespread as that for sequencing technologies, especially when the acquisition, understanding, and application of metabolic approaches in wheat populations remain more difficult and even arguably underutilized. In this review, we briefly introduce the techniques used in the acquisition of metabolomics data and their utility in large-scale identification of functional candidate genes. Considerable progress has been made in delivering improved varieties, suggesting that the inclusion of information concerning these metabolites and genes and metabolic pathways enables a more explicit understanding of phenotypic traits and, as such, this procedure could serve as an -omics-informed roadmap for executing similar improvement strategies in wheat and other species.

Designing Future Crops: Genomics-Assisted Breeding Comes of Age

Over the past decade, genomics-assisted breeding (GAB) has been instrumental in harnessing the potential of modern genome resources and characterizing and exploiting allelic variation for germplasm enhancement and cultivar development. Sustaining GAB in the future (GAB 2.0) will rely upon a suite of new approaches that fast-track targeted manipulation of allelic variation for creating novel diversity and facilitate their rapid and efficient incorporation in crop improvement programs. Genomic breeding strategies that optimize crop genomes with accumulation of beneficial alleles and purging of deleterious alleles will be indispensable for designing future crops. In coming decades, GAB 2.0 is expected to play a crucial role in breeding more climate-smart crop cultivars with higher nutritional value in a cost-effective and timely manner.

Unlocking big data doubled the accuracy in predicting the grain yield in hybrid wheat

The potential of big data to support businesses has been demonstrated in financial services, manufacturing, and telecommunications. Here, we report on efforts to enter a new data era in plant breeding by collecting genomic and phenotypic information from 12,858 wheat genotypes representing 6575 single-cross hybrids and 6283 inbred lines that were evaluated in six experimental series for yield in field trials encompassing ~125,000 plots. Integrating data resulted in twofold higher prediction ability compared with cases in which hybrid performance was predicted across individual experimental series. Our results suggest that combining data across breeding programs is a particularly appropriate strategy to exploit the potential of big data for predictive plant breeding. This paradigm shift can contribute to increasing yield and resilience, which is needed to feed the growing world population.

Genomic prediction of hybrid crops allows disentangling dominance and epistasis

We revisited, in a genomic context, the theory of hybrid genetic evaluation models of hybrid crosses of pure lines, as the current practice is largely based on infinitesimal model assumptions. Expressions for covariances between hybrids due to additive substitution effects and dominance and epistatic deviations were analytically derived. Using dense markers in a GBLUP analysis, it is possible to split specific combining ability into dominance and across-groups epistatic deviations, and to split general combining ability (GCA) into within-line additive effects and within-line additive by additive (and higher order) epistatic deviations. We analyzed a publicly available maize data set of Dent × Flint hybrids using our new model (called GCA-model) up to additive by additive epistasis. To model higher order interactions within GCAs, we also fitted "residual genetic" line effects. Our new GCA-model was compared with another genomic model which assumes a uniquely defined effect of genes across origins. Most variation in hybrids is accounted by GCA. Variances due to dominance and epistasis have similar magnitudes. Models based on defining effects either differently or identically across heterotic groups resulted in similar predictive abilities for hybrids. The currently used model inflates the estimated additive genetic variance. This is not important for hybrid predictions but has consequences for the breeding scheme-e.g. overestimation of the genetic gain within heterotic group. Therefore, we recommend using GCA-model, which is appropriate for genomic prediction and variance component estimation in hybrid crops using genomic data, and whose results can be practically interpreted and used for breeding purposes.

On the use of GBLUP and its extension for GWAS with additive and epistatic effects

Genomic best linear unbiased prediction (GBLUP) is the most widely used model for genome-wide predictions. Interestingly, it is also possible to perform genome-wide association studies (GWAS) based on GBLUP. Although the estimated marker effects in GBLUP are shrunken and the conventional test based on such effects has low power, it was observed that a modified test statistic can be produced and the result of test was identical to a standard GWAS model. Later, a mathematical proof was given for the special case that there is no fixed covariate in GBLUP. Since then, the new approach has been called “GWAS by GBLUP”. Nevertheless, covariates such as environmental and subpopulation effects are very common in GBLUP. Thus, it is necessary to confirm the equivalence in the general case. Recently, the concept was generalized to GWAS for epistatic effects and the new approach was termed rapid epistatic mixed-model association analysis (REMMA) because it greatly improved the computational efficiency. However, the relationship between REMMA and the standard GWAS model has not been investigated. In this study, we first provided a general mathematical proof of the equivalence between “GWAS by GBLUP” and the standard GWAS model for additive effects. Then, we compared REMMA with the standard GWAS model for epistatic effects by a theoretical investigation and by empirical data analyses. We hypothesized that the similarity of the two models is influenced by the relative contribution of additive and epistatic effects to the phenotypic variance, which was verified by empirical and simulation studies.

Multi-omics-based prediction of hybrid performance in canola

Complementing or replacing genetic markers with transcriptomic data and use of reproducing kernel Hilbert space regression based on Gaussian kernels increases hybrid prediction accuracies for complex agronomic traits in canola. In plant breeding, hybrids gained particular importance due to heterosis, the superior performance of offspring compared to their inbred parents. Since the development of new top performing hybrids requires labour-intensive and costly breeding programmes, including testing of large numbers of experimental hybrids, the prediction of hybrid performance is of utmost interest to plant breeders. In this study, we tested the effectiveness of hybrid prediction models in spring-type oilseed rape (Brassica napus L./canola) employing different omics profiles, individually and in combination. To this end, a population of 950 F₁ hybrids was evaluated for seed yield and six other agronomically relevant traits in commercial field trials at several locations throughout Europe. A subset of these hybrids was also evaluated in a climatized glasshouse regarding early biomass production. For each of the 477 parental rapeseed lines, 13,201 single nucleotide polymorphisms (SNPs), 154 primary metabolites, and 19,479 transcripts were determined and used as predictive variables. Both, SNP markers and transcripts, effectively predict hybrid performance using (genomic) best linear unbiased prediction models (gBLUP). Compared to models using pure genetic markers, models incorporating transcriptome data resulted in significantly higher prediction accuracies for five out of seven agronomic traits, indicating that transcripts carry important information beyond genomic data. Notably, reproducing kernel Hilbert space regression based on Gaussian kernels significantly exceeded the predictive abilities of gBLUP models for six of the seven agronomic traits, demonstrating its potential for implementation in future canola breeding programmes.

Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge

We propose a novel Bayesian approach that robustifies genomic modeling by leveraging expert knowledge (EK) through prior distributions. The central component is the hierarchical decomposition of phenotypic variation into additive and nonadditive genetic variation, which leads to an intuitive model parameterization that can be visualized as a tree. The edges of the tree represent ratios of variances, for example broad-sense heritability, which are quantities for which EK is natural to exist. Penalized complexity priors are defined for all edges of the tree in a bottom-up procedure that respects the model structure and incorporates EK through all levels. We investigate models with different sources of variation and compare the performance of different priors implementing varying amounts of EK in the context of plant breeding. A simulation study shows that the proposed priors implementing EK improve the robustness of genomic modeling and the selection of the genetically best individuals in a breeding program. We observe this improvement in both variety selection on genetic values and parent selection on additive values; the variety selection benefited the most. In a real case study, EK increases phenotype prediction accuracy for cases in which the standard maximum likelihood approach did not find optimal estimates for the variance components. Finally, we discuss the importance of EK priors for genomic modeling and breeding, and point to future research areas of easy-to-use and parsimonious priors in genomic modeling.

Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge

We propose a novel Bayesian approach that robustifies genomic modeling by leveraging expert knowledge (EK) through prior distributions. The central component is the hierarchical decomposition of phenotypic variation into additive and nonadditive genetic variation, which leads to an intuitive model parameterization that can be visualized as a tree. The edges of the tree represent ratios of variances, for example broad-sense heritability, which are quantities for which EK is natural to exist. Penalized complexity priors are defined for all edges of the tree in a bottom-up procedure that respects the model structure and incorporates EK through all levels. We investigate models with different sources of variation and compare the performance of different priors implementing varying amounts of EK in the context of plant breeding. A simulation study shows that the proposed priors implementing EK improve the robustness of genomic modeling and the selection of the genetically best individuals in a breeding program. We observe this improvement in both variety selection on genetic values and parent selection on additive values; the variety selection benefited the most. In a real case study, EK increases phenotype prediction accuracy for cases in which the standard maximum likelihood approach did not find optimal estimates for the variance components. Finally, we discuss the importance of EK priors for genomic modeling and breeding, and point to future research areas of easy-to-use and parsimonious priors in genomic modeling.

Heterosis and Hybrid Crop Breeding: A Multidisciplinary Review

Although hybrid crop varieties are among the most popular agricultural innovations, the rationale for hybrid crop breeding is sometimes misunderstood. Hybrid breeding is slower and more resource-intensive than inbred breeding, but it allows systematic improvement of a population by recurrent selection and exploitation of heterosis simultaneously. Inbred parental lines can identically reproduce both themselves and their F1 progeny indefinitely, whereas outbred lines cannot, so uniform outbred lines must be bred indirectly through their inbred parents to harness heterosis. Heterosis is an expected consequence of whole-genome non-additive effects at the population level over evolutionary time. Understanding heterosis from the perspective of molecular genetic mechanisms alone may be elusive, because heterosis is likely an emergent property of populations. Hybrid breeding is a process of recurrent population improvement to maximize hybrid performance. Hybrid breeding is not maximization of heterosis per se, nor testing random combinations of individuals to find an exceptional hybrid, nor using heterosis in place of population improvement. Though there are methods to harness heterosis other than hybrid breeding, such as use of open-pollinated varieties or clonal propagation, they are not currently suitable for all crops or production environments. The use of genomic selection can decrease cycle time and costs in hybrid breeding, particularly by rapidly establishing heterotic pools, reducing testcrossing, and limiting the loss of genetic variance. Open questions in optimal use of genomic selection in hybrid crop breeding programs remain, such as how to choose founders of heterotic pools, the importance of dominance effects in genomic prediction, the necessary frequency of updating the training set with phenotypic information, and how to maintain genetic variance and prevent fixation of deleterious alleles.

Optimization of training sets for genomic prediction of early-stage single crosses in maize

Training population optimization algorithms are useful for efficiently training genomic prediction models for single-cross performance, especially if the population is extended beyond only realized crosses to all possible single crosses. Genomic prediction of single-cross performance could allow effective evaluation of all possible single crosses between all inbreds developed in a hybrid breeding program. The objectives of the present study were to investigate the effect of different levels of relatedness on genomic predictive ability of single crosses, evaluate the usefulness of deterministic formula to forecast prediction accuracy in advance, and determine the potential for TRS optimization based on prediction error variance (PEVmean) and coefficient of determination (CDmean) criteria. We used 481 single crosses made by crossing 89 random recombinant inbred lines (RILs) belonging to the Iowa stiff stalk synthetic group with 103 random RILs belonging to the non-stiff stalk synthetic heterotic group. As expected, predictive ability was enhanced by ensuring close relationships between TRSs and target sets, even when TRS sizes were smaller. We found that designing a TRS based on PEVmean or CDmean criteria is useful for increasing the efficiency of genomic prediction of maize single crosses. We went further and extended the sampling space from that of all observed single crosses to all possible single crosses, providing a much larger genetic space within which to design a training population. Using all possible single crosses increased the advantage of the PEVmean and CDmean methods based on expected prediction accuracy. This finding suggests that it may be worthwhile using an optimization algorithm to select a training population from all possible single crosses to maximize efficiency in training accurate models for hybrid genomic prediction.

Genetic Dissection of Hybrid Performance and Heterosis for Yield-Related Traits in Maize

Heterosis contributes a big proportion to hybrid performance in maize, especially for grain yield. It is attractive to explore the underlying genetic architecture of hybrid performance and heterosis. Considering its complexity, different from former mapping method, we developed a series of linear mixed models incorporating multiple polygenic covariance structures to quantify the contribution of each genetic component (additive, dominance, additive-by-additive, additive-by-dominance, and dominance-by-dominance) to hybrid performance and midparent heterosis variation and to identify significant additive and non-additive (dominance and epistatic) quantitative trait loci (QTL). Here, we developed a North Carolina II population by crossing 339 recombinant inbred lines with two elite lines (Chang7-2 and Mo17), resulting in two populations of hybrids signed as Chang7-2 × recombinant inbred lines and Mo17 × recombinant inbred lines, respectively. The results of a path analysis showed that kernel number per row and hundred grain weight contributed the most to the variation of grain yield. The heritability of midparent heterosis for 10 investigated traits ranged from 0.27 to 0.81. For the 10 traits, 21 main (additive and dominance) QTL for hybrid performance and 17 dominance QTL for midparent heterosis were identified in the pooled hybrid populations with two overlapping QTL. Several of the identified QTL showed pleiotropic effects. Significant epistatic QTL were also identified and were shown to play an important role in ear height variation. Genomic selection was used to assess the influence of QTL on prediction accuracy and to explore the strategy of heterosis utilization in maize breeding. Results showed that treating significant single nucleotide polymorphisms as fixed effects in the linear mixed model could improve the prediction accuracy under prediction schemes 2 and 3. In conclusion, the different analyses all substantiated the different genetic architecture of hybrid performance and midparent heterosis in maize. Dominance contributes the highest proportion to heterosis, especially for grain yield, however, epistasis contributes the highest proportion to hybrid performance of grain yield.