Integrated multi-omics analyses and genome-wide association studies reveal prime candidate genes of metabolic and vegetative growth variation in canola.

Genome-wide association studies (GWAS) identified thousands of genetic loci associated with complex plant traits, including many traits of agronomical importance. However, functional interpretation of GWAS results remains challenging because of large candidate regions due to linkage disequilibrium. High-throughput omics technologies, such as genomics, transcriptomics, proteomics and metabolomics open new avenues for integrative systems biological analyses and help to nominate systems information supported (prime) candidate genes. In the present study, we capitalise on a diverse canola population with 477 spring-type lines which was previously analysed by high-throughput phenotyping of growth-related traits and by RNA sequencing and metabolite profiling for multi-omics-based hybrid performance prediction. We deepened the phenotypic data analysis, now providing 123 time-resolved image-based traits, to gain insight into the complex relations during early vegetative growth and reanalysed the transcriptome data based on the latest Darmor-bzh v10 genome assembly. Genome-wide association testing revealed 61 298 robust quantitative trait loci (QTL) including 187 metabolite QTL, 56814 expression QTL and 4297 phenotypic QTL, many clustered in pronounced hotspots. Combining information about QTL colocalisation across omics layers and correlations between omics features allowed us to discover prime candidate genes for metabolic and vegetative growth variation. Prioritised candidate genes for early biomass accumulation include A06p05760.1_BnaDAR (PIAL1), A10p16280.1_BnaDAR, C07p48260.1_BnaDAR (PRL1) and C07p48510.1_BnaDAR (CLPR4). Moreover, we observed unequal effects of the Brassica A and C subgenomes on early biomass production.

Can biochemical traits bridge the gap between genomics and plant performance? A study in rice under drought.

The possibility of introducing metabolic/biochemical phenotyping to complement genomics-based predictions in breeding pipelines has been considered for years. Here we examine to what extent and under what environmental conditions metabolic/biochemical traits can effectively contribute to understanding and predicting plant performance. In this study, multivariable statistical models based on flag leaf central metabolism and oxidative stress status were used to predict grain yield (GY) performance for 271 indica rice (Oryza sativa) accessions grown in the field under well-watered and reproductive stage drought conditions. The resulting models displayed significantly higher predictability than multivariable models based on genomic data for the prediction of GY under drought (Q2 = 0.54-0.56 versus 0.35) and for stress-induced GY loss (Q2 = 0.59-0.64 versus 0.03-0.06). Models based on the combined datasets showed predictabilities similar to metabolic/biochemical-based models alone. In contrast to genetic markers, models with enzyme activities and metabolite values also quantitatively integrated the effect of physiological differences such as plant height on GY. The models highlighted antioxidant enzymes of the ascorbate-glutathione cycle and a lipid oxidation stress marker as important predictors of rice GY stability under drought at the reproductive stage, and these stress-related variables were more predictive than leaf central metabolites. These findings provide evidence that metabolic/biochemical traits can integrate dynamic cellular and physiological responses to the environment and can help bridge the gap between the genome and the phenome of crops as predictors of GY performance under drought.

Downy mildew resistance is genetically mediated by prophylactic production of phenylpropanoids in hop.

Downy mildew in hop (Humulus lupulus L.) is caused by Pseudoperonospora humuli and generates significant losses in quality and yield. To identify the biochemical processes that confer natural downy mildew resistance (DMR), a metabolome- and genome-wide association study was performed. Inoculation of a high density genotyped F1 hop population (n = 192) with the obligate biotrophic oomycete P. humuli led to variation in both the levels of thousands of specialized metabolites and DMR. We observed that metabolites of almost all major phytochemical classes were induced 48 hr after inoculation. But only a small number of metabolites were found to be correlated with DMR and these were enriched with phenylpropanoids. These metabolites were also correlated with DMR when measured from the non-infected control set. A genome-wide association study revealed co-localization of the major DMR loci and the phenylpropanoid pathway markers indicating that the major contribution to resistance is mediated by these metabolites in a heritable manner. The application of three putative prophylactic phenylpropanoids led to a reduced degree of leaf infection in susceptible genotypes, confirming their protective activity either directly or as precursors of active compounds.

Multi-omics-based prediction of hybrid performance in canola.

KEY MESSAGE: 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 F1 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.

Age-dependent loss of seed viability is associated with increased lipid oxidation and hydrolysis.

The accumulation of reactive oxygen species has been associated with a loss of seed viability. Therefore, we have investigated the germination ability of a range of seed stocks, including two wheat collections and one barley collection that had been dry-aged for 5-40 years. Metabolite profiling analysis revealed that the accumulation of glycerol was negatively correlated with the ability to germinate in all seed sets. Furthermore, lipid degradation products such as glycerol phosphates and galactose were accumulated in some seed sets. A quantitative analysis of nonoxidized and oxidized lipids was performed in the wheat seed set that showed the greatest variation in germination. This analysis revealed that the levels of fully acylated and nonoxidized storage lipids like triacylglycerols and structural lipids like phospho- and galactolipids were decreasing. Moreover, the abundance of oxidized variants and hydrolysed products such as mono-/diacylglycerols, lysophospholipids, and fatty acids accumulated as viability decreased. The proportional formation of oxidized and nonoxidized fatty acids provides evidence for an enzymatic hydrolysis of specifically oxidized lipids in dry seeds. The results link reactive oxygen species with lipid oxidation, structural damage, and death in long-term aged seeds.

Biomarkers for grain yield stability in rice under drought stress.

Crop yield stability requires an attenuation of the reduction of yield losses caused by environmental stresses such as drought. Using a combination of metabolomics and high-throughput colorimetric assays, we analysed central metabolism and oxidative stress status in the flag leaf of 292 indica rice (Oryza sativa) accessions. Plants were grown in the field and were, at the reproductive stage, exposed to either well-watered or drought conditions to identify the metabolic processes associated with drought-induced grain yield loss. Photorespiration, protein degradation, and nitrogen recycling were the main processes involved in the drought-induced leaf metabolic reprogramming. Molecular markers of drought tolerance and sensitivity in terms of grain yield were identified using a multivariate model based on the values of the metabolites and enzyme activities across the population. The model highlights the central role of the ascorbate-glutathione cycle, particularly dehydroascorbate reductase, in minimizing drought-induced grain yield loss. In contrast, malondialdehyde was an accurate biomarker for grain yield loss, suggesting that drought-induced lipid peroxidation is the major constraint under these conditions. These findings highlight new breeding targets for improved rice grain yield stability under drought.

Vacuolar processing enzyme 4 contributes to maternal control of grain size in barley by executing programmed cell death in the pericarp.

The angiosperm embryo and endosperm are limited in space because they grow inside maternal seed tissues. The elimination of cell layers of the maternal seed coat by programmed cell death (PCD) could provide space and nutrition to the filial organs. Using the barley (Hordeum vulgare L.) seed as a model, we elucidated the role of vacuolar processing enzyme 4 (VPE4) in cereals by using an RNAi approach and targeting the enzymatic properties of the recombinant protein. A comparative characterization of transgenic versus wild-type plants included transcriptional and metabolic profiling, flow cytometry, histology and nuclear magnetic imaging of grains. The recombinant VPE4 protein exhibited legumain and caspase-1 properties in vitro. Pericarp disintegration was delayed in the transgenic grains. Although the VPE4 gene and enzymatic activity was decreased in the early developing pericarp, storage capacity and the size of the endosperm and embryo were reduced in the mature VPE4-repressed grains. The persistence of the pericarp in the VPE4-affected grains constrains endosperm and embryo growth and leads to transcriptional reprogramming, perturbations in signalling and adjustments in metabolism. We conclude that VPE4 expression executes PCD in the pericarp, which is required for later endosperm filling, and argue for a role of PCD in maternal control of seed size in cereals.

Structure Annotation and Quantification of Wheat Seed Oxidized Lipids by High-Resolution LC-MS/MS.

Lipid oxidation is a process ubiquitous in life, but the direct and comprehensive analysis of oxidized lipids has been limited by available analytical methods. We applied high-resolution liquid chromatography-mass spectrometry (LC-MS) and tandem mass spectrometry (MS/MS) to quantify oxidized lipids (glycerides, fatty acids, phospholipids, lysophospholipids, and galactolipids) and implemented a platform-independent high-throughput-amenable analysis pipeline for the high-confidence annotation and acyl composition analysis of oxidized lipids. Lipid contents of 90 different naturally aged wheat (Triticum aestivum) seed stocks were quantified in an untargeted high-resolution LC-MS experiment, resulting in 18,556 quantitative mass-to-charge ratio features. In a posthoc liquid chromatography-tandem mass spectrometry experiment, high-resolution MS/MS spectra (5 mD accuracy) were recorded for 8,957 out of 12,080 putatively monoisotopic features of the LC-MS data set. A total of 353 nonoxidized and 559 oxidized lipids with up to four additional oxygen atoms were annotated based on the accurate mass recordings (1.5 ppm tolerance) of the LC-MS data set and filtering procedures. MS/MS spectra available for 828 of these annotations were analyzed by translating experimentally known fragmentation rules of lipids into the fragmentation of oxidized lipids. This led to the identification of 259 nonoxidized and 365 oxidized lipids by both accurate mass and MS/MS spectra and to the determination of acyl compositions for 221 nonoxidized and 295 oxidized lipids. Analysis of 15-year aged wheat seeds revealed increased lipid oxidation and hydrolysis in seeds stored in ambient versus cold conditions.

A naturally occurring promoter polymorphism of the Arabidopsis FUM2 gene causes expression variation, and is associated with metabolic and growth traits.

Fumarate and malate are known intermediates of the TCA cycle, a mitochondrial metabolic pathway generating NADH for respiration. Arabidopsis thaliana and other Brassicaceae contain an additional cytosolic fumarase (FUM2) that functions in carbon assimilation and nitrogen use. Here, we report the identification of a hitherto unknown FUM2 promoter insertion/deletion (InDel) polymorphism found between the Col-0 and C24 accessions, which also divides a large number of Arabidopsis accessions carrying either the Col-0 or the C24 allele. The polymorphism consists of two stretches of 2.1 and 3.8 kb, which are both absent from the promotor region of Col-0 FUM2. By analysing mutants as well as mapping and natural populations with contrasting FUM2 alleles, the promotor insertion was linked to reduced FUM2 mRNA expression, reduced fumarase activity and reduced fumarate/malate ratio in leaves. In a large population of 174 natural accessions, the polymorphism was also found to be associated with the fumarate/malate ratio, malate and fumarate levels, and with dry weight at 15 days after sowing (DAS). The association with biomass production was confirmed in an even larger (251) accession population for dry weight at 22 DAS. The dominant Col-0 allele that results in increased fumarate/malate ratios and enhanced biomass production is predominantly found in central/eastern European accessions, whereas the C24 type allele is prevalent on the Iberian Peninsula, west of the Rhine and in the British Isles. Our findings support the role of FUM2 in diurnal carbon storage, and point to a growth advantage of accessions carrying the FUM2 Col-0 allele.

Genetic dissection of metabolite variation in Arabidopsis seeds: evidence for mQTL hotspots and a master regulatory locus of seed metabolism.

To gain insight into genetic factors controlling seed metabolic composition and its relationship to major seed properties, an Arabidopsis recombinant inbred line (RIL) population, derived from accessions Col-0 and C24, was studied using an MS-based metabolic profiling approach. Relative intensities of 311 polar primary metabolites were used to identify associated genomic loci and to elucidate their interactions by quantitative trait locus (QTL) mapping. A total of 786 metabolic QTLs (mQTLs) were unequally distributed across the genome, forming several hotspots. For the branched-chain amino acid leucine, mQTLs and candidate genes were elucidated in detail. Correlation studies displayed links between metabolite levels, seed protein content, and seed weight. Principal component analysis revealed a clustering of samples, with PC1 mapping to a region on the short arm of chromosome IV. The overlap of this region with mQTL hotspots indicates the presence of a potential master regulatory locus of seed metabolism. As a result of database queries, a series of candidate regulatory genes, including bZIP10, were identified within this region. Depending on the search conditions, metabolic pathway-derived candidate genes for 40-61% of tested mQTLs could be determined, providing an extensive basis for further identification and characterization of hitherto unknown genes causal for natural variation of Arabidopsis seed metabolism.

Down-regulation of the sucrose transporters HvSUT1 and HvSUT2 affects sucrose homeostasis along its delivery path in barley grains.

Sucrose transport and partitioning are crucial for seed filling. While many plasma-membrane-localised sucrose transporters (SUT1 family members) have been analysed in seeds, the functions of vacuolar SUT2 members are still obscure. In barley grains, expression of HvSUT1 and HvSUT2 overlap temporally and spatially, suggesting concerted functions to regulate sucrose homeostasis. Using HvSUT2-RNAi plants, we found that grains were also deficient in HvSUT1 expression and seemingly sucrose-limited during mid-to-late grain filling. Transgenic endosperms accumulated less starch and dry weight, although overall sucrose and hexose contents were higher. Comprehensive transcript and metabolite profiling revealed that genes related to glycolysis, the tricarboxylic acid cycle, starch and amino acid synthesis, grain maturation, and abscisic acid signalling were down-regulated together with most glycolytic intermediates and amino acids. Sucrose was increased along the sucrose delivery route in the nucellar projection, the endosperm transfer cells, and the starchy endosperm, indicating that suppressed transporter activity diminished sucrose efflux from vacuoles, which generated sugar deficiency in the cytoplasm. Thus, endosperm vacuoles may buffer sucrose concentrations to regulate homeostasis at grain filling. Transcriptional changes revealed that limited endosperm sucrose initiated sugar starvation responses, such as sugar recycling from starch, hemicelluloses and celluloses together with vacuolar protein degradation, thereby supporting formation of nucleotide sugars. Barley endosperm cells can thus suppress certain pathways to retrieve resources to maintain essential cell functions.

Metabolic and transcriptional transitions in barley glumes reveal a role as transitory resource buffers during endosperm filling.

During grain filling in barley (Hordeum vulgare L. cv. Barke) reserves are remobilized from vegetative organs. Glumes represent the vegetative tissues closest to grains, senesce late, and are involved in the conversion of assimilates. To analyse glume development and metabolism related to grain filling, parallel transcript and metabolite profiling in glumes and endosperm were performed, showing that glume metabolism and development adjusts to changing grain demands, reflected by specific signatures of metabolite and transcript abundances. Before high endosperm sink strength is established by storage product accumulation, glumes form early, intermediary sink organs, shifting then to remobilizing and exporting source organs. Metabolic and transcriptional transitions occur at two phases: first, at the onset of endosperm filling, as a consequence of endosperm sink activity and assimilate depletion in endosperm and vascular tissues; second, at late grain filling, by developmental ageing and senescence. Regulation of and transition between phases are probably governed by specific NAC and WRKY transcription factors, and both abscisic and jasmonic acid, and are accompanied by changed expression of specific nitrogen transporters. Expression and metabolite profiling suggest glume-specific mechanisms of assimilate conversion and translocation. In summary, grain filling and endosperm sink strength coordinate phase changes in glumes via metabolic, hormonal, and transcriptional control. This study provides a comprehensive view of barley glume development and metabolism, and identifies candidate genes and associated pathways, potentially important for breeding improved grain traits.

Transcription factor sensor system for parallel quantification of metabolites on-chip.

Steadily growing demands for identification and quantification of cellular metabolites in higher throughput have brought a need for new analytical technologies. Here, we developed a synthetic biological sensor system for quantifying metabolites from biological cell samples. For this, bacterial transcription factors were exploited, which bind to or dissociate from regulatory DNA elements in response to physiological changes in the cellular metabolite concentration range. Representatively, the bacterial pyruvate dehydrogenase (PdhR), trehalose (TreR), and l-arginine (ArgR) repressor proteins were functionalized to detect pyruvate, trehalose-6-phosphate (T6P), and arginine concentration in solution. For each transcription factor the mutual binding behavior between metabolite and DNA, their working range, and othogonality were determined. High-throughput, parallel processing, and automation were achieved through integration of the metabolic sensor system on a microfluidic large-scale integration (mLSI) chip platform. To demonstrate the functionality of the integrated metabolic sensor system, we measured diurnal concentration changes of pyruvate and the plant signaling molecule T6P within cell etxracts of Arabidopsis thaliana rosettes. The transcription factor sensor system is of generic nature and extendable on the microfluidic chip.

Differentiation of endosperm transfer cells of barley: a comprehensive analysis at the micro-scale.

Barley endosperm cells differentiate into transfer cells (ETCs) opposite the nucellar projection. To comprehensively analyse ETC differentiation, laser microdissection-based transcript and metabolite profiles were obtained from laser microdissected tissues and cell morphology was analysed. Flange-like secondary-wall ingrowths appeared between 5 and 7 days after pollination within the three outermost cell layers. Gene expression analysis indicated that ethylene-signalling pathways initiate ETC morphology. This is accompanied by gene activity related to cell shape control and vesicle transport, with abundant mitochondria and endomembrane structures. Gene expression analyses indicate predominant formation of hemicelluloses, glucuronoxylans and arabinoxylans, and transient formation of callose, together with proline and 4-hydroxyproline biosynthesis. Activation of the methylation cycle is probably required for biosynthesis of phospholipids, pectins and ethylene. Membrane microdomains involving sterols/sphingolipids and remorins are potentially involved in ETC development. The transcriptional activity of assimilate and micronutrient transporters suggests ETCs as the main uptake organs of solutes into the endosperm. Accordingly, the endosperm grows maximally after ETCs are fully developed. Up-regulated gene expression related to amino acid catabolism, C:N balances, carbohydrate oxidation, mitochondrial activity and starch degradation meets high demands for respiratory energy and carbohydrates, required for cell proliferation and wall synthesis. At 10 days after pollination, ETCs undergo further differentiation, potentially initiated by abscisic acid, and metabolism is reprogrammed as shown by activated storage and stress-related processes. Overall, the data provide a comprehensive view of barley ETC differentiation and development, and identify candidate genes and associated pathways.

A tyrosine aminotransferase involved in tocopherol synthesis in Arabidopsis.

The metabolic function of the predicted Arabidopsis tyrosine aminotransferase (TAT) encoded by the At5g53970 gene was studied using two independent knock-out mutants. Gas chromatography-mass spectrometry based metabolic profiling revealed a specific increase in tyrosine levels, supporting the proposed function of At5g53970 as a tyrosine-specific aminotransferase not involved in tyrosine biosynthesis, but rather in utilization of tyrosine for other metabolic pathways. The TAT activity of the At5g53970-encoded protein was verified by complementation of the Escherichia coli tyrosine auxotrophic mutant DL39, and in vitro activity of recombinantly expressed and purified At5g53970 was found to be specific for tyrosine. To investigate the physiological role of At5g53970, the consequences of reduction in tyrosine utilization on metabolic pathways having tyrosine as a substrate were analysed. We found that tocopherols were substantially reduced in the mutants and we conclude that At5g53970 encodes a TAT important for the synthesis of tocopherols in Arabidopsis.

The production of male-sterile wheat plants through split barnase expression is promoted by the insertion of introns and flexible peptide linkers.

The successful use of transgenic plants depends on the strong and stable expression of the heterologous genes. In this study, three introns (PSK7-i1 and PSK7-i3 from Petunia and UBQ10-i1 from Arabidopsis) were tested for their ability to enhance the tapetum-specific expression of a split barnase transgene. We also analyzed the effects of introducing multiple copies of flexible peptide linkers that bridged the fusion domains of the assembled protein. The barnase fragments were assembled into a functional cytotoxin via intein-mediated trans-splicing, thus leading to male sterility through pollen ablation. A total of 14 constructs carrying different combinations of introns and peptide linkers were transformed into wheat plants. The resulting populations (between 41 and 301 independent plants for each construct) were assayed for trait formation. Depending on which construct was used, there was an increase of up to fivefold in the proportion of plants exhibiting male sterility compared to the populations harboring unmodified constructs. Furthermore, the average barnase copy number in the plants displaying male sterility could be reduced. The metabolic profiles of male-sterile transgenic plants and non-transgenic plants were compared using gas chromatography-mass spectrometry. The profiles generated from leaf tissues displayed no differences, thus corroborating the anther specificity of barnase expression. The technical advances achieved in this study may be a valuable contribution for future improvement of transgenic crop systems.

Corrigendum: Facing energy limitations - approaches to increase basil (Ocimum basilicum L.) growth and quality by different increasing light intensities emitted by a broadband LED light spectrum (400-780 nm).

[This corrects the article DOI: 10.3389/fpls.2022.1055352.].

Facing energy limitations - approaches to increase basil (Ocimum basilicum L.) growth and quality by different increasing light intensities emitted by a broadband LED light spectrum (400-780 nm).

Based on the current trend towards broad-bandwidth LED light spectra for basil productions in multi-tiered controlled-environment horticulture, a recently developed white broad-bandwidth LED light spectrum (400-780 nm) including far-red wavelengths with elevated red and blue light fractions was employed to cultivate basil. Four Ocimum basilicum L. cultivars (cv. Anise, cv. Cinnamon, cv. Dark Opal and cv. Thai Magic) were exposed to two different rising light intensity conditions (ILow and IHigh). In dependence of the individual cultivar-specific plant height increase over time, basil cultivars were exposed to light intensities increasing from ~ 100 to ~ 200 µmol m-2 s-1 under ILow, and from 200 to 400 µmol m-2 s-1 under IHigh (due to the exponential light intensity increases with decreasing proximity to the LED light fixtures). Within the first experiment, basils' morphological developments, biomass yields and time to marketability under both light conditions were investigated and the energy consumptions were determined to calculate the basils' light use efficiencies. In detail, cultivar-dependent differences in plant height, leaf and branch pair developments over time are described. In comparison to the ILow light conditions, IHigh resulted in accelerated developments and greater yields of all basil cultivars and expedited their marketability by 3-5 days. However, exposure to light intensities above ~ 300 µmol m-2 s-1 induced light avoidance responses in the green-leafed basil cultivars cv. Anise, cv. Cinnamon and cv. Thai Magic. In contrast, ILow resulted in consumer-preferred visual qualities and greater biomass efficiencies of the green-leafed basil cultivars and are discussed as a result of their ability to adapt well to low light conditions. Contrarily to the green-leafed cultivars, purple-leafed cv. Dark Opal developed insufficiently under ILow, but remained light-tolerant under IHigh, which is related to its high anthocyanin contents. In a second experiment, cultivars' volatile organic compound (VOC) contents and compositions over time were investigated. While VOC contents per gram of leaf dry matter gradually decreased in purple-leafed cv. Dark Opal between seedling stage to marketability, their contents gradually increased in the green cultivars. Regardless of the light treatment applied, cultivar-specific VOC compositions changed tremendously in a developmental stage-dependent manner.

The carrot monoterpene synthase gene cluster on chromosome 4 harbours genes encoding flavour-associated sabinene synthases.

In plants, low molecular weight terpenes produced by terpene synthases (TPS) contribute to multiple ecologically and economically important traits. The present study investigates a carrot terpene synthase gene cluster on chromosome 4 associated with volatile monoterpene production. Two carrot mutants, yellow and cola, which are contrasting in the content of low molecular weight terpenes, were crossed to develop an F2 mapping population. The mapping analysis revealed overlapping QTLs on chromosome 4 for sabinene, α-thujene, α-terpinene, γ-terpinene, terpinen-4-ol and 4-carene. The genomic region of this locus includes a cluster of five terpene synthase genes (DcTPS04, DcTPS26, DcTPS27, DcTPS54 and DcTPS55). DcTPS04 and DcTPS54 displayed genotype- and tissue-specific variation in gene expression. Based on the QTL mapping results and the gene expression patterns, DcTPS04 and DcTPS54 were selected for functional characterization. In vitro enzyme assays showed that DcTPS54 is a single-product enzyme catalysing the formation of sabinene, whereas DcTPS04 is a multiple-product terpene synthase producing α-terpineol as a major product and four additional products including sabinene, ß-limonene, ß-pinene and myrcene. Furthermore, we developed a functional molecular marker that could discriminate carrot genotypes with different sabinene content in a set of 85 accessions.