miR164-MhNAC1 regulates apple root nitrogen uptake under low nitrogen stress.

Nitrogen is an essential nutrient for plant growth and serves as a signaling molecule to regulate gene expression inducing physiological, growth and developmental responses. An excess or deficiency of nitrogen may have adverse effects on plants. Studying nitrogen uptake will help us understand the molecular mechanisms of utilization for targeted molecular breeding. Here, we identified and functionally validated an NAC (NAM-ATAF1/2-CUC2) transcription factor based on the transcriptomes of two apple rootstocks with different nitrogen uptake efficiency. NAC1, a target gene of miR164, directly regulates the expression of the high-affinity nitrate transporter (MhNRT2.4) and citric acid transporter (MhMATE), affecting root nitrogen uptake. To examine the role of MhNAC1 in nitrogen uptake, we produced transgenic lines that overexpressed or silenced MhNAC1. Silencing MhNAC1 promoted nitrogen uptake and citric acid secretion in roots, and enhanced plant tolerance to low nitrogen conditions, while overexpression of MhNAC1 or silencing miR164 had the opposite effect. This study not only revealed the role of the miR164-MhNAC1 module in nitrogen uptake in apple rootstocks but also confirmed that citric acid secretion in roots affected nitrogen uptake, which provides a research basis for efficient nitrogen utilization and molecular breeding in apple.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing.

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Constitutive expression of apple endo-POLYGALACTURONASE1 in fruit induces early maturation, alters skin structure and accelerates softening.

During fruit ripening, polygalacturonases (PGs) are key contributors to the softening process in many species. Apple is a crisp fruit that normally exhibits only minor changes to cell walls and limited fruit softening. Here, we explore the effects of PG overexpression during fruit development using transgenic apple lines overexpressing the ripening-related endo-POLYGALACTURONASE1 gene. MdPG1-overexpressing (PGox) fruit displayed early maturation/ripening with black seeds, conversion of starch to sugars and ethylene production occurring by 80 days after pollination (DAP). PGox fruit exhibited a striking, white-skinned phenotype that was evident from 60 DAP and most likely resulted from increased air spaces and separation of cells in the hypodermis due to degradation of the middle lamellae. Irregularities in the integrity of the epidermis and cuticle were also observed. By 120 DAP, PGox fruit cracked and showed lenticel-associated russeting. Increased cuticular permeability was associated with microcracks in the cuticle around lenticels and was correlated with reduced cortical firmness at all time points and extensive post-harvest water loss from the fruit, resulting in premature shrivelling. Transcriptomic analysis suggested that early maturation was associated with upregulation of genes involved in stress responses, and overexpression of MdPG1 also altered the expression of genes involved in cell wall metabolism (e.g. ß-galactosidase, MD15G1221000) and ethylene biosynthesis (e.g. ACC synthase, MD14G1111500). The results show that upregulation of PG not only has dramatic effects on the structure of the fruit outer cell layers, indirectly affecting water status and turgor, but also has unexpected consequences for fruit development.

Genotype and phenotype data standardization, utilization and integration in the big data era for agricultural sciences.

Large-scale genotype and phenotype data have been increasingly generated to identify genetic markers, understand gene function and evolution and facilitate genomic selection. These datasets hold immense value for both current and future studies, as they are vital for crop breeding, yield improvement and overall agricultural sustainability. However, integrating these datasets from heterogeneous sources presents significant challenges and hinders their effective utilization. We established the Genotype-Phenotype Working Group in November 2021 as a part of the AgBioData Consortium (https://www.agbiodata.org) to review current data types and resources that support archiving, analysis and visualization of genotype and phenotype data to understand the needs and challenges of the plant genomic research community. For 2021-22, we identified different types of datasets and examined metadata annotations related to experimental design/methods/sample collection, etc. Furthermore, we thoroughly reviewed publicly funded repositories for raw and processed data as well as secondary databases and knowledgebases that enable the integration of heterogeneous data in the context of the genome browser, pathway networks and tissue-specific gene expression. Based on our survey, we recommend a need for (i) additional infrastructural support for archiving many new data types, (ii) development of community standards for data annotation and formatting, (iii) resources for biocuration and (iv) analysis and visualization tools to connect genotype data with phenotype data to enhance knowledge synthesis and to foster translational research. Although this paper only covers the data and resources relevant to the plant research community, we expect that similar issues and needs are shared by researchers working on animals. Database URL: https://www.agbiodata.org.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes.

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

A multiplexed plant-animal SNP array for selective breeding and species conservation applications.

Reliable and high-throughput genotyping platforms are of immense importance for identifying and dissecting genomic regions controlling important phenotypes, supporting selection processes in breeding programs, and managing wild populations and germplasm collections. Amongst available genotyping tools, single nucleotide polymorphism arrays have been shown to be comparatively easy to use and generate highly accurate genotypic data. Single-species arrays are the most commonly used type so far; however, some multi-species arrays have been developed for closely related species that share single nucleotide polymorphism markers, exploiting inter-species cross-amplification. In this study, the suitability of a multiplexed plant-animal single nucleotide polymorphism array, including both closely and distantly related species, was explored. The performance of the single nucleotide polymorphism array across species for diverse applications, ranging from intra-species diversity assessments to parentage analysis, was assessed. Moreover, the value of genotyping pooled DNA of distantly related species on the single nucleotide polymorphism array as a technique to further reduce costs was evaluated. Single nucleotide polymorphism performance was generally high, and species-specific single nucleotide polymorphisms proved suitable for diverse applications. The multi-species single nucleotide polymorphism array approach reported here could be transferred to other species to achieve cost savings resulting from the increased throughput when several projects use the same array, and the pooling technique adds another highly promising advancement to additionally decrease genotyping costs by half.

Dihydrochalcone glycoside biosynthesis in Malus is regulated by two MYB-like transcription factors and is required for seed development.

Dihydrochalcones (DHCs) including phlorizin (phloretin 2'-O-glucoside) and its positional isomer trilobatin (phloretin 4'-O-glucoside) are the most abundant phenylpropanoids in apple (Malus spp.). Transcriptional regulation of DHC production is poorly understood despite their importance in insect- and pathogen-plant interactions in human physiology research and in pharmaceuticals. In this study, segregation in hybrid populations and bulked segregant analysis showed that the synthesis of phlorizin and trilobatin in Malus leaves are both single-gene-controlled traits. Promoter sequences of PGT1 and PGT2, two glycosyltransferase genes involved in DHC glycoside synthesis, were shown to discriminate Malus with different DHC glycoside patterns. Differential PGT1 and PGT2 promoter activities determined DHC glycoside accumulation patterns between genotypes. Two transcription factors containing MYB-like DNA-binding domains were then shown to control DHC glycoside patterns in different tissues, with PRR2L mainly expressed in leaf, fruit, flower, stem, and seed while MYB8L mainly expressed in stem and root. Further hybridizations between specific genotypes demonstrated an absolute requirement for DHC glycoside production in Malus during seed development which explains why no Malus spp. with a null DHC chemotype have been reported.

Comparative genomics uncovered differences between clinical and environmental populations of Vibrio parahaemolyticus in New Zealand.

Vibrio parahaemolyticus has been identified as an emerging human pathogen worldwide with cases undergoing a global expansion over recent decades in phase with climate change. New Zealand had remained free of outbreaks until 2019, but different outbreaks have been reported consecutively since then. To provide new insights into the recent emergence of cases associated with outbreak clones over recent years, a comparative genomic study was carried out using a selection of clinical (mostly outbreak) and environmental isolates of V. parahaemolyticus obtained in New Zealand between 1973 and 2021. Among 151 isolates of clinical (n=60) and environmental (n=91) origin, 47 sequence types (STs) were identified, including 31 novel STs. The population of environmental isolates generated 30 novel STs, whereas only 1 novel ST (ST2658) was identified among the population of clinical isolates. The novel clinical ST was a single-locus variant of the pandemic ST36 strain, indicating further evolution of this pandemic strain. The environmental isolates exhibited a significant genetic heterogeneity compared to the clinical isolates. The whole-genome phylogeny separated the population of clinical isolates from their environmental counterparts, clearly indicating their distant genetic relatedness. In addition to differences in ancestral profiles and genetic relatedness, these two groups of isolates exhibited a profound difference in their virulence profiles. While the entire population of clinical isolates harboured the thermostable direct haemolysin (tdh) and/or the thermostable-related haemolysin (trh), only a few isolates of environmental origin possessed the same virulence genes. In contrast to tdh and trh, adhesin-encoding genes, vpadF and MSHA, showed a significantly (P<0.001) greater association with the environmental isolates compared to the clinical isolates. The effectors, VopQ, VPA0450 and VopS, which belong to T3SS1, were ubiquitous, being present in each isolate regardless of its origin. The effectors VopC and VopA, which belong to T3SS2, were rarely detected in any of the examined isolates. Our data indicate that the clinical and environmental isolates of V. parahaemolyticus from New Zealand differ in their population structures, ancestral profiles, genetic relatedness and virulence profiles. In addition, we identified numerous unique non-synonymous single-nucleotide polymorphisms (nsSNPs) in adhesins and effectors, exclusively associated with the clinical isolates tested, which may suggest a possible role of these mutations in the overall virulence of the clinical isolates.

The hops (Humulus lupulus) genome contains a mid-sized terpene synthase family that shows wide functional and allelic diversity.

BACKGROUND: Hops (Humulus lupulus L.) are a dioecious climbing perennial, with the dried mature "cones" (strobili) of the pistillate/female inflorescences being widely used as both a bittering agent and to enhance the flavour of beer. The glandular trichomes of the bract and bracteole flowering structures of the cones produce an abundance of secondary metabolites, such as terpenoids, bitter acids and prenylated phenolics depending on plant genetics, developmental stage and environment. More knowledge is required on the functional and allelic diversity of terpene synthase (TPS) genes responsible for the biosynthesis of volatile terpenes to assist in flavour-directed hop breeding. RESULTS: Major volatile terpene compounds were identified using gas chromatography-mass spectrometry (GC-MS) in the ripe cones of twenty-one hop cultivars grown in New Zealand. All cultivars produced the monoterpene ß-myrcene and the sesquiterpenes α-humulene and ß-caryophyllene, but the quantities varied broadly. Other terpenes were found in large quantities in only a smaller subset of cultivars, e.g. ß-farnesene (in seven cultivars) and α-pinene (in four). In four contrasting cultivars (Wakatu™, Wai-iti™, Nelson Sauvin™, and 'Nugget'), terpene production during cone development was investigated in detail, with concentrations of some of the major terpenes increasing up to 1000-fold during development and reaching maximal levels from 50-60 days after flowering. Utilising the published H. lupulus genome, 87 putative full-length and partial terpene synthase genes were identified. Alleles corresponding to seven TPS genes were amplified from ripe cone cDNA from multiple cultivars and subsequently functionally characterised by transient expression in planta. Alleles of the previously characterised HlSTS1 produced humulene/caryophyllene as the major terpenes. HlRLS alleles produced (R)-(-)-linalool, whilst alleles of two sesquiterpene synthase genes, HlAFS1 and HlAFS2 produced α-farnesene. Alleles of HlMTS1, HlMTS2 and HlTPS1 were inactive in all the hop cultivars studied. CONCLUSIONS: Alleles of four TPS genes were identified and shown to produce key aroma volatiles in ripe hop cones. Multiple expressed but inactive TPS alleles were also identified, suggesting that extensive loss-of-function has occurred during domestication and breeding of hops. Our results can be used to develop hop cultivars with novel/improved terpene profiles using marker-assisted breeding strategies to select for, or against, specific TPS alleles.

Multi-omics analysis reveals the mechanism of bHLH130 responding to low-nitrogen stress of apple rootstock.

Nitrogen is critical for plant growth and development. With the increase of nitrogen fertilizer application, nitrogen use efficiency decreases, resulting in wasted resources. In apple (Malus domestica) rootstocks, the potential molecular mechanism for improving nitrogen uptake efficiency to alleviate low-nitrogen stress remains unclear. We utilized multi-omics approaches to investigate the mechanism of nitrogen uptake in two apple rootstocks with different responses to nitrogen stress, Malus hupehensis and Malus sieversii. Under low-nitrogen stress, Malus sieversii showed higher efficiency in nitrogen uptake. Multi-omics analysis revealed substantial differences in the expression of genes involved in flavonoid and lignin synthesis pathways between the two materials, which were related to the corresponding metabolites. We discovered that basic helix-loop-helix 130 (bHLH130) transcription factor was highly negatively associated with the flavonoid biosynthetic pathway. bHLH130 may directly bind to the chalcone synthase gene (CHS) promoter and inhibit its expression. Overexpressing CHS increased flavonoid accumulation and nitrogen uptake. Inhibiting bHLH130 increased flavonoid biosynthesis while decreasing lignin accumulation, thus improving nitrogen uptake efficiency. These findings revealed the molecular mechanism by which bHLH130 regulates flavonoid and lignin biosyntheses in apple rootstocks under low-nitrogen stress.

Strigolactone regulates adventitious root formation via the MdSMXL7-MdWRKY6-MdBRC1 signaling cascade in apple.

Propagation through stem cuttings is a popular method worldwide for species such as fruit tree rootstocks and forest trees. Adventitious root (AR) formation from stem cuttings is crucial for effective and successful clonal propagation of apple rootstocks. Strigolactones (SLs) are newly identified hormones involved in AR formation. However, the regulatory mechanisms underpinning this process remain elusive. In the present study, weighted gene co-expression network analysis, as well as rooting assays using stable transgenic apple materials, revealed that MdBRC1 served as a key gene in the inhibition of AR formation by SLs. We have demonstrated that MdSMXL7 and MdWRKY6 synergistically regulated MdBRC1 expression, depending on the interactions of MdSMXL7 and MdWRKY6 at the protein level downstream of SLs as well as the direct promoter binding on MdBRC1 by MdWRKY6. Furthermore, biochemical studies and genetic analysis revealed that MdBRC1 inhibited AR formation by triggering the expression of MdGH3.1 in a transcriptional activation pathway. Finally, the present study not only proposes a component, MdWRKY6, that enables MdSMXL7 to regulate MdBRC1 during the process of SL-controlled AR formation in apple, but also provides prospective target genes to enhance AR formation capacity using CRISPR (i.e. clustered regularly interspaced short palindromic repeats) technology, particularly in woody plants.

GWAS provides new insights into the genetic mechanisms of phytochemicals production and red skin colour in apple.

Understanding the genetic architecture of apple phytochemicals, and their interplay with conventional selection traits, is critical for the development of new apple cultivars with enhanced health benefits. Apple accessions (n = 344) used for this genome-wide association study (GWAS) represented the wide diversity of metabolic profiles in the domesticated and wild Malus genepools. Fruit samples were phenotyped for 34 metabolites, including a stable vitamin C glycoside "ascorbic acid 2-ß-glucoside" (AA-2ßG), and the accessions were genotyped using the Apple 20 K SNP Array. Several fruit quality traits, including red skin over-colour (OCOL), were also assessed. Wild Malus accessions showed at least 2-fold higher average content of several metabolites (e.g. ascorbic acid, chlorogenic acid, phloridzin, and trilobatin) than Malus domestica accessions. Several new genomic regions and potential candidate genes underpinning the genetic diversity of apple phytochemicals were identified. The percentage of phenotypic variance explained by the best SNP ranged between 3% and 21% for the different metabolites. Novel association signals for OCOL in the syntenic regions on chromosomes 13 and 16 suggested that whole genome duplication has played a role in the evolution of apple red skin colour. Genetic correlations between phytochemicals and sensory traits were moderate. This study will assist in the selection of Malus accessions with specific phytochemical profiles to establish innovative genomics-based breeding strategies for the development of apple cultivars with enhanced nutritional value.

High-density linkage map construction in an autotetraploid blueberry population and detection of quantitative trait loci for anthocyanin content.

Highbush blueberry (Vaccinium corymbosum, 2n = 4x = 48) is the most cultivated type of blueberry, both in New Zealand and overseas. Its perceived nutritional value is conferred by phytonutrients, particularly anthocyanins. Identifying the genetic mechanisms that control the biosynthesis of these metabolites would enable faster development of cultivars with improved fruit qualities. Here, we used recently released tools for genetic mapping in autotetraploids to build a high-density linkage map in highbush blueberry and to detect quantitative trait loci (QTLs) for fruit anthocyanin content. Genotyping was performed by target sequencing, with ∼18,000 single nucleotide polymorphism (SNP) markers being mapped into 12 phased linkage groups (LGs). Fruits were harvested when ripe for two seasons and analyzed with high-performance liquid chromatography-mass spectrometry (HPLC-MS): 25 different anthocyanin compounds were identified and quantified. Two major QTLs that were stable across years were discovered, one on LG2 and one on LG4, and the underlying candidate genes were identified. Interestingly, the presence of anthocyanins containing acylated sugars appeared to be under strong genetic control. Information gained in this study will enable the design of molecular markers for marker-assisted selection and will help build a better understanding of the genetic control of anthocyanin biosynthesis in this crop.

Elevating fruit carotenoid content in apple (Malus x domestica Borkh).

Carotenoid compounds accumulate to confer coloration to plant tissues and have some established health benefits in humans. These pigments have antioxidant properties and are precursors of vitamin A, which is important for human vision. Apple is widely consumed globally, but most commercial apple cultivars have low fruit carotenoid content because these pigments accumulate mostly in the fruit skin rather than the flesh (the majority of the edible portion). Although carotenoids accumulate in the early stages of fruit development, much of this carotenoid is lost by fruit maturity as a result of low biosynthetic rate, rapid turnover of compounds and/or lack of storage capacity in these tissues. Improving apple fruit carotenoid content through traditional breeding or genetic technologies, will take a long time because of the extended juvenile phase of the trees and limited germplasm diversity within many commercial breeding programs. This process, however, can be accelerated by fundamental understanding of the apple carotenoid biosynthetic pathway and the mechanisms controlling the metabolic steps. The availability of a well annotated apple genome sequence has led to the identification of apple carotenoid gene families and potential transcription factors. This is an important step since the knowledge could be used to elevate carotenoid content either through breeding or genetic transformation techniques. Here, we provide an overview of carotenogenesis in apple and outline the methods employed to improve the carotenoid content of this horticultural crop.

The Coordinated Action of MYB Activators and Repressors Controls Proanthocyanidin and Anthocyanin Biosynthesis in Vaccinium.

Vaccinium berries are regarded as "superfoods" owing to their high concentrations of anthocyanins, flavonoid metabolites that provide pigmentation and positively affect human health. Anthocyanin localization differs between the fruit of cultivated highbush blueberry (V. corymbosum) and wild bilberry (V. myrtillus), with the latter having deep red flesh coloration. Analysis of comparative transcriptomics across a developmental series of blueberry and bilberry fruit skin and flesh identified candidate anthocyanin regulators responsible for this distinction. This included multiple activator and repressor transcription factors (TFs) that correlated strongly with anthocyanin production and had minimal expression in blueberry (non-pigmented) flesh. R2R3 MYB TFs appeared key to the presence and absence of anthocyanin-based pigmentation; MYBA1 and MYBPA1.1 co-activated the pathway while MYBC2.1 repressed it. Transient overexpression of MYBA1 in Nicotiana benthamiana strongly induced anthocyanins, but this was substantially reduced when co-infiltrated with MYBC2.1. Co-infiltration of MYBC2.1 with MYBA1 also reduced activation of DFR and UFGT, key anthocyanin biosynthesis genes, in promoter activation studies. We demonstrated that these TFs operate within a regulatory hierarchy where MYBA1 activated the promoters of MYBC2.1 and bHLH2. Stable overexpression of VcMYBA1 in blueberry elevated anthocyanin content in transgenic plants, indicating that MYBA1 is sufficient to upregulate the TF module and activate the pathway. Our findings identify TF activators and repressors that are hierarchically regulated by SG6 MYBA1, and fine-tune anthocyanin production in Vaccinium. The lack of this TF module in blueberry flesh results in an absence of anthocyanins.

Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.).

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

The genome of low-chill Chinese plum "Sanyueli" (Prunus salicina Lindl.) provides insights into the regulation of the chilling requirement of flower buds.

Chinese plum (Prunus salicina Lindl.) is a stone fruit that belongs to the Prunus genus and plays an important role in the global production of plum. In this study, we report the genome sequence of the Chinese plum "Sanyueli", which is known to have a low-chill requirement for flower bud break. The assembled genome size was 282.38 Mb, with a contig N50 of 1.37 Mb. Over 99% of the assembly was anchored to eight pseudochromosomes, with a scaffold N50 of 34.46 Mb. A total of 29,708 protein-coding genes were predicted from the genome and 46.85% (132.32 Mb) of the genome was annotated as repetitive sequence. Bud dormancy is influenced by chilling requirement in plum and partly controlled by DORMANCY ASSOCIATED MADS-box (DAM) genes. Six tandemly arrayed PsDAM genes were identified in the assembled genome. Sequence analysis of PsDAM6 in "Sanyueli" revealed the presence of large insertions in the intron and exon regions. Transcriptome analysis indicated that the expression of PsDAM6 in the dormant flower buds of "Sanyueli" was significantly lower than that in the dormant flower buds of the high chill requiring "Furongli" plum. In addition, PsDAM6 expression was repressed by chilling treatment. The genome sequence of "Sanyueli" plum provides a valuable resource for elucidating the molecular mechanisms responsible for the regulation of chilling requirements, and it is also useful for the identification of the genes involved in the control of other important agronomic traits and molecular breeding in plum.

Hierarchical regulation of MYBPA1 by anthocyanin- and proanthocyanidin-related MYB proteins is conserved in Vaccinium species.

Members of the Vaccinium genus bear fruits rich in anthocyanins, a class of red-purple flavonoid pigments that provide human health benefits, although the localization and concentrations of anthocyanins differ between species: blueberry (V. corymbosum) has white flesh, while bilberry (V. myrtillus) has red flesh. Comparative transcriptomics between blueberry and bilberry revealed that MYBPA1.1 and MYBA1 strongly correlated with the presence of anthocyanins, but were absent or weakly expressed in blueberry flesh. MYBPA1.1 had a biphasic expression profile, correlating with both proanthocyanidin biosynthesis early during fruit development and anthocyanin biosynthesis during berry ripening. MYBPA1.1 was unable to induce anthocyanin or proanthocyanidin accumulation in Nicotiana benthamiana, but activated promoters of flavonoid biosynthesis genes. The MYBPA1.1 promoter is directly activated by MYBA1 and MYBPA2 proteins, which regulate anthocyanins and proanthocyanidins, respectively. Our findings suggest that the lack of VcMYBA1 expression in blueberry flesh results in an absence of VcMYBPA1.1 expression, which are both required for anthocyanin regulation. In contrast, VmMYBA1 is well expressed in bilberry flesh, up-regulating VmMYBPA1.1, allowing coordinated regulation of flavonoid biosynthesis genes and anthocyanin accumulation. The hierarchal model described here for Vaccinium may also occur in a wider group of plants as a means to co-regulate different branches of the flavonoid pathway.