Evolution and functional divergence of glycosyltransferase genes shaped the quality and cold tolerance of tea plants.

Plant glycosyltransferases (UGTs) play a key role in plant growth and metabolism. Here, we examined the evolutionary landscape among UGTs in 28 fully sequenced species from early algae to angiosperms. Our findings revealed a distinctive expansion and contraction of UGTs in the G and H groups in tea (Camellia sinensis), respectively. Whole-genome duplication and tandem duplication events jointly drove the massive expansion of UGTs, and the interplay of natural and artificial selection has resulted in marked functional divergence within the G group of the sinensis-type tea population. In Cluster II of group G, differences in substrate selection (e.g., Abscisic Acid) of the enzymes encoded by UGT genes led to their functional diversification, and these genes influence tolerance to abiotic stresses such as low temperature and drought via different modes of positive and negative regulation, respectively. UGTs in Cluster III of the G group have diverse aroma substrate preferences, which contributes a diverse aroma spectrum of the sinensis-type tea population. All Cluster III genes respond to low-temperature stress, whereas UGTs within Cluster III-1, shaped by artificial selection, are unresponsive to drought. This suggests that artificial selection of tea plants focused on improving quality and cold tolerance as primary targets.

Kinetics of Secoisolariciresinol Glucosyltransferase LuUGT74S1 and Its Mutants.

The lignan secoisolariciresinol (SECO) diglucoside (SDG) is a phytoestrogen with diverse effects. LuUGT74S1 glucosylates SECO to SDG, whereby only small amounts of the monoglucoside SMG are formed intermediately, which exhibit increased activity. To identify critical amino acids that are important for enzymatic activity and the SMG/SDG ratio, 3D structural modeling and docking, as well as site-directed mutation studies, were performed. Enzyme assays with ten mutants revealed that four of them had identical kinetic data to LuUGT74S1, while three showed reduced and one increased catalytic efficiency kcat/Km. S82F and E189L substitutions resulted in the complete absence of activity. A17 and Q136 are crucial for the conversion of SMG to SDG as A17S and Q136F mutants exhibited the highest SMG/SDG ratios of 0.7 and 0.4. Kinetic analyses show that diglucosylation is an essentially irreversible reaction, while monoglycosylation is kinetically favored. The results lay the foundation for the biotechnological production of SMG.

2,4-Dihydroxybenzoic Acid, a Novel SA Derivative, Controls Plant Immunity via UGT95B17-Mediated Glucosylation: A Case Study in Camellia Sinensis.

The plant hormone salicylic acid (SA) plays critical roles in plant innate immunity. Several SA derivatives and associated modification are identified, whereas the range and modes of action of SA-related metabolites remain elusive. Here, the study discovered 2,4-dihydroxybenzoic acid (2,4-DHBA) and its glycosylated form as native SA derivatives in plants whose accumulation is largely induced by SA application and Ps. camelliae-sinensis (Pcs) infection. CsSH1, a 4/5-hydroxylase, catalyzes the hydroxylation of SA to 2,4-DHBA, and UDP-glucosyltransferase UGT95B17 catalyzes the formation of 2,4-DHBA glucoside. Down-regulation reduced the accumulation of 2,4-DHBA glucosides and enhanced the sensitivity of tea plants to Pcs. Conversely, overexpression of UGT95B17 increased plant disease resistance. The exogenous application of 2,4-DHBA and 2,5-DHBA, as well as the accumulation of DHBA and plant resistance comparison, indicate that 2,4-DHBA functions as a potentially bioactive molecule and is stored mainly as a glucose conjugate in tea plants, differs from the mechanism described in Arabidopsis. When 2,4-DHBA is applied exogenously, UGT95B17-silenced tea plants accumulated more 2,4-DHBA than SA and showed induced resistance to Pcs infection. These results indicate that 2,4-DHBA glucosylation positively regulates disease resistance and highlight the role of 2,4-DHBA as potentially bioactive molecule in the establishment of basal resistance in tea plants.

UGT89AC1-mediated quercetin glucosylation is induced upon herbivore damage and enhances Camellia sinensis resistance to insect feeding.

Quercetin is a key flavonol in tea plants (Camellia sinensis (L.) O. Kuntze) with various health benefits, and it often occurs in the form of glucosides. The roles of quercetin and its glucosylated forms in plant defense are generally not well-studied, and remain unknown in the defense of tea. Here, we found higher contents of quercetin glucosides and a decline of the aglucone upon Ectropis grisescens (E. grisescens) infestation of tea. Nine UGTs were strongly induced, among which UGT89AC1 exhibited the highest activity toward quercetin in vitro and in vivo. The mass of E. grisescens larvae that fed on plants with repressed UGT89AC1 or varieties with lower levels of UGT89AC1 was significantly lower than that of larvae fed on controls. Artificial diet supplemented with quercetin glucoside also reduced the larval growth rate, whereas artificial diet supplemented with free quercetin had no significant effect on larval growth. UGT89AC1 was located in both the cytoplasm and nucleus, and its expression was modulated by JA, JA-ILE, and MeJA. These findings demonstrate that quercetin glucosylation serves a defensive role in tea against herbivory. Our results also provide novel insights into the ecological relevance of flavonoid glycosides under biotic stress in plants.

Distribution of and Temporal Variation in Volatiles in Tea (Camellia sinensis) Flowers during the Opening Stages.

Tea (Camellia sinensis) flowers emit a large amount of volatiles that attract pollinators. However, few studies have characterized temporal and spatial variation in tea floral volatiles. To investigate the distribution of volatiles within tea flowers and their variation among opening stages, volatile components from different parts of tea flowers and different opening stages were collected by headspace solid-phase microextraction and analyzed by gas chromatography-mass spectrometry. A total of 51 volatile compounds of eight chemical classes were identified in the tea flowers. Volatile compounds were most abundant in tea flowers of the Shuchazao cultivar. Acetophenone, 1-phenylethanol, 2-phenylethanol, and benzyl alcohol were the most abundant volatiles. Terpenes were common in the sepals, and benzoids were common in the stamens. The fatty acid derivatives were mainly distributed in the pistils and receptacles and were less abundant in the petals, sepals, and stamens. During the opening phase of tea flowers, the volatile content increased 12-fold, which mainly stemmed from the increase in benzoids. These results enhance our understanding of the formation of volatiles in tea flowers.

A geraniol synthase regulates plant defense via alternative splicing in tea plants.

Geraniol is an important contributor to the pleasant floral scent of tea products and one of the most abundant aroma compounds in tea plants; however, its biosynthesis and physiological function in response to stress in tea plants remain unclear. The proteins encoded by the full-length terpene synthase (CsTPS1) and its alternative splicing isoform (CsTPS1-AS) could catalyze the formation of geraniol when GPP was used as a substrate in vitro, whereas the expression of CsTPS1-AS was only significantly induced by Colletotrichum gloeosporioides and Neopestalotiopsis sp. infection. Silencing of CsTPS1 and CsTPS1-AS resulted in a significant decrease of geraniol content in tea plants. The geraniol content and disease resistance of tea plants were compared when CsTPS1 and CsTPS1-AS were silenced. Down-regulation of the expression of CsTPS1-AS reduced the accumulation of geraniol, and the silenced tea plants exhibited greater susceptibility to pathogen infection than control plants. However, there was no significant difference observed in the geraniol content and pathogen resistance between CsTPS1-silenced plants and control plants in the tea plants infected with two pathogens. Further analysis showed that silencing of CsTPS1-AS led to a decrease in the expression of the defense-related genes PR1 and PR2 and SA pathway-related genes in tea plants, which increased the susceptibility of tea plants to pathogens infections. Both in vitro and in vivo results indicated that CsTPS1 is involved in the regulation of geraniol formation and plant defense via alternative splicing in tea plants. The results of this study provide new insights into geraniol biosynthesis and highlight the role of monoterpene synthases in modulating plant disease resistance via alternative splicing.

Volatile compound-mediated plant-plant interactions under stress with the tea plant as a model.

Plants respond to environmental stimuli via the release of volatile organic compounds (VOCs), and neighboring plants constantly monitor and respond to these VOCs with great sensitivity and discrimination. This sensing can trigger increased plant fitness and reduce future plant damage through the priming of their own defenses. The defense mechanism in neighboring plants can either be induced by activation of the regulatory or transcriptional machinery, or it can be delayed by the absorption and storage of VOCs for the generation of an appropriate response later. Despite much research, many key questions remain on the role of VOCs in interplant communication and plant fitness. Here we review recent research on the VOCs induced by biotic (i.e. insects and pathogens) and abiotic (i.e. cold, drought, and salt) stresses, and elucidate the biosynthesis of stress-induced VOCs in tea plants. Our focus is on the role of stress-induced VOCs in complex ecological environments. Particularly, the roles of VOCs under abiotic stress are highlighted. Finally, we discuss pertinent questions and future research directions for advancing our understanding of plant interactions via VOCs.

(Z)-3-Hexenol integrates drought and cold stress signaling by activating abscisic acid glucosylation in tea plants.

Cold and drought stresses severely limit crop production and can occur simultaneously. Although some transcription factors and hormones have been characterized in plants subjected each stress, the role of metabolites, especially volatiles, in response to cold and drought stress exposure is rarely studied due to lack of suitable models. Here, we established a model for studying the role of volatiles in tea (Camellia sinensis) plants experiencing cold and drought stresses simultaneously. Using this model, we showed that volatiles induced by cold stress promote drought tolerance in tea plants by mediating reactive oxygen species and stomatal conductance. Needle trap microextraction combined with GC-MS identified the volatiles involved in the crosstalk and showed that cold-induced (Z)-3-hexenol improved the drought tolerance of tea plants. In addition, silencing C. sinensis alcohol dehydrogenase 2 (CsADH2) led to reduced (Z)-3-hexenol production and significantly reduced drought tolerance in response to simultaneous cold and drought stress. Transcriptome and metabolite analyses, together with plant hormone comparison and abscisic acid (ABA) biosynthesis pathway inhibition experiments, further confirmed the roles of ABA in (Z)-3-hexenol-induced drought tolerance of tea plants. (Z)-3-Hexenol application and gene silencing results supported the hypothesis that (Z)-3-hexenol plays a role in the integration of cold and drought tolerance by stimulating the dual-function glucosyltransferase UGT85A53, thereby altering ABA homeostasis in tea plants. Overall, we present a model for studying the roles of metabolites in plants under multiple stresses and reveal the roles of volatiles in integrating cold and drought stresses in plants.

Genome-wide identification of UDP-glycosyltransferases in the tea plant (Camellia sinensis) and their biochemical and physiological functions.

Tea (Camellia sinensis) has been an immensely important commercially grown crop for decades. This is due to the presence of essential nutrients and plant secondary metabolites that exhibit beneficial health effects. UDP-glycosyltransferases (UGTs) play an important role in the diversity of such secondary metabolites by catalysing the transfer of an activated sugar donor to acceptor molecules, and thereby creating a huge variety of glycoconjugates. Only in recent years, thanks to the sequencing of the tea plant genome, have there been increased efforts to characterise the UGTs in C. sinensis to gain an understanding of their physiological role and biotechnological potential. Based on the conserved plant secondary product glycosyltransferase (PSPG) motif and the catalytically active histidine in the active site, UGTs of family 1 in C. sinensis are identified here, and shown to cluster into 21 groups in a phylogenetic tree. Building on this, our current understanding of recently characterised C. sinensis UGTs (CsUGTs) is highlighted and a discussion on future perspectives made.

Acceptors and Effectors Alter Substrate Inhibition Kinetics of a Plant Glucosyltransferase NbUGT72AY1 and Its Mutants.

One of the main obstacles in biocatalysis is the substrate inhibition (SI) of enzymes that play important roles in biosynthesis and metabolic regulation in organisms. The promiscuous glycosyltransferase UGT72AY1 from Nicotiana benthamiana is strongly substrate-inhibited by hydroxycoumarins (inhibitory constant Ki < 20 µM), but only weakly inhibited when monolignols are glucosylated (Ki > 1000 µM). Apocarotenoid effectors reduce the inherent UDP-glucose glucohydrolase activity of the enzyme and attenuate the SI by scopoletin derivatives, which could also be achieved by mutations. Here, we studied the kinetic profiles of different phenols and used the substrate analog vanillin, which has shown atypical Michaelis-Menten kinetics in previous studies, to examine the effects of different ligands and mutations on the SI of NbUGT72AY1. Coumarins had no effect on enzymatic activity, whereas apocarotenoids and fatty acids strongly affected SI kinetics by increasing the inhibition constant Ki. Only the F87I mutant and a chimeric version of the enzyme showed weak SI with the substrate vanillin, but all mutants exhibited mild SI when sinapaldehyde was used as an acceptor. In contrast, stearic acid reduced the transferase activity of the mutants to varying degrees. The results not only confirm the multi-substrate functionality of NbUGT72AY1, but also reveal that the enzymatic activity of this protein can be fine-tuned by external metabolites such as apocarotenoids and fatty acids that affect SI. Since these signals are generated during plant cell destruction, NbUGT72AY1 likely plays an important role in plant defense by participating in the production of lignin in the cell wall and providing direct protection through the formation of toxic phytoalexins.

Abscisic Acid Synthesis and Signaling during the Ripening of Raspberry (Rubus idaeus 'Heritage') Fruit.

The raspberry (Rubus idaeus L.) fruit is characterized by its richness in functional molecules and high nutritional value, but the high rate of fruit softening limits its quality during postharvest. Raspberry drupelets have a particular ripening regulation, depending partially on the effect of ethylene produced from the receptacle. However, the possible role of abscisic acid (ABA) in the modulation of quality parameters during the ripening of raspberry is unclear. This study characterized the fruit quality-associated parameters and hormonal contents during fruit development in two seasons. The quality parameters showed typical changes during ripening: a drastic loss of firmness, increase in soluble solids content, loss of acidity, and turning to a red color from the large green stage to fully ripe fruit in both seasons. A significant increase in the ABA content was observed during the ripening of drupelets and receptacles, with the higher content in the receptacle of ripe and overripe stages compared to the large green stage. Moreover, identification of ABA biosynthesis-(9-cis-epoxycarotenoid dioxygenase/NCED) and ABA receptor-related genes (PYRs-like receptors) showed three genes encoding RiNCEDs and nine genes for RiPYLs. The expression level of these genes increased from the large green stage to the full-ripe stage, specifically characterized by a higher expression of RiNCED1 in the receptacle tissue. This study reports a consistent concomitant increase in the ABA content and the expression of RiNCED1, RiPYL1, and RiPYL8 during the ripening of the raspberry fruit, thus supporting the role for ABA signaling in drupelets.

Apocarotenoids are allosteric effectors of a dimeric plant glycosyltransferase involved in defense and lignin formation.

Glycosyltransferases are nature's versatile tools to tailor the functionalities of proteins, carbohydrates, lipids, and small molecules by transferring sugars. Prominent substrates are hydroxycoumarins such as scopoletin, which serve as natural plant protection agents. Similarly, C13-apocarotenoids, which are oxidative degradation products of carotenoids/xanthophylls, protect plants by repelling pests and attracting pest predators. We show that C13-apocarotenoids interact with the plant glycosyltransferase NbUGT72AY1 and induce conformational changes in the enzyme catalytic center ultimately reducing its inherent UDP-α-d-glucose glucohydrolase activity and increasing its catalytic activity for productive hydroxycoumarin substrates. By contrast, C13-apocarotenoids show no effect on the catalytic activity toward monolignol lignin precursors, which are competitive substrates. In vivo studies in tobacco plants (Nicotiana benthamiana) confirmed increased glycosylation activity upon apocarotenoid supplementation. Thus, hydroxycoumarins and apocarotenoids represent specialized damage-associated molecular patterns, as they each provide precise information about the plant compartments damaged by pathogen attack. The molecular basis for the C13-apocarotenoid-mediated interplay of two plant protective mechanisms and their function as allosteric enhancers opens up potential applications of the natural products in agriculture and pharmaceutical industry.

QTL analysis for ascorbic acid content in strawberry fruit reveals a complex genetic architecture and association with GDP-L-galactose phosphorylase.

Strawberry (Fragaria × ananassa) fruits are an excellent source of L-ascorbic acid (AsA), a powerful antioxidant for plants and humans. Identifying the genetic components underlying AsA accumulation is crucial for enhancing strawberry nutritional quality. Here, we unravel the genetic architecture of AsA accumulation using an F1 population derived from parental lines 'Candonga' and 'Senga Sengana', adapted to distinct Southern and Northern European areas. To account for environmental effects, the F1 and parental lines were grown and phenotyped in five locations across Europe (France, Germany, Italy, Poland and Spain). Fruit AsA content displayed normal distribution typical of quantitative traits and ranged five-fold, with significant differences among genotypes and environments. AsA content in each country and the average in all of them was used in combination with 6,974 markers for quantitative trait locus (QTL) analysis. Environmentally stable QTLs for AsA content were detected in linkage group (LG) 3A, LG 5A, LG 5B, LG 6B and LG 7C. Candidate genes were identified within stable QTL intervals and expression analysis in lines with contrasting AsA content suggested that GDP-L-Galactose Phosphorylase FaGGP(3A), and the chloroplast-located AsA transporter gene FaPHT4;4(7C) might be the underlying genetic factors for QTLs on LG 3A and 7C, respectively. We show that recessive alleles of FaGGP(3A) inherited from both parental lines increase fruit AsA content. Furthermore, expression of FaGGP(3A) was two-fold higher in lines with high AsA. Markers here identified represent a useful resource for efficient selection of new strawberry cultivars with increased AsA content.

Subfunctionalization of a monolignol to a phytoalexin glucosyltransferase is accompanied by substrate inhibition.

Uridine diphosphate-dependent glycosyltransferases (UGTs) mediate the glycosylation of plant metabolites, thereby altering their physicochemical properties and bioactivities. Plants possess numerous UGT genes, with the encoded enzymes often glycosylating multiple substrates and some exhibiting substrate inhibition kinetics, but the biological function and molecular basis of these phenomena are not fully understood. The promiscuous monolignol/phytoalexin glycosyltransferase NbUGT72AY1 exhibits substrate inhibition (Ki) at 4 µM scopoletin, whereas the highly homologous monolignol StUGT72AY2 is inhibited at 190 µM. We therefore used hydrogen/deuterium exchange mass spectrometry and structure-based mutational analyses of both proteins and introduced NbUGT72AY1 residues into StUGT72AY2 and vice versa to study promiscuity and substrate inhibition of UGTs. A single F87I and chimeric mutant of NbUGT72AY1 showed significantly reduced scopoletin substrate inhibition, whereas its monolignol glycosylation activity was almost unaffected. Reverse mutations in StUGT72AY2 resulted in increased scopoletin glycosylation, leading to enhanced promiscuity, which was accompanied by substrate inhibition. Studies of 3D structures identified open and closed UGT conformers, allowing visualization of the dynamics of conformational changes that occur during catalysis. Previously postulated substrate access tunnels likely serve as drainage channels. The results suggest a two-site model in which the second substrate molecule binds near the catalytic site and blocks product release. Mutational studies showed that minor changes in amino acid sequence can enhance the promiscuity of the enzyme and add new capabilities such as substrate inhibition without affecting existing functions. The proposed subfunctionalization mechanism of expanded promiscuity may play a role in enzyme evolution and highlights the importance of promiscuous enzymes in providing new functions.

Comparative Physicochemical and Biochemical Characterization of Small-Molecule Glucosides.

Glycosylation of small molecules can significantly improve their physicochemical and biological properties. Only recently, decisive improvements in the biotechnological production of small-molecule glucosides (SMGs) have resulted in a large number of these compounds now being commercially available. In this study, we have analyzed a number of physical, chemical, and biological parameters of 31 SMGs, including solubility, stability, melting and pyrolysis points, partition coefficient log P, minimum inhibitory concentration against Escherichia coli (MIC), and enzymatic degradability. The properties such as water solubility, pH stability, and MICs of the glycosides were strongly dependent on the structures of the respective aglycones, which is why the SMG clustered according to their aglycones in most cases. Phenolic and furanone glucosides were readily hydrolyzed by saliva and skin microflora, whereas monoterpenol glycosides were poorer substrates for the enzymes involved. The results of this comparative analysis of SMGs provide valuable information for elucidating the biological functions of SMGs and the future technological applications of these useful natural products.

Single-cell transcriptome atlas reveals developmental trajectories and a novel metabolic pathway of catechin esters in tea leaves.

The tea plant is an economically important woody beverage crop. The unique taste of tea is evoked by certain metabolites, especially catechin esters, whereas their precise formation mechanism in different cell types remains unclear. Here, a fast protoplast isolation method was established and the transcriptional profiles of 16 977 single cells from 1st and 3rd leaves were investigated. We first identified 79 marker genes based on six isolated tissues and constructed a transcriptome atlas, mapped developmental trajectories and further delineated the distribution of different cell types during leaf differentiation and genes associated with cell fate transformation. Interestingly, eight differently expressed genes were found to co-exist at four branch points. Genes involved in the biosynthesis of certain metabolites showed cell- and development-specific characteristics. An unexpected catechin ester glycosyltransferase was characterized for the first time in plants by a gene co-expression network in mesophyll cells. Thus, the first single-cell transcriptional landscape in woody crop leave was reported and a novel metabolism pathway of catechin esters in plants was discovered.

Strawberry fruit FanCXE1 carboxylesterase is involved in the catabolism of volatile esters during the ripening process.

Volatile compounds produced during ripening of strawberry are key determinants of fruit quality and consumer preference. Strawberry volatiles are largely esters which are synthesized by alcohol acyltransferases (AATs) and degraded by carboxylesterases (CXEs). Although CXE activity can have a marked influence on volatile contents in ripe strawberry fruits, CXE function and regulation in them are poorly known. Here, we report the biochemical and functional characterization of the fruit receptacle-specific and ripening-related carboxylesterase FanCXE1. The expression of the corresponding gene was found to be antagonistically regulated by auxins and abscisic acid, key hormones that regulate fruit growth and ripening in strawberry. In vitro, FanCXE1 was able to hydrolyze artificial ester substrates similar to those produced by ripe strawberry fruits. Transient suppression of the FanCXE1 gene by RNAi resulted in an increase of important volatile esters such as methyl hexanoate, methyl butanoate and ethyl hexanoate as well as a decrease of the alcohols hexenol and linanool. The results of this work enhance our understanding of the molecular basis for volatile syntheses and facilitate production of better flavored strawberry fruits by introduction of the relevant alleles into common cultivars.

Identification of allergenomic signatures in allergic and well-tolerated apple genotypes using LC-MS/MS.

The apple fruit (Malus domestica L. Borkh) is one of the most popular fruits worldwide. Beyond their beneficial properties, apples contain proteins that trigger allergic reactions in susceptible consumers. Mal d1 to d4 are allergens present in a variety of different isoforms in apples. In this study, we used proteomics to quantify all four Mal d proteins in 52 apple genotypes with varying allergenic potentials. A total of 195, 17, 14, and 18 peptides were found to be related to Mal d1, d2, d3, and d4 proteins, respectively of which 25 different Mal d proteins could be unambiguously identified. The allergenic potential of the Mal d isoforms was characterized by comparing the isoform abundance with the allergenic score of genotypes from oral challenge tests. The detected Mal d peptides presumably have different IgE binding properties and could be used as potential molecular markers to discriminate between hypoallergenic and hyperallergenic cultivars.

High Resolution Quantitative Trait Locus Mapping and Whole Genome Sequencing Enable the Design of an Anthocyanidin Reductase-Specific Homoeo-Allelic Marker for Fruit Colour Improvement in Octoploid Strawberry (Fragaria × ananassa).

Fruit colour is central to the sensorial and nutritional quality of strawberry fruit and is therefore a major target in breeding programmes of the octoploid cultivated strawberry (Fragaria × ananassa). The red colour of the fruit is caused by the accumulation of anthocyanins, which are water-soluble flavonoids. To facilitate molecular breeding, here we have mapped with high resolution fruit colour quantitative trait loci (QTLs) (COLOUR, scored visually as in selection programmes) and associated flavonoid metabolic QTLs (5 anthocyanins compounds together with 8 flavonols and flavan-3-ols) to specific subgenomes of cultivated strawberry. Two main colour-related QTLs were located on the LG3A linkage group (Fragaria vesca subgenome). Genetic mapping, transcriptome analysis and whole genome sequencing enabled the detection of a homoeo-allelic variant of ANTHOCYANIDIN REDUCTASE (ANR) underlying the major male M3A COLOUR and pelargonidin-3-glucoside (PgGs) QTLs (up to ∼20% of explained variance). Consistent with previously published functional studies, ANR transcript abundance was inversely related with PgGs content in contrasted progeny individuals. Genetic segregation analyses further indicated that a molecular marker designed using an 18 bp deletion found in the 5'UTR of the candidate ANR homoeo-allelic variant is effective in identifying genotypes with intense red fruit colour. Our study provides insights into the genetic and molecular control of colour-related traits in strawberry and further defines a genetic marker for marker-assisted selection of new strawberry varieties with improved colour. The QTLs detected and the underlying candidate genes are different from those described to date, emphasising the importance of screening a wide diversity of genetic resources in strawberry.