Genome-wide identification and characterisation of R2R3-MYB genes in sugar beet (Beta vulgaris)
© Stracke et al.; licensee BioMed Central Ltd. 2014
Received: 26 May 2014
Accepted: 17 September 2014
Published: 25 September 2014
The R2R3-MYB genes comprise one of the largest transcription factor gene families in plants, playing regulatory roles in plant-specific developmental processes, metabolite accumulation and defense responses. Although genome-wide analysis of this gene family has been carried out in some species, the R2R3-MYB genes in Beta vulgaris ssp. vulgaris (sugar beet) as the first sequenced member of the order Caryophyllales, have not been analysed heretofore.
We present a comprehensive, genome-wide analysis of the MYB genes from Beta vulgaris ssp. vulgaris (sugar beet) which is the first species of the order Caryophyllales with a sequenced genome. A total of 70 R2R3-MYB genes as well as genes encoding three other classes of MYB proteins containing multiple MYB repeats were identified and characterised with respect to structure and chromosomal organisation. Also, organ specific expression patterns were determined from RNA-seq data. The R2R3-MYB genes were functionally categorised which led to the identification of a sugar beet-specific clade with an atypical amino acid composition in the R3 domain, putatively encoding betalain regulators. The functional classification was verified by experimental confirmation of the prediction that the R2R3-MYB gene Bv_iogq encodes a flavonol regulator.
This study provides the first step towards cloning and functional dissection of the role of MYB transcription factor genes in the nutritionally and evolutionarily interesting species B. vulgaris. In addition, it describes the flavonol regulator BvMYB12, being the first sugar beet R2R3-MYB with an experimentally proven function.
Transcriptional control of gene expression influences almost all biological processes in eukaryotic cells or organisms. Transcription factors perform this function, alone or complexed with other proteins, by activating or repressing (or both) the recruitment of RNA polymerase to specific genes. The large number and diversity of transcription factors is related to their substantial regulatory complexity .
MYB proteins are widely distributed in all eukaryotic organisms and constitute one of the largest transcription factor families in the plant kingdom. MYB proteins are defined by a highly conserved MYB DNA-binding domain, mostly located at the N-terminus, generally consisting of up to four imperfect amino acid sequence repeats (R) of about 52 amino acids, each forming three alpha–helices . The second and third helices of each repeat build a helix–turn–helix (HTH) structure with three regularly spaced tryptophan (or hydrophobic) residues, forming a hydrophobic core . The third helix of each repeat is the DNA recognition helix that makes direct contact with DNA . During DNA contact, two MYB repeats are closely packed in the major groove, so that the two recognition helices bind cooperatively to the specific DNA recognition sequence motif.
MYB proteins can be divided into different classes depending on the number of adjacent repeats (one, two, three or four). The three repeats of the prototypic MYB protein c-Myb  are referred to as R1, R2 and R3, and repeats from other MYB proteins are named according to their similarity. Plant R1R2R3-type MYB (MYB3R) proteins have been proposed to play divergent roles in cell cycle control ,, similar to the functions of their animal homologs.
Most plant MYB genes encode R2R3-MYB class proteins, containing two repeats ,, which are thought to have evolved from an R1R2R3-MYB gene ancestor, by the loss of the sequences encoding the R1 repeat and subsequent expansion of the gene family -. R2R3-MYB transcription factors have a modular structure, with the N-terminal MYB domain as DNA-binding domain and an activation or repression domain usually located at the highly variable C-terminus. Components for the establishment of protein-protein interactions with other components of the eukaryotic transcriptional machinery have been detected in the N-terminal module -.
Based on the conservation of the MYB domain and of common amino acid motifs in the C-terminal domains, R2R3-MYB proteins have been divided into several subgroups which often group proteins with functional relationship. The reliability of the subgroups defined on the basis of phylogenetic analysis is also supported by additional criteria, such as the gene structure and the presence and position of introns . Most of these subgroups, defined first for the proteins of A. thaliana,,, are also present, and are sometimes expanded, in other higher plants. Comparative phylogenetic studies have identified new R2R3-MYB subgroups in other plant species for which there are no representatives in A. thaliana (e.g. in rice, poplar and grapevine), suggesting that these proteins might have specialised functions which were either lost in A. thaliana or were acquired after divergence from the last common ancestor -.
As initially described in the first plant MYB gene family review , the expansion of the plant-specific R2R3-MYB gene family is thought to be correlated with the increase in complexity of plants, particularly in Angiosperms. Consequently, the functions of R2R3-MYB genes are likely associated with regulating plant-specific processes including primary and secondary metabolism, developmental processes, cell fate and identity and responses to biotic and abiotic stresses ,,.
With the growing number of fully sequenced plant genomes, the identification of R2R3-MYB genes has increased in recent times. Based on their well conserved MYB domains, R2R3-MYB gene families have been annotated genome-wide in A. thaliana (126 members) , Zea mays (157 members) , Oryza sativa (102 members) , Vitis vinifera (117 members) , Populus trichocarpa (192 members) , Glycine max (244 members) , Cucumis sativus (55 members)  and Malus x domestica (222 members) . Given the potential roles of R2R3-MYB proteins in the regulation of gene expression, secondary metabolism, and responses to environmental stresses, and that Beta vulgaris ssp. vulgaris (order Caryophyllales) is the first non-rosid, non-asterid eudicot for which the genome has been sequenced , it is of interest to achieve a complete identification and classification of MYB genes in this species with respect to the number, chromosome locations, phylogenetic relationships, conserved motifs as well as expression patterns. Particularly, since sugar beet is an important crop of the temperate climates as a source for bioethanol as well as animal feed and provides nearly 30% of the worlds annual sugar production .
In the present study, we describe the R2R3-MYB gene family by means of in silico analysis of the B. vulgaris genome sequence, in order to predict protein domain architectures, and to assess the extent of conservation and divergence between B. vulgaris and A. thaliana gene families, thus leading to a functional classification of the sugar beet MYB genes on the basis of phylogenetic analyses. Furthermore, RNA-seq data was used to analyse expression in different B. vulgaris organs and to compare expression patterns of closely grouped co-orthologs. To validate the functional classification, a candidate gene was chosen for cDNA isolation and subsequent functional analysis by transient transactivation assays and complementation of an orthologous A. thaliana mutant. We identified the R2R3-MYB gene Bv_iogq activating two flavonol biosynthesis enzyme promoters and complementing the flavonol-deficient myb11 myb12 myb111 mutant, and thus encoding a functional flavonol biosynthesis regulator. Our findings provide the first step towards further investigations on the biological and molecular functions of MYB transcription factors in the economically and evolutionarily interesting species B. vulgaris.
Results and discussion
The annotated genome sequence of B. vulgaris has recently become available. It has been obtained from the double haploid breeding line KWS2320 . The sequence has been assigned to nine chromosomes and B. vulgaris was predicted to contain 27,421 protein-coding genes (RefBeet) in 567 Mb from which 85% are chromosomally assigned.
Identification and genomic distribution of B. vulgaris R2R3-MYB genes
List of annotated MYB genes with two or more repeats in the B. vulgaris ssp. vulgaris (KWS2320) genome
Position on pseudochr.
Landmark MYB in clade
Protein length [aa]
The total number of identified MYB genes was, compared to other plant species, low in B. vulgaris. Even if some MYB genes may have been missed due to gaps in the reference sequence, this does not adequately explain the small number. High numbers of MYB genes in a species are mainly attributed to ancestral whole genome duplication events as known for A. thaliana, O. sativa, P. trichocarpa, G. max and M. x domestica-. The absence of a recent lineage-specific whole genome duplication event in B. vulgaris is further substantiated by the detection of only 70 R2R3-MYB genes, because a lack of this duplication event can easily explain the small number of MYB genes in this species. This interpretation is in accordance with the findings in cucumber, where the number of R2R3-MYB genes has been reported to be 55 .
We further determined physically linked sister BvMYB gene pairs along the nine chromosomes (Figure 1, marked with vertical black bars), which form clusters and may have evolved from local intrachromosomal duplication events that result in tandem arrangement of the duplicated gene. Three gene pairs have been identified: one on chromosome 7 consisting of the closely related genes Bv_ahzu and Bv_eztu, a second on the top of chromosome 8 with Bv_dxny and Bv_zguf and a third on an unlinked scaffold (0254.scaffold00675) constituted of Bv_dani and Bv_sjwa. The BvMYB genes of the two latter pairs were physically located near to each other without intervening annotated genes between. In total, about 5% (6 of 75) of BvMYBs were involved in tandem duplication, which is the same value as reported for MYB genes in soybean . Moreover, an incomplete gene pair was observed on chromosome 2, where a solitary typically R2R3-MYB "third exon" containing sequences encoding a part of a R3 repeat and the C-terminal region was found about 18.7 kb downstream of Bv_jkkr showing 88% identity on cDNA- and 82% identity on deduced protein level to the third exon of the near Bv_jkkr.
Gene structure analysis revealed that most BvR2R3-MYB genes (53 of 70, 76%) follow the previously reported rule of having two introns and three exons, and display the highly conserved splicing arrangement that has also been reported for other plant species ,,. Eleven BvR2R3-MYB genes (16%) have one intron and two exons and four (6%) were one exon genes. Only two BvR2R3-MYB genes have more than three exons: Bv_mxck with four exons and Bv_zeqy with twelve exons (Table 1). The complex exon-intron structure of Bv_zeqy is known from its A. thaliana orthologs AtMYB88 and AtMYB124/FOUR LIPS (FLP) containing ten and eleven exons, respectively, and more than the typically zero to two introns in the MYB domain coding sequences ,. This supports their close evolutionary relationship, but also indicates the conservation of this intron pattern in evolution since the split of the Caryophyllales from the precursor of rosids and asterids.
Sequence features of the MYB domains
In addition to the highly conserved tryptophan residues, we observed amino acid residues that are conserved in all B. vulgaris R2R3-MYBs: D11, C43, R46 in the R2 repeat and E64, G76, R89 and T90 in the R3 repeat. Further highly conserved residues of the B. vulgaris R2R3-MYB domains are: G4, E10, L14, G35, R38, K41, R44, N49, L51 and P53 in the R2 repeat and I82, A83, N92, K95 and N96 in the R3 repeat (Figure 2). These highly conserved amino acid residues are mainly located in the third helix and the turn of the helix-turn-helix (HTH) motif, which is in good accordance with the findings in other plant species.
Phylogenetic analysis of the B. vulgaris MYB family
To explore the putative function of the predicted B. vulgaris MYBs, we assigned them to functional clades known from A. thaliana, which was chosen because most of our knowledge about plant MYB genes has been obtained from studies of this major plant model. As known from similar studies, most MYB proteins sharing similar functions cluster in the same phylogenetic clades, suggesting that most closely-related MYBs could recognise similar target genes and possess redundant, overlapping, and/or cooperative functions.
We also observed species-specific clades and clades containing B. vulgaris or A. thaliana MYB proteins together with landmark MYBs from other species. It should be noted that we use the term "species-specific" in the context of the current set of species with sequenced genomes. Significantly more genome sequences would be required to resolve the presence and absence of genes or clades at the genus or family level. As indicated also from other studies, the observation of species-specific clades may be taken as a hint that MYB genes may have been acquired or been lost in a single species, during the following divergence from the most recent common ancestor. For example, members of the clade C2 (subgroup S12, with HIGH ALIPHATIC GLUCOSINOLATE1 (AtHAG1), HIGH INDOLIC GLUCOSINOLATE1 (AtHIG1), ALTERED TRYPTOPHAN REGULATION1 (AtATR1)) have been identified as glucosinolate biosynthesis regulators -. No BvMYBs were grouped within this clade, containing members which are predominantly present in plants of the glucosinolate compounds accumulating Brassicaceae family. A previous study indicated that this clade was derived from a duplication event before Arabidopsis diverged from Brassica . Another MYB clade without B. vulgaris orthologs was C20 (S15), that include the landmark R2R3-MYBs AtMYB0/GLABRA1 and AtMYB66/WEREWOLF (Figure 3), which are known to be involved in epidermal cell development leading to the formation of trichomes and root hairs -. Similar observations have been made in the non-rosids maize and soybean ,, both not containing C20 (S15) orthologs, while the rosids grape and poplar do ,. As GLABRA1-like MYB genes have been hypothesised to have been acquired in rosids after the rosid-asterids division ,, the absence of BvMYBs in this clade is consistent with this hypothesis, since Caryophyllales branched off before the separation of asterids and rosids. Trichome formation in B. vulgaris thus could be regulated by genes of the evolutionary older MIXTA clade C1 (S9)  whose members are also known to play a role in multicellular trichome formation . Therefore, Bv_dcmm as the only sugar beet gene in this clade, is a candidate to encode a trichome development regulator. Further clades without BvMYBs are C5 (S10) containing the A. thaliana proteins AtMYB9, AtMYB107 with unknown functional assignment, C16 (S5) containing R2R3-MYB landmark anthocyanin regulators from monocots and C22 (S6) including the landmark R2R3-MYB factors PRODUCTION OF ANTHOCYANIN PIGMENT1 (AtPAP1), ANTHOCYANIN2 (PhAN2) and ROSEA1 (AmROSEA1), which are known to regulate anthocyanin biosynthesis in many species -. The lack of BvMYBs in the two latter clades fits to the observation, that plants of the genus Beta do not produce anthocyanin pigments . The Caryophyllales is the single order in the plant kingdom that contains taxa that have replaced anthocyanins with chemically distinct but functionally identical red and yellow pigments - the indole-derived betalains, named from Beta vulgaris, from which betalains were first extracted . Although betalains have functions analogous to those of anthocyanins as pigments, anthocyanins and betalains are mutually exclusive pigments in plants .
Three clades do not contain any A. thaliana MYB (Figure 3). Two of them contain landmark MYBs and B. vulgaris MYBs: C18 (S5) and C15, functionally assigned to proanthocyanidin regulation and repression of flavonoid biosynthesis, respectively. One clade contains only BvMYBs. Clade C21, constituted of the two BvMYBs (Bv_ralf and Bv_jkkr), was found close to the "anthocyanin" regulator representing clade C22 (S6) in the phylogenetic trees (Figure 3, Additional file 3). This clade could be described as a lineage-specific expansion in B. vulgaris, reflecting a species-specific adaptation. Two classical, linked beet pigment loci, RED (R) and YELLOW (Y) ,, are known to influence the production of betalains in beet. Recently, the R locus was shown to encode a cytochrome P450 (CYP76AD1, Bv_ucyh in KWS2320) . It has been hypothesised that the betalain pathway may have co-opted the anthocyanin regulators because both pigment types are produced in a similar temporal and spatial pattern . Thus, a R2R3-MYB-type transcriptional activator, homologues to the anthocyanin pigment regulators, is thought to control betalain biosynthesis through regulation of the biosynthetic enzymes R (CYP76AD1) and DODA (4,5-DOPA-dioxygenase) ,. The two C21 BvMYB genes Bv_jkkr and Bv_ralf are both located on chromosome 2 close to the R locus, as indicated by the chromosome-based gene designations (Bv2g027795_jkkr, Bv2g029890_ucyh, Bv2g030925_ralf). As the R and Y loci are known to be linked at a genetically distance of about 7.5 cM ,, this makes the R2R3-MYB genes Bv_ralf and Bv_jkkr candidates for encoding the Y regulator. A closer inspection of the MYB domains of Bv_ralf and Bv_jkkr, and comparison to anthocyanin regulator MYBs, revealed some interesting features of C21 clade MYBs which could cause the separation of C21 and C22 MYBs (Figure 3). C21 MYBs do not contain the bHLH-binding consensus motif [D/E]Lx2[R/K]x3Lx6Lx3R  found in all bHLH-interacting R2R3-MYBs, suggesting that Bv_ralf and Bv_jkkr, in contrast to the C22 anthocyanin regulator MYBs, do not interact with bHLH proteins. A key amino acid residue in the R3 repeat, identified by Heppel et al. in separating anthocyanin regulators (A89) from proanthocyanidin regulators (G89), is represented by isoleucine (I89) in C21 MYBs. Furthermore, one of the amino acid residues in the third helix of R2 repeat, known to be directly involved in DNA binding (position 44), differs from those found at this position in anthocyanin regulators (Additional file 1), maybe indicating a different target promoter specificity. In this direction also the atypical cysteine at the highly conserved position 98 in the MYB domain of Bv_jkkr may be important as discussed below.
A motif search in the BvMYB proteins for the above mentioned bHLH-interaction motif  identified seven BvMYBs containing this motif and thus putatively interacting with bHLH proteins. These seven BvMYBs were all functionally assigned to clades containing (potentially) known bHLH-interacting R2R3-MYBs: Bv_crae, Bv_swwi and Bv_dwki were assigned to the proanthocyanidin regulator clade C17 (S5), Bv_ihfg to clade C15 containing the negative flavonoid regulator FaMYB1 , Bv_cjuq to clade C19 containing the general phenylpropanoid pathway regulator VvMYB5a  and Bv_gjwr and Bv_iquc to the "repressors" clade C14 (S4).
Six B. vulgaris MYB proteins, Bv_uksi, Bv_josh, Bv_cfqe Bv_ijmc, Bv_swwi and Bv_pgya did not cluster in any of the identified clades or subgroups or showed ambiguous placements between the different phylogenetic trees, implying that these BvMYB proteins may have specialised roles that were acquired or expanded in B. vulgaris during the process of genome evolution.
Basically, higher numbers of AtMYBs than their BvMYB orthologs were found in most clades, indicating that they were duplicated after the divergence of the two lineages. For example, the phylogeny for C28 (S18) and C36 (S21) included only three BvMYBs and eight AtMYBs, or C25 (S13) with one BvMYB and four AtMYBs. By contrast, three BvMYBs and two AtMYB were found in clade C26. As mentioned above, the higher number of MYB genes in A. thaliana is presumably mainly due to an ancestral duplication of the entire genome and subsequent rearrangements, followed by gene loss and extensive local gene duplications taken place in the evolution of Arabidopsis .
Expression profiles for B. vulgaris MYB genes in different developmental stages
In large transcription factor families functional redundancy is not unusual. Thus, a particular transcription factor has often to be studied and characterised in the context of the whole family. In this context, the gene expression pattern can provide important clues for gene function. We used genome-mapped global RNA-seq Illumina reads from the B. vulgaris reference genotype KWS2320  to analyse the expression of the 75 BvMYB genes in different organs and developmental stages: seedlings, taproot (conical white fleshy root), young and old leaves, inflorescences and seeds. Filtered RNA-seq reads were aligned to the genome reference sequence and the number of mapped reads per annotated transcript were quantified and statistically compared across the analysed samples giving normalised RNA-seq read values which are given in Additional file 3. In common with other transcription factor genes, many of the BvMYBs exhibited low transcript abundance levels, as determined by the RNA-seq analysis.
Some clustered BvMYB genes (Figure 1), which are often considered as paralogs (created by a duplication event within the genome), showed similar expression profiles, while other clustered BvMYB genes did not. Bv_ahzs and Bv_eztu (both in the development-related clade 26), are mainly expressed in inflorescences and seeds with only low (or no) expression in the other analysed organs. The two clustered genes Bv_dxny and Bv_zguf (both in the defence-related clade 7) showed different expression profiles with Bv_zguf mainly expressed in seedlings and seeds and Bv_dxny further expressed in leaves. These results indicate, that Bv_ahzs and Bv_eztu could be functionally redundant genes, while Bv_dxny and Bv_zguf could be (partly) involved in various aspects of defence processes.
Functional analysis of the R2R3-MYB gene Bv_iogq
In the context of functional analysis we selected Bv_iogq to start with. Bv_iogq is the only BvMYB in clade C13 (S7) clustering together with landmark MYBs having been implicated as flavonol- and phlobaphene-specific activators of flavonoid biosynthesis during plant development. This especially includes the Z. mays factor ZmP, controlling the accumulation of red phlobaphene pigments in pericarps , the V. vinifera flavonol regulator VvMYBF1  and the A. thaliana PRODUCTION OF FLAVONOL GLYCOSIDES (PFG) family constituted of AtMYB12/AtPFG1, AtMYB11/AtPFG2 and AtMYB111/AtPFG3 . AtMYB12 was shown to regulate CHALCONE SYNTHASE (CHS) expression by binding the MYB-recognition element MRE CHS . The three A. thaliana PFG-family factors were further found to have overlapping functions in flavonol pathway-specific gene regulation but to play organ-specific roles. Concordantly, A. thaliana seedlings of the pfg triple mutant (myb11 myb12 myb111) do not form flavonols under standard growth conditions, whereas the accumulation of other flavonoids is not affected .
We tested if Bv_iogq codes for a functional flavonol regulator. The Bv_iogq open reading frame (ORF) encodes a predicted protein with 387 amino acid residues, with a calculated molecular mass of 43.6 kD. Analysis of the deduced amino acid sequence revealed that Bv_iogq contains both motifs which have been described to be characteristic for flavonol regulators. The subgroup S7 motif (GRTxRSxMK)  was found to be present in the C-terminus of Bv_iogq with one amino acid substitution (GRTSRWAMQ), and the SG7-2 motif ([W/x][L/x]LS; ) was found at the very C-terminal end of Bv_iogq.
Since flavonol glycosides are known to accumulate in B. vulgaris leaves ,, and our expression data indicated Bv_iogq to be expressed in young leaves (Figure 4), this organ was chosen as source to clone the cDNA of the putative flavonol biosynthesis regulator Bv_iogq.
Bv_iogq activates promoters of the flavonoid pathway structural genes
Complementation of the flavonol-deficient phenotype of A. thaliana multiple pfg mutant
The genome-wide identification, chromosomal organisation, functional classification and expression analyses of B. vulgaris MYB genes provide an overall insight into this transcription factor-encoding gene family and their potential involvement in growth and development processes. This study provides the first step towards cloning and functional dissection to uncover the role of MYB genes in this economically and evolutionarily interesting species, as representative for the order Caryophyllales. The functional classification was successfully verified by experimental confirmation of the prediction that the R2R3-MYB gene Bv_iogq encodes a flavonol regulator, thus making BvMYB12 the first sugar beet R2R3-MYB with an experimentally proven function.
Database search for MYB protein coding genes in the B. vulgaris genome
An initial search for MYB protein coding genes in B. vulgaris was performed using a consensus R2R3-MYB DNA binding domain sequence  as protein query in TBLASTN  searches on the beet genome sequence (RefBeet, http://www.genomforschung.uni-bielefeld.de/en/projects/annobeet). To confirm the obtained amino acid sequences, the putative MYB sequences were manually analysed for the presence of an intact MYB domain. All B. vulgaris MYB proteins were inspected to ensure that the putative gene models contained two or more (multiple) MYB repeats, and that the gene models mapped to unique loci in the genome. Redundant sequences were discarded from the data set to obtain unique BvMYB genes. The identified BvMYB genes were matched with the automatically annotated MYB genes from the study of Dohm et al. (). The open reading frames of the identified BvMYB have been verified by mapping, visualisation and manually inspection of combined RNA-seq data. Additionally, manually identified BvMYB genes received a unique immutable four-letter-ID and a gene code describing the position on pseudochromosomes or unanchored scaffolds. Multiple FASTA files with protein sequences used in this work, CDS and CDSi (CDS plus introns) sequences of BvMYBs are found as Additional files 4, 5 and 6.
Sequence conservation analysis
To analyse the sequence features of the MYB domain of B. vulgaris R2R3-MYB proteins, the sequences of R2 and R3 MYB repeats of 70 BvR2R3-MYB proteins were aligned with ClustalOmega (1.2.0) (http://www.ebi.ac.uk/Tools/msa/clustalo/) using the default parameters. The sequence logos for R2 and R3 MYB repeats were produced from the multiple alignment files by the web based application WebLogo (http://weblogo.berkeley.edu/logo.cgi)  using default settings.
133 A. thaliana MYB protein sequences were obtained from TAIR (http://www.arabidopsis.org/). 51 well-known landmark plant R2R3-MYB protein sequences were collected from GenBank at the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/). We additionally considered atypical multiple MYB-repeat proteins from A. thaliana in the phylogenetic analysis to determine orthologs in the B. vulgaris genome: AtCDC5, AtMYB4R1 and five AtMYB3R. Phylogenetic trees were constructed from ClustalW aligned full-length MYB proteins (75 BvMYBs, 133 AtMYBs and 51 plant landmark MYBs) using MEGA5.2  with full-length proteins and default settings. Two statistical methods were used for tree generation: (i) neighbor-joining (NJ) method with poison correction and bootstrap analysis with 1000 replicates, (ii) maximum parsimony (MP) method with 1000 bootstrap replications. Classification of the B. vulgaris MYBs was performed according to their phylogenetic relationships with their corresponding A. thaliana and landmark MYB proteins.
Expression analysis from RNA-sequencing data
RNA-seq raw data from inflorescences, taproots, old leaves, seeds and seedlings of the B. vulgaris double haploid reference line KWS2320 have been produced in the work of Dohm et al. and are deposited at the NCBI Short Read Archive (SRA) under the accession numbers SRX287608-SRX287615. RNA-seq raw data from RNA isolated from young leaves (third and forth leaf without midrip) have been generated by paired-end sequencing (2× 101 bases) on an Illumina HiSeq1500 instrument and are deposited together with a detailed documentation at the SRA under the accession number SRX647324. Raw RNA-seq reads were trimmed using Trimmomatic-0.22  and the options ILLUMINACLIP, LEADING:3, TRAILING:3, SLIDINGWINDOW:4:15 as well as MINLEN:36. For Tophat  mappings versus the RefBeet sequence, only trimmed paired reads generated by Trimmomatic-0.22 were used (Additional file 7). For mapping of long read sets (>55 bases read length) and short read sets (≤55 bases read length) the mate inner distance was set to 175 and 150, respectively and the option solexa1.3-quals was switched on. Mapping results are displayed in Additional file 7. Resulting BAM files were quality filtered according a minimum mapping quality of five and converted using SAMtools . SAM files were sorted afterwards (Unix command sort -s -k 1,1). Reads per gene were counted based on manually corrected gene structures stored in gff3 format and the sorted SAM files by htseq-count (option t = CDS) that is part of HTSeq version 0.5.4 (http://www-huber.embl.de/users/anders/HTSeq/doc/overview.html). Single count tables for each organ created by HTSeq were joined using a customised Perl script and subsequently imported in R by means of the DESeq package . Data design, sample assignment and normalisation were performed as described previously . The correction factors of different data sets were computed as 4.02, 1.25, 3.53, 0.11, 1.47, 0.41 for seedling, taproot, young leaf, old leaf, inflorescence and seed, respectively.
Hierachical clustering and heatmap visualisation of RNA-seq data
Centroid-linkage hierachical clustering of log2 transformed normalised RNA-seq data was performed using Cluster 3.0 software  (http://bonsai.hgc.jp/~mdehoon/software/cluster/software.htm). The hierachical clustering results were visualised using Java TreeView  (http://jtreeview.sourceforge.net).
Plant material and growth conditions
B. vulgaris KWS2320 seeds were surface-sterilised with 70% ethanol and sown on sterile, wet cotton wool in a covered plastic beaker. Seedlings were grown for one week at 22°C giving 16 hours light per day. All other plant material was harvested from greenhouse/soil-grown B. vulgaris plants.
RNA isolation and cDNA synthesis
RNA was isolated from 50 mg plant material using the NucleoSpin ® RNA Plant kit (Macherey-Nagel) according to the manufacturers instruction. cDNA was generated from 1 μg total RNA as described previously .
A cDNA fragment corresponding to the full-length ORF of Bv_iogq (1164 bp) was obtained by PCR using attB recombination sites containing sense RS1050 (5'-attB1-ccATGGGGAGAGCACCGTGTTGCGAAAAG-3') and antisense RS1051 (5'-attB2-gTTATGAAAGTAGCCAATCAAGCATAGC-3') primers with Phusion® polymerase (New England BioLabs) on cDNAs prepared from B. vulgaris KWS2320 young leaves as template. The resulting PCR product was introduced into the GATEWAY® vector pDONRzeo using BP clonase (Invitrogen), giving a Bv_iogq entry clone (RSt810). The Bv_iogq cDNA sequence was submitted to EMBL/GenBank; accession number KJ707238.
Complementation experiments with Bv_iogq
The Bv_iogq entry clone was recombined by LR-reaction (Invitrogen) into the binary destination vector pLEELA  giving a 2×35S:Bv_iogq construct (RSt815). The resulting plasmid was used in a complementation experiment transforming the flavonol-deficient A. thaliana triple mutant myb11 myb12 myb111 (NASC N9815) via Agrobacterium tumefaciens (GV3101::pM90RK)  according to the floral dip protocol. Flavonol-containing methanolic extracts of T2 seedlings of different transgenic lines were analysed by high performance thin layer chromatography (HPTLC) and diphenylboric acid 2-aminoethyl ester (DPBA)-staining as described elsewhere .
Transient protoplast assays
Growth of A. thaliana At7 cell culture, protoplast isolation and transfection experiments for the detection of transient expression were performed as described by Mehrtens et al.. In the co-transfection experiments, a total of 25 μg of pre-mixed plasmid DNA was transfected, consisting of 10 μg of reporter plasmid (transcriptional GUS fusion), 10 μg of effector plasmid and 5 μg of the luciferase (LUC) transfection control (4xproUBI:LUC) . Protoplasts were incubated for 20 h at 26°C in the dark before LUC and GUS enzyme activities were measured. Specific GUS activity is given in pmol 4-methylumbelliferone (MU) per mg protein per min. The error bars indicate the standard deviation of four GUS values determined for the respective reporter/effector combination. The 35S:Bv_iogq effector plasmid (RSt816) was generated by LR-reaction (Invitrogen) of the Bv_iogq Entry clone (RSt810) with the destination vector pBTdest . A 1064 bp promoter fragment of proBvCHS_phae was obtained by PCR on genomic B. vulgaris KWS2320 DNA using attB recombination sites containing primers BP60 (5'-attB1-CCACAAGAACATCTTGTAAGAGC-3') and BP55 (5'-attB2-CAACTCGTGAGTATGAAGAAATATG-3'). The 1137 bp promoter fragment of proBvFLS_ckiq was gained with the attB primers RS1128 (5'-attB1-TATACCTTGTAATATCACTTAAATCT-3') and RS1129 (5'-attB2-CCATTGAATATCAACTCAGGATTTTC-3'). The resulting PCR products were introduced via BP and LR reactions into the GATEWAY® destination vector pDISCO  giving the proBvCHS_phae-GUS (DH33) and proBvFLS_ckiq-GUS (RSt830) reporter constructs. The other effector- and reporter constructs have been described elsewere ,,.
Availability of supporting data
Phylogenetic data (trees and the data used to generate them) have been deposited in TreeBASE respository and is available under the URL http://purl.org/phylo/treebase/phylows/study/TB2:S16360.
RNA-seq data from inflorescences, taproots, young leaves, old leaves, seeds and seedlings of the B. vulgaris line KWS2320 are deposited at the NCBI Short Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra/) under the accession numbers SRX287608-SRX287615 and SRX647324.
We are grateful to Melanie Kuhlmann for excellent technical assistance. We thank Sebastian Packheiser who contributed to cloning and plant transformation. This work was supported by the Ministry of Education and Science (BMBF) grant “AnnoBeet: Annotation des Genoms der Zuckerrübe unter Berücksichtigung von Genfunktionen und struktureller Variabilität für Nutzung von Genomdaten in der Pflanzenbiotechnologie" (FKZ 0315962 A). A further "thanks" to our industrial partners KWS Saat AG and Syngenta Seeds GmbH. We acknowledge support for the Article Processing Charge by the German Research Foundation and the Open Access Publication Fund of Bielefeld University Library.
- Riechmann JL, Heard J, Martin G, Reuber L, Jiang CZ, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, Creelman R, Pilgrim M, Broun P, Zhang JZ, Ghandehari D, Sherman BK, Yu CL: Arabidopsis transcription factors: Genome-wide comparative analysis among eukaryotes. Science. 2000, 290 (5499): 2105-2110. 10.1126/science.290.5499.2105.View ArticlePubMedGoogle Scholar
- Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L: MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15: 573-581. 10.1016/j.tplants.2010.06.005.View ArticlePubMedGoogle Scholar
- Ogata K, Kanei-Ishii C, Sasaki M, Hatanaka H, Nagadoi A, Enari M, Nakamura H, Nishimura Y, Ishii S, Sarai A: The cavity in the hydrophobic core of Myb DNA-binding domain is reserved for DNA recognition and trans-activation. Nat Struct Biol. 1996, 3 (2): 178-187. 10.1038/nsb0296-178.View ArticlePubMedGoogle Scholar
- Jia L, Clegg MT, Jiang T: Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes. Plant Physiol. 2004, 134 (2): 575-585. 10.1104/pp.103.027201.PubMed CentralView ArticlePubMedGoogle Scholar
- Klempnauer KH, Gonda TJ, Bishop JM: Nucleotide sequence of the retroviral leukemia gene v-myb and its cellular progenitor c-myb: the architecture of a transduced oncogene. Cell. 1982, 31: 453-463. 10.1016/0092-8674(82)90138-6.View ArticlePubMedGoogle Scholar
- Ito M: Conservation and diversification of three-repeat Myb transcription factors in plants. J Plant Res. 2005, 118 (1): 61-69. 10.1007/s10265-005-0192-8.View ArticlePubMedGoogle Scholar
- Haga N, Kato K, Murase M, Araki S, Kubo M, Demura T, Suzuki K, Muller I, Voss U, Jurgens G, Ito M: R1R2R3-Myb proteins positively regulate cytokinesis through activation of KNOLLE transcription in Arabidopsis thaliana. Development. 2007, 134 (6): 1101-1110. 10.1242/dev.02801.View ArticlePubMedGoogle Scholar
- Jin H, Martin C: Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol. 1999, 41 (5): 577-585. 10.1023/A:1006319732410.View ArticlePubMedGoogle Scholar
- Rosinski JA, Atchley WR: Molecular evolution of the Myb family of transcription factors: evidence for polyphyletic origin. J Mol Evol. 1998, 46 (1): 74-83. 10.1007/PL00006285.View ArticlePubMedGoogle Scholar
- Braun EL, Grotewold E: Newly discovered plant c-myb-like genes rewrite the evolution of the plant myb gene family. Plant Physiol. 1999, 121: 21-24. 10.1104/pp.121.1.21.PubMed CentralView ArticlePubMedGoogle Scholar
- Kranz H, Scholz K, Weisshaar B: c-MYB oncogene-like genes encoding three MYB repeats occur in all major plant lineages. Plant J. 2000, 21 (2): 231-235. 10.1046/j.1365-313x.2000.00666.x.View ArticlePubMedGoogle Scholar
- Dias AP, Braun EL, McMullen MD, Grotewold E: Recently duplicated maize R2R3 Myb genes provide evidence for distinct mechanisms of evolutionary divergence after duplication. Plant Physiol. 2003, 131 (2): 610-620. 10.1104/pp.012047.PubMed CentralView ArticlePubMedGoogle Scholar
- Zimmermann IM, Heim MA, Weisshaar B, Uhrig JF: Comprehensive identification of Arabidopsis thaliana MYB transcription factors interacting with R/B-like BHLH proteins. Plant J. 2004, 40: 22-34. 10.1111/j.1365-313X.2004.02183.x.View ArticlePubMedGoogle Scholar
- Song S, Qi T, Huang H, Ren Q, Wu D, Chang C, Peng W, Liu Y, Peng J, Xie D: The Jasmonate-ZIM domain proteins interact with the R2R3-MYB transcription factors MYB21 and MYB24 to affect Jasmonate-regulated stamen development in Arabidopsis. Plant Cell. 2011, 23 (3): 1000-1013. 10.1105/tpc.111.083089.PubMed CentralView ArticlePubMedGoogle Scholar
- Du H, Yang SS, Liang Z, Feng BR, Liu L, Huang YB, Tang YX: Genome-wide analysis of the MYB transcription factor superfamily in soybean. BMC Plant Biol. 2012, 12 (106): 106-10.1186/1471-2229-12-106.PubMed CentralView ArticlePubMedGoogle Scholar
- Kranz HD, Denekamp M, Greco R, Jin H, Leyva A, Meissner RC, Petroni K, Urzainqui A, Bevan M, Martin C, Smeekens S, Tonelli C, Paz-Ares J, Weisshaar B: Towards functional characterisation of the members of the R2R3-MYB gene family from Arabidopsis thaliana. Plant J. 1998, 16: 263-276. 10.1046/j.1365-313x.1998.00278.x.View ArticlePubMedGoogle Scholar
- Stracke R, Werber M, Weisshaar B: The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001, 4: 447-456. 10.1016/S1369-5266(00)00199-0.View ArticlePubMedGoogle Scholar
- Jiang C, Gu X, Peterson T: Identification of conserved gene structures and carboxy-terminal motifs in the Myb gene family of Arabidopsis and Oryza sativa L. ssp. indica. Genome Biol. 2004, 5 (7): R46-10.1186/gb-2004-5-7-r46.PubMed CentralView ArticlePubMedGoogle Scholar
- Matus JT, Aquea F, Arce-Johnson P: Analysis of the grape MYB R2R3 subfamily reveals expanded wine quality-related clades and conserved gene structure organization across Vitis and Arabidopsis genomes. BMC Plant Biol 2008, 8(83).Google Scholar
- Wilkins O, Nahal H, Foong J, Provart NJ, Campbell MM: Expansion and diversification of the Populus R2R3-MYB family of transcription factors. Plant Physiol. 2009, 149 (2): 981-993. 10.1104/pp.108.132795.PubMed CentralView ArticlePubMedGoogle Scholar
- Martin C, Paz-Ares J: MYB transcription factors in plants. Trends Genet. 1997, 13 (2): 67-73. 10.1016/S0168-9525(96)10049-4.View ArticlePubMedGoogle Scholar
- Du H, Feng BR, Yang SS, Huang YB, Tang YX: The R2R3-MYB transcription factor gene family in maize. PLoS One. 2012, 7 (6): e37463-10.1371/journal.pone.0037463.PubMed CentralView ArticlePubMedGoogle Scholar
- Yanhui C, Xiaoyuan Y, Kun H, Meihua L, Jigang L, Zhaofeng G, Zhiqiang L, Yunfei Z, Xiaoxiao W, Xiaoming Q, Yunping S, Li Z, Xiaohui D, Jingchu L, Xing-Wang D, Zhangliang C, Hongya G, Li-Jia Q: The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006, 60 (1): 107-124. 10.1007/s11103-005-2910-y.View ArticlePubMedGoogle Scholar
- Li Q, Zhang C, Li J, Wang L, Ren Z: Genome-wide identification and characterization of R2R3MYB family in Cucumis sativus. PLoS One. 2012, 7 (10): e47576-10.1371/journal.pone.0047576.PubMed CentralView ArticlePubMedGoogle Scholar
- Cao ZH, Zhang SZ, Wang RK, Zhang RF, Hao YJ: Genome wide analysis of the apple MYB transcription factor family allows the identification of MdoMYB121 gene confering abiotic stress tolerance in plants. PLoS One. 2013, 8 (7): e69955-10.1371/journal.pone.0069955.PubMed CentralView ArticlePubMedGoogle Scholar
- Dohm JC, Minoche AE, Holtgrawe D, Capella-Gutierrez S, Zakrzewski F, Tafer H, Rupp O, Sorensen TR, Stracke R, Reinhardt R, Goesmann A, Kraft T, Schulz B, Stadler PF, Schmidt T, Gabaldon T, Lehrach H, Weisshaar B, Himmelbauer H: The genome of the recently domesticated crop plant sugar beet (Beta vulgaris). Nature. 2014, 505 (7484): 546-549. 10.1038/nature12817.View ArticlePubMedGoogle Scholar
- The Arabidopsis Genome Initiative: Analysis of the genome sequence of the flowering plant Arabidopsis thaliana. Nature. 2000, 408 (6814): 796-815. 10.1038/35048692.View ArticleGoogle Scholar
- Goff SA, Ricke D, Lan T-H, Presting G, Wang R, Dunn M, Glazebrook J, Sessions A, Oeller P, Varma H, Hadley D, Hutchison D, Martin C, Katagiri F, Lange BM, Moughamer T, Xia Y, Budworth P, Zhong J, Miguel T, Paszkowski U, Zhang S, Colbert M, Sun W-l, Chen L, Cooper B, Park S, Wood TC, Mao L, Quail P, et al: A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002, 296 (5565): 92-100. 10.1126/science.1068275.View ArticlePubMedGoogle Scholar
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, et al: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604. 10.1126/science.1128691.View ArticlePubMedGoogle Scholar
- Schmutz J, Cannon SB, Schlueter J, Ma J, Mitros T, Nelson W, Hyten DL, Song Q, Thelen JJ, Cheng J, Xu D, Hellsten U, May GD, Yu Y, Sakurai T, Umezawa T, Bhattacharyya MK, Sandhu D, Valliyodan B, Lindquist E, Peto M, Grant D, Shu S, Goodstein D, Barry K, Futrell-Griggs M, Abernathy B, Du J, Tian Z, Zhu L, et al: Genome sequence of the palaeopolyploid soybean. Nature. 2010, 463 (7278): 178-183. 10.1038/nature08670.View ArticlePubMedGoogle Scholar
- Velasco R, Zharkikh A, Affourtit J, Dhingra A, Cestaro A, Kalyanaraman A, Fontana P, Bhatnagar SK, Troggio M, Pruss D, Salvi S, Pindo M, Baldi P, Castelletti S, Cavaiuolo M, Coppola G, Costa F, Cova V, Dal Ri A, Goremykin V, Komjanc M, Longhi S, Magnago P, Malacarne G, Malnoy M, Micheletti D, Moretto M, Perazzolli M, Si-Ammour A, Vezzulli S, et al: The genome of the domesticated apple (Malus x domestica Borkh.). Nature Genetics. 2010, 42 (10): 833-839. 10.1038/ng.654.View ArticlePubMedGoogle Scholar
- Dohm JC, Lange C, Holtgräwe D, Sörensen TR, Borchardt D, Schulz B, Lehrach H, Weisshaar B, Himmelbauer H: Palaeohexaploid ancestry for Caryophyllales inferred from extensive gene-based physical and genetic mapping of the sugar beet genome (Beta vulgaris). Plant J. 2012, 70 (3): 528-540. 10.1111/j.1365-313X.2011.04898.x.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28 (10): 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMedGoogle Scholar
- Celenza JL, Quiel JA, Smolen GA, Merrikh H, Silvestro AR, Normanly J, Bender J: The Arabidopsis ATR1 Myb transcription factor controls indolic glucosinolate homeostasis. Plant Physiol. 2005, 137 (1): 253-262. 10.1104/pp.104.054395.PubMed CentralView ArticlePubMedGoogle Scholar
- Sonderby IE, Hansen BG, Bjarnholt N, Ticconi C, Halkier BA, Kliebenstein DJ: A systems biology approach identifies a R2R3 MYB gene subfamily with distinct and overlapping functions in regulation of aliphatic glucosinolates. PLoS One. 2007, 2 (12): e1322-10.1371/journal.pone.0001322.PubMed CentralView ArticlePubMedGoogle Scholar
- Gigolashvili T, Berger B, Mock HP, Muller C, Weisshaar B, Flugge UI: The transcription factor HIG1/MYB51 regulates indolic glucosinolate biosynthesis in Arabidopsis thaliana. Plant J. 2007, 50 (5): 886-901. 10.1111/j.1365-313X.2007.03099.x.View ArticlePubMedGoogle Scholar
- Gigolashvili T, Engqvist M, Yatusevich R, Muller C, Flugge UI: HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. New Phytol. 2008, 177 (3): 627-642. 10.1111/j.1469-8137.2007.02295.x.View ArticlePubMedGoogle Scholar
- Oppenheimer DG, Herman PL, Sivakumaran S, Esch J, Marks MD: A myb gene required for leaf trichome differentiation in Arabidopsis is expressed in stipules. Cell. 1991, 67: 483-493. 10.1016/0092-8674(91)90523-2.View ArticlePubMedGoogle Scholar
- Lee MM, Schiefelbein J: WEREWOLF, a MYB-related protein in Arabidopsis, is a position-dependent regulator of epidermal cell patterning. Cell. 1999, 24 (5): 473-483. 10.1016/S0092-8674(00)81536-6.View ArticleGoogle Scholar
- Kirik V, Lee MM, Wester K, Herrmann U, Zheng Z, Oppenheimer D, Schiefelbein J, Hulskamp M: Functional diversification of MYB23 and GL1 genes in trichome morphogenesis and initiation. Development. 2005, 132 (7): 1477-1485. 10.1242/dev.01708.View ArticlePubMedGoogle Scholar
- Serna L, Martin C: Trichomes: different regulatory networks lead to convergent structures. Trends Plant Sci. 2006, 11 (6): 274-280. 10.1016/j.tplants.2006.04.008.View ArticlePubMedGoogle Scholar
- Brockington SF, Alvarez-Fernandez R, Landis JB, Alcorn K, Walker RH, Thomas MM, Hileman LC, Glover BJ: Evolutionary analysis of the MIXTA gene family highlights potential targets for the study of cellular differentiation. Mol Biol Evol. 2013, 30 (3): 526-540. 10.1093/molbev/mss260.View ArticlePubMedGoogle Scholar
- Perez-Rodriguez M, Jaffe FW, Butelli E, Glover BJ, Martin C: Development of three different cell types is associated with the activity of a specific MYB transcription factor in the ventral petal of Antirrhinum majus flowers. Development. 2005, 132 (2): 359-370. 10.1242/dev.01584.View ArticlePubMedGoogle Scholar
- Quattrocchio F, Wing J, van der Woude K, Souer E, De Vetten N, Mol J, Koes R: Molecular analysis of the anthocyanin2 gene of Petunia and its role in the evolution of flower color. Plant Cell. 1999, 11: 1433-1444. 10.1105/tpc.11.8.1433.PubMed CentralView ArticlePubMedGoogle Scholar
- Borevitz JO, Xia YJ, Blount J, Dixon RA, Lamb C: Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000, 12 (12): 2383-2393. 10.1105/tpc.12.12.2383.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwinn K, Venail J, Shang Y, Mackay S, Alm V, Butelli E, Oyama R, Bailey P, Davies K, Martin C: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonzalez A, Zhao M, Leavitt JM, Lloyd AM: 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.View ArticlePubMedGoogle Scholar
- Kimler L, Mears J, Mabry TJ, Rösler H: On the question of the mutual exclusiveness of batalains and anthocyanins. Taxon. 1970, 19 (6): 875-878. 10.2307/1218301.View ArticleGoogle Scholar
- Schmidt OT, Schönleben W: Zur Kenntnis der Farbstoffe der Roten Rübe, II. Zeitschrift für Naturforschung. 1957, 12b: 262-263.Google Scholar
- Stafford HA: Anthocyanins and betalains: evolution of the mutually exclusive pathways. Plant Sci. 1994, 101: 91-98. 10.1016/0168-9452(94)90244-5.View ArticleGoogle Scholar
- Keller W: Inheritance of some major color types in beets. J Agric Res. 1936, 52: 27-38.Google Scholar
- Goldman IL, Austin D: Linkage among the R, Y and BI loci in table beet. Theor Appl Genet. 2000, 100: 337-343. 10.1007/s001220050044.View ArticleGoogle Scholar
- Hatlestad GJ, Sunnadeniya RM, Akhavan NA, Gonzalez A, Goldman IL, McGrath JM, Lloyd AM: The beet R locus encodes a new cytochrome P450 required for red betalain production. Nat Genet. 2012, 44 (7): 816-820. 10.1038/ng.2297.View ArticlePubMedGoogle Scholar
- Christinet L, Burdet FX, Zaiko M, Hinz UZ, Zrÿd JP: Characterization and functional identification of a novel plant 4,5-extradiol dioxygenase involved in betalain pigment biosynthesis in Portulaca grandiflora. Plant Physiol. 2004, 134 (1): 265-274. 10.1104/pp.103.031914.PubMed CentralView ArticlePubMedGoogle Scholar
- Heppel SC, Jaffé FW, Takos AM, Schellmann S, Rausch T, Walker AR, Bogs J: Identification of key amino acids for the evolution of promoter target specificity of anthocyanin and proanthocyanidin regulating MYB factors. Plant Mol Biol. 2013, 82 (4-5): 457-471. 10.1007/s11103-013-0074-8.View ArticlePubMedGoogle Scholar
- Aharoni A, De Vos CHR, Wein M, Sun Z, Greco R, Kroon A, Mol JNM, O'Connell AP: 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.View ArticlePubMedGoogle Scholar
- Deluc L, Barrieu F, Marchive C, Lauvergeat V, Decendit A, Richard T, Carde JP, Merillon JM, Hamdi S: Characterization of a grapevine R2R3-MYB transcription factor that regulates the phenylpropanoid pathway. Plant Physiol. 2006, 140 (2): 499-511. 10.1104/pp.105.067231.PubMed CentralView ArticlePubMedGoogle Scholar
- Grotewold E, Drummond BJ, Bowen B, Peterson T: The myb-homologous P gene controls phlobaphene pigmentation in maize floral organs by directly activating a flavonoid biosynthetic gene subset. Cell. 1994, 76 (3): 543-553. 10.1016/0092-8674(94)90117-1.View ArticlePubMedGoogle Scholar
- Czemmel S, Stracke R, Weisshaar B, Cordon N, Harris NN, Walker AR, Robinson SP, Bogs J: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Stracke R, Ishihara H, Huep G, Barsch A, Mehrtens F, Niehaus K, Weisshaar B: Differential regulation of closely related R2R3-MYB transcription factors controls flavonol accumulation in different parts of the Arabidopsis thaliana seedling. Plant J. 2007, 50 (4): 660-677. 10.1111/j.1365-313X.2007.03078.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Mehrtens F, Kranz H, Bednarek P, Weisshaar B: 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Gardner RL, Kerst AF, Wilson DM, Payne MG: Beta vulgaris L.: The characterization of three polyphenols isolated from the leaves. Phytochemistry. 1967, 6: 417-422. 10.1016/S0031-9422(00)86299-8.View ArticleGoogle Scholar
- Richardson M: Flavonols and C-glycosylflavonoids of the Caryophyllales. Biochem System Ecol. 1978, 6: 283-286. 10.1016/0305-1978(78)90046-7.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410. 10.1016/S0022-2836(05)80360-2.View ArticlePubMedGoogle Scholar
- Crooks GE, Hon G, Chandonia J-M, Brenner SE: WebLogo: A Sequence Logo Generator. Genome Res. 2004, 14: 1188-1190. 10.1101/gr.849004.PubMed CentralView ArticlePubMedGoogle Scholar
- Lohse M, Bolger AM, Nagel A, Fernie AR, Lunn JE, Stitt M, Usadel B: RobiNA: a user-friendly, integrated software solution for RNA-Seq-based transcriptomics. Nucleic Acids Res. 2012, 40 (Web Server issue): W622-W627. 10.1093/nar/gks540.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14 (4): R36-10.1186/gb-2013-14-4-r36.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- De Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20 (9): 1453-1454. 10.1093/bioinformatics/bth078.View ArticlePubMedGoogle Scholar
- Saldanha AJ: Java Treeview-extensible visualization of microarray data. Bioinformatics. 2004, 20 (17): 3246-3248. 10.1093/bioinformatics/bth349.View ArticlePubMedGoogle Scholar
- Jakoby M, Wang HY, Reidt W, Weisshaar B, Bauer P: FRU (BHLH029) is required for induction of iron mobilization genes in Arabidopsis thaliana. FEBS Lett. 2004, 577 (3): 528-534. 10.1016/j.febslet.2004.10.062.View ArticlePubMedGoogle Scholar
- Koncz C, Schell J: The promoter of TL-DNA gene 5 controls the tissue-specific expression of chimaeric genes carried by a novel type of Agrobacterium binary vector. Mol Gen Genet. 1986, 204: 383-396. 10.1007/BF00331014.View ArticleGoogle Scholar
- Sprenger-Haussels M, Weisshaar B: Transactivation properties of parsley proline rich bZIP transcription factors. Plant J. 2000, 22 (1): 1-8. 10.1046/j.1365-313x.2000.00687.x.View ArticlePubMedGoogle Scholar
- Baudry A, Heim MA, Dubreucq B, Caboche M, Weisshaar B, Lepiniec L: 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.View ArticlePubMedGoogle Scholar
- Stracke R, Jahns O, Keck M, Tohge T, Niehaus K, Fernie AR, Weisshaar B: Analysis of production of flavonol glycosides-dependent flavonol glycoside accumulation in Arabidopsis thaliana plants reveals MYB11-, MYB12- and MYB111-independent flavonol glycoside accumulation. New Phytol. 2010, 188 (4): 985-1000. 10.1111/j.1469-8137.2010.03421.x.View ArticlePubMedGoogle Scholar
- Hartmann U, Valentine WJ, Christie JM, Hays J, Jenkins GI, Weisshaar B: Identification of UV/blue light-response elements in the Arabidopsis thaliana chalcone synthase promoter using a homologous protoplast transient expression system. Plant Mol Biol. 1998, 36: 741-754. 10.1023/A:1005921914384.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.