GWA Mapping of Anthocyanin Accumulation Reveals Balancing Selection of MYB90 in Arabidopsis thaliana

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.

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Genome-wide analysis of R2R3-MYB transcription factors in Boehmeria nivea (L.) gaudich revealed potential cadmium tolerance and anthocyanin biosynthesis genes

Gene family, especially MYB as one of the largest transcription factor family in plants, the study of its subfunctional characteristics is a key step in the study of plant gene function. The sequencing of ramie genome provides a good opportunity to study the organization and evolutionary characters of the ramie MYB gene at the whole genome level. In this study, a total of 105 BnGR2R3-MYB genes were identified from ramie genome and subsequently grouped into 35 subfamilies according to phylogeny divergence and sequences similarity. Chromosomal localization, gene structure, synteny analysis, gene duplication, promoter analysis, molecular characteristics and subcellular localization were accomplished using several bioinformatics tools. Collinearity analysis showed that the segmental and tandem duplication events is the dominant form of the gene family expansion, and duplications prominent in distal telomeric regions. Highest syntenic relationship was obtained between BnGR2R3-MYB genes and that of Apocynum venetum (88). Furthermore, transcriptomic data and phylogenetic analysis revealed that BnGMYB60, BnGMYB79/80 and BnGMYB70 might inhibit the biosynthesis of anthocyanins, and UPLC-QTOF-MS data further supported the results. qPCR and phylogenetic analysis revealed that the six genes (BnGMYB9, BnGMYB10, BnGMYB12, BnGMYB28, BnGMYB41, and BnGMYB78) were cadmium stress responsive genes. Especially, the expression of BnGMYB10/12/41 in roots, stems and leaves all increased more than 10-fold after cadmium stress, and in addition they may interact with key genes regulating flavonoid biosynthesis. Thus, a potential link between cadmium stress response and flavonoid synthesis was identified through protein interaction network analysis. The study thus provided significant information into MYB regulatory genes in ramie and may serve as a foundation for genetic enhancement and increased productivity.

An R3-MYB repressor, BnCPC forms a feedback regulation with MBW complex to modulate anthocyanin biosynthesis in Brassica napus

BACKGROUND

Anthocyanins are metabolites of phenylpropanoid pathway, and involves in diverse processes of plant development and adaptation, which are regulated by the MYB-bHLH-WD40 (MBW) protein complexes. Many R2R3-MYB activators have been well characterized, but the MYB repressors in anthocyanin biosynthesis were recognized recently, which are also important in modulating phenylpropanoid metabolism in plants. The regulatory mechanism of anthocyanin biosynthesis in oil crop Brassica napus remains to be revealed.

RESULTS

In this study, we identified an anthocyanin repressor BnCPC in B. napus. BnCPC encoded a typical R3-MYB protein containing a conserved [D/E]Lx2[R/K]x3Lx6Lx3R motif for interaction with bHLH proteins. Overexpression of BnCPC in B. napus inhibited anthocyanin accumulation, especially under anthocyanin inducible conditions. Protein-protein interaction and dual-luciferase assays confirmed that BnCPC could compete with BnPAP1 to interact with bHLHs (BnTT8 and BnEGL3), and repress the expression of anthocyanin biosynthetic genes (e.g., BnDFR) that activated by MBW complexes. Moreover, we found BnCPC inhibited the MBW complex-induced BnCPC activity.

CONCLUSIONS

Overall, this research demonstrated that BnCPC repressed anthocyanin biosynthesis by affecting the formation of MBW complex, and formed a feedback loop to regulate anthocyanin accumulation in B. napus.

Effects of light quality, photoperiod, CO2 concentration, and air temperature on chlorogenic acid and rutin accumulation in young lettuce plants

Environmental stimuli modulate plant metabolite accumulation, facilitating adaptation to stressful conditions. In this study, the effects of blue and red light, photoperiod, CO2 concentration, and air temperature on the chlorogenic acid (CGA) and rutin contents of lettuce (Lactuca sativa L.) were evaluated. Under continuous blue light and a high CO2 concentration (1000 ppm), the CGA level increased. The increased expression of phenylalanine ammonia-lyase (PAL) and activity of its product were correlated with high expression of cinnamate 4-hydroxylase (C4H) and coumarate 3-hydroxylase (C3H). Furthermore, changes in PAL activity altered the CGA content in lettuce exposed to the three environmental factors, blue light, continuous lighting and high CO2 concentration. In addition, the expression levels of genes related to flavonoid biosynthesis increased in accordance with the promotion of CGA accumulation by the environmental factors. Under continuous blue light, 400 ppm CO2 promoted rutin accumulation to a greater degree compared to 1000 ppm CO2, by downregulating DFR expression. Low air temperature induced CGA accumulation in lettuce grown under continuous blue light and 1000 ppm CO2. Therefore, light quality, photoperiod, CO2 concentration, and air temperature exert synergistic effects on the CGA and rutin contents of lettuce by modulating activity in the corresponding biosynthesis pathways.

Comparison of Metabolome and Transcriptome of Flavonoid Biosynthesis in Two Colors of Coreopsis tinctoria Nutt

Coreopsis tinctoria Nutt. (C. tinctoria) has a long history of application and high economic and medicinal value. Flavonoids, the main active components of C. tinctoria, are widely studied in pharmacology and food development. However, the flavonoid biosynthesis pathway in C. tinctoria is unclear. In this study, we comprehensively compared the transcriptomes and metabolite profiles of two colors of C. tinctoria flowers (LS and JS) at different flowering stages. A total of 165 flavonoids (46 flavonoids, 42 flavonols, 22 anthocyanins, 18 chalcones, 12 dihydroflavonols, nine isoflavones, eight dihydroflavonoids, six flavanols, and two tannins) were identified in LS and JS at different flowering stages. Thirty-three metabolites (11 anthocyanins, 11 flavonols, seven flavonoids, two dihydroflavonols, one dihydroflavone, and one chalcone) were found to be statistically significantly different in the LS vs. JS groups. LS flowers accumulated higher levels of 10 anthocyanins (seven cyanidins and three pelargonidins) than JS flowers. Furthermore, candidate genes related to the regulation of flavonoid and anthocyanin synthesis were identified and included 28 structural genes (especially F3H, Cluster-28756.299649, and 3GT, Cluster-28756.230942) in LS and JS, six key differentially expressed transcription factors (especially MYB90a, Cluster-28756.143139) in LS and JS, and 17 other regulators (mainly including transporter proteins and others) in LS. Our results provide valuable information for further studies on the mechanism underlying flavonoid biosynthesis in C. tinctoria.

Incipient diploidization of the medicinal plant Perilla within 10,000 years

Perilla is a young allotetraploid Lamiaceae species widely used in East Asia as herb and oil plant. Here, we report the high-quality, chromosome-scale genomes of the tetraploid (Perilla frutescens) and the AA diploid progenitor (Perilla citriodora). Comparative analyses suggest post Neolithic allotetraploidization within 10,000 years, and nucleotide mutation in tetraploid is 10% more than in diploid, both of which are dominated by G:C → A:T transitions. Incipient diploidization is characterized by balanced swaps of homeologous segments, and subsequent homeologous exchanges are enriched towards telomeres, with excess of replacements of AA genes by fractionated BB homeologs. Population analyses suggest that the crispa lines are close to the nascent tetraploid, and involvement of acyl-CoA: lysophosphatidylcholine acyltransferase gene for high α-linolenic acid content of seed oil is revealed by GWAS. These resources and findings provide insights into incipient diploidization and basis for breeding improvement of this medicinal plant.

Integrated Metabolomic and Transcriptomic Analysis Reveals the Flavonoid Regulatory Network by Eutrema EsMYB90

Flavonoids are representative secondary metabolites with different metabolic functions in plants. Previous study found that ectopic expression of EsMYB90 from Eutremasalsugineum could strongly increase anthocyanin content in transgenic tobacco via regulating the expression of anthocyanin biosynthesis genes. In the present research, metabolome analysis showed that there existed 130 significantly differential metabolites, of which 23 metabolites enhanced more than 1000 times in EsMYB90 transgenic tobacco leaves relative to the control, and the top 10 of the increased metabolites included caffeic acid, cyanidin O-syringic acid, myricetin and naringin. A total of 50 markedly differential flavonoids including flavones (14), flavonols (13), flavone C-glycosides (9), flavanones (7), catechin derivatives (5), anthocyanins (1) and isoflavone (1) were identified, of which 46 metabolites were at a significantly enhanced level. Integrated analysis of metabolome and transcriptome revealed that ectopic expression of EsMYB90 in transgenic tobacco leaves is highly associated with the prominent up-regulation of 16 flavonoid metabolites and the corresponding 42 flavonoid biosynthesis structure genes in phenylpropanoid/flavonoid pathways. Dual luciferase assay documented that EsMYB90 strongly activated the transcription of NtANS and NtDFR genes via improving their promoter activity in transiently expressed tobacco leaves, suggesting that EsMYB90 functions as a key regulator on anthocyanin and flavonoid biosynthesis. Taken together, the crucial regulatory role of EsMYB90 on enhancing many flavonoid metabolite levels is clearly demonstrated via modulating flavonoid biosynthesis gene expression in the leaves of transgenic tobacco, which extends our understanding of the regulating mechanism of MYB transcription factor in the phenylpropanoid/flavonoid pathways and provides a new clue and tool for further investigation and genetic engineering of flavonoid metabolism in plants.

A genome-wide association study uncovers a critical role of the RsPAP2 gene in red-skinned Raphanus sativus L

Radish (Raphanus sativus L.) taproot contains high concentrations of flavonoids, including anthocyanins (ATCs), in red-skinned genotypes. However, little information on the genetic regulation of ATC biosynthesis in radish is available. A genome-wide association study of radish red skin color was conducted using whole-genome sequencing data derived from 179 radish genotypes. The R2R3-MYB transcription factor production of anthocyanin pigment 2 (PAP2) gene was found in the region associated with a leading SNP located on chromosome 2. The amino acid sequence encoded by the RsPAP2 gene was different from those of the other published RsMYB genes responsible for the red skin color of radish. The overexpression of the RsPAP2 gene resulted in ATC accumulation in Arabidopsis and radish, which was accompanied by the upregulation of several ATC-related structural genes. RsPAP2 was found to bind the RsUFGT and RsTT8 promoters, as shown by a dual-luciferase reporter system and a yeast one-hybrid assay. The promoter activities of the RsANS, RsCHI, RsPAL, and RsUFGT genes could be strongly activated by coinfiltration with RsPAP2 and RsTT8. These findings showed the effectiveness of GWAS in identifying candidate genes in radish and demonstrated that RsPAP2 could (either directly or together with its cofactor RsTT8) regulate the transcript levels of ATC-related genes to promote ATC biosynthesis, facilitating the genetic enhancement of ATC contents and other related traits in radish.

Transcriptomic profiling of purple broccoli reveals light-induced anthocyanin biosynthetic signaling and structural genes

Purple Broccoli (Brassica oleracea L. var italica) attracts growing attention as a functional food. Its purple coloration is due to high anthocyanin amounts. Light represents a critical parameter affecting anthocyanins biosynthesis. In this study, 'Purple Broccoli', a light-responding pigmentation cultivar, was assessed for exploring the mechanism underlying light-induced anthocyanin biosynthesis by RNA-Seq. Cyanidin, delphinidin and malvidin derivatives were detected in broccoli head samples. Shading assays and RNA-seq analysis identified the flower head as more critical organ compared with leaves. Anthocyanin levels were assessed at 0, 7 and 11 days, respectively, with further valuation by RNA-seq under head-shading and light conditions. RNA sequences were de novo assembled into 50,329 unigenes, of which 38,701 were annotated against four public protein databases. Cluster analysis demonstrated that anthocyanin/phenylpropanoid biosynthesis, photosynthesis, and flavonoid biosynthesis in cluster 8 were the main metabolic pathways regulated by light and had showed associations with flower head growth. A total of 2,400 unigenes showed differential expression between the light and head-shading groups in cluster 8, including 650 co-expressed, 373 specifically expressed under shading conditions and 1,377 specifically expressed under normal light. Digital gene expression (DGE) analysis demonstrated that light perception and the signal transducers CRY3 and HY5 may control anthocyanin accumulation. Following shading, 15 structural genes involved in anthocyanin biosynthesis were downregulated, including PAL, C4H, 4CL, CHS, CHI, F3H and DFR. Moreover, six BoMYB genes (BoMYB6-1, BoMYB6-2, BoMYB6-3, BoMYB6-4, BoMYBL2-1 and BoMYBL2-2) and three BobHLH genes (BoTT8_5-1, BoTT8_5-2 and BoEGL5-3) were critical transcription factors controlling anthocyanin accumulation under light conditions. Based on these data, a light-associated anthocyanin biosynthesis pathway in Broccoli was proposed. This information could help improve broccoli properties, providing novel insights into the molecular mechanisms underpinning light-associated anthocyanin production in purple vegetables.

Identification of the Eutrema salsugineum EsMYB90 gene important for anthocyanin biosynthesis

BACKGROUND

Anthocyanins contribute to coloration and antioxidation effects in different plant tissues. MYB transcription factors have been demonstrated to be a key regulator for anthocyanin synthesis in many plants. However, little information was available about the MYB genes in the halophyte species Eutrema salsugineum.

RESULT

Here we report the identification of an important anthocyanin biosynthesis regulator EsMYB90 from Eutrema salsugineum, which is a halophyte tolerant to multiple abiotic stresses. Our phylogenetic and localization analyses supported that EsMYB90 is an R2R3 type of MYB transcription factor. Ectopic expression of EsMYB90 in tobacco and Arabidopsis enhanced pigmentation and anthocyanin accumulation in various organs. The transcriptome analysis revealed that 42 genes upregulated by EsMYB90 in 35S:EsMYB90 tobacco transgenic plants are required for anthocyanin biosynthesis. Moreover, our qRT-PCR results showed that EsMYB90 promoted expression of early (PAL, CHS, and CHI) and late (DFR, ANS, and UFGT) anthocyanin biosynthesis genes in stems, leaves, and flowers of 35S:EsMYB90 tobacco transgenic plants.

CONCLUSIONS

Our results indicated that EsMYB90 is a MYB transcription factor, which regulates anthocyanin biosynthesis genes to control anthocyanin biosynthesis. Our work provides a new tool to enhance anthocyanin production in various plants.

Identification of the Eutrema salsugineum EsMYB90 gene important for anthocyanin biosynthesis

BACKGROUND

Anthocyanins contribute to coloration and antioxidation effects in different plant tissues. MYB transcription factors have been demonstrated to be a key regulator for anthocyanin synthesis in many plants. However, little information was available about the MYB genes in the halophyte species Eutrema salsugineum.

RESULT

Here we report the identification of an important anthocyanin biosynthesis regulator EsMYB90 from Eutrema salsugineum, which is a halophyte tolerant to multiple abiotic stresses. Our phylogenetic and localization analyses supported that EsMYB90 is an R2R3 type of MYB transcription factor. Ectopic expression of EsMYB90 in tobacco and Arabidopsis enhanced pigmentation and anthocyanin accumulation in various organs. The transcriptome analysis revealed that 42 genes upregulated by EsMYB90 in 35S:EsMYB90 tobacco transgenic plants are required for anthocyanin biosynthesis. Moreover, our qRT-PCR results showed that EsMYB90 promoted expression of early (PAL, CHS, and CHI) and late (DFR, ANS, and UFGT) anthocyanin biosynthesis genes in stems, leaves, and flowers of 35S:EsMYB90 tobacco transgenic plants.

CONCLUSIONS

Our results indicated that EsMYB90 is a MYB transcription factor, which regulates anthocyanin biosynthesis genes to control anthocyanin biosynthesis. Our work provides a new tool to enhance anthocyanin production in various plants.

Overexpression of the trehalose-6-phosphate phosphatase family gene AtTPPF improves the drought tolerance of Arabidopsis thaliana

BACKGROUND

Trehalose-6-phosphate phosphatases (TPPs), which are encoded by members of the TPP gene family, can improve the drought tolerance of plants. However, the molecular mechanisms underlying the dynamic regulation of TPP genes during drought stress remain unclear. In this study, we explored the function of an Arabidopsis TPP gene by conducting comparative analyses of a loss-of-function mutant and overexpression lines.

RESULTS

The loss-of-function mutation of Arabidopsis thaliana TPPF, a member of the TPP gene family, resulted in a drought-sensitive phenotype, while a line overexpressing TPPF showed significantly increased drought tolerance and trehalose accumulation. Compared with wild-type plants, tppf1 mutants accumulated more H₂O₂ under drought, while AtTPPF-overexpressing plants accumulated less H₂O₂ under drought. Overexpression of AtTPPF led to increased contents of trehalose, sucrose, and total soluble sugars under drought conditions; these compounds may play a role in scavenging reactive oxygen species. Yeast one-hybrid and luciferase activity assays revealed that DREB1A could bind to the DRE/CRT element within the AtTPPF promoter and activate the expression of AtTPPF. A transcriptome analysis of the TPPF-overexpressing plants revealed that the expression levels of drought-repressed genes involved in electron transport activity and cell wall modification were upregulated, while those of stress-related transcription factors related to water deprivation were downregulated. These results indicate that, as well as its involvement in regulating trehalose and soluble sugars, AtTPPF is involved in regulating the transcription of stress-responsive genes.

CONCLUSION

AtTPPF functions in regulating levels of trehalose, reactive oxygen species, and sucrose levels during drought stress, and the expression of AtTPPF is activated by DREB1A in Arabidopsis. These findings shed light on the molecular mechanism by which AtTPPF regulates the response to drought stress.

MdMYBL2 helps regulate cytokinin-induced anthocyanin biosynthesis in red-fleshed apple (Malus sieversii f. niedzwetzkyana) callus

Anthocyanin biosynthesis is induced by cytokinins, and is regulated by MYB transcription factors. However, the underlying molecular mechanisms have not been fully characterised. In the present study, red-fleshed apple callus were induced from the leaves of an R6/R6 homozygous line, which was the hybrid offspring of Malus sieversii f. niedzwetzkyana and 'Fuji'. We analysed the callus anthocyanin contents in response to different cytokinin concentrations. We observed that cytokinin treatments upregulated the expression of anthocyanin structural genes MdDFR and MdUFGT and transcription factor genes MdMYB10 and MdbHLH3. Additionally, the expression of MdMYBL2, which encodes the bHLH and EAR motifs, was inhibited by cytokinin treatments. The MdMYBL2-overexpressing callus had lower anthocyanin contents than the wild-type controls. We noted that the expression levels of anthocyanin biosynthesis structural genes MdDFR and MdUFGT and transcription factor genes MdMYB10 and MdbHLH3 were strongly suppressed in the transgenic callus. Subsequent yeast two-hybrid, bimolecular fluorescence complementation, and pull-down assays indicated that MdMYBL2 interacts with MdbHLH3, which may influence the expression of anthocyanin biosynthesis-related genes. Our findings may provide new insights into how MYB transcription factors influence the cytokinin-regulated anthocyanin biosynthesis in red-fleshed apples.

Discovering variation of secondary metabolite diversity and its relationship with disease resistance in Cornus florida L

Understanding intraspecific relationships between genetic and functional diversity is a major goal in the field of evolutionary biology and is important for conserving biodiversity. Linking intraspecific molecular patterns of plants to ecological pressures and trait variation remains difficult due to environment-driven plasticity. Next-generation sequencing, untargeted liquid chromatography-mass spectrometry (LC-MS) profiling, and interdisciplinary approaches integrating population genomics, metabolomics, and community ecology permit novel strategies to tackle this problem. We analyzed six natural populations of the disease-threatened Cornus florida L. from distinct ecological regions using genotype-by-sequencing markers and LC-MS-based untargeted metabolite profiling. We tested the hypothesis that higher genetic diversity in C. florida yielded higher chemical diversity and less disease susceptibility (screening hypothesis), and we also determined whether genetically similar subpopulations were similar in chemical composition. Most importantly, we identified metabolites that were associated with candidate loci or were predictive biomarkers of healthy or diseased plants after controlling for environment. Subpopulation clustering patterns based on genetic or chemical distances were largely congruent. While differences in genetic diversity were small among subpopulations, we did observe notable similarities in patterns between subpopulation averages of rarefied-allelic and chemical richness. More specifically, we found that the most abundant compound of a correlated group of putative terpenoid glycosides and derivatives was correlated with tree health when considering chemodiversity. Random forest biomarker and genomewide association tests suggested that this putative iridoid glucoside and other closely associated chemical features were correlated to SNPs under selection.

Mapping the Arabidopsis Metabolic Landscape by Untargeted Metabolomics at Different Environmental Conditions

Metabolic genome-wide association studies (mGWAS), whereupon metabolite levels are regarded as traits, can help unravel the genetic basis of metabolic networks. A total of 309 Arabidopsis accessions were grown under two independent environmental conditions (control and stress) and subjected to untargeted LC-MS-based metabolomic profiling; levels of the obtained hydrophilic metabolites were used in GWAS. Our two-condition-based GWAS for more than 3000 semi-polar metabolites resulted in the detection of 123 highly resolved metabolite quantitative trait loci (p ≤ 1.0E-08), 24.39% of which were environment-specific. Interestingly, differently from natural variation in Arabidopsis primary metabolites, which tends to be controlled by a large number of small-effect loci, we found several major large-effect loci alongside a vast number of small-effect loci controlling variation of secondary metabolites. The two-condition-based GWAS was followed by integration with network-derived metabolite-transcript correlations using a time-course stress experiment. Through this integrative approach, we selected 70 key candidate associations between structural genes and metabolites, and experimentally validated eight novel associations, two of them showing differential genetic regulation in the two environments studied. We demonstrate the power of combining large-scale untargeted metabolomics-based GWAS with time-course-derived networks both performed under different abiotic environments for identifying metabolite-gene associations, providing novel global insights into the metabolic landscape of Arabidopsis.