Can metabolic prediction be an alternative to genomic prediction in barley?

Multi-omics assists genomic prediction of maize yield with machine learning approaches

With the improvement of high-throughput technologies in recent years, large multi-dimensional plant omics data have been produced, and big-data-driven yield prediction research has received increasing attention. Machine learning offers promising computational and analytical solutions to interpret the biological meaning of large amounts of data in crops. In this study, we utilized multi-omics datasets from 156 maize recombinant inbred lines, containing 2496 single nucleotide polymorphisms (SNPs), 46 image traits (i-traits) from 16 developmental stages obtained through an automatic phenotyping platform, and 133 primary metabolites. Based on benchmark tests with different types of prediction models, some machine learning methods, such as Partial Least Squares (PLS), Random Forest (RF), and Gaussian process with Radial basis function kernel (GaussprRadial), achieved better prediction for maize yield, albeit slight difference for method preferences among i-traits, genomic, and metabolic data. We found that better yield prediction may be caused by various capabilities in ranking and filtering data features, which is found to be linked with biological meaning such as photosynthesis-related or kernel development-related regulations. Finally, by integrating multiple omics data with the RF machine learning approach, we can further improve the prediction accuracy of grain yield from 0.32 to 0.43. Our research provides new ideas for the application of plant omics data and artificial intelligence approaches to facilitate crop genetic improvements.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-024-01454-z.

Underutilized legumes: nutrient status and advanced breeding approaches for qualitative and quantitative enhancement

Underutilized/orphan legumes provide food and nutritional security to resource-poor rural populations during periods of drought and extreme hunger, thus, saving millions of lives. The Leguminaceae, which is the third largest flowering plant family, has approximately 650 genera and 20,000 species and are distributed globally. There are various protein-rich accessible and edible legumes, such as soybean, cowpea, and others; nevertheless, their consumption rate is far higher than production, owing to ever-increasing demand. The growing global urge to switch from an animal-based protein diet to a vegetarian-based protein diet has also accelerated their demand. In this context, underutilized legumes offer significant potential for food security, nutritional requirements, and agricultural development. Many of the known legumes like Mucuna spp., Canavalia spp., Sesbania spp., Phaseolus spp., and others are reported to contain comparable amounts of protein, essential amino acids, polyunsaturated fatty acids (PUFAs), dietary fiber, essential minerals and vitamins along with other bioactive compounds. Keeping this in mind, the current review focuses on the potential of discovering underutilized legumes as a source of food, feed and pharmaceutically valuable chemicals, in order to provide baseline data for addressing malnutrition-related problems and sustaining pulse needs across the globe. There is a scarcity of information about underutilized legumes and is restricted to specific geographical zones with local or traditional significance. Around 700 genera and 20,000 species remain for domestication, improvement, and mainstreaming. Significant efforts in research, breeding, and development are required to transform existing local landraces of carefully selected, promising crops into types with broad adaptability and economic viability. Different breeding efforts and the use of biotechnological methods such as micro-propagation, molecular markers research and genetic transformation for the development of underutilized crops are offered to popularize lesser-known legume crops and help farmers diversify their agricultural systems and boost their profitability.

Phenomic selection in wheat breeding: prediction of the genotype-by-environment interaction in multi-environment breeding trials

Phenomic prediction of wheat grain yield and heading date in different multi-environmental trial scenarios is accurate. Modelling the genotype-by-environment interaction effect using phenomic data is a potentially low-cost complement to genomic prediction. The performance of wheat cultivars in multi-environmental trials (MET) is difficult to predict because of the genotype-by-environment interactions (G × E). Phenomic selection is supposed to be efficient for modelling the G × E effect because it accounts for non-additive effects. Here, phenomic data are near-infrared (NIR) spectra obtained from plant material. While phenomic selection has recently been shown to accurately predict wheat grain yield in single environments, its accuracy needs to be investigated for MET. We used four datasets from two winter wheat breeding programs to test and compare the predictive abilities of phenomic and genomic models for grain yield and heading date in different MET scenarios. We also compared different methods to model the G × E using different covariance matrices based on spectra. On average, phenomic and genomic prediction abilities are similar in all different MET scenarios. Better predictive abilities were obtained when G × E effects were modelled with NIR spectra than without them, and it was better to use all the spectra of all genotypes in all environments for modelling the G × E. To facilitate the implementation of phenomic prediction, we tested MET designs where the NIR spectra were measured only on the genotype-environment combinations phenotyped for the target trait. Missing spectra were predicted with a weighted multivariate ridge regression. Intermediate predictive abilities for grain yield were obtained in a sparse testing scenario and for new genotypes, which shows that phenomic selection is an efficient and practicable prediction method for dealing with G × E.

Integration of DNA Methylation and Transcriptome Data Improves Complex Trait Prediction in Hordeum vulgare

Whole-genome multi-omics profiles contain valuable information for the characterization and prediction of complex traits in plants. In this study, we evaluate multi-omics models to predict four complex traits in barley (Hordeum vulgare); grain yield, thousand kernel weight, protein content, and nitrogen uptake. Genomic, transcriptomic, and DNA methylation data were obtained from 75 spring barley lines tested in the RadiMax semi-field phenomics facility under control and water-scarce treatment. By integrating multi-omics data at genomic, transcriptomic, and DNA methylation regulatory levels, a higher proportion of phenotypic variance was explained (0.72–0.91) than with genomic models alone (0.55–0.86). The correlation between predictions and phenotypes varied from 0.17–0.28 for control plants and 0.23–0.37 for water-scarce plants, and the increase in accuracy was significant for nitrogen uptake and protein content compared to models using genomic information alone. Adding transcriptomic and DNA methylation information to the prediction models explained more of the phenotypic variance attributed to the environment in grain yield and nitrogen uptake. It furthermore explained more of the non-additive genetic effects for thousand kernel weight and protein content. Our results show the feasibility of multi-omics prediction for complex traits in barley.

Leveraging a graft collection to develop metabolome-based trait prediction for the selection of tomato rootstocks with enhanced salt tolerance

Grafting has been demonstrated to significantly enhance the salt tolerance of crops. However, breeding efforts to develop enhanced graft combinations are hindered by knowledge-gaps as to how rootstocks mediate scion-response to salt stress. We grafted the scion of cultivated M82 onto rootstocks of 254 tomato accessions and explored the morphological and metabolic responses of grafts under saline conditions (EC = 20 dS m^-1) as compared to self-grafted M82 (SG-M82). Correlation analysis and Least Absolute Shrinkage and Selection Operator were performed to address the association between morphological diversification and metabolic perturbation. We demonstrate that grafting the same variety onto different rootstocks resulted in scion phenotypic heterogeneity and emphasized the productivity efficiency of M82 irrespective of the rootstock. Spectrophotometric analysis to test lipid oxidation showed largest variability of malondialdehyde (MDA) equivalents across the population, while the least responsive trait was the ratio of fruit fresh weight to total fresh weight (FFW/TFW). Generally, grafts showed greater values for the traits measured than SG-M82, except for branch number and wild race-originated rootstocks; the latter were associated with smaller scion growth parameters. Highly responsive and correlated metabolites were identified across the graft collection including malate, citrate, and aspartate, and their variance was partly related to rootstock origin. A group of six metabolites that consistently characterized exceptional graft response was observed, consisting of sorbose, galactose, sucrose, fructose, myo-inositol, and proline. The correlation analysis and predictive modelling, integrating phenotype- and leaf metabolite data, suggest a potential predictive relation between a set of leaf metabolites and yield-related traits.

Improvement of prediction ability by integrating multi-omic datasets in barley

Background

Genomic prediction (GP) based on single nucleotide polymorphisms (SNP) has become a broadly used tool to increase the gain of selection in plant breeding. However, using predictors that are biologically closer to the phenotypes such as transcriptome and metabolome may increase the prediction ability in GP. The objectives of this study were to (i) assess the prediction ability for three yield-related phenotypic traits using different omic datasets as single predictors compared to a SNP array, where these omic datasets included different types of sequence variants (full-SV, deleterious-dSV, and tolerant-tSV), different types of transcriptome (expression presence/absence variation-ePAV, gene expression-GE, and transcript expression-TE) sampled from two tissues, leaf and seedling, and metabolites (M); (ii) investigate the improvement in prediction ability when combining multiple omic datasets information to predict phenotypic variation in barley breeding programs; (iii) explore the predictive performance when using SV, GE, and ePAV from simulated 3’end mRNA sequencing of different lengths as predictors.

Results

The prediction ability from genomic best linear unbiased prediction (GBLUP) for the three traits using dSV information was higher than when using tSV, all SV information, or the SNP array. Any predictors from the transcriptome (GE, TE, as well as ePAV) and metabolome provided higher prediction abilities compared to the SNP array and SV on average across the three traits. In addition, some (di)-similarity existed between different omic datasets, and therefore provided complementary biological perspectives to phenotypic variation. Optimal combining the information of dSV, TE, ePAV, as well as metabolites into GP models could improve the prediction ability over that of the single predictors alone.

Conclusions

The use of integrated omic datasets in GP model is highly recommended. Furthermore, we evaluated a cost-effective approach generating 3’end mRNA sequencing with transcriptome data extracted from seedling without losing prediction ability in comparison to the full-length mRNA sequencing, paving the path for the use of such prediction methods in commercial breeding programs.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08337-7).

Recent applications of metabolomics in plant breeding

Metabolites play a central role in maintaining organismal life and in defining crop phenotypes, such as nutritional value, fragrance, color, and stress resistance. Among the 'omes' in biology, the metabolome is the closest to the phenotype. Consequently, metabolomics has been applied to crop improvement as method for monitoring changes in chemical compositions, clarifying the mechanisms underlying cellular functions, discovering markers and diagnostics, and phenotyping for mQTL, mGWAS, and metabolite-genome predictions. In this review, 359 reports of the most recent applications of metabolomics to plant breeding-related studies were examined. In addition to the major crops, more than 160 other crops including rare medicinal plants were considered. One bottleneck associated with using metabolomics is the wide array of instruments that are used to obtain data and the ambiguity associated with metabolite identification and quantification. To further the application of metabolomics to plant breeding, the features and perspectives of the technology are discussed.

Metabolomic selection for enhanced fruit flavor

Consumers often regard heirloom fruit varieties grown in the garden as more flavorful than commercial varieties purchased at the grocery store. While plant breeders have historically focused on improving producer-orientated traits such as yield, consumer-oriented traits such as flavor have regularly been neglected. This is, in part, due to the difficulty associated with measuring the sensory perceptions of flavor. Here, we combine fruit chemical and consumer sensory panel information to train machine learning models that can predict how flavorful a fruit will be from its chemistry. By increasing the throughput of flavor evaluations, these models will help plant breeders to integrate flavor earlier in the breeding pipeline and aid in the design of varieties with exceptional flavor profiles.

Although they are staple foods in cuisines globally, many commercial fruit varieties have become progressively less flavorful over time. Due to the cost and difficulty associated with flavor phenotyping, breeding programs have long been challenged in selecting for this complex trait. To address this issue, we leveraged targeted metabolomics of diverse tomato and blueberry accessions and their corresponding consumer panel ratings to create statistical and machine learning models that can predict sensory perceptions of fruit flavor. Using these models, a breeding program can assess flavor ratings for a large number of genotypes, previously limited by the low throughput of consumer sensory panels. The ability to predict consumer ratings of liking, sweet, sour, umami, and flavor intensity was evaluated by a 10-fold cross-validation, and the accuracies of 18 different models were assessed. The prediction accuracies were high for most attributes and ranged from 0.87 for sourness intensity in blueberry using XGBoost to 0.46 for overall liking in tomato using linear regression. Further, the best-performing models were used to infer the flavor compounds (sugars, acids, and volatiles) that contribute most to each flavor attribute. We found that the variance decomposition of overall liking score estimates that 42% and 56% of the variance was explained by volatile organic compounds in tomato and blueberry, respectively. We expect that these models will enable an earlier incorporation of flavor as breeding targets and encourage selection and release of more flavorful fruit varieties.

Gain time to adapt: How sorghum acquires tolerance to salinity

Salinity is a global environmental threat to agricultural production and food security around the world. To delineate salt-induced damage from adaption events we analysed a pair of sorghum genotypes which are contrasting in their response to salt stress with respect to physiological, cellular, metabolomic, and transcriptional responses. We find that the salt-tolerant genotype Della can delay the transfer of sodium from the root to the shoot, more swiftly deploy accumulation of proline and antioxidants in the leaves and transfer more sucrose to the root as compared to its susceptible counterpart Razinieh. Instead Razinieh shows metabolic indicators for a higher extent photorespiration under salt stress. Following sodium accumulation by a fluorescent dye in the different regions of the root, we find that Della can sequester sodium in the vacuoles of the distal elongation zone. The timing of the adaptive responses in Della leaves indicates a rapid systemic signal from the roots that is travelling faster than sodium itself. We arrive at a model where resistance and susceptibility are mainly a matter of temporal patterns in signalling.

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.

Prediction of complex phenotypes using the Drosophila melanogaster metabolome

Understanding the genotype-phenotype map and how variation at different levels of biological organization is associated are central topics in modern biology. Fast developments in sequencing technologies and other molecular omic tools enable researchers to obtain detailed information on variation at DNA level and on intermediate endophenotypes, such as RNA, proteins and metabolites. This can facilitate our understanding of the link between genotypes and molecular and functional organismal phenotypes. Here, we use the Drosophila melanogaster Genetic Reference Panel and nuclear magnetic resonance (NMR) metabolomics to investigate the ability of the metabolome to predict organismal phenotypes. We performed NMR metabolomics on four replicate pools of male flies from each of 170 different isogenic lines. Our results show that metabolite profiles are variable among the investigated lines and that this variation is highly heritable. Second, we identify genes associated with metabolome variation. Third, using the metabolome gave better prediction accuracies than genomic information for four of five quantitative traits analyzed. Our comprehensive characterization of population-scale diversity of metabolomes and its genetic basis illustrates that metabolites have large potential as predictors of organismal phenotypes. This finding is of great importance, e.g., in human medicine, evolutionary biology and animal and plant breeding.

Genome-wide association study on metabolite accumulation in a wild barley NAM population reveals natural variation in sugar metabolism

Metabolites play a key role in plants as they are routing plant developmental processes and are involved in biotic and abiotic stress responses. Their analysis can offer important information on the underlying processes. Regarding plant breeding, metabolite concentrations can be used as biomarkers instead of or in addition to genetic markers to predict important phenotypic traits (metabolic prediction). In this study, we applied a genome-wide association study (GWAS) in a wild barley nested association mapping (NAM) population to identify metabolic quantitative trait loci (mQTL). A set of approximately 130 metabolites, measured at early and late sampling dates, was analysed. For four metabolites from the early and six metabolites from the late sampling date significant mQTL (grouped as 19 mQTL for the early and 25 mQTL for the late sampling date) were found. Interestingly, all of those metabolites could be classified as sugars. Sugars are known to be involved in signalling, plant growth and plant development. Sugar-related genes, encoding mainly sugar transporters, have been identified as candidate genes for most of the mQTL. Moreover, several of them co-localized with known flowering time genes like Ppd-H1, HvELF3, Vrn-H1, Vrn-H2 and Vrn-H3, hinting on the known role of sugars in flowering. Furthermore, numerous disease resistance-related genes were detected, pointing to the signalling function of sugars in plant resistance. An mQTL on chromosome 1H in the region of 13 Mbp to 20 Mbp stood out, that alone explained up to 65% of the phenotypic variation of a single metabolite. Analysis of family-specific effects within the diverse NAM population showed the available natural genetic variation regarding sugar metabolites due to different wild alleles. The study represents a step towards a better understanding of the genetic components of metabolite accumulation, especially sugars, thereby linking them to biological functions in barley.

Fungal natural products galaxy: Biochemistry and molecular genetics toward blockbuster drugs discovery

Secondary metabolites synthesized by fungi have become a precious source of inspiration for the design of novel drugs. Indeed, fungi are prolific producers of fascinating, diverse, structurally complex, and low-molecular-mass natural products with high therapeutic leads, such as novel antimicrobial compounds, anticancer compounds, immunosuppressive agents, among others. Given that these microorganisms possess the extraordinary capacity to secrete diverse chemical scaffolds, they have been highly exploited by the giant pharma companies to generate small molecules. This has been made possible because the isolation of metabolites from fungal natural sources is feasible and surpasses the organic synthesis of compounds, which otherwise remains a significant bottleneck in the drug discovery process. Here in this comprehensive review, we have discussed recent studies on different fungi (pathogenic, non-pathogenic, commensal, and endophytic/symbiotic) from different habitats (terrestrial and marines), the specialized metabolites they biosynthesize, and the drugs derived from these specialized metabolites. Moreover, we have unveiled the logic behind the biosynthesis of vital chemical scaffolds, such as NRPS, PKS, PKS-NRPS hybrid, RiPPS, terpenoids, indole alkaloids, and their genetic mechanisms. Besides, we have provided a glimpse of the concept behind mycotoxins, virulence factor, and host immune response based on fungal infections.

Machine learning approaches for crop improvement: Leveraging phenotypic and genotypic big data

Highly efficient and accurate selection of elite genotypes can lead to dramatic shortening of the breeding cycle in major crops relevant for sustaining present demands for food, feed, and fuel. In contrast to classical approaches that emphasize the need for resource-intensive phenotyping at all stages of artificial selection, genomic selection dramatically reduces the need for phenotyping. Genomic selection relies on advances in machine learning and the availability of genotyping data to predict agronomically relevant phenotypic traits. Here we provide a systematic review of machine learning approaches applied for genomic selection of single and multiple traits in major crops in the past decade. We emphasize the need to gather data on intermediate phenotypes, e.g. metabolite, protein, and gene expression levels, along with developments of modeling techniques that can lead to further improvements of genomic selection. In addition, we provide a critical view of factors that affect genomic selection, with attention to transferability of models between different environments. Finally, we highlight the future aspects of integrating high-throughput molecular phenotypic data from omics technologies with biological networks for crop improvement.

Merging Genomics and Transcriptomics for Predicting Fusarium Head Blight Resistance in Wheat

Genomic selection with genome-wide distributed molecular markers has evolved into a well-implemented tool in many breeding programs during the last decade. The resistance against Fusarium head blight (FHB) in wheat is probably one of the most thoroughly studied systems within this framework. Aside from the genome, other biological strata like the transcriptome have likewise shown some potential in predictive breeding strategies but have not yet been investigated for the FHB-wheat pathosystem. The aims of this study were thus to compare the potential of genomic with transcriptomic prediction, and to assess the merit of blending incomplete transcriptomic with complete genomic data by the single-step method. A substantial advantage of gene expression data over molecular markers has been observed for the prediction of FHB resistance in the studied diversity panel of breeding lines and released cultivars. An increase in prediction ability was likewise found for the single-step predictions, although this can mostly be attributed to an increased accuracy among the RNA-sequenced genotypes. The usage of transcriptomics can thus be seen as a complement to already established predictive breeding pipelines with pedigree and genomic data, particularly when more cost-efficient multiplexing techniques for RNA-sequencing will become more accessible in the future.

Genomic Dissection of Peduncle Morphology in Barley through Nested Association Mapping

Straw biomass and stability are crucial for stable yields. Moreover, straw harbors the potential to serve as a valuable raw material for bio-economic processes. The peduncle is the top part of the last shoot internode and carries the spike. This study investigates the genetic control of barley peduncle morphology. Therefore, 1411 BC1S3 lines of the nested association mapping (NAM) population “Halle Exotic Barley 25” (HEB-25), generated by crossing the spring barley elite cultivar Barke with an assortment of 25 exotic barley accessions, were used. Applying 50k Illumina Infinium iSelect SNP genotyping yielded new insights and a better understanding of the quantitative trait loci (QTL) involved in controlling the peduncle diameter traits, we found the total thickness of peduncle tissues and the area of the peduncle cross-section. We identified three major QTL regions on chromosomes 2H and 3H mainly impacting the traits. Remarkably, the exotic allele at the QTL on chromosome 3H improved all three traits investigated in this work. Introgressing this QTL in elite cultivars might facilitate to adjust peduncle morphology for improved plant stability or enlarged straw biomass production independent of flowering time and without detrimental effects on grain yield.