Metabolome Analysis of Multi-Connected Biparental Chromosome Segment Substitution Line Populations

Beyond pathways: Accelerated flavonoids candidate identification and novel exploration of enzymatic properties using combined mapping populations of wheat

Although forward-genetics-metabolomics methods such as mGWAS and mQTL have proven effective in providing myriad loci affecting metabolite contents, they are somehow constrained by their respective constitutional flaws such as the hidden population structure for GWAS and insufficient recombinant rate for QTL. Here, the combination of mGWAS and mQTL was performed, conveying an improved statistical power to investigate the flavonoid pathways in common wheat. A total of 941 and 289 loci were, respectively, generated from mGWAS and mQTL, within which 13 of them were co-mapped using both approaches. Subsequently, the mGWAS or mQTL outputs alone and their combination were, respectively, utilized to delineate the metabolic routes. Using this approach, we identified two MYB transcription factor encoding genes and five structural genes, and the flavonoid pathway in wheat was accordingly updated. Moreover, we have discovered some rare-activity-exhibiting flavonoid glycosyl- and methyl-transferases, which may possess unique biological significance, and harnessing these novel catalytic capabilities provides potentially new breeding directions. Collectively, we propose our survey illustrates that the forward-genetics-metabolomics approaches including multiple populations with high density markers could be more frequently applied for delineating metabolic pathways in common wheat, which will ultimately contribute to metabolomics-assisted wheat crop improvement.

Metabolite-based genome-wide association studies enable the dissection of the genetic bases of flavonoids, betaine and spermidine in wolfberry (Lycium)

Wolfberry is a plant with medicinal and food values. However, its bioactive ingredients and the corresponding genetic bases have not been determined. Here, we de novo generated a chromosome-level genome assembly for wolfberry, yielding a genome sequence of ~1.77 Gb with contig N50 of 50.55 Mb and 39 224 predicted gene models. A variation map, using 307 re-sequenced accessions, was called based on this genome assembly. Furthermore, the fruit metabolome of these accessions was profiled using 563 annotated metabolites, which separated Lycium barbarum L. and non-L. barbarum L. The flavonoids, coumarins, alkaloids and nicotinic acid contents were higher in the former than in the latter. A metabolite-based genome-wide association study mapped 156 164 significant single nucleotide polymorphisms corresponding to 340 metabolites. This included 19 219 unique lead single nucleotide polymorphisms in 1517 significant association loci, of which three metabolites, flavonoids, betaine and spermidine, were highlighted. Two candidate genes, LbUGT (evm.TU.chr07.2692) and LbCHS (evm.TU.chr07.2738), with non-synonymous mutations, were associated with the flavonoids content. LbCHS is a structural gene that interacts with a nearby MYB transcription factor (evm.TU.chr07.2726) both in L. barbarum and L. ruthenicum. Thus, these three genes might be involved in the biosynthesis/metabolism of flavonoids. LbSSADH (evm.TU.chr09.627) was identified as possibly participating in betaine biosynthesis/metabolism. Four lycibarbarspermidines (E-G and O) were identified, and only the lycibarbarspermidines O content was higher in L. barbarum varieties than in non-L. barbarum varieties. The evm.TU.chr07.2680 gene associated with lycibarbarspermidines O was annotated as an acetyl-CoA-benzylalcohol acetyltransferase, suggesting that it is a candidate gene for spermidine biosynthesis. These results provide novel insights into the specific metabolite profile of non-L. barbarum L. and the genetic bases of flavonoids, betaine and spermidine biosynthesis/metabolism.

Dynamic metabolite QTL analyses provide novel biochemical insights into kernel development and nutritional quality improvement in common wheat

Despite recent advances in crop metabolomics, the genetic control and molecular basis of the wheat kernel metabolome at different developmental stages remain largely unknown. Here, we performed widely targeted metabolite profiling of kernels from three developmental stages (grain-filling kernels [FKs], mature kernels [MKs], and germinating kernels [GKs]) using a population of 159 recombinant inbred lines. We detected 625 annotated metabolites and mapped 3173, 3143, and 2644 metabolite quantitative trait loci (mQTLs) in FKs, MKs, and GKs, respectively. Only 52 mQTLs were mapped at all three stages, indicating the high stage specificity of the wheat kernel metabolome. Four candidate genes were functionally validated by in vitro enzymatic reactions and/or transgenic approaches in wheat, three of which mediated the tricin metabolic pathway. Metabolite flux efficiencies within the tricin pathway were evaluated, and superior candidate haplotypes were identified, comprehensively delineating the tricin metabolism pathway in wheat. Finally, additional wheat metabolic pathways were re-constructed by updating them to incorporate the 177 candidate genes identified in this study. Our work provides new information on variations in the wheat kernel metabolome and important molecular resources for improvement of wheat nutritional quality.

Plant biochemical genetics in the multiomics era

Our understanding of plant biology has been revolutionized by modern genetics and biochemistry. However, biochemical genetics can be traced back to the foundation of Mendelian genetics; indeed, one of Mendel's milestone discoveries of seven characteristics of pea plants later came to be ascribed to a mutation in a starch branching enzyme. Here, we review both current and historical strategies for the elucidation of plant metabolic pathways and the genes that encode their component enzymes and regulators. We use this historical review to discuss a range of classical genetic phenomena including epistasis, canalization, and heterosis as viewed through the lens of contemporary high-throughput data obtained via the array of approaches currently adopted in multiomics studies.

Metabolic-GWAS provides insights into genetic architecture of seed metabolome in buckwheat

BACKGROUND

Buckwheat (Fagopyrum spp.), belonging to the Polygonaceae family, is an ancient pseudo-cereal with high nutritional and nutraceutical properties. Buckwheat proteins are gluten-free and show balanced amino acid and micronutrient profiles, with higher content of health-promoting bioactive flavonoids that make it a golden crop of the future. Plant metabolome is increasingly gaining importance as a crucial component to understand the connection between plant physiology and environment and as a potential link between the genome and phenome. However, the genetic architecture governing the metabolome and thus, the phenome is not well understood. Here, we aim to obtain a deeper insight into the genetic architecture of seed metabolome in buckwheat by integrating high throughput metabolomics and genotyping-by-sequencing applying an array of bioinformatics tools for data analysis.

RESULTS

High throughput metabolomic analysis identified 24 metabolites in seed endosperm of 130 diverse buckwheat genotypes. The genotyping-by-sequencing (GBS) of these genotypes revealed 3,728,028 SNPs. The Genome Association and Prediction Integrated Tool (GAPIT) assisted in the identification of 27 SNPs/QTLs linked to 18 metabolites. Candidate genes were identified near 100 Kb of QTLs, providing insights into several metabolic and biosynthetic pathways.

CONCLUSIONS

We established the metabolome inventory of 130 germplasm lines of buckwheat, identified QTLs through marker trait association and positions of potential candidate genes. This will pave the way for future dissection of complex economic traits in buckwheat.

Natural variation of respiration-related traits in plants

Plant respiration is one of the greatest global metabolic fluxes, but rates of respiration vary massively both within different cell types as well as between different individuals and different species. Whilst this is well known, few studies have detailed population-level variation of respiration until recently. The last 20 years have seen a renaissance in studies of natural variance. In this review, we describe how experimental breeding populations and collections of large populations of accessions can be used to determine the genetic architecture of plant traits. We further detail how these approaches have been used to study the rate of respiration per se as well as traits that are intimately associated with respiration. The review highlights specific breakthroughs in these areas but also concludes that the approach should be more widely adopted in the study of respiration per se as opposed to the more frequently studied respiration-related traits.

The pathway of melatonin biosynthesis in common wheat (Triticum aestivum)

Melatonin (Mel) is a multifunctional biomolecule found in both animals and plants. In plants, the biosynthesis of Mel from tryptophan (Trp) has been delineated to comprise of four consecutive reactions. However, while the genes encoding these enzymes in rice are well characterized no systematic evaluation of the overall pathway has, as yet, been published for wheat. In the current study, the relative contents of six Mel-pathway-intermediates including Trp, tryptamine (Trm), serotonin (Ser), 5-methoxy tryptamine (5M-Trm), N-acetyl serotonin (NAS) and Mel, were determined in 24 independent tissues spanning the lifetime of wheat. These studies indicated that Trp was the most abundant among the six metabolites, followed by Trm and Ser. Next, the candidate genes expressing key enzymes involved in the Mel pathway were explored by means of metabolite-based genome-wide association study (mGWAS), wherein two TDC genes, a T5H gene and one SNAT gene were identified as being important for the accumulation of Mel pathway metabolites. Moreover, a 463-bp insertion within the T5H gene was discovered that may be responsible for variation in Ser content. Finally, a ASMT gene was found via sequence alignment against its rice homolog. Validations of these candidate genes were performed by in vitro enzymatic reactions using proteins purified following recombinant expression in Escherichia coli, transient gene expression in tobacco, and transgenic approaches in wheat. Our results thus provide the first comprehensive investigation into the Mel pathway metabolites, and a swift candidate gene identification via forward-genetics strategies, in common wheat.

Detecting and pyramiding target QTL for plant- and grain-related traits via chromosomal segment substitution line of rice

INTRODUCTION

Plant height and grain length are important agronomic traits in rice, exhibiting a strong effect on plant architecture and grain quality of rice varieties.

METHODS

Methods: A novel rice chromosomal segment substitution line (CSSL), i.e., CSSL-Z1357, with significantly increased plant height (PH) and grain length (GL) was identified from CSSLs constructed by using Nipponbare as a receptor and a restorer line Xihui 18 as a donor. Seven agronomic traits of PH, PL, GL, GW, GPP, SPP, and TGW were phenotyped, and REML implemented in HPMIXED of SAS were used to detect the QTL for these traits. Secondary CSSLs were screened out via marker-assisted selection (MAS) to estimate the additive and epistatic effects of detected QTLs, evaluating the potential utilization of pyramiding the target QTLs for yield and quality improvement of rice varieties.

RESULTS AND DISCUSSION

Results and Discussion: CSSL-Z1357 carried nine segments from Xihui 18 with an average segment length of 4.13 Mb. The results show that the long grain of CSSL-Z1357 was caused by the increased number of surface cells and the length of the inner glume. Thirteen quantitative trait loci were identified via the F2 population of Nipponbare/CSSL-Z1357, including three each for GL (qGL-3, qGL-6, and qGL-7) and PH (qPH-1, qPH-7, and qPH-12I), among which qGL-3 increased GL by 0.23 mm with synergistic allele from CSSL-Z1357. Additionally, three single (S1 to S3), two double (D1, D2), and one triple segment (T1) substitution lines were developed in F3 via MAS. Results show that pyramiding the segments from Chr.3 (qGL-3 and qPH-3), Chr.6 (qGL-6 and qPH-6), and Chr.7 (Null and qPH-7) tended to result in better phenotype of increased GL and PH and decreased grain width, providing a potential basis for enhancing grain yield and quality in rice breeding.

Twenty years of rice genomics research: From sequencing and functional genomics to quantitative genomics

Since the completion of the rice genome sequencing project in 2005, we have entered the era of rice genomics, which is still in its ascendancy. Rice genomics studies can be classified into three stages: structural genomics, functional genomics, and quantitative genomics. Structural genomics refers primarily to genome sequencing for the construction of a complete map of rice genome sequence. This is fundamental for rice genetics and molecular biology research. Functional genomics aims to decode the functions of rice genes. Quantitative genomics is large-scale sequence- and statistics-based research to define the quantitative traits and genetic features of rice populations. Rice genomics has been a transformative influence on rice biological research and contributes significantly to rice breeding, making rice a good model plant for studying crop sciences.

Fine Mapping of Two Interacting Loci for Transmission Ratio Distortion in Rice (Oryza sativa L.)

Transmission ratio distortion (TRD) denotes the observed allelic or genotypic frequency deviation from the expected Mendelian segregation ratios in the offspring of a heterozygote. TRD can severely hamper gene flow between and within rice species. Here, we report the fine mapping and characterization of two loci (TRD4.1 and TRD4.2) for TRD using large F₂ segregating populations, which are derived from rice chromosome segment substitution lines, each containing a particular genomic segment introduced from the japonica cultivar Nipponbare (NIP) into the indica cultivar Zhenshan (ZS97). The two loci exhibited a preferential transmission of ZS97 alleles in the derived progeny. Reciprocal crossing experiments using near-isogenic lines harboring three different alleles at TRD4.1 suggest that the gene causes male gametic selection. Moreover, the transmission bias of TRD4.2 was diminished in heterozygotes when they carried homozygous TRD4.1 ^ZS97. This indicates an epistatic interaction between these two loci. TRD4.2 was mapped into a 35-kb region encompassing one candidate gene that is specifically expressed in the reproductive organs in rice. These findings broaden the understanding of the genetic mechanisms of TRD and offer an approach to overcome the barrier of gene flow between the subspecies in rice, thus facilitating rice improvement by introgression breeding.

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.

Integration of Dual Stress Transcriptomes and Major QTLs from a Pair of Genotypes Contrasting for Drought and Chronic Nitrogen Starvation Identifies Key Stress Responsive Genes in Rice

We report here the genome-wide changes resulting from low N (N-W+), low water (N+W-)) and dual stresses (N-W-) in root and shoot tissues of two rice genotypes, namely, IR 64 (IR64) and Nagina 22 (N22), and their association with the QTLs for nitrogen use efficiency. For all the root parameters, except for root length under N-W+, N22 performed better than IR64. Chlorophyll a, b and carotenoid content were higher in IR64 under N+W+ treatment and N-W+ and N+W- stresses; however, under dual stress, N22 had higher chlorophyll b content. While nitrite reductase, glutamate synthase (GS) and citrate synthase assays showed better specific activity in IR64, glutamate dehydrogenase showed better specific activity in N22 under dual stress (N-W-); the other N and C assimilating enzymes showed similar but low specific activities in both the genotypes. A total of 8926 differentially expressed genes (DEGs) were identified compared to optimal (N+W+) condition from across all treatments. While 1174, 698 and 903 DEGs in IR64 roots and 1197, 187 and 781 in N22 roots were identified, nearly double the number of DEGs were found in the shoot tissues; 3357, 1006 and 4005 in IR64 and 4004, 990 and 2143 in N22, under N-W+, N+W- and N-W- treatments, respectively. IR64 and N22 showed differential expression in 15 and 11 N-transporter genes respectively, under one or more stress treatments, out of which four showed differential expression also in N+W- condition. The negative regulators of N- stress, e.g., NIGT1, OsACTPK1 and OsBT were downregulated in IR64 while in N22, OsBT was not downregulated. Overall, N22 performed better under dual stress conditions owing to its better root architecture, chlorophyll and porphyrin synthesis and oxidative stress management. We identified 12 QTLs for seed and straw N content using 253 recombinant inbred lines derived from IR64 and N22 and a 5K SNP array. The QTL hotspot region on chromosome 6 comprised of 61 genes, of which, five were DEGs encoding for UDP-glucuronosyltransferase, serine threonine kinase, anthocyanidin 3-O-glucosyltransferase, and nitrate induced proteins. The DEGs, QTLs and candidate genes reported in this study can serve as a major resource for both rice improvement and functional biology.

Genomics and breeding innovations for enhancing genetic gain for climate resilience and nutrition traits

KEY MESSAGE

Integrating genomics technologies and breeding methods to tweak core parameters of the breeder's equation could accelerate delivery of climate-resilient and nutrient rich crops for future food security. Accelerating genetic gain in crop improvement programs with respect to climate resilience and nutrition traits, and the realization of the improved gain in farmers' fields require integration of several approaches. This article focuses on innovative approaches to address core components of the breeder's equation. A prerequisite to enhancing genetic variance (σ2g) is the identification or creation of favorable alleles/haplotypes and their deployment for improving key traits. Novel alleles for new and existing target traits need to be accessed and added to the breeding population while maintaining genetic diversity. Selection intensity (i) in the breeding program can be improved by testing a larger population size, enabled by the statistical designs with minimal replications and high-throughput phenotyping. Selection priorities and criteria to select appropriate portion of the population too assume an important role. The most important component of breeder's equation is heritability (h2). Heritability estimates depend on several factors including the size and the type of population and the statistical methods. The present article starts with a brief discussion on the potential ways to enhance σ2g in the population. We highlight statistical methods and experimental designs that could improve trait heritability estimation. We also offer a perspective on reducing the breeding cycle time (t), which could be achieved through the selection of appropriate parents, optimizing the breeding scheme, rapid fixation of target alleles, and combining speed breeding with breeding programs to optimize trials for release. Finally, we summarize knowledge from multiple disciplines for enhancing genetic gains for climate resilience and nutritional traits.

Understanding Omics Driven Plant Improvement and de novo Crop Domestication: Some Examples

In the current era, one of biggest challenges is to shorten the breeding cycle for rapid generation of a new crop variety having high yield capacity, disease resistance, high nutrient content, etc. Advances in the "-omics" technology have revolutionized the discovery of genes and bio-molecules with remarkable precision, resulting in significant development of plant-focused metabolic databases and resources. Metabolomics has been widely used in several model plants and crop species to examine metabolic drift and changes in metabolic composition during various developmental stages and in response to stimuli. Over the last few decades, these efforts have resulted in a significantly improved understanding of the metabolic pathways of plants through identification of several unknown intermediates. This has assisted in developing several new metabolically engineered important crops with desirable agronomic traits, and has facilitated the de novo domestication of new crops for sustainable agriculture and food security. In this review, we discuss how "omics" technologies, particularly metabolomics, has enhanced our understanding of important traits and allowed speedy domestication of novel crop plants.

Integrating multi-omics data for crop improvement

Our agricultural systems are now in urgent need to secure food for a growing world population. To meet this challenge, we need a better characterization of plant genetic and phenotypic diversity. The combination of genomics, transcriptomics and metabolomics enables a deeper understanding of the mechanisms underlying the complex architecture of many phenotypic traits of agricultural relevance. We review the recent advances in plant genomics to see how these can be integrated with broad molecular profiling approaches to improve our understanding of plant phenotypic variation and inform crop breeding strategies.

Promoting Human Nutrition and Health through Plant Metabolomics: Current Status and Challenges

Simple Summary

This review summarizes the status, applications, and challenges of plant metabolomics in the context of crop breeding, food quality and safety, and human nutrition and health. It also highlights the importance of plant metabolomics in elucidating biochemical and genetic bases of traits associated with nutritive and healthy beneficial foods and other plant products to secure food supply, to ensure food quality, to protect humans from malnutrition and other diseases. Meanwhile, this review calls for comprehensive collaborations to accelerate relevant researches and applications in the context of human nutrition and health.

Abstract

Plant metabolomics plays important roles in both basic and applied studies regarding all aspects of plant development and stress responses. With the improvement of living standards, people need high quality and safe food supplies. Thus, understanding the pathways involved in the biosynthesis of nutritionally and healthily associated metabolites in plants and the responses to plant-derived biohazards in humans is of equal importance to meet people’s needs. For each, metabolomics has a vital role to play, which is discussed in detail in this review. In addition, the core elements of plant metabolomics are highlighted, researches on metabolomics-based crop improvement for nutrition and safety are summarized, metabolomics studies on plant natural products including traditional Chinese medicine (TCM) for health promotion are briefly presented. Challenges are discussed and future perspectives of metabolomics as one of the most important tools to promote human nutrition and health are proposed.

Metabolite-based genome-wide association study enables dissection of the flavonoid decoration pathway of wheat kernels

The marriage of metabolomic approaches with genetic design has proven a powerful tool in dissecting diversity in the metabolome and has additionally enhanced our understanding of complex traits. That said, such studies have rarely been carried out in wheat. In this study, we detected 805 metabolites from wheat kernels and profiled their relative contents among 182 wheat accessions, conducting a metabolite-based genome-wide association study (mGWAS) utilizing 14 646 previously described polymorphic SNP markers. A total of 1098 mGWAS associations were detected with large effects, within which 26 candidate genes were tentatively designated for 42 loci. Enzymatic assay of two candidates indicated they could catalyse glucosylation and subsequent malonylation of various flavonoids and thereby the major flavonoid decoration pathway of wheat kernel was dissected. Moreover, numerous high-confidence genes associated with metabolite contents have been provided, as well as more subdivided metabolite networks which are yet to be explored within our data. These combined efforts presented the first step towards realizing metabolomics-associated breeding of wheat.

5Gs for crop genetic improvement

Here we propose a 5G breeding approach for bringing much-needed disruptive changes to crop improvement. These 5Gs are Genome assembly, Germplasm characterization, Gene function identification, Genomic breeding (GB), and Gene editing (GE). In our view, it is important to have genome assemblies available for each crop and a deep collection of germplasm characterized at sequencing and agronomic levels for identification of marker-trait associations and superior haplotypes. Systems biology and sequencing-based mapping approaches can be used to identify genes involved in pathways leading to the expression of a trait, thereby providing diagnostic markers for target traits. These genes, markers, haplotypes, and genome-wide sequencing data may be utilized in GB and GE methodologies in combination with a rapid cycle breeding strategy.

Decoding Plant-Environment Interactions That Influence Crop Agronomic Traits

To ensure food security in the face of increasing global demand due to population growth and progressive urbanization, it will be crucial to integrate emerging technologies in multiple disciplines to accelerate overall throughput of gene discovery and crop breeding. Plant agronomic traits often appear during the plants' later growth stages due to the cumulative effects of their lifetime interactions with the environment. Therefore, decoding plant-environment interactions by elucidating plants' temporal physiological responses to environmental changes throughout their lifespans will facilitate the identification of genetic and environmental factors, timing and pathways that influence complex end-point agronomic traits, such as yield. Here, we discuss the expected role of the life-course approach to monitoring plant and crop health status in improving crop productivity by enhancing the understanding of plant-environment interactions. We review recent advances in analytical technologies for monitoring health status in plants based on multi-omics analyses and strategies for integrating heterogeneous datasets from multiple omics areas to identify informative factors associated with traits of interest. In addition, we showcase emerging phenomics techniques that enable the noninvasive and continuous monitoring of plant growth by various means, including three-dimensional phenotyping, plant root phenotyping, implantable/injectable sensors and affordable phenotyping devices. Finally, we present an integrated review of analytical technologies and applications for monitoring plant growth, developed across disciplines, such as plant science, data science and sensors and Internet-of-things technologies, to improve plant productivity.

Computational approaches to unravel the pathways and evolution of specialized metabolism

Specialized metabolites serve as a chemical arsenal that protects plants from abiotic stress, pathogens, and herbivores, and they are an essential component of our nutrition and medicine. Despite their importance, we are still at the beginning of unravelling biosynthetic pathways that produce these compounds, which is needed to produce more resilient and nutritious crops, expand our inventory of useful biomolecules, and give valuable insights into plant evolution. This review focuses on the evolution of specialized metabolism in the plant kingdom and the state-of-the-art approaches used to identify the biosynthetic pathways of these useful compounds.

Comparative analysis of metabolome of rice seeds at three developmental stages using a recombinant inbred line population

Plants are considered an important food and nutrition source for humans. Despite advances in plant seed metabolomics, knowledge about the genetic and molecular bases of rice seed metabolomes at different developmental stages is still limited. Here, using Zhenshan 97 (ZS97) and Minghui 63 (MH63), we performed a widely targeted metabolic profiling in seeds during grain filling, mature seeds and germinating seeds. The diversity between MH63 and ZS97 was characterized in terms of the content of metabolites and the metabolic shifting across developmental stages. Taking advantage of the ultra-high-density genetic map of a population of 210 recombinant inbred lines (RILs) derived from a cross between ZS97 and MH63, we identified 4681 putative metabolic quantitative trait loci (mQTLs) in seeds across the three stages. Further analysis of the mQTLs for the codetected metabolites across the three stages revealed that the genetic regulation of metabolite accumulation was closely related to developmental stage. Using in silico analyses, we characterized 35 candidate genes responsible for 30 structurally identified or annotated compounds, among which LOC_Os07g04970 and LOC_Os06g03990 were identified to be responsible for feruloylserotonin and l-asparagine content variation across populations, respectively. Metabolite-agronomic trait association and colocation between mQTLs and phenotypic quantitative trait loci (pQTLs) revealed the complexity of the metabolite-agronomic trait relationship and the corresponding genetic basis.