Metabolomic analysis coupled with advanced genetic populations represents a powerful tool with which to investigate the plant metabolome. However, genetic analyses of the rice (Oryza sativa) metabolome have been conducted mainly using natural accessions or a single biparental population. Here, the flag leaves from three interconnected chromosome segment substitution line populations with a common recurrent genetic background were used to dissect rice metabolic diversity. We effectively used multiple interconnected biparental populations, constructed by introducing genomic segments into Zhenshan 97 from ACC10 (A/Z), Minghui 63 (M/Z), and Nipponbare (N/Z), to map metabolic quantitative trait loci (mQTL). A total of 1,587 mQTL were generated, of which 684, 479, and 722 were obtained from the A/Z, M/Z, and N/Z chromosome segment substitution line populations, respectively, and we designated 99 candidate genes for 367 mQTL. In addition, 1,001 mQTL were generated specifically from joint linkage analysis with 25 candidate genes assigned. Several of these candidates were validated, such as LOC_Os07g01020 for the in vivo content of pyridoxine and its derivative and LOC_Os04g25980 for cis-zeatin glucosyltransferase activity. We propose a novel biosynthetic pathway for O-methylapigenin C-pentoside and demonstrated that LOC_Os04g11970 encodes a component of this pathway through fine-mapping. We postulate that the methylated apigenin may confer plant disease resistance. This study demonstrates the power of using multiple interconnected populations to generate a large number of veritable mQTL. The combined results are discussed in the context of functional metabolomics and the possible features of assigned candidates underlying respective metabolites.

Rice (Oryza sativa) produces phenolamides and diterpenoids as major phytoalexins. Although the biosynthetic pathways of phenolamides and diterpenoids in plants have been revealed, knowledge of their accumulation regulatory mechanisms remains limited, and, in particular, no co-regulatory factor has been identified to date. Here, using a combined co-expression and evolutionary analysis, we identified the basic helix–loop–helix (bHLH) transcription factor OsbHLH5 as a positive bifunctional regulator of phenolamide and diterpenoid biosynthesis in rice. Metabolomic analysis revealed that OsbHLH5 significantly increased the content of phenolamides (such as feruloyl tryptamine (Fer-Trm) and p-coumaroyl tyramine (Cou-Tyr)) and diterpenoid phytoalexins (such as momilactones A, momilactones B) in the overexpression lines, while their content was reduced in the OsbHLH5 knockout lines. Gene expression and dual-luciferase assays revealed that OsbHLH5 activates phenolamide biosynthetic genes (including putrescine hydroxycinnamoyltransferase 3 (OsPHT3), tyramine hydroxycinnamoyltransferases 1/2 (OsTHT1/2), and tryptamine benzoyltransferase 2 (OsTBT2)) as well as diterpenoid biosynthetic genes (including copalyl diphosphate synthase 4 (OsCPS4) and kaurene synthase-like 4/7/10/11 (OsKSL4/7/10/11)). Furthermore, we have demonstrated that OsbHLH5 is induced by jasmonic acid (JA), while pathogen inoculation assays indicated that the overexpression of OsbHLH5 in transgenic rice plants leads to enhanced resistance to Xanthomonas oryzae pv. oryzae (Xoo). Overall, we have identified a positive regulator of phenolamide and diterpenoid biosynthesis and have demonstrated that biotic stress induces phytoalexin accumulation partly in an OsbHLH5-dependent manner, providing new insights into the metabolic interactions involved in pathogen response and offering valuable gene resources for the development, through genetic improvement, of new rice varieties that are resistant to diseases.