==== Front Plant Physiol Plant Physiol plphys Plant Physiology 0032-0889 1532-2548 Oxford University Press US 37129240 10.1093/plphys/kiad250 kiad250 Update AcademicSubjects/SCI01270 AcademicSubjects/SCI01280 AcademicSubjects/SCI02286 AcademicSubjects/SCI02287 AcademicSubjects/SCI02288 Plphys/50 Vaccinium as a comparative system for understanding of complex flavonoid accumulation profiles and regulation in fruit https://orcid.org/0000-0002-8579-529X Albert Nick W New Cultivar Innovation, The New Zealand Institute for Plant and Food Research Limited (PFR), Palmerston North, 4472, New Zealand https://orcid.org/0000-0002-0032-8247 Iorizzo Massimo Plants for Human Health Institute, North Carolina State University, Kannapolis, NC 28081, USA Department of Horticultural Science, North Carolina State University, Raleigh, NC 27607, USA https://orcid.org/0009-0009-4979-1652 Mengist Molla F Plants for Human Health Institute, North Carolina State University, Kannapolis, NC 28081, USA Department of Horticultural Science, North Carolina State University, Raleigh, NC 27607, USA https://orcid.org/0000-0001-9224-056X Montanari Sara New Cultivar Innovation, PFR, Motueka, 7198, New Zealand https://orcid.org/0000-0001-8840-3565 Zalapa Juan USDA-ARS, Vegetable Crops Research Unit, Department of Horticulture, University of Wisconsin-Madison, Madison, WI 53706, USA https://orcid.org/0000-0002-3187-1461 Maule Andrew Plants for Human Health Institute, North Carolina State University, Kannapolis, NC 28081, USA Department of Horticultural Science, North Carolina State University, Raleigh, NC 27607, USA https://orcid.org/0000-0001-6836-3041 Edger Patrick P Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA MSU AgBioResearch, Michigan State University, East Lansing, MI 48824, USA https://orcid.org/0000-0002-0974-364X Yocca Alan E Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA HudsonAlpha Institute for Biotechnology, 601 Genome Way, Huntsville, 35806, USA https://orcid.org/0000-0001-9238-9647 Platts Adrian E Department of Horticulture, Michigan State University, East Lansing, MI 48824, USA HudsonAlpha Institute for Biotechnology, 601 Genome Way, Huntsville, 35806, USA https://orcid.org/0000-0002-3321-7471 Pucker Boas Institute of Plant Biology & BRICS, TU Braunschweig, Braunschweig, 38106, Germany https://orcid.org/0000-0002-1732-3688 Espley Richard V PFR, Auckland, 1142, New Zealand Author for correspondence: nick.albert@plantandfood.co.nz (N.W.A.), richard.espley@plantandfood.co.nz (R.V.E.) The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (https://academic.oup.com/plphys/pages/General-Instructions) are: Richard V. Espley (Richard.espley@plantandfood.co.nz), Nick W. Albert (nick.albert@plantandfood.co.nz). Conflict of interest statement. None declared. 7 2023 02 5 2023 02 5 2023 192 3 16961710 04 4 2023 06 4 2023 11 5 2023 © The Author(s) 2023. Published by Oxford University Press on behalf of American Society of Plant Biologists. 2023 https://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract The genus Vaccinium L. (Ericaceae) contains premium berryfruit crops, including blueberry, cranberry, bilberry, and lingonberry. Consumption of Vaccinium berries is strongly associated with various potential health benefits, many of which are attributed to the relatively high concentrations of flavonoids, including the anthocyanins that provide the attractive red and blue berry colors. Because these phytochemicals are increasingly appealing to consumers, they have become a crop breeding target. There has been substantial recent progress in Vaccinium genomics and genetics together with new functional data on the transcriptional regulation of flavonoids. This is helping to unravel the developmental control of flavonoids and identify genetic regions and genes that can be selected for to further improve Vaccinium crops and advance our understanding of flavonoid regulation and biosynthesis across a broader range of fruit crops. In this update we consider the recent progress in understanding flavonoid regulation in fruit crops, using Vaccinium as an example and highlighting the significant gains in both genomic tools and functional analysis. ==== Body pmcAdvances New Vaccinium genomic resources link genetic and physical maps Recent publication of blueberry, cranberry, and bilberry genome sequences Shared community resource hub at Genome Database for Vaccinium (www.vaccinium.org) New QTLs for anthocyanin content identified in blueberry and cranberry A network of MYB transcription factors has been identified to co-ordinately regulate flavonoid accumulation to help more broadly understand flavonoid control Increased data linking flavonoid pathway gene transcription with flavonoid chemistry, applicable to many crops Introduction Flavonoid biosynthesis is one of the best characterized specialized metabolite pathways in plants, owing to a long history of research on mutants that are phenotypically conspicuous. Examples include petunias (Petunia hybrida) and snapdragons (Antirrhinum majus) with altered flower colors, maize (Zea mays) with altered cob colors, and Arabidopsis with pale seed coats (transparent testa/tt mutants). Studies from these and other systems have helped build an understanding of the flavonoid biosynthetic pathway (Fig. 1A; (reviewed in Tohge et al. 2017; Albert et al. 2022). However, we are still learning about the accumulation of flavonoids; as they increasingly become a breeding target for nutritional enhancement in crops, so the need for genetic markers becomes more important. Berry crops, such as blueberry, often contain high concentrations of diverse flavonoids, making them useful to study the biosynthesis and control of flavonoids. The recent generation of Vaccinium genome assemblies, transcriptome assemblies, and metabolomes opened the opportunity to study flavonoid production in a range of fruit crops. Figure 1. Biosynthesis, inheritance, and regulation of flavonoid traits in Vaccinium.A) The biosynthetic pathway producing the major flavonoid classes reported in Vaccinium species: flavonols, PAs, and anthocyanins. Dihydroflavonols can vary in hydroxylation on the B-ring: dihydrokaemferol (DHK: 4ʹ), dihydroquercetin (DHQ: 3ʹ,4ʹ), and dihydromyricetin (DHM: 3ʹ,4ʹ,5ʹ). The hydroxylation pattern is retained by subsequent compounds in the pathway, corresponding to the 3 major anthocyanidins pelargonidin (Pel), cyanidin (Cy), and delphinidin (Del), and methylated derivatives peonidin (Peo), petunidin (Pet), and malvidin (Mal). LDOX and ANS can refer to the same gene, although some species, including Vaccinium, possess distinct LDOX and ANS genes. B) Anchored QTLs for flavonoids and color detected in blueberry and cranberry. Genomic QTL locations for anthocyanin acylation and glycosylation mapped in blueberry chr2 and chr4, respectively, across 2 mapping populations. QTL regions in blue represent those mapped in NxHP (Montanari et al. 2022), and QTLs in red are those mapped in DSxJ (Mengist et al. 2022b). Markers spanning 2-LOD intervals were used to anchor these QTLs. The sections in a grey shaded box indicate the genomic regions in chr2 (chr2: 8,005,747–11,180,545 bp) and chr4 (chr4: 55,208,167–60,978,559) where the regions spanning 2-LOD intervals for the same trait (e.g. acylated anthocyanin) overlap across the 2 mapping populations/studies. Genomic locations of QTLs for TAcy, PAC, color variation, and color intensity mapped in cranberry across 3 mapping populations and 2 genetic studies (Diaz-Garcia et al. 2018, 2021). QTL regions in blue represent those mapped in CNJ04–02, QTLs in red represent QTLs mapped in CNJ02–01, and QTLs in green represent QTLs mapped in GRYG. Markers spanning the 1.5-LOD interval were used to anchor these QTLs. The grey box sections indicate the genomic regions in chr3 (chr3: 8,700,846–12,176,781 bp, and chr3:16,922,077–20,907,119), where the region spanning the 1.5-LOD intervals for the same or related traits (e.g. TAcy and color intensity) overlap across the populations/studies. C) Flavonoid profiles change during bilberry, blueberry, and cranberry fruit development. The shift from PAs during early, immature stages are linked with SG5/MYBPA2 genes, and anthocyanin accumulation at ripening is linked with SG6/MYBA genes. MYBPA1 has a biphasic expression profile, correlating with both PAs and anthocyanins. D) Anthocyanins and PAs are regulated by MBW TF complexes, which activate flavonoid biosynthetic gene expression. The target genes regulated by MBW complexes are primarily determined by the type of R2R3-MYB. The exact stoichiometry is unknown, but MBW complexes can contain 2 MYBs via dimerized bHLH proteins. E) MBW complexes containing MYBA and MYBPA2 proteins directly regulate MYBPA1 genes, which may then be reinforced by MYBPA1 itself. The genus Vaccinium L. (Ericaceae) comprises approximately 450 diverse species, including the following major commercial crops: northern highbush blueberry (V. corymbosum L.), rabbiteye blueberry (V. virgatum, formerly V. ashei), lowbush blueberry (V. angustifolium), bilberry (V. myrtillus L.), cranberry (V. macrocarpon), and lingonberry (V. vitis-idaea L.) (Luby et al. 1991; Retamales and Hancock 2018). The berries of many wild Vaccinium species were traditionally used as foods and medicines; however, they remained undomesticated until the 20th century, when selection and breeding efforts began (Ehlenfeldt 2008; Vorsa and Zalapa 2019). There is now an increasingly sophisticated approach to Vaccinium crop improvement (Edger et al. 2022). Vaccinium ploidy spans diploid (2×), tetraploid (4×), and hexaploid (6×) species, although their ancestry and taxonomy are not yet fully resolved (Lyrene et al. 2003). Recent progress developing genetic resources for Vaccinium includes genome assemblies for blueberry, bilberry, and cranberry (Colle et al. 2019; Diaz-Garcia et al. 2021; Wu et al. 2021; Mengist et al. 2022a ) and a centralized genomics hub for genomic and transcriptomic sequences: Genome Database for Vaccinium (www.vaccinium.org). Global production and consumption of blueberries and other Vaccinium berries is growing, benefiting from an increasing body of research demonstrating their potential health promotion (Norberto et al. 2013). This translates to consumer behaviors, with health benefits being a key driver for blueberry-purchasing decisions (Gilbert et al. 2014). The health benefits of Vaccinium berries are attributed to their phytochemical content, including various flavonoid and other phenolic compounds. For instance, Vaccinium berry consumption has been associated with numerous health benefits: alleviation of inflammation, fever, digestive and urinary tract disorders; acting as an antiseptic (Tundis et al. 2021); association with a reduced risk of cardiovascular disease, some cancers, and type 2 diabetes; and benefits to neuroprotection, cognition, and vision (Ma et al. 2018; Kalt et al. 2020). Although much of this research has focused on blueberry and cranberry (Nemzer et al. 2022), there is also increasing evidence for the health benefits of other Vaccinium species, including bilberry (Sharma and Lee 2022) and lingonberry (Kowalska 2021). Flavonoids have diverse physiological roles in planta (reviewed in Albert et al. 2022). Flavonols provide protection against UV-B radiation and abiotic and biotic stresses and have different roles according to tissue type and developmental stage. In vegetative tissues, anthocyanins help to screen excess light in stressed tissues and possibly act as antioxidants, protecting cells against reactive oxygen species–induced oxidative damage. Anthocyanins perform important roles for animal signaling, providing visual cues to attract pollinators and coloring fruit upon ripening to attract seed distributers. Another class of flavonoids, proanthocyanidins (PAs), are astringent compounds that deter herbivory. They are abundant in Vaccinium vegetative tissues (Hokkanen et al. 2009) but also important in fruit by deterring consumption before fruit (and seed) maturation. Vaccinium berries contain high concentrations of anthocyanins, which are generally restricted to the fruit skin, but in some species, like bilberry, they can also be present in the flesh (Riihinen et al. 2008). The anthocyanins in Vaccinium include various glycosides of cyanidin, delphinidin, malvidin, peonidin, and petunidin, and the concentrations and moiety decorations can vary within and between species (Kalt et al. 1999; Veberic et al. 2015). Flavonoid production appears spatially orchestrated within the berry as it develops, with the specialized metabolism of flavonoids following a pattern as the fruit matures (Dare et al. 2022). This suggests synchronized transcriptional control of the pathway. Although there is good evidence that anthocyanin production is developmentally coordinated, the type and concentration can vary considerably with geographic location and environmental conditions (Zoratti et al. 2015; Karppinen et al. 2016). Vaccinium berries also contain an impressive variety of other flavonoids and other notable phenylpropanoids. In blueberry, these are largely hydroxycinnamic acids, stilbenes, PAs, and flavonols, all being at their highest abundance early in fruit development and decreasing in concentration as anthocyanin production rises (Dare et al. 2022). It is quite possible that a variety of phytochemicals contributes to the overall bioactivity of Vaccinium berries. For example, Carey et al. (2013) showed that blueberry extract had greater efficacy than pure anthocyanin or stilbenes in protecting against stress-induced inflammation. They suggested synergistic effects between anthocyanins, stilbenes, and other secondary metabolites present in blueberries may be responsible for the increased efficacy of extracts compared with pure compounds. Indeed, synergistic effects of anthocyanin and stilbenes have been demonstrated for an inflammatory bowel disease model, using transgenic tomato extracts producing anthocyanins, flavonols, and stilbenes, or all of these, showing significant improvement in reversing disease. This approach allows the contribution of different classes of phenolic metabolites to be assessed in different combinations in a well-controlled manner. Similar effects were also observed comparing white and red grape skins, which contain different quantities of anthocyanins and stilbenes (Scarano et al. 2017). This review distils recent progress in understanding flavonoid regulation using Vaccinium as an example of where recent significant gains in both genomic tools and functional analysis have been made. There is strong momentum within the Vaccinium research community to further progress to drive increasingly sophisticated breeding approaches but also to develop knowledge that will helpus understand flavonoid regulation in other fruit crops. Flavonoid biosynthesis and chemical diversity in Vaccinium Most flavonoids are 15-carbon structures comprised of 3 rings called the A-, B-, and C-rings. The first committed step toward flavonoid biosynthesis is catalyzed by chalcone synthase (CHS), often encoded by multigene families of differentially regulated genes. Although conversion of naringenin chalcone into naringenin occurs spontaneously in vitro, this reaction is catalyzed by chalcone isomerase (CHI) in planta. The flavonoid backbone is then hydroxylated by flavanone 3-hydroxylase (F3H) and may be hydroxylated further by flavonoid 3ʹhydroxylase (F3ʹH) and flavonoid 3ʹ5ʹhydroxylase (F3ʹ5ʹH). The 3 resulting dihydroflavonols, which differ in the number of hydroxyls on the B-ring, progress into the pathway toward branches for flavonols or anthocyanins/PAs. The metabolic flux into the 2 competing branches depends on substrate preferences of dihydroflavonol 4-reductase (DFR) and flavonol synthase (FLS). In bilberry, flavonols accumulate from early stages of fruit development, primarily in the skin (Dare et al. 2022). A similar pattern was observed in strawberry (Fragaria vesca), where flavonols dominate during early developmental stages of fruits, whereas anthocyanin accumulate at late stages (Miosic et al. 2014). In blueberry and bilberry fruit, these are primarily quercetin glycosides with traces of myricetin glycosides (or methylated derivatives), although kaempferol-based flavonols are detected in leaves (Riihinen et al. 2008). FLS genes were assigned to 3 lineages in the Vaccinium species. One lineage shows perfect conservation of all amino acids that have been reported to be functionally important (Pucker et al. 2020). The other 2 lineages differ at position S225 and E295 from the Arabidopsis FLS1, respectively. It is possible that one of these substitutions altered the substrate specificity of FLS. Spatio-temporal differences in gene expression of these FLS genes could explain the flavonol pattern observed in leaves and fruits. Methylated derivatives of quercetin and myricetin are also produced, but this occurs later in fruit development when anthocyanin biosynthesis occurs, coinciding with flavonoid O-methyl transferase expression (Günther et al. 2020). Cranberry has a different flavonol profile, with significant concentrations of myricetin-glycosides in addition to quercetin-glycosides (Vvedenskaya et al. 2004; Singh et al. 2009). Flavonoid biosynthesis genes are generally conserved across plant species (Supplemental Table S1). Differences affect the flavone biosynthesis, which relies on the P450 cytochrome flavone synthase II (FNSII) in most plant species. However, no FNSII ortholog was discovered in Arabidopsis, kiwifruit (Actinidia chinensis), cranberry, or bilberry. Only species of the Apiaceae are known to produce flavones through the F3H-derived flavone synthase I. Another important difference between species is the hydroxylation pattern of anthocyanins. F3′5'H is required for the formation of delphinidins and the resulting blue pigmentation. The F3′5'H enzyme functionality evolved multiple times independently from the F3′H lineage. Likely functional F3′5'H candidates were discovered in grape (Vitis vinifera), kiwifruit, and the 3 Vaccinium species. This functionality is missing in Arabidopsis, apple (Malus domestica), and strawberry. Arabidopsis does not harbor leucoanthocyanidin reductase (LAR), which leads to lower diversity of PAs compared with other species. Orthologs of 3AT and 3MAT were not detected in most of the species, which suggests that these might be lineage-specific functions. A parallel or convergent evolution of this functionality remains possible. Copy number differences between species seem to be tandem duplications that form gene arrays. The anthocyanin color is primarily determined by the hydroxylation state of the B-ring: pelargonidin/orange-red (4ʹ), cyanidin/red-magenta (3ʹ4ʹ), and delphinidin/purple-blue (3ʹ4ʹ5ʹ). Pelargonidin-based anthocyanins are not reported in Vaccinium species, suggesting that, as observed in other species (Johnson et al. 2001), DFR has strong substrate preference for dihydroquercetin (3ʹ4ʹ) or dihydromyricetin (3ʹ4ʹ5ʹ) and not dihydrokaemferol (4ʹ). However, an inspection of the DFR sequences in the Vaccinium species revealed N at the critical position corresponding to 134 in the reference sequence analyzed by Johnson et al. 2001, which suggests that these enzymes should be able to accept all 3 dihydroflavonols. No DFR sequence with the dihydrokaemperol excluding D at this position was detected. It is possible that additional amino acid residues are restricting the accepted substrate or that a strong hydroxylase activity restricts the available substrate to dihydroquercetin and dihydromyrecitin. Blueberries and bilberries have complex anthocyanin profiles consisting of cyanidin- and delphinidin-based anthocyanins (and methylated forms: peonidin, malvidin, petunidin) compared with red-fruited cranberries, lingonberries, and some other wild Vaccinium species, which produce cyanidin/peonidin. The lack of delphinidin-based anthocyanins in these red-fruited species could be due to the absence of F3ʹ5ʹH activity. This is similar to apple and strawberry, which do not have F3′5′H activity due to a lack of the corresponding gene (Han et al. 2010). Differences regarding the F3′5′H activity have also been reported for grape cultivars and are also largely due to transcriptional regulation (Falginella et al. 2010). However, cranberry is reported to produce trihydroxylated flavonols (myricetin) in fruit (Vvedenskaya et al. 2004), and the cranberry genome has an intact F3ʹ5ʹH gene (vmacro08261-RA) as well as other apparently pseudogenized copies. Thus, the lack of delphinidin may be because the timing of F3ʹ5ʹH activity does not coincide with the expression of other anthocyanin biosynthesis genes. Multiple genes encode F3ʹ5ʹH in blueberries, but it is unclear if these are differentially regulated for producing different flavonoid classes (e.g. flavonols vs anthocyanins). Formation of anthocyanidins, common to anthocyanin and PA pathways, is catalyzed by anthocyanidin synthase (ANS) or leucoanthocyanidin dioxygenase (LDOX). In some species these refer to a single gene (e.g. Arabidopsis), but in others, including blueberry (Lafferty et al. 2022a), they can be phylogenetically separated into distinct clades, with ANS genes regulated for anthocyanin production and LDOX genes for PAs (Jun et al. 2018). Glycosylation of the anthocyanidin, e.g. by UDP-glucose:flavonoid-3-O-glucosyltransferase (UFGT), forms anthocyanin pigments. In blueberries and cranberry, the major flavonoid glycosylation patterns are typically galactosides and arabinosides, with lower concentrations of glucosides (Prior et al. 2001). In contrast, V. oxycoccus (small cranberry) contains a greater proportion of glucosides than galactosides (Vorsa and Polashock 2005), suggesting a great variety of anthocyanin decorations may exist within nondomesticated Vaccinium species. Additional modifications/decorations of the anthocyanin molecules can generate even greater diversity of anthocyanin structures through glycosylation in other positions (5GT) or acylation (acyl transferase). Anthocyanins with 3-glucose 6-acetyl modifications are found in some highbush blueberries cultivars (e.g. “Nui”), but in other Vaccinium taxa these are present at low concentrations (e.g. lowbush blueberry) or not detectable (rabbiteye blueberry, cranberry) (Prior et al. 2001; Günther et al. 2020). These anthocyanin modifications alter solubility, stability, and the ability of anthocyanins to form supra-molecular complexes with other flavonoids and phenolics (e.g. co-pigmentation with flavonols). They also alter bioaccessibility and bioavailability when consumed (He et al. 2005; Mengist et al. 2020). After synthesis at the endoplasmatic reticulum, anthocyanins are transported to the central vacuole for long-term storage. Two different models for anthocyanin transport have been proposed: (1) vesicle-mediated or (2) direct transport through the cytoplasm and across the tonoplast (reviewed by Pucker and Selmar 2022). The latter is thought to require the activity of a specific glutathione-S-transferase and additional transporters [e.g. including multi-drug and toxin extrusion transporter (Zhao et al. 2011) and ATP-binding cassette subfamily C transporter (Goodman et al. 2004)]. The transport of anthocyanins is poorly understood and requires additional research. For PA biosynthesis, leucoanthocyanidins can be converted into (+)-catechin by LAR, which can either go on to polymerize into PAs or, as demonstrated by Jun et al. (2018), be converted into anthocyanidin (e.g. Cy) by LDOX. Anthocyanidins are then reduced by anthocyanidin reductase (ANR) into (−)-epicatechin. Blueberry and other Vaccinium species have both ANS and LDOX genes, which are differentially expressed during anthocyanin or PA biosynthesis, respectively (Lafferty et al. 2022b). This not only provides redundant activity, but it may also contribute to metabolite complexity, as proposed for Medicago truncatula by Jun et al. (2018). The PA monomers can react to form oligomers and polymers, forming PAs, although this aspect is not fully understood. Transport to the vacuole requires specific proton pumps (e.g. AtTT13) and MATE transporters (Zhao and Dixon 2009; Appelhagen et al. 2015). Vaccinium fruit tend to have complex flavonoid profiles. This includes mixtures of methylated flavonoids and varying types of glycosylation and acylation. Some of the genes responsible for modifying and decorating the flavonoids are not yet identified. Flavonoid glycosyl transferases and other decorating enzymes belong to very large gene families, and their functionality is difficult to infer from sequences (Yonekura-Sakakibara and Saito 2014; Bontpart et al. 2018). Flavonoids are not the only bioactive compounds in Vaccinium fruit, and it is likely synergistic actions with other phenolics [e.g. chlorogenic acid (CHA) and stilbenes] occur to confer health benefits when consumed. Quantitative genetics can make a significant contribution to developing novel crops with improved health benefits by identifying loci associated with particular flavonoid species and phenolic compounds to assist breeding efforts. Flavonoid genetics in blueberry and cranberry With the recent development/adoption of advanced genotyping platforms and genomic resources in Vaccinium, improvements in understanding the genetic mechanisms controlling flavonoid accumulation have been made for blueberry and cranberry. In this review, we summarize data from available flavonoid studies in these 2 important crops, emphasizing quantitative trait loci (QTLs) that were stable across years and/or detected across multiple genetic backgrounds, and we anchor them to the physical map. Regions spanning these QTLs are the most likely targets for studies on flavonoid genetics in these crops. Anthocyanin QTLs in blueberry For blueberries, 2 studies focused on dissecting the genetic mechanism controlling CHA content and anthocyanin content and composition (Mengist et al. 2022b; Montanari et al. 2022). Mengist et al. (2022b) used a mapping population (DSxJ), representing highbush blueberry cultivars “Draper” and “Jewel,” and identified 180 QTLs for total and individual anthocyanin content, relative anthocyanin composition, and CHA (Supplemental Table S2). The study highlighted identification of QTLs for CHA, total anthocyanin content, and the conjugations of anthocyanin with the different sugar moieties and acylation. These QTLs were stable across the years and explained a large fraction (up to 80%) of the phenotypic variation. QTLs for CHA and anthocyanin acylation overlapped in the same region of chr2, and haplotype analysis indicated that these QTLs were located in different haplotypes, which implies they could be independently selected. Montanari et al. (2022) used a different mapping population (NxHB) representing 2 different highbush cultivars, “Nui” and “Hortblue Petite.” Consistent with the results of Mengist et al. (2022b), they mapped stable and strong effect QTLs on chr2 and chr4 that regulated the concentration of acylated anthocyanins and glucoside-based anthocyanins, respectively. In this review, we anchored the 2 clusters of QTLs detected on chrs 2 and 4 in the 2 studies to the W85 v2 genome (Mengist et al. 2022b), showing that they overlap (Fig. 1B). The QTLs controlling acylation on chr 2 span a common 3.2-Mb-long region between positions 8.0 and 11.2 Mb, and the QTLs linked to glycosylation mapped to chr4 overlap over a 4.8-2Mb-long region between positions 56.2 and 61.0 Mb (Fig. 1B). Co-mapping of these QTLs across different genetic backgrounds represents the first level of QTL validation. Initial efforts to mine these regions for candidate genes by Mengist et al. (2022b) and Montanari et al. (2022) identified a putative UDP-glucosyl transferase gene (Vcev1_p0.Chr04.11988), 2 putative DFR genes (Vcev1_p0.Chr4.11916, Vcev1_p0.Chr4.11918), and a MYB transcription factor (TF) (Vcev1_p0.Chr4.11990) as candidate underpinning the chr4 QTL. Two BADH acyltransferases genes (Vcev1_p0.Chr02.03371 and Vcev1_p0.Chr02.03383), a MYB TF (Vcev1_p0.Chr2.03364), and a DFR gene (Vcev1_p0.Chr2.03323), were identified as candidates underlying the chr2 QTL. Candidacy of the 2 BADH-acyltransferase genes in chr2 and the UDP-glucosyl transferase in chr4 was also supported by gene expression analysis by Mengist et al. (2022b). Anthocyanin and proanthocyanin QTLs in cranberry To date, only one study has examined the genetics of flavonoid content in cranberry (Diaz-Garcia et al. 2018). This focussed on total anthocyanin content (TAcy) and PA content (PAC), evaluated using spectrophotometric methods, with fruit color intensity and variation measured by digital image analysis. QTL analysis was performed using an integrated map representing 3 mapping populations (CNJ02–01, CNJ04–02, and GRYG) that enabled direct comparison of the locations of the QTLs detected across the 3 families (Supplemental Table S3). Several QTLs for TAcy and/or PAC overlapped with QTLs for color variation and/or intensity, which indicated that these traits are correlated (Diaz-Garcia et al. 2018). Seven regions in chrs 1, 3, 6, 9, and 12 harbored QTLs associated with the same/correlated traits (e.g. TAcy and color intensity) across 2 or 3 populations, suggesting a stable effect across genetic backgrounds. However, only 2 of these 7 regions contained QTLs that explained a considerable proportion (>5%) of phenotypic variation. In particular, one region on chr3 of the GRYG and CNJ02–01 populations harbored overlapping QTLs for color intensity, color variation, and PAC that explained up to 19% of phenotypic variation. The second hotspot was also identified on chr3 in the CNJ02–01 and CNJ04–02 populations, and it included QTLs for all the traits (PAC, TAcy, color variation, color intensity) that explained up to 39.5% of the phenotypic variation. In 2021, Diaz-Garcia et al. (2021) reanalyzed the same TAcy data using an improved linkage map of population CNJ02–01 and confirmed a strong effect of QTL in the same chr3 region that explained up to 51.2% of phenotypic variance. Here, anchoring of the overlapping QTLs on chr3 to the “Stevens’ genome identified 2 regions spanning 3.5 and 4 Mb. Within the first region, located between 8.7 and 12.2 Mb on chr 3, they identified a cluster of 3 R2R3 MYB110 TF genes (vmacro18045:10437125-10439365, vmacro18044:10493992-10499550, and vmacro18043:10561041-10564609) as candidate for the control of TAcy in cranberry. Two of these genes, vmacro18044 and vmacro18045, exhibited extremely high similarity with blueberry MYBA TF, which is known to directly activate anthocyanin biosynthesis (Plunkett et al. 2018). A third unannotated MYBA1 gene was identified in cranberry within this locus (Wu et al. 2021). In the second region, which spanned position 16.9–20.9 Mb on chr 3, a DFR gene was annotated as a potential candidate underlying this QTL. Additionally, an analysis of the anthocyanin content in a diverse cranberry germplasm suggested that genetic diversity for a range of anthocyanin traits exists. This includes variation in the content of cyanidin and its methylated derivative, peonidin, and variation in glycosylation (galactosides/arabinosides/glucosides) (Vorsa et al. 2003). QTLs associated with these traits have not yet been reported for cranberry. The studies reported here for blueberry and cranberry add to a list of works on the genetic control of flavonoid content in fruit crops, mostly in species belonging to the family Rosaceae, such as strawberry (Labadie et al. 2020, Pott et al. 2020), peach (Prunus persica; Abdelghafar et al. 2020), Asian plum (Prunus salicina; Valderrama-Soto et al. 2021; Battistoni et al. 2022), sweet cherry (Prunus avium; Calle et al. 2021), and apple (Kumar et al. 2022). Such studies probably have been boosted by the recent improvements in flavonoid identification and quantification methodologies and the increased interest around their beneficial health properties. All these studies reported on the quantitative nature of flavonoid content and highlighted how this trait is affected by the environment. Nevertheless, major QTLs and chromosome hotspots were detected in all species and were stable across seasons/environments, suggesting that a small number of genes have a strong effect on the accumulation of flavonoids. This is important because it opens the door to the application of marker-assisted selection for flavonoid content in fruit crops. Of all the different classes of flavonoids, anthocyanins have received the greatest attention in fruit crops. However, most studies investigated the genetic control of anthocyanin content but not their composition. On the contrary, both Mengist et al. (2022b) and Montanari et al. (2022) showed that most genetic variation for blueberry anthocyanin is associated with glycosylation and acylation of the molecule, and genetic contribution for total anthocyanin is limited relative to their composition. An analogous study performed by Sun et al. (2020) in grape found QTLs linked to both the methylation and acylation of the anthocyanin molecules in the skin of grape berries. Analyzing the specific composition of anthocyanins in other fruit crops, including cranberry, is therefore important to fully understand the genetic regulation of these compounds. In all these fruit crops, including blueberry and cranberry, genes encoding MYB TFs were identified in the major QTL regions, confirming their well-known role in anthocyanin regulation. However, other candidate genes that appeared of particular interest in blueberry are the BADH-acyltransferases and the UDP-glucosyl transferases, which have only been reported in peach (Abdelghafar et al. 2020) and Citrus (Mou et al. 2021). It is, however, important to note that genes involved in anthocyanin acylation and glycosylation belong to very large gene families, for which the function is not well conserved across species, making them difficult to identify. Continuing to leverage genetic studies is likely to be the most effective way to identify these gene. Because work on flavonoid genetics in Vaccinium is still in its relative infancy, the next obvious step would be to expand association analyses to flavonoid classes other than anthocyanins. Looking at the variation of a large number of flavonoid compounds and their biochemical relationships can provide insights into their regulatory gene networks and contribute to a better understanding of their biosynthetic pathways, as was shown in strawberry by Pott et al. (2020). Flavonoid regulation The regulation of anthocyanin and proanthocyanin biosynthesis is well characterized. The flavonoid biosynthesis genes are transcriptionally regulated, and flavonoids accumulate in the tissues where coordinated expression of the biosynthesis and transport genes for a particular metabolite occurs. This is mediated by TFs, particularly those of the MYB family. The R2R3-MYB genes can be phylogenetically clustered into subgroups (SG), which generally have conserved functions among flowering plants; e.g. those belonging to SG5 regulate PAs, SG6 regulates anthocyanins, and SG7 regulates flavonols or flavones (Stracke et al. 2001; Albert and Allan 2021). These insights have largely come from genetic models such as maize, petunia, Antirrhinum, and Arabidopsis (reviewed in Albert et al. 2022), but it is becoming clear flavonoid regulation is more nuanced, particularly in crops with complex flavonoid profiles that change throughout different developmental stages (Fig. 1C). The MYB-bHLH-WDR complex In addition to pathway-specific MYBs, basic helix-loop-helix (bHLH) and WD-repeat (WDR) TFs are essential for regulating anthocyanins and PAs. These TFs can form MYB-bHLH-WDR (MBW) complexes (Fig. 1D), which bind to the promoters of target biosynthesis genes, activating transcription and ultimately resulting in metabolite accumulation (Baudry et al. 2004; Albert et al. 2014; Xu et al. 2015). The bHLHs involved belong to SG IIIf but can be further separated into bHLH-1 and bHLH-2 subtypes, corresponding to AtEGL3/PhJAF13/AmDel and AtTT8/PhAN1/AmIncI, respectively (Feller et al. 2011; Albert et al. 2021). These proteins share some activity but also have nonredundant functions, varying between species; for example, TT8 (bHLH-2) is essential for PA regulation in Arabidopsis seed coats but has redundant activity for anthocyanin regulation with the bHLH-1 proteins (EGL3 and GL3) (Nesi et al. 2000; Baudry et al. 2006; Gonzalez et al. 2008). In petunia and other Solanaceae species, AN1 is essential for anthocyanin regulation, vacuolar acidification, and seed coat formation (presumably PAs) despite expressing the functional bHLH1 gene, JAF13 (Montefiori et al. 2015; Spelt et al. 2000, 2002). Hierarchical regulation of bHLH2 genes by MBW complexes is now well established in anthocyanin and PA regulation and may be initiated by MBW complexes containing bHLH-1 proteins (Baudry et al. 2006; Albert et al. 2014); it is conserved in Vaccinium (Lafferty et al. 2022b). This provides a feed-forward mechanism to enhance responses by ensuring sufficient abundance of key proteins and can also expand the expression domain of bHLH2 genes into additional tissues (Albert et al. 2021). Other conserved aspects of MBW activity, such as repression by SG4 MYB and R3-MYB proteins (Albert et al. 2014), have been identified and characterized in Vaccinium species for anthocyanin and PA regulation (Lafferty et al. 2022b). Given the shared roles for WDR and bHLH components within MBW complexes, it is primarily the R2R3-MYB (activator) that determines the specificity of the genes and pathways regulated. Flavonol regulation Flavonols help protect plants against abiotic stress, acting as UV-B sunscreens and scavengers of reactive oxygen species. Consequently, their biosynthesis is typically promoted by light (including UVA and UVB) in vegetative tissues, immature flowers, or developing fruit, particularly in exposed skin tissues. Their production is regulated by SG7 MYBs, such as AtMYB12 in Arabidopsis and VvMYBF1 in grape (Mehrtens et al. 2005; Czemmel et al. 2009), but also by bZIP TFs like HY5 (Hartmann et al. 2005), functioning independently of MBW complexes. These directly activate the promoters of genes necessary for flavonol production (CHS, CHI, F3H, and FLS). SG7 MYBs have been identified and characterized in numerous fruit crops, including apple (MsMYB22;Wang et al. 2017), pear (Pyrus bretschneideri)PbMYB112b;Zhai et al. 2019), strawberry (Martínez-Rivas et al. 2023), and recently VmMYBF in bilberry (Karppinen et al. 2021). Genetic variants of these regulators can confer quantitative variation in flavonol content, such as pink-fruited tomato (yellow locus: Slmyb12; Ballester et al. 2010), pepper (Capsicum annum: CaMYB12-like; Wu et al. 2023), and pear (Zhai et al. 2019). There has been recent progress in understanding flavonol regulation, but, given the wide variation in flavonol concentration and the influence of environment, it is likely that there are other regulatory pathways and genes yet to be discovered. Anthocyanin regulation The SG6 MYBs are central for regulating anthocyanins, forming transcription factor complexes with bHLH and WDR proteins (MBW) that activate genes for anthocyanin biosynthesis and transport. These genes often exist in small gene families, with distinct expression profiles, and it is common to have genes linked with vegetative pigmentation (stress/light responsive) and those with more specific roles in flower pigmentation or fruit color (Schwinn et al. 2006; Albert et al. 2011; Matus et al. 2017). Although initially characterized in model plant species, significant insights have also come from fruit crops. This is because variants of SG6 MYB genes can give rise to fruit with different colors or patterns, which have been selected for or against. Thus, SG6 MYB alleles have been identified in numerous crops that are responsible for determining anthocyanin content or patterning (Paauw et al. 2019; reviewed in Albert and Allan 2021). Despite the well-established role of SG6 MYBs in controlling anthocyanins, these remained elusive for some time in Vaccinium species compared with other fruit species. Early studies identified MYBPA1 genes as the candidate primary regulator rather than a typical SG6 MYB, based on expression patterns linked with pigmentation. For example, silencing a bilberry homologue of the MADS-box ripening regulator FRUITFUL resulted in a loss of MYBPA1 expression and anthocyanins (Jaakola et al. 2010). Similarly, albino mutants of bilberry and bog bilberry (V. uligonosum) showed loss of MYBPA1 expression during berry ripening (Primetta et al. 2015; Zorenc et al. 2017). The MYBPA1-type proteins belong to a distinct subgroup of R2R3-MYB originally identified in grape (VvMYBPA1) and absent from Arabidopsis. VvMYBPA1 was initially considered a PA regulator because it could activate the expression of the PA-biosynthesis genes ANR and LAR required to make catechin and epicatechin monomers, and it could complement the Arabidopsis tt2 mutant (Bogs et al. 2007). VvMYBPA1 expression peaks early during fruit development when PA biosynthesis occurs and again at véraison in the skin, where anthocyanin accumulation occurs (Bogs et al. 2007). Yet VvMYBPA1 was unable to activate UFGT (Bogs et al. 2007), a key step for anthocyanin biosynthesis in grape that is regulated by the VvMYBA proteins (SG6) (Walker et al. 2007; Matus et al. 2017). Similar temporal patterns of metabolite accumulation are shared between grape and blueberry/bilberry (Czemmel et al. 2012; Primetta et al. 2015; Karppinen et al. 2016) (Fig. 1B), suggesting that both MYBPA1 and SG6/MYBA genes may be necessary for activating anthocyanin production in these species, in contrast to simpler model species. The SG6 MYB genes were identified in blueberry by taking a targeted approach. VcMYBA1 was isolated from blueberry fruit tissues and shown to be highly expressed in the fruit skin during ripening, when anthocyanin biosynthesis occurs (Plunkett et al. 2018). Furthermore, it was able to complement a defined SG6 MYB mutant (Antirrhinum, roseadorsea) and induce anthocyanin accumulation when infiltrated into Nicotiana benthamiana leaves. MYBA1 transgenes activated the promoters of several flavonoid genes, including UFGT, a critical biosynthesis step regulated specifically for anthocyanins in Vaccinium. Subsequently, SG6 MYBs from other Vaccinium species were identified and characterized. In bilberry, VmMYBA1 is associated with ripening-associated anthocyanin accumulation in fruit and is downregulated in 2 wild mutants that produce albino fruit (Karppinen et al. 2021). Furthermore, treating mature green bilberry or blueberry fruit with abscisic acid rapidly accelerated ripening and anthocyanin accumulation (Karppinen et al. 2018; Oh et al. 2018) and rapidly induced VmMYBA1 expression in bilberry (Karppinen et al. 2021). A MYBA1 orthologue from wufanshu (V. bracteatum) has been characterized with even stronger anthocyanin biosynthesis activation properties than blueberry VcMYBA (Zhang et al. 2021). Analysis of the recent genome assemblies for tetraploid highbush blueberry (Colle et al. 2019), cranberry (Diaz-Garcia et al. 2021), and bilberry found this locus is complex, containing numerous MYBA copies, pseudogenes, and repeat elements (Wu et al. 2021). The location of the MYBA locus is syntenic between cranberry, blueberry, and bilberry, but considerable variation in the locus structure exists between species, suggesting gene duplication events have occurred since species divergence. These gene copies can be separated into 2 clades, “MYBA1” and “MYBA2” (Wu et al. 2021), suggesting existence of 2 ancestral genes with different promoter sequences that determine pigmentation patterns. In species studied to date, the MYBA1 genes are expressed in fruit during ripening, while MYBA2 copies may contribute to regulating pigmentation in other tissues, for example, light-regulated pigmentation in vegetative or floral tissues (Karppinen et al. 2021). The identification of the MYBA genes in Vaccinium species raises questions of how anthocyanin regulation is balanced with the production of other types of flavonoids during fruit development and how the regulatory mechanisms might overlap. PA regulation In Vaccinium species PAs are a major class of flavonoids that accumulate throughout vegetative and reproductive tissues. Recent studies in bilberry and blueberry have identified several distinct MYB classes and gene families associated with PA regulation (Karppinen et al. 2021; Lafferty et al. 2022b). In bilberry this includes at least 5 SG5/MYBPA2 genes, 2 MYBPA1 genes, 2 MYB5 genes, and a PAR/MYBPA3 gene (Karppinen et al. 2021). MYB5 contributes to seedcoat formation and seed mucilage in Arabidopsis and M. truncatula (Gonzalez et al. 2009; Li et al. 2009; Liu et al. 2014) and vacuole acidification in petunia and citrus (Quattrocchio et al. 2006; Zhang et al. 2019). MYB5 expression can enhance PA accumulation in M. truncatula, petunia, and tobacco (Deluc et al. 2008; Cavallini et al. 2014) but seems to target a smaller subset of PA biosynthetic genes than SG5 MYBs (Liu et al. 2014). Liu et al (2014) showed that MBW complexes containing MtMYB14 (SG5) and MtMYB5 exist and act synergistically. The PAR/MYBPA3 gene also contributes to PA regulation in M. truncatula (Verdier et al. 2012), and the bilberry homologue induced production of the PA precursor flavan-3-ol (gallocatechin) monomers when expressed in tobacco (Karppinen et al. 2021). The number of classes of MYBs involved in PA regulation and the lack of mutants for these (natural or induced) in grape or Vaccinium makes it difficult to establish the contribution of each regulator type. There is clearly some redundancy between these proteins for some target genes (e.g. MYBPA1 and MYBPA2 both regulate DFR and ANR), but they are not functionally equivalent, suggesting additional target genes for transport or polymerisation. Why so many MYBs? Despite studies in grape, persimmon (Diospyros kaki), and poplar (Populus species) (Bogs et al. 2007; Akagi et al. 2010; James et al. 2017), the role of MYBPA1 proteins for flavonoid regulation has been unclear because their gene expression and target genes are associated with both PA and anthocyanin biosynthesis. However, recent evidence from Vaccinium helps to clarify their role as dual regulators, contributing to both pathways, reconciling findings from other fruit species. Silencing MYBPA1 in bilberry resulted in a visible loss of anthocyanins and reduced expression of CHS, F3ʹ5ʹH, ANS, and LAR1b, demonstrating that MYBPA1 makes an essential contribution to anthocyanin regulation (Karppinen et al. 2021). Blueberry VcMYBPA1.1 could not induce anthocyanin or PA accumulation in tobacco leaves alone, but a modified version fused to a strong transcriptional repressor (MYBPA1.1-CREST) was a potent inhibitor when coexpressed with VcMYBA1 (of anthocyanin) or VcMYBPA2.2/VcMYBPA2.3 (of PAs). Promoter activation assays again showed MYBPA1.1 had strong activation activity against the promoters of CHS, F3ʹ5ʹH, and ANR, moderate activity against the DFR and ANS, but none against UFGT or the MATE8 anthocyanin transporter (Lafferty et al. 2022a). These findings support a role for regulating key steps within the flavonoid pathway despite lacking the ability to regulate all the steps required for metabolite production and transport. How is the activity of multiple MYBs coordinated for metabolite production? The answer lies in hierarchy: some TFs activate others. Hierarchical regulation of MYB repressor and bHLH2 genes is well established for anthocyanin and PA regulation (Albert et al. 2014), but regulation of MYB activators by MBW complexes is a recent addition. It was recently shown that MYBPA1 genes are directly regulated by MBW complexes containing MYBPA2 or MYBA proteins (Lafferty et al. 2022b) (Fig. 1D). This key finding reconciles the expression profiles of MYBPA1 in grape and Vaccinium species, with peak expression corresponding with PA regulation in immature fruit and a second peak when anthocyanin occurs (Bogs et al. 2007) (Fig. 1C). It also explains why grape hairy root cultures overexpressing VvMYBPA2 have increased VvMYBPA1 expression (Terrier et al. 2009). In apple, overexpressing MdMYBPA1 enhances PA accumulation, and MdMYBPA2 and MdbHLH33 were shown to activate MdMYBPA1 promoter (Wang et al. 2018). Anthocyanin accumulation was also enhanced by MdMYBPA1 overexpression with low temperature when MdMYB10 (SG6 MYB) was induced (Wang et al. 2018), suggesting that cooperative regulation of core flavonoid genes by MYBPA1 proteins with pathway specific MYBs is common in fruit. SG5 MYBs may also contribute to anthocyanin regulation and modification. A recent study in strawberry characterized FaMYB123 (SG5) as a regulator of anthocyanin biosynthesis genes, including ANS, and genes for flavonoid decoration and acylation (glycosyl- and malonyl-transferases) but not those for PA biosynthesis (Martínez-Rivas et al. 2023). The activity was not redundant with the primary anthocyanin regulator, FaMYB10 (SG6), suggesting these also cooperatively act to regulate the production of more specialized anthocyanin types during fruit ripening. In blueberry, acylated anthocyanin species are restricted to late stages of fruit development when acyl transferases are expressed (Günther et al. 2020). It is not clear how these are regulated, but the lack of these more decorated forms of anthocyanins in vegetative tissues (Plunkett et al 2018) suggests additional regulators may be involved. Regulation of F3′5′H—a special target? The presence of flavonoids with 3ʹ4ʹ5ʹ hydroxylation (e.g. delphinidin-based anthocyanins, myricetin-based flavonols) is sporadic throughout angiosperm families, with F3ʹ5ʹH thought to have evolved numerous times from ancestral F3ʹH genes (Seitz et al. 2015). This may be reflected by differences in regulation of F3ʹ5ʹH between groups of plants. For example, within the Solanaceae, F3ʹ5ʹH genes are regulated by SG6 MYBs (MBW) (Albert et al. 2014; Butelli et al. 2021), whereas in other plants activity is restricted to particular tissues, developmental stages, or conditions. In grape, differential regulation of F3ʹ5ʹH genes occurs between different SG6 MYBs: VviMYBA1 strongly activates several F3ʹ5ʹH genes, but VviMYBA6 and VviMYBA7 have weak to no activity. This results in different anthocyanin profiles between berries (higher delphinidin:cyanidin ratio) and vegetative tissues because these MYBs have distinct spatial and temporal expression patterns, VviMYBA1 is the key determinant of fruit pigmentation, and VviMYBA6/7 regulate anthocyanin in vegetative tissues (Matus et al. 2017). In blueberry, vegetative tissues accumulate only simple cyanidin-glycosides instead of the complex mixtures of cyanidin, delphinidin, petunidin, and malvidin-based anthocyanins present in fruit (Plunkett et al. 2018). This suggests a lack of F3ʹ5ʹH activity that may be even more restricted than in grape. However, it remains unclear how this difference in anthocyanin profiles arises, given different observations and data. Overexpression of VcMYBA1 in blueberry plants resulted in red vegetative pigmentation, suggesting cyanidin accumulation (Lafferty et al. 2022a). This is somewhat surprising, given these plants express MYBPA1, MYBPA2 in leaves, associated with PA accumulation, in addition to the MYBA1 transgene (Lafferty et al. 2022a)—3 coexpressed MYBs that were each shown to activate the blueberry F3ʹ5ʹH in promoter activation assays (Lafferty et al. 2022b). And finally, cranberry has F3ʹ5ʹH activity associated with flavonol production during early stages of fruit development but lacks delphinidin-based anthocyanins (Vvedenskaya et al. 2004). Understanding how this occurs is of interest, given the role of F3ʹ5ʹH for generating diverse flavonoid compounds with altered properties, including color, stability, and health attributes. It is possible that distinct F3ʹ5ʹH genes are differentially regulated by flavonol versus anthocyanin MYBs, or perhaps it is regulated by a distinct regulator, as occurs in kiwifruit (Actinidia spp.) (Peng et al. 2019). Spatial and temporal patterning Within Vaccinium, many species have anthocyanin restricted to the berry skin/peel, but there are notable red-fleshed exceptions. The best characterized is bilberry, but the red-flesh trait is shared with other members within Vaccinium section Myrtillus and sporadically arises within other parts of the genus (e.g. deerberry, V. staminium). These differences are exquisitely captured by mass-spectrometry imaging (Dare et al. 2022), which showed the spatial distribution of flavonoids as they are produced at different developmental stages. Most striking is the advance of anthocyanin patterning in line with maturity, with a strong correlation with sugars and, ultimately, how all this apparently coordinated metabolite production can be defined by chemical signatures that progress from the skin into the flesh. How is this patterning established? In other species, altered expression of the SG6 MYB genes are central to generating new pigmentation patterns (reviewed in Albert et al. 2022). This can involve cis regulatory changes, such as promoter variants (e.g. apple; Espley et al. 2009) or transposon activation/enhancement (e.g. strawberry; Castillejo et al. 2020), or trans changes, altering the expression of upstream regulators of the MYB gene (e.g. red-flesh peach; Zhou et al. 2015; Hara-Kitagawa et al. 2020). In blueberry, MYBA1 transcripts are highly abundant in the skin tissues at ripening and present in flesh tissues in trace amounts, and in bilberry, MYBA1 transcripts are also expressed in the flesh, although still less than in skin (Lafferty et al. 2022a). This argues against red-flesh in bilberry being a loss of anthocyanin repressors, such as SG4 R2R3-MYBs, R3-MYB repressors, or SPLs, which antagonize the activity of MBW complexes (Aharoni et al. 2001; Gou et al. 2011; Albert et al. 2014; Cao et al. 2017). It is unknown if the genetic basis of red flesh in bilberry is linked to the MYBA locus itself or to another region, because it is not polymorphic for flesh color, preventing QTL mapping. Insights may come from inter-specific hybrids between deerberry and blueberry, which show that red flesh from deerberry is a dominant (or semi-dominant) trait (Lyrene 2021). If mapping populations can be generated, this would help address the basis of this trait in deerberry, which may be a shared mechanism in other red-fleshed Vaccinium species. Small-RNAs are increasingly recognized for their roles regulating the expression of key TFs that control flavonoid metabolism. This includes miR828 and miR858, which target SG5 R2R3-MYB and MYBPA1 genes involved in PA regulation in many plants and commonly have an inverse expression pattern to the transcript abundance for the target MYB gene in fruit (e.g. apple, kiwifruit, persimmon) (Li et al. 2019; Yang et al. 2020; Zhang et al. 2022). A similar pattern occurs in blueberry, with SG5/MYBPA2 transcript accumulation inverse to miR858b during blueberry fruit development, with degraded transcript detected at late stages when intact mRNAs are not abundant (Li et al. 2018). The pattern of PA accumulation in blueberry (and other fruits) closely follows the expression profile of SG5/MYBPA2 genes (Li et al. 2018; Lafferty et al. 2022a); thus, miRNA regulation may determine the temporal expression pattern of MYBPA2 genes during berry development. Genetic variation A pangenome for Vaccinium was recently developed (Yocca 2022). A single reference genome sequence does not contain all genes within a species because genes may be present in every individual (core) or absent in at least one individual (auxiliary; Golicz et al. 2020). The Vaccinium pangenome was mined to assess gene presence across dozens of blueberry genotypes (Supplemental Fig. S1). Although most flavonoid biosynthesis genes and transcription factors are retained across all genotypes, some are missing in a few genotypes and others are only present in a single genome. Variation in retention of these genes may lead to cross-accession differences in anthocyanin accumulation or structure, which must be accounted for in future breeding decisions. Future directions The recent surge in genomic resources and research has driven progress in understanding the production and regulation of the key flavonoids in Vaccinium. This now extends to functional data on the transcriptional regulation of the flavonoid pathway and the development of genetic markers. It also emphasizes the need to match resources across the range of Vaccinium species, allowing us to understand how, among other quality traits, breeders can develop improved cultivars and exploit the potential flavonoid-derived health benefits. Clearly, there is need for improved genetic resources in other Vaccinium species to help drive the breeding and selection of key traits, including the flavonoids that are so intrinsically linked with health properties. This need extends to the establishment of common gardens with sufficient germplasm resources for GxE-type experiments. Finally, the construction of Vaccinium pangenomes will help our understanding of wild relatives and the potential for introgression of valuable traits for further crop improvement. The rich diversity of flavonoids in Vaccinium species, coupled with both the newly acquired and developing genetic resources, will provide new tools to understand flavonoid production, which can be applied across a range of fruit crops. These tools can build on those developed in grape, strawberry, apple, and other crops to answer questions such as what governs the pathway flux toward di- or trihydroxylated anthocyanins or the decoration and acylation of flavonoids as well as the factors that control flavonoid accumulation in different fruit tissues (flesh and skin). Supplementary Material kiad250_Supplementary_Data Click here for additional data file. Acknowledgments We apologize to colleagues whose work could not be cited owing to space limitations. We thank PFR colleagues Drs Kevin Davies and Rebecca Bloomer for helpful advice on the manuscript. N.W.A. and R.V.E. were supported by The New Zealand Government Ministry of Business, Innovation and Employment Endeavour Fund “Filling the Void” (contract C11X1704). M.I. and M.M.F. were supported by the United States Department of Agriculture National Institute of Food and Agriculture, Hatch project 1008691, and the National Institute of Food and Agriculture, United States Department of Agriculture, under award number 2019–51181–30015, project “VacciniumCAP: Leveraging genetic and genomic resources to enable development of blueberry and cranberry cultivars with improved fruit quality attributes.” Supplementary data Supplemental Fig. S1 . Presence-absence variation (PAV) of anthocyanin biosynthesis genes across northern highbush (NHB; A) and southern highbush (SHB; B) blueberry genotypes. Light blue indicates the gene is both syntenic and present. Supplemental Table S1 . Summary of flavonoid biosynthetic gene predictions from Vaccinium species and Arabidopsis, apple, strawberry, and grape. Predicted peptide sequences of Vitis vinifera PN40024 v4.1, Actinidia chinensis Red5 v1.0, Malus domestica v1.1, Fragaria vesca v4.0.a2, Vaccinium caesariense W85-20 v2.0, Vaccinium myrtillus isolate NK2018 v1.0, and Vaccinium macrocarpon cv Stevens v1.0 were retrieved. KIPEs v3.2.2 was applied with the flavonoid bait sequence collection v3.0 to identify the anthocyanin biosynthesis genes in these datasets. Proteins involved in anthocyanin decoration and transport were additionally checked through phylogenetic trees based on MAFFT v7 alignments and constructed with FastTree v2.1 (Price et al. 2010; Katoh and Standley 2013). Supplemental Table S2 . Summary of QTL for chlorogenic acid and anthocyanins concentration, anthocyanin composition mapped in blueberry genetic studies. Data were extracted from Mengist et al. (2022b) and Montanari et al. (2022). Supplemental Table S3 . Summary of quantitative trait loci for TAcy, PAC, color variation, and color intensity mapped in cranberry genetic studies. Data were extracted from Diaz-Garcia et al. (2018) and Diaz-Garcia et al. (2021). ==== Refs References Abdelghafar A , da Silva LingeC, OkieWR, GasicK. Mapping QTLs for phytochemical compounds and fruit quality in peach. Mol Breed. 2020:40 (3 ):1–18. 10.1007/s11032-020-01114-y Aharoni A , De VosCR, WeinM, SunZ, GrecoR, KroonA, MolJN, O'ConnellAP. The strawberry FaMYB1 transcription factor suppresses anthocyanin and flavonol accumulation in transgenic tobacco. Plant J. 2001:28 (3 ):319–332. 10.1046/j.1365-313X.2001.01154.x 11722774 Akagi T , IkegamiA, YonemoriK. Dkmyb2 wound-induced transcription factor of persimmon (Diospyros kaki thunb.), contributes to proanthocyanidin regulation. Planta. 2010:232 (5 ):1045–1059. 10.1007/s00425-010-1241-7 20690029 Albert NW , AllanAC. MYB Genes involved in domestication and crop improvement. In: Roberts JA, editor. Annual plant reviews. 2021. p. 199–242. 10.1002/9781119312994.apr0767 Albert NW , ButelliE, MossSMA, PiazzaP, WaiteCN, SchwinnKE, DaviesKM, MartinC. Discrete bHLH transcription factors play functionally overlapping roles in pigmentation patterning in flowers of Antirrhinum majus. New Phytologist. 2021:231 (2 ):849–863. 10.1111/nph.17142 33616943 Albert NW , DaviesKM, LewisDH, ZhangH, MontefioriM, BrendoliseC, BoaseMR, NgoH, JamesonPE, SchwinnKE. A conserved network of transcriptional activators and repressors regulates anthocyanin pigmentation in eudicots. Plant Cell. 2014:26 (3 ):962–980. 10.1105/tpc.113.122069 24642943 Albert NW , LaffertyDJ, MossSMA, DaviesKM. Flavonoids—flowers, fruit, forage and the future. Journal of the Royal Society of New Zealand. 2022:53 (3 ): 304–331 10.1080/03036758.2022.2034654 Albert NW , LewisDH, ZhangH, SchwinnKE, JamesonPE, DaviesKM. Members of an R2R3-MYB transcription factor family in Petunia are developmentally and environmentally regulated to control complex floral and vegetative pigmentation patterning. Plant J. 2011:65 (5 ):771–784. 10.1111/j.1365-313X.2010.04465.x 21235651 Appelhagen I , NordholtN, SeidelT, SpeltK, KoesR, QuattrochioF, SagasserM, WeisshaarB. TRANSPARENT TESTA 13 is a tonoplast P3A-ATPase required for vacuolar deposition of proanthocyanidins in Arabidopsis thaliana seeds. Plant J. 2015:82 (5 ):840–849. 10.1111/tpj.12854 25891958 Ballester AR , MolthoffJ, de VosR, HekkertBTL, OrzaezD, Fernaݩndez-MorenoJP, TripodiP, GrandilloS, MartinC, HeldensJ, et al Biochemical and molecular analysis of pink tomatoes: deregulated expression of the gene encoding transcription factor SlMYB12 leads to pink tomato fruit color. Plant Physiol. 2010:152 (1 ):71–84. 10.1104/pp.109.147322 19906891 Battistoni B , SalazarJ, VegaW, Valderrama-SotoD, Jiménez-MuñozP, Sepúlveda-GonzálezA, AhumadaS, ChoI, GardanaCS, MoralesH. An upgraded, highly saturated linkage map of Japanese plum (prunus salicina lindl.), and identification of a new major locus controlling the flavan-3-ol composition in fruits. Front Plant Sci. 2022:13 : 805744. 10.3389/fpls.2022.805744 35310655 Baudry A , CabocheM, LepiniecL. TT8 Controls its own expression in a feedback regulation involving TTG1 and homologous MYB and bHLH factors, allowing a strong and cell-specific accumulation of flavonoids in Arabidopsis thaliana. Plant J. 2006:46 (5 ):768–779. 10.1111/j.1365-313X.2006.02733.x 16709193 Baudry A , HeimMA, DubreucqB, CabocheM, WeisshaarB, LepiniecL. TT2, TT8, and TTG1 synergistically specify the expression of BANYULS and proanthocyanidin biosynthesis in Arabidopsis thaliana. Plant J. 2004:39 (3 ):366–380. 10.1111/j.1365-313X.2004.02138.x 15255866 Bogs J , JafféFW, TakosAM, WalkerAR, RobinsonSP. The grapevine transcription factor VvMYBPA1 regulates proanthocyanidin synthesis during fruit development. Plant Physiol. 2007:143 (3 ):1347–1361. 10.1104/pp.106.093203 17208963 Bontpart T , FerreroM, KhaterF, MarlinT, VialetS, Vallverdù-QueraltA, PinasseauL, AgeorgesA, CheynierV, TerrierN. Focus on putative serine carboxypeptidase-like acyltransferases in grapevine. Plant Physiol Biochem. 2018:130 :356–366. 10.1016/j.plaphy.2018.07.023 30055344 Butelli E , BullingK, HillL, MartinC. Beyond purple tomatoes: combined strategies targeting anthocyanins to generate crimson, magenta, and indigo fruit. Horticulturae. 2021:7 (9 ):327. 10.3390/horticulturae7090327 Calle A , SerradillaMJ, WünschA. QTL Mapping of phenolic compounds and fruit color in sweet cherry using a 6 + 9 K SNP array genetic map. Sci Hortic. 2021:280 :109900. 10.1016/j.scienta.2021.109900 Cao X , QiuZ, WangX, Van GiangT, LiuX, WangJ, WangX, GaoJ, GuoY, DuY. A putative R3 MYB repressor is the candidate gene underlying atroviolacium, a locus for anthocyanin pigmentation in tomato fruit. J Exp Bot. 2017:68 (21–22 ):5745–5758. 10.1093/jxb/erx382 29186488 Carey AN , FisherDR, RimandoAM, GomesSM, BielinskiDF, Shukitt-HaleB. Stilbenes and anthocyanins reduce stress signaling in BV-2 mouse microglia. J Agric Food Chem. 2013:61 (25 ):5979–5986. 10.1021/jf400342g 23731018 Castillejo C , WaurichV, WagnerH, RamosR, OizaN, MuñozP, TriviñoJC, CaruanaJ, LiuZ, CoboN, et al Allelic variation of MYB10 is the major force controlling natural variation in skin and flesh color in strawberry (Fragaria spp.) fruit. Plant Cell. 2020:32 (12 ):3723–3749. 10.1105/tpc.20.00474 33004617 Cavallini E , ZenoniS, FinezzoL, GuzzoF, ZamboniA, AvesaniL, TornielliGB. Functional diversification of grapevine MYB5a and MYB5b in the control of flavonoid biosynthesis in a petunia anthocyanin regulatory mutant. Plant Cell Physiol. 2014:55 (3 ):517–534. 10.1093/pcp/pct190 24363289 Colle M , LeisnerCP, WaiCM, OuS, BirdKA, WangJ, WisecaverJH, YoccaAE, AlgerEI, TangH, et al Haplotype-phased genome and evolution of phytonutrient pathways of tetraploid blueberry. GigaScience. 2019:8 (3 ):giz012. 10.1093/gigascience/giz1012 30715294 Czemmel S , HeppelSC, BogsJ. R2r3 MYB transcription factors: key regulators of the flavonoid biosynthetic pathway in grapevine. Protoplasma. 2012:249 (S2 ):109–118. 10.1007/s00709-012-0380-z Czemmel S , StrackeR, WeisshaarB, CordonN, HarrisNN, WalkerAR, RobinsonSP, BogsJ. The grapevine R2R3-MYB transcription factor VvMYBF1 regulates flavonol synthesis in developing grape berries. Plant Physiol. 2009:151 (3 ):1513–1530. 10.1104/pp.109.142059 19741049 Dare AP , GüntherCS, GreyAC, GuoG, DemaraisNJ, CordinerS, McGhieTK, BoldinghH, HuntM, DengC, et al Resolving the developmental distribution patterns of polyphenols and related primary metabolites in bilberry (Vaccinium myrtillus) fruit. Food Chem. 2022:374 :131703. 10.1016/j.foodchem.2021.131703 34902814 Deluc L , BogsJ, WalkerAR, FerrierT, DecenditA, MerillonJ-M, RobinsonSP, BarrieuF. The transcription factor VvMYB5b contributes to the regulation of anthocyanin and proanthocyanidin biosynthesis in developing grape berries. Plant Physiol. 2008:147 (4 ):2041–2053. 10.1104/pp.108.118919 18539781 Diaz-Garcia L , Garcia-OrtegaLF, González-RodríguezM, DelayeL, IorizzoM, ZalapaJ. Chromosome-level genome assembly of the American cranberry (Vaccinium macrocarpon ait.) and its wild relative Vaccinium microcarpum. Front Plant Sci. 2021:12 :633310. 10.3389/fpls.2021.633310 33643360 Diaz-Garcia L , SchlautmanB, Covarrubias-PazaranG, MauleA, Johnson-CicaleseJ, GrygleskiE, VorsaN, ZalapaJ. Massive phenotyping of multiple cranberry populations reveals novel QTLs for fruit anthocyanin content and other important chemical traits. Mol Genet Genomics. 2018:293 (6 ):1379–1392. 10.1007/s00438-018-1464-z 29967963 Edger PP , IorizzoM, BassilNV, BenevenutoJ, FerrãoLFV, GiongoL, HummerK, LawasLMF, LeisnerCP, LiC, et al There and back again; historical perspective and future directions for Vaccinium breeding and research studies. Hortic Res. 2022:9 :uhac083. 10.1093/hr/uhac083 35611183 Ehlenfeldt M . Domestication of the highbush blueberry at Whitesbog, 1911–1916. In: IX International Vaccinium Symposium 810, Whitesbog, NJ, 147–152. 2008. Espley RV , BrendoliseC, ChagneD, Kutty-AmmaS, GreenS, VolzR, PutterillJ, SchoutenHJ, GardinerSE, HellensRP. Multiple repeats of a promoter segment causes transcription factor autoregulation in red apples. Plant Cell. 2009:21 (1 ):168–183. 10.1105/tpc.108.059329 19151225 Falginella L , CastellarinSD, TestolinR, GambettaGA, MorganteM, Di GasperoG. Expansion and subfunctionalisation of flavonoid 3′, 5′-hydroxylases in the grapevine lineage. BMC Genomics. 2010:1 :1–8. 10.1186/1471-2164-11-562 Feller A , MacHemerK, BraunEL, GrotewoldE. Evolutionary and comparative analysis of MYB and bHLH plant transcription factors. Plant J. 2011:66 (1 ):94–116. 10.1111/j.1365-313X.2010.04459.x 21443626 Gilbert JL , OlmsteadJW, ColquhounTA, LevinLA, ClarkDG, MoskowitzHR. Consumer-assisted selection of blueberry fruit quality traits. HortScience horts. 2014:49 (7 ):864–873. 10.21273/HORTSCI.49.7.864 Golicz AA , BayerPE, BhallaPL, BatleyJ, EdwardsD. Pangenomics comes of age: from bacteria to plant and animal applications. Trends Genet. 2020:36 (2 ):132–145. 10.1016/j.tig.2019.11.006 31882191 Gonzalez A , MendenhallJ, HuoY, LloydA. TTG1 complex MYBs, MYB5 and TT2, control outer seed coat differentiation. Dev Biol. 2009:325 (2 ):412–421. 10.1016/j.ydbio.2008.10.005 18992236 Gonzalez A , ZhaoM, LeavittJM, LloydAM. Regulation of the anthocyanin biosynthetic pathway by the TTG1/bHLH/myb transcriptional complex in Arabidopsis seedlings. Plant J. 2008:53 (5 ):814–827. 10.1111/j.1365-313X.2007.03373.x 18036197 Goodman CD , CasatiP, WalbotV. A multidrug resistance–associated protein involved in anthocyanin transport in Zea mays. Plant Cell. 2004:16 (7 ):1812–1826. 10.1105/tpc.022574 15208386 Gou J-Y , FelippesFF, LiuC-J, WeigelD, WangJ-W. Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor. Plant Cell. 2011:23 (4 ):1512–1522. 10.1105/tpc.111.084525 21487097 Günther CS , DareAP, McGhieTK, DengC, LaffertyDJ, PlunkettBJ, GriersonERP, TurnerJL, JaakolaL, AlbertNW, et al Spatiotemporal modulation of flavonoid metabolism in blueberries. Front Plant Sci. 2020:11 :545. 10.3389/fpls.2020.00545 32477384 Han Y , VimolmangkangS, Soria-GuerraRE, Rosales-MendozaS, ZhengD, LyginAV, KorbanSS. Ectopic expression of apple F3′H genes contributes to anthocyanin accumulation in the Arabidopsis tt7 mutant grown under nitrogen stress. Plant Physiol. 2010:153 (2 ):806–820. 10.1104/pp.109.152801 20357139 Hara-Kitagawa M , UnokiY, HiharaS, OdaK. Development of simple PCR-based DNA marker for the red-fleshed trait of a blood peach ‘tenshin-suimitsuto’. Mol Breed. 2020:40 (1 ):5. 10.1007/s11032-019-1068-2 Hartmann U , SagasserM, MehrtensF, StrackeR, WeisshaarB. Differential combinatorial interactions of cis-acting elements recognized by R2R3-MYB, BZIP, and BHLH factors control light-responsive and tissue-specific activation of phenylpropanoid biosynthesis genes. Plant Mol Biol. 2005:57 (2 ):155–171. 10.1007/s11103-004-6910-0 15821875 He J , MagnusonBA, GiustiMM. Analysis of anthocyanins in rat intestinal contents - impact of anthocyanin chemical structure on fecal excretion. J Agric Food Chem. 2005:53 (8 ):2859–2866. 10.1021/jf0479923 15826031 Hokkanen J , MattilaS, JaakolaL, PirttiläAM, TolonenA. Identification of phenolic compounds from lingonberry (Vaccinium vitis-idaea L.), bilberry (Vaccinium myrtillus L.) and hybrid bilberry (Vaccinium x intermedium ruthe L.) leaves. J Agric Food Chem. 2009:57 (20 ):9437–9447. 10.1021/jf9022542 19788243 Jaakola L , PooleM, JonesMO, Kamarainen-KarppinenT, KoskimakiJJ, HohtolaA, HaggmanH, FraserPD, ManningK, KingGJ, et al A SQUAMOSA MADS box gene involved in the regulation of anthocyanin accumulation in bilberry fruits. Plant Physiol. 2010:153 (4 ):1619–1629. 10.1104/pp.110.158279 20566708 James AM , MaDW, MellwayR, GesellA, YoshidaK, WalkerV, TranL, StewartD, ReicheltM, SuvantoJ, et al Poplar MYB115 and MYB134 transcription factors regulate proanthocyanidin synthesis and structure. Plant Physiol. 2017:174 (1 ):154–171. 10.1104/pp.16.01962 28348066 Johnson ET , RyuS, YiH, ShinB, CheongH, ChoiG. Alteration of a single amino acid changes the substrate specificity of dihydroflavonol 4-reductase. Plant J. 2001:25 (3 ):325–333. 10.1046/j.1365-313x.2001.00962.x 11208024 Jun JH , XiaoX, RaoX, DixonRA. Proanthocyanidin subunit composition determined by functionally diverged dioxygenases. Nat Plants. 2018:4 (12 ):1034–1043. 10.1038/s41477-018-0292-9 30478357 Kalt W , CassidyA, HowardLR, KrikorianR, StullAJ, TremblayF, Zamora-RosR. Recent research on the health benefits of blueberries and their anthocyanins. Adv Nutr. 2020:11 (2 ):224–236. 10.1093/advances/nmz065 31329250 Kalt W , McDonaldJ, RickerR, LuX. Anthocyanin content and profile within and among blueberry species. Can J Plant Sci. 1999:79 (4 ):617–623. 10.4141/P99-009 Karppinen K , LaffertyDJ, AlbertNW, MikkolaN, McGhieT, AllanAC, AfzalBM, HäggmanH, EspleyRV, JaakolaL. MYBA And MYBPA transcription factors co-regulate anthocyanin biosynthesis in blue-colored berries. New Phytologist. 2021:232 (3 ):1350–1367. 10.1111/nph.17669 34351627 Karppinen K , TegelbergP, HäggmanH, JaakolaL. Abscisic acid regulates anthocyanin biosynthesis and gene expression associated with cell wall modification in ripening bilberry (Vaccinium myrtillus L.) fruits. Front Plant Sci. 2018:9 :1259. 10.3389/fpls.2018.01259 30210522 Karppinen K , ZorattiL, NguyenquynhN, HäggmanH, JaakolaL. On the developmental and environmental regulation of secondary metabolism in Vaccinium spp. berries. Front Plant Sci. 2016:7 :655. 10.3389/fpls.2016.00655 27242856 Katoh K , StandleyDM. MAFFT Multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013:30 (4 ):772–780. 10.1093/molbev/mst010 23329690 Kowalska K . Lingonberry (Vaccinium vitis-idaea L.) fruit as a source of bioactive compounds with health-promoting effects—a review. Int J Mol Sci. 2021:22 (10 ):5126. 10.3390/ijms22105126 34066191 Kumar S , MolloyC, HuntM, DengCH, WiedowC, AndreC, DareA, McGhieT. GWAS Provides new insights into the genetic mechanisms of phytochemicals production and red skin color in apple. Hortic Res. 2022:9 :uhac218. 10.1093/hr/uhac218 36479587 Labadie M , VallinG, PetitA, RingL, HoffmannT, GastonA, PotierA, SchwabW, RothanC, DenoyesB. Metabolite quantitative trait loci for flavonoids provide new insights into the genetic architecture of strawberry (Fragaria × ananassa) fruit quality. J Agric Food Chem. 2020:68 (25 ):6927–6939. 10.1021/acs.jafc.0c01855 32469530 Lafferty DJ , EspleyRV, DengCH, DareAP, GüntherCS, JaakolaL, KarppinenK, BoaseMR, WangL, LuoH, et al The coordinated action of MYB activators and repressors controls proanthocyanidin and anthocyanin biosynthesis in Vaccinium. Front Plant Sci. 2022a:13 :910155. 10.3389/fpls.2022.910155 35812927 Lafferty DJ , EspleyRV, DengCH, GüntherCS, PlunkettB, TurnerJL, JaakolaL, KarppinenK, AllanAC, AlbertNW. Hierarchical regulation of MYBPA1 by anthocyanin- and proanthocyanidin-related MYB proteins is conserved in Vaccinium species. J Exp Bot. 2022b:73 (5 ):1344–1356. 10.1093/jxb/erab460 34664645 Li Y , CuiW, WangR, LinM, ZhongY, SunL, QiX, FangJ. MicroRNA858-mediated regulation of anthocyanin biosynthesis in kiwifruit (Actinidia arguta) based on small RNA sequencing. PLoS One. 2019:14 :e0217480. 10.1371/journal.pone.0217480 31120996 Li SF , MillikenON, PhamH, SeyitR, NapoliR, PrestonJ, KoltunowAM, ParishRW. The Arabidopsis MYB5 transcription factor regulates mucilage synthesis, seed coat development, and trichome morphogenesis. Plant Cell. 2009:21 (1 ):72–89. 10.1105/tpc.108.063503 19136646 Li G , WangY, LouX, LiH, ZhangC. Identification of blueberry miRNAs and their targets based on high-throughput sequencing and degradome analyses. Int J Mol Sci. 2018:19 (4 ):983. 10.3390/ijms19040983 29587414 Liu C , JunJH, DixonRA. MYB5 And MYB14 play pivotal roles in seed coat polymer biosynthesis in Medicago truncatula. Plant Physiol. 2014:165 (4 ):1424–1439. 10.1104/pp.114.241877 24948832 Luby JJ , BallingtonJR, DraperAD, PliszkaK, AustinME. Blueberries and cranberries (Vaccinium). Genetic Resources of Temperate Fruit and Nut Crops. 1991:290 :393–458. doi: 10.17660/ActaHortic.1991.290.9 Lyrene PM . Breeding cultivars from blueberry × deerberry hybrids: progress and prospects. HortScience. 2021:56 (4 ):439–446. 10.21273/HORTSCI15619-20 Lyrene PM , VorsaN, BallingtonJR. Polyploidy and sexual polyploidization in the genus Vaccinium. Euphytica. 2003:133 (1 ):27–36. 10.1023/A:1025608408727 Ma L , SunZ, ZengY, LuoM, YangJ. Molecular mechanism and health role of functional ingredients in blueberry for chronic disease in human beings. Int J Mol Sci. 2018:19 (9 ):2785. 10.3390/ijms19092785 30223619 Martínez-Rivas FJ , Blanco-PortalesR, SerratosaMP, Ric-VarasP, Guerrero-SánchezV, Medina-PucheL, MoyanoL, MercadoJA, AlseekhS, CaballeroJL, et al Molina-Hidalgo FJ FaMYB123 interacts with FabHLH3 to regulate the late steps of anthocyanin and flavonol biosynthesis during ripening. Plant J. 2023:114 (3 ):683–698. 10.1111/tpj.16166 36840368 Matus JT , CavalliniE, LoyolaR, HöllJ, FinezzoL, Dal SantoS, VialetS, CommissoM, RomanF, SchubertA, et al A group of grapevine MYBA transcription factors located in chromosome 14 control anthocyanin synthesis in vegetative organs with different specificities compared with the berry color locus. Plant J. 2017:91 (2 ):220–236. 10.1111/tpj.13558 28370629 Mehrtens F , KranzH, BednarekP, WeisshaarB. The Arabidopsis transcription factor MYB12 is a flavonol-specific regulator of phenylpropanoid biosynthesis. Plant Physiol. 2005:138 (2 ):1083–1096. 10.1104/pp.104.058032 15923334 Mengist MF , BostanH, De PaolaD, TeresiSJ, PlattsAE, CremonaG, QiX, MackeyT, BassilNV, AshrafiH, et al Autopolyploid inheritance and a heterozygous reciprocal translocation shape chromosome genetic behavior in tetraploid blueberry (Vaccinium corymbosum). New Phytologist. 2022a:237 (3 ):1024–1039. 10.1111/nph.18428 35962608 Mengist MF , BurtchH, DebeloH, PottorffM, BostanH, NunnC, CorbinS, KayCD, BassilN, HummerK, et al Development of a genetic framework to improve the efficiency of bioactive delivery from blueberry. Sci Rep. 2020:10 (1 ):17311. 10.1038/s41598-020-74280-w 33057109 Mengist MF , GraceMH, MackeyT, MunozB, PuckerB, BassilN, LubyC, FerruzziM, LilaMA, IorizzoM. Dissecting the genetic basis of bioactive metabolites and fruit quality traits in blueberries (Vaccinium corymbosum L.). Front Plant Sci. 2022b:13 :964656. 10.3389/fpls.2022.964656 36119607 Miosic S , ThillJ, MilosevicM, GoschC, PoberS, MolitorC, EjazS, RompelA, StichK, HalbwirthH. Dihydroflavonol 4-reductase genes encode enzymes with contrasting substrate specificity and show divergent gene expression profiles in Fragaria species. PLoS One. 2014:9 (11 ):e112707. 10.1371/journal.pone.0112707 25393679 Montanari S , ThomsonS, CordinerS, GüntherCS, MillerP, DengCH, McGhieT, KnäbelM, FosterT, TurnerJ, et al High-density linkage map construction in an autotetraploid blueberry population and detection of quantitative trait loci for anthocyanin content. Front Plant Sci. 2022:13 :965397. 10.3389/fpls.2022.965397 36247546 Montefiori M , BrendoliseC, DareAP, Lin-WangK, DaviesKM, HellensRP, AllanAC. In the Solanaceae, a hierarchy of bHLHs confer distinct target specificity to the anthocyanin regulatory complex. J Exp Bot. 2015:66 (5 ):1427–1436. 10.1093/jxb/eru494 25628328 Mou J , ZhangZ, QiuH, LuY, ZhuX, FanZ, ZhangQ, YeJ, FernieAR, ChengY, et al Multiomics-based dissection of citrus flavonoid metabolism using a Citrus reticulata × Poncirus trifoliata population. Hortic Res. 2021:8 :56. 10.1038/s41438-021-00472-8 33642588 Nemzer BV , Al-TaherF, YashinA, RevelskyI, YashinY. Cranberry: chemical composition, antioxidant activity and impact on human health: overview. Molecules. 2022:27 (5 ):1503. 10.3390/molecules27051503 35268605 Nesi N , DebeaujonI, JondC, PelletierG, CabocheM, LepiniecL. The TT8 gene encodes a basic helix-loop-helix domain protein required for expression of DFR and BAN genes in Arabidopsis siliques. Plant Cell. 2000:12 (10 ):1863–1878. 10.1105/tpc.12.10.1863 11041882 Norberto S , SilvaS, MeirelesM, FariaA, PintadoM, CalhauC. Blueberry anthocyanins in health promotion: a metabolic overview. J Funct Foods. 2013:5 (4 ):1518–1528. 10.1016/j.jff.2013.08.015 Oh HD , YuDJ, ChungSW, CheaS, LeeHJ. Abscisic acid stimulates anthocyanin accumulation in ‘Jersey’ highbush blueberry fruits during ripening. Food Chem. 2018:244 :403–407. 10.1016/j.foodchem.2017.10.051 29120800 Paauw M , KoesR, QuattrocchioFM. Alteration of flavonoid pigmentation patterns during domestication of food crops. J Exp Bot. 2019:70 (15 ):3719–3735. 10.1093/jxb/erz141 30949670 Peng Y , Lin-WangK, CooneyJM, WangT, EspleyRV, AllanAC. Differential regulation of the anthocyanin profile in purple kiwifruit (Actinidia species). Hortic Res. 2019:6 (1 ):3. 10.1038/s41438-018-0076-4 30622721 Plunkett BJ , EspleyRV, DareAP, WarrenBAW, GriersonERP, CordinerS, TurnerJL, AllanAC, AlbertNW, DaviesKM, et al MYBA From blueberry (Vaccinium section Cyanococcus) is a subgroup 6 type R2R3MYB transcription factor that activates anthocyanin production. Front Plant Sci. 2018:9 :1300. 10.3389/fpls.2018.01300 30254656 Pott DM , VallarinoJG, Cruz-RusE, WillmitzerL, Sánchez-SevillaJF, AmayaI, OsorioS. Genetic analysis of phenylpropanoids and antioxidant capacity in strawberry fruit reveals mQTL hotspots and candidate genes. Sci Rep. 2020:10 (1 ):1–15. 10.1038/s41598-020-76946-x 31913322 Price MN , DehalPS, ArkinAP. Fasttree 2–approximately maximum-likelihood trees for large alignments. PloS One. 2010:5 (3 ):e9490. 10.1371/journal.pone.0009490 20224823 Primetta AK , KarppinenK, RiihinenKR, JaakolaL. Metabolic and molecular analyses of white mutant Vaccinium berries show down-regulation of MYBPA1-type R2R3 MYB regulatory factor. Planta. 2015:242 (3 ):631–643. 10.1007/s00425-015-2363-8 26168981 Prior RL , LazarusSA, CaoG, MuccitelliH, HammerstoneJF. Identification of procyanidins and anthocyanins in blueberries and cranberries (Vaccinium spp.) using high-performance liquid chromatography/mass spectrometry. J Agric Food Chem. 2001:49 (3 ):1270–1276. 10.1021/jf001211q 11312849 Pucker B , ReiherF, SchilbertHM. Automatic identification of players in the flavonoid biosynthesis with application on the biomedicinal plant Croton tiglium. Plants. 2020:9 (9 ):1103. 10.3390/plants9091103 32867203 Pucker B , SelmarD. Biochemistry and molecular basis of intracellular flavonoid transport in plants. Plants. 2022:11 (7 ):963. 10.3390/plants11070963 35406945 Quattrocchio F , VerweijW, KroonA, SpeltC, MolJ, KoesR. PH4 Of petunia is an R2R3 MYB protein that activates vacuolar acidification through interactions with basic-helix-loop-helix transcription factors of the anthocyanin pathway. Plant Cell. 2006:18 (5 ):1274–1291. 10.1105/tpc.105.034041 16603655 Retamales J , HancockJ. Blueberry taxonomy and breeding. Blueberries. 2018:2 :18–60. 10.1079/9781780647265.0018 Riihinen K , JaakolaL, KärenlampiS, HohtolaA. Organ-specific distribution of phenolic compounds in bilberry (Vaccinium myrtillus) and ‘Northblue’ blueberry (Vaccinium corymbosum x V. angustifolium). Food Chem. 2008:110 (1 ):156–160. 10.1016/j.foodchem.2008.01.057 26050178 Scarano A , ButelliE, De SantisS, CavalcantiE, HillL, De AngelisM, GiovinazzoG, ChieppaM, MartinC, SantinoA. Combined dietary anthocyanins, flavonols, and stilbenoids alleviate inflammatory bowel disease symptoms in mice. Front Nutr. 2017:4 :75. 10.3389/fnut.2017.00075 29473042 Schwinn K , VenailJ, ShangY, MackayS, AlmV, ButelliE, OyamaR, BaileyP, DaviesK, MartinC. A small family of MYB-regulatory genes controls floral pigmentation intensity and patterning in the genus Antirrhinum. Plant Cell. 2006:18 (4 ):831–851. 10.1105/tpc.105.039255 16531495 Seitz C , AmeresS, SchlangenK, ForkmannG, HalbwirthH. Multiple evolution of flavonoid 3′,5′-hydroxylase. Planta. 2015:242 (3 ):561–573. 10.1007/s00425-015-2293-5 25916309 Sharma A , LeeHJ. Anti-inflammatory activity of bilberry (Vaccinium myrtillus L.). Curr Issues Mol Biol. 2022:44 (10 ):4570–4583. 10.3390/cimb44100313 36286028 Singh AP , WilsonT, KalkAJ, CheongJ, VorsaN. Isolation of specific cranberry flavonoids for biological activity assessment. Food Chem. 2009:116 (4 ):963–968. 10.1016/j.foodchem.2009.03.062 20161027 Spelt C , QuattrocchioF, MolJNM, KoesR. Anthocyanin1 of petunia encodes a basic helix-loop-helix protein that directly activates transcription of structural anthocyanin genes. Plant Cell. 2000:12 (9 ):1619–1631. 10.1105/tpc.12.9.1619 11006336 Spelt C , QuattrocchioF, MolJ, KoesR. ANTHOCYANIN1 Of petunia controls pigment synthesis, vacuolar pH, and seed coat development by genetically distinct mechanisms. Plant Cell. 2002:14 (9 ):2121–2135. 10.1105/tpc.003772 12215510 Stracke R , WerberM, WeisshaarB. The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001:4 (5 ):447–456. 10.1016/S1369-5266(00)00199-0 11597504 Sun L , LiS, JiangJ, TangX, FanX, ZhangY, LiuJ, LiuC. New quantitative trait locus (QTLs) and candidate genes associated with the grape berry color trait identified based on a high-density genetic map. BMC Plant Biol. 2020:20 (1 ):302. 10.1186/s12870-020-02517-x 32605636 Terrier N , TorregrosaL, AgeorgesA, VialetS, VerrièsC, CheynierV, RomieuC. Ectopic expression of VvMybPA2 promotes proanthocyanidin biosynthesis in grapevine and suggests additional targets in the pathway. Plant Physiol. 2009:149 (2 ):1028–1041. 10.1104/pp.108.131862 19098092 Tohge T , de SouzaLP, FernieAR. Current understanding of the pathways of flavonoid biosynthesis in model and crop plants. J Exp Bot. 2017:68 (15 ):4013–4028. 10.1093/jxb/erx177 28922752 Tundis R , TenutaMC, LoizzoMR, BonesiM, FinettiF, TrabalziniL, DeguinB. Vaccinium species (Ericaceae): from chemical composition to bio-functional activities. Applied Sciences. 2021:11 (12 ):5655. 10.3390/app11125655 Valderrama-Soto D , SalazarJ, Sepúlveda-GonzálezA, Silva-AndradeC, GardanaC, MoralesH, BattistoniB, Jiménez-MuñozP, GonzálezM, Peña-NeiraÁ. Detection of quantitative trait loci controlling the content of phenolic compounds in an Asian plum (Prunus salicina L.) F1 population. Front Plant Sci. 2021:12 :679059. 10.3389/fpls.2021.679059 34305977 Veberic R , SlatnarA, BizjakJ, StamparF, Mikulic-PetkovsekM. Anthocyanin composition of different wild and cultivated berry species. LWT-Food Sci Technol. 2015:60 (1 ):509–517. 10.1016/j.lwt.2014.08.033 Verdier J , ZhaoJ, Torres-JerezI, GeS, LiuC, HeX, MysoreKS, DixonRA, UdvardiMK. MtPAR MYB transcription factor acts as an on switch for proanthocyanidin biosynthesis in Medicago truncatula. Proc Natl Acad Sci USA. 2012:109 (5 ):1766–1771. 10.1073/pnas.1120916109 22307644 Vorsa N , PolashockJJ. Alteration of anthocyanin glycosylation in cranberry through interspecific hybridization. J Am Soc Hort Sci. 2005:130 (5 ):711–715. 10.21273/JASHS.130.5.711 Vorsa N , PolashockJ, CunninghamD, RoderickR. Genetic inferences and breeding implications from analysis of cranberry germplasm anthocyanin profiles. J Am Soc Hort Sci. 2003:128 (5 ):691–697. 10.21273/JASHS.128.5.0691 Vorsa N , ZalapaJ. Domestication, genetics, and genomics of the American cranberry. Plant Breed Rev. 2019:43 :279–315. 10.1002/9781119616801.ch8 Vvedenskaya IO , RosenRT, GuidoJE, RussellDJ, MillsKA, VorsaN. Characterization of flavonols in cranberry (Vaccinium macrocarpon) powder. J Agric Food Chem. 2004:52 (2 ):188–195. 10.1021/jf034970s 14733493 Walker AR , LeeE, BogsJ, McDavidDA, ThomasMR, RobinsonSP. White grapes arose through the mutation of two similar and adjacent regulatory genes. Plant J. 2007:49 (5 ):772–785. 10.1111/j.1365-313X.2006.02997.x 17316172 Wang N , QuC, JiangS, ChenZ, XuH, FangH, SuM, ZhangJ, WangY, LiuW, et al The proanthocyanidin-specific transcription factor MdMYBPA1 initiates anthocyanin synthesis under low-temperature conditions in red-fleshed apples. Plant J. 2018:96 (1 ):39–55. 10.1111/tpj.14013 29978604 Wang N , XuH, JiangS, ZhangZ, LuN, QiuH, QuC, WangY, WuS, ChenX. MYB12 And MYB22 play essential roles in proanthocyanidin and flavonol synthesis in red-fleshed apple (Malus sieversii f. niedzwetzkyana). Plant J. 2017:90 (2 ):276–292. 10.1111/tpj.13487 28107780 Wu C , DengC, HilarioE, AlbertNW, LaffertyD, GriersonERP, PlunkettBJ, ElboroughC, SaeiA, GuntherCS, et al A chromosome-scale assembly of the bilberry genome identifies a complex locus controlling berry anthocyanin composition. Mol Ecol Resour. 2021:345–360. 10.1111/1755-0998.13467 34260155 Wu Y , Popovsky-SaridS, TikunovY, BorovskyY, BaruchK, VisserRG, ParanI, BovyA. CaMYB12-like underlies a major QTL for flavonoid content in pepper (Capsicum annuum) fruit. New Phytologist. 2023:237 (6 ):2255–2267. 10.1111/nph.18693 36545937 Xu W , DubosC, LepiniecL. Transcriptional control of flavonoid biosynthesis by MYB-bHLH-WDR complexes. Trends Plant Sci. 2015:20 (3 ):176–185. 10.1016/j.tplants.2014.12.001 25577424 Yang S , ZhangM, XuL, LuoZ, ZhangQ. Mir858b inhibits proanthocyanidin accumulation by the repression of DkMYB19 and DkMYB20 in persimmon. Front Plant Sci. 2020:11 : 576378. 10.3389/fpls.2020.57637 Yocca AE . Leveraging angiosperm pangenomics to understand genome evolution. Ph.D. thesis. East Lansing, MI, Michigan State University. 2022. Yonekura-Sakakibara K , SaitoK. Function, structure, and evolution of flavonoid glycosyltransferases in plants. In: Romani A, Lattanzio V, Quideau S, editors. Recent advances in polyphenol research. 2014. p. 61–82. 10.1002/9781118329634.ch3 Zhai R , ZhaoY, WuM, YangJ, LiX, LiuH, WuT, LiangF, YangC, WangZ, et al The MYB transcription factor PbMYB12b positively regulates flavonol biosynthesis in pear fruit. BMC Plant Biol. 2019:19 (1 ):85. 10.1186/s12870-019-1687-0 30791875 Zhang Y-L , Lin-WangK, AlbertNW, ElboroughC, EspleyRV, AndreCM, FangZ-Z. Identification of a strong anthocyanin activator, VbMYBA, from berries of Vaccinium bracteatum thunb. Front Plant Sci. 2021:12 :697212. 10.3389/fpls.2021.697212 34938303 Zhang B , YangH-J, QuD, ZhuZ-Z, YangY-Z, ZhaoZ-Y. The MdBBX22–miR858–MdMYB9/11/12 module regulates proanthocyanidin biosynthesis in apple peel. Plant Biotechnol J. 2022:20 (9 ):1683–1700. 10.1111/pbi.13839 35527510 Zhang Y , YeJ, LiuC, XuQ, LongL, DengX. Citrus PH4–noemi regulatory complex is involved in proanthocyanidin biosynthesis via a positive feedback loop. J Exp Bot. 2019:71 (4 ):1306–1321. 10.1093/jxb/erz506 Zhao J , DixonRA. MATE Transporters facilitate vacuolar uptake of epicatechin 3′-O-glucoside for proanthocyanidin biosynthesis in medicago truncatula and arabidopsis. Plant Cell. 2009:21 (8 ):2323–2340. 10.1105/tpc.109.067819 19684242 Zhao J , HuhmanD, ShadleG, HeXZ, SumnerLW, TangY, DixonRA. MATE2 Mediates vacuolar sequestration of flavonoid glycosides and glycoside malonates in Medicago truncatula. Plant Cell. 2011:23 (4 ):1536–1555. 10.1105/tpc.110.080804 21467581 Zhou H , Lin-WangK, WangH, GuC, DareAP, EspleyRV, HeH, AllanAC, HanY. Molecular genetics of blood-fleshed peach reveals activation of anthocyanin biosynthesis by NAC transcription factors. Plant J. 2015:82 (1 ):105–121. 10.1111/tpj.12792 25688923 Zoratti L , JaakolaL, HäggmanH, GiongoL. Anthocyanin profile in berries of wild and cultivated Vaccinium spp. Along altitudinal gradients in the Alps. J Agric Food Chem. 2015:63 (39 ):8641–8650. 10.1021/acs.jafc.5b02833 26373665 Zorenc Z , VebericR, SlatnarA, KoronD, MiosicS, ChenM-H, Haselmair-GoschC, HalbwirthH, Mikulic-PetkovsekM. A wild ‘albino’ bilberry (Vaccinium myrtillus L.) from Slovenia shows three bottlenecks in the anthocyanin pathway and significant differences in the expression of several regulatory genes compared to the common blue berry type. PLoS One. 2017:12 (12 ):e0190246. 10.1371/journal.pone.0190246 29272302