- Research article
- Open Access
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 Biology volume 8, Article number: 83 (2008)
The MYB superfamily constitutes the most abundant group of transcription factors described in plants. Members control processes such as epidermal cell differentiation, stomatal aperture, flavonoid synthesis, cold and drought tolerance and pathogen resistance. No genome-wide characterization of this family has been conducted in a woody species such as grapevine. In addition, previous analysis of the recently released grape genome sequence suggested expansion events of several gene families involved in wine quality.
We describe and classify 108 members of the grape R2R3 MYB gene subfamily in terms of their genomic gene structures and similarity to their putative Arabidopsis thaliana orthologues. Seven gene models were derived and analyzed in terms of gene expression and their DNA binding domain structures. Despite low overall sequence homology in the C-terminus of all proteins, even in those with similar functions across Arabidopsis and Vitis, highly conserved motif sequences and exon lengths were found. The grape epidermal cell fate clade is expanded when compared with the Arabidopsis and rice MYB subfamilies. Two anthocyanin MYBA related clusters were identified in chromosomes 2 and 14, one of which includes the previously described grape colour locus. Tannin related loci were also detected with eight candidate homologues in chromosomes 4, 9 and 11.
This genome wide transcription factor analysis in Vitis suggests that clade-specific grape R2R3 MYB genes are expanded while other MYB genes could be well conserved compared to Arabidopsis. MYB gene abundance, homology and orientation within particular loci also suggests that expanded MYB clades conferring quality attributes of grapes and wines, such as colour and astringency, could possess redundant, overlapping and cooperative functions.
The MYB superfamily exemplifies how both conserved and divergent domains are present within a transcription factor family. MYB genes are exclusive to eukaryotes [1, 2] and although in animals their function is restricted to the control of the cell division and differentiation , they have diverse functions in plants .
MYB genes in plants are modulated by diverse hormones  and participate in key processes related to epidermal cell destiny  (which includes pigmentation  and formation of trichomes ), seed development , response to drought  and cold [11, 12], pathogen-disease resistance [13–15] stomatal movements [16, 17], phytochrome A-dependent light-sensing responses  and sucrose related responses , among many other functions.
Different MYB families have evolved after duplications of their DNA binding domains . MYB DNA binding domains are 100–160 residues in length, depending on the number of imperfect repeats (named R) in the N-terminal region. From the different classes identified, the R2R3 subfamily is the most abundant in plants [20–22]. Each repeat adopts a helix-helix-turn-helix structure which interacts with regulatory elements in the promoter, while the C-terminal region is responsible for establishing protein-protein interactions with other components of the eukaryotic transcriptional machinery. Based on their well conserved DNA-binding domains, R2R3 MYB genes have been annotated in both monocotyledonous and dicotyledonous genomes. In Arabidopsis thaliana, there are an estimated 126 R2R3 subfamily members , although only some of them have been functionally characterised.
Despite MYB genes have also been identified in woody species such as poplar, grape and apple, no genome-wide characterization of this family has been conducted in such species. In grapes, the MYB genes characterised to date are all involved in the control of flavonoid synthesis (MYBA , MYB5a (MYBCS-1) , MYB5b , and MYBPA1 ). Since quality parameters in grapes and wines such as colour, bitterness and astringency are determined by flavonoid compounds accumulating in fruit tissues, these studies have focussed on anthocyanin and proanthocyanidin synthesis and their regulation via MYB factors. In species such as Arabidopsis, petunia and maize, it has been determined that MYB proteins can act together with bHLH and WDR proteins in a transcriptional complex capable of regulating flavonoid synthesis and other processes concerning epidermal cell identity [27, 28]. The affinity of MYB proteins for the DNA binding sites in each promoter can be influenced by the interaction with these factors.
The annotated genome sequence of Vitis vinifera has recently become available  after an Italian-French effort to obtain a highly homozygous Pinot noir (PN40024) genome . In addition, an Italian consortium sequenced a highly heterozygous Pinot noir genotype (ENTAV115) [31, 32]. These two initiatives provided a genomic platform for studying this fruit crop. To date, its sequence has been assigned to 19 chromosomes  or linkage groups , the same number as in Poplar (Populus trichocarpa) . While the PN40024 genotype is predicted to contain 30,434 genes in 487 Mb (46% of the genome; ), the ENTAV115 clone is predicted to have 29,585 genes within 531 Mb . Although some supercontigs are still being correctly assembled and oriented in each chromosome, the sequence shows that 41.4% of the genome is composed of transposable elements mostly located in introns, explaining its high variability and heterozygocity . Paralogous regions in the grape genome arose by gene duplication and are present in clusters, forcing gene expansion and diversification. The best example of this event is the terpene synthase gene family, with 30–40 members in Arabidopsis and almost 90 members in grape . To date, information regarding expansion events of the MYB gene superfamily in Vitis is absent.
In this work, we describe the putative R2R3 MYB gene subfamily by means of in silico analysis of the grape genome sequence, in order to predict protein domain architectures, and to assess the extent of conservation and divergence between grape and Arabidopsis gene families. Candidate genes were chosen for isolation and their expressions were tested in different grape organs to compare expression patterns of closely grouping co-orthologues. This extended analysis is the first to be performed on a dicotyledonous species other than Arabidopsis. Possible gene retention, loss and expansion processes of MYB genes are discussed.
Results and Discussion
Identification and chromosomal distribution of grape R2R3 MYB genes
The genome sequence of the homozygous PN40024 genotype of Vitis vinifera cv. Pinot noir  was searched for MYB gene models. These are predicted from combining ab initio models together with V. vinifera complementary DNA sequences (from EST databases) and alignments of gene/protein models from other species . A consensus R2R3 MYB DNA binding domain sequence (Additional file 1) was used in a BLATsearch, which maintains an index of the entire genome in memory and finds sequences of 20 amino acids or more with 80% or higher similarity in the genome. The MYB superfamily has been defined as the most abundant transcription factor family in plants , with at least 198 members in Arabidopsis and 183 in rice . In grapevine, a total of 279 MYB genes were previously estimated , including members of the R1R2R3, R2R3, atypical MYB and MYB-related gene subfamilies. The Arabidopsis and rice R2R3 subfamilies are composed of 126 and 109 members, respectively , while in grape 108 candidate genes were identified in this study. These were further analysed in terms of their putative protein sequences, exon lengths and their positions and orientations in each chromosome (Additional file 2).
This analysis revealed that grape MYB genes are distributed in almost all chromosomes, except for chromosome (chr) 10 (Figure 1). However, thirteen MYB gene models have not yet been assigned to a particular chromosome and it is possible that some of them are positioned in chr10. Chromosomes with more abundant MYB gene models are 1, 3 and 5 in Arabidopsis and 2, 4 and 14 in the Vitis genome. These were subsequently analysed for gene cluster expansion events.
Exon/intron organisation of Vitis and Arabidopsis MYB families
Exon structures are highly conserved when the grape and Arabidopsis R2R3 MYB subfamilies are compared (Figures 2 and 3). Exons 1 and 2 code for almost the entire R2R3 DNA binding domain, although this pattern differs in complex multiexonic genes (e.g. AtMYB88 and 124). As shown in Figure 2A, exon 1 appears to be the most restricted in length, while exons 2 and 3 are more variable (<100 bp – >1000 bp). Presence of a fourth and fifth exon is exclusive to some particular genes. Despite this variability, the modal lengths of the first two exons are very similar (exon 1, 133 bp; exon 2, 130 bp) and highly conserved (exon 1, 38.7% occurrence; exon 2, 63.2% occurrence). These occurrence values are very similar in Arabidopsis, suggesting that MYB DNA binding domains could be partially conserved because exons coding for this domain have restricted lengths during the evolution of plant species. It is interesting to note that the exon lengths observed for the R2R3 MYB subfamily coincide with the average exon length predicted from models of other grape gene families (130 bp, ).
Although exon 3 is the most diverse in size, MYB subfamilies from both species are similarly distributed when the first three exon lengths are considered (Figure 2B). Exon 3 codes for the last region of the R3 repeat and for the C-region of the protein. Changes both in length and sequence of this exon could have generated functional divergence between MYB homologues within and between species, leading to different functional motifs and domains .
In order to gain further insights into exon conservation, the split positions of the R2R3 repeats within the DNA binding domain were classified in the grape and Arabidopsis gene models (Figure 3; note that this classification does not consider the final number of exons). It is possible to distinguish four groups in Arabidopsis, the first three of which are also present in grape (Figure 3A). All groups differ in the organisation of exons determining the R2R3 domain. In general, the R2R3 domain consists of 106 amino acids. The N-terminal region before this domain may vary within each group. Groups I and II have their R2 repeat split between exons 1 and 2, but in Group I, the R3 repeat is divided between exons 2 and 3 (Figure 3B; note that the last exon in Groups I, II and III may have additional introns depending on the gene). In Group III, the first exon codes the complete R2 repeat but R3 is split between in exons 1 and 2. Group IV consists only of 1 exon coding for the complete MYB protein: members of this group have not been identified to date in grapes. AtMYB88 and AtMYB124 are not included in any group of the Arabidopsis plot because of their more complex exon/intron organization (10–11 exons, respectively).
Seven Arabidopsis MYB genes and the grape gene model GSVIVT16947001 could not be positioned in these groups. Together with MYB118 and 106, all share a longer N-terminal region before the start of the R2 repeat (between 52 and 213 residues). Regarding their exon boundaries, the R2 and R3 repeats are organised differently and can not be classified by these means. The other Arabidopsis MYBs and their grape gene homologues are maintained in the same respective groups. Despite the occurrence of different duplication and diversification events in these two species, MYB exons are highly conserved. Length and sequence restrictions in the first two exons could help maintain functionality and achieve sequence specificity of the R2R3 domain, while variations in the exons coding for the C-terminal regions could facilitate the gaining of new or cooperative functions when new domains are included.
Updated functional clustering of the Arabidopsis R2R3 MYB Subfamily
Taking the complete Arabidopsis R2R3 MYB subfamily (126 members, Additional file 3), a phylogenetic tree was constructed to update new functional clades which could also be present in the grape subfamily (Figure 4). Most genes sharing similar functions were clustered in the same phylogenetic clades and subclades (functional clades), suggesting that most closely-related MYBs in the same species could recognise similar target genes and possess redundant, overlapping, and/or cooperative functions.
The resulting tree is validated as it shows the same subgroups (1–17) observed in phylogenetic trees constructed by Kranz et al  and Yanhui et al . New functions have been recently identified, including MYBs involved in glucosinolate biosynthesis [37–41] and hormone induced responses. These new functions are also clustered within this phylogenetic analysis (Figure 4), although a particularly polyphyletic feature within this Subfamily is observed: flower development. Flower organs were conceived as an adaptive process late in plant evolution. Several mechanisms involving cell cycle regulation, secondary metabolite accumulation and hormonal control have been related to reproductive male organ development. MYB genes regulating these events have been recently identified (AtMYB103 [42, 43], AtMYB26 , AtMYB24 , AtMYB21 , AtMYB32  and AtMYB101, AtMYB65 and AtMYB33 [48, 49]). Although they cluster in different clades within the MYB subfamily and their target genes seem to be different, all ultimately regulate anther development.
The R2R3 domain provides both the MYB nuclear localisation signals  and regulatory specificity . C-terminal regions often possess protein-protein interaction or DNA-repression motifs. Such is the case of the C2 motif repressor clade. Members of this group possess the C1 and C2 motifs, known to participate in bHLH interactions  and promoter repression, respectively (e.g. CH4 promoter for AtMYB4;  and RD29A promoter for AtMYB8; ).
Among flavonoids, flavonols are known to have functions such as filtering UV radiation and auxin retention and transport (reviewed in ). They also have pharmaceutical properties when consumed by animals. Recently, their synthesis has been shown to be regulated by three closely related MYBs (AtMYB11, 12 and 111 ). These MYBs are found together in the 'production of flavonol glycosides' clade (Figure 4) and act at different times and tissues throughout plant development . This new clade appears to share an ancestral MYB gene both with the 'C2-repressor motif' clade and the 'epidermal cell fate' clade (which includes the anthocyanin and trichome related subclades). These clades regulate different processes but share up to some extent the control of phenolic compound synthesis.
Phylogenetic analysis of the grape R2R3 MYB Subfamily
The same alignment protein weight matrices and phylogeny algorithms used for the construction of the Arabidopsis tree were employed for the integrated analysis of the Arabidopsis, rice and Vitis MYB subfamilies (Additional files 4 and 5). Grape protein models with a complete R2R3 domain (91 sequences), seven characterised grape proteins and 126 Arabidopsis homologues were included in the first phylogeny. The remaining grape models generally possessed an incomplete R2 repeat, so were not included in this analysis.
Sequence alignment and phylogenetic tree analysis using both the single R2R3 DNA-binding domain and the full predicted protein sequences were produced (Additional file 4A and 4B, respectively). Outside the DNA binding domain, MYB proteins had the greatest sequence divergence. Although a full sequence alignment could be highly sensitive to the dissimilarities in the C-terminal regions outside the conserved motifs, both trees are extremely similar regarding topology and bootstrap values. Possible orthologues between both species appear outside and inside the functional clades previously described for Arabidopsis (Additional files 4 and 5).
In some cases, it was easy to identify two putative orthologues in Arabidopsis and Vitis, since they were grouped in pairs within a clade. In other situations, grape and Arabidopsis homologues were clustered by species within a clade (e.g. proteins from the C2-repressor motif clade, with the exception of MYB4). Their functional characterisation will help prove which genes are orthologues within each clade.
None of the grape gene models were grouped within the Arabidopsis 'glucosinolate' clade (Additional file 4). Several lines of evidence suggests that this clade resulted from a β-type duplication  which evolved independently in the order Brassicales as an adaptive response to herbivory [54, 55], explaining its absence in rice and grape.
For several transcription factor families, linear correlations of abundance have been identified in Arabidopsis and Vitis . In Arabidopsis, the MYB superfamily consists of 198 R1R2R3, R2R3, MYB related and atypical MYB genes , in contrast to the 279 estimated MYB superfamily members in Vitis. In addition, because of the anatomical and physiological differences between grapes and Arabidopsis, it can be suggested that some clades could be differentially expanded when comparing Arabidopsis and Vitis R2R3 subfamilies. Both inside and outside of some of these functional clades, grape gene expansions were observed as gene pairs or clusters (highlighted with asterisks in Additional file 4). Such is the case of the Vitis-specific MYB5a and MYB5b, the MYB4 gene pair within the C2-repressor motif clade, MYB101 ('GAMYB like/miRNA regulated' clade) and MYB 38/RAX1 ('axillary meristem' clade). Outside the functional clades, expanded MYBs are represented by dichotomous pairs of grape MYB15, 17, 35, 67 and 102 homologues, whose functions have yet to be determined. In the case of the 'production of flavonol glycosides' clade, one member (possibly MYB11) may be absent in grape. This could be a result of the loss of a clade member or a mis-annotation of the Vitis genome. It is possible that new R2R3 MYB genes could be identified in the future as annotations improve. For instance, AtMYB88 and AtMYB124 grape homologues, which belong to the 'guard cell division restriction' clade, have yet to be found, possibly because their complex exon/intron structure (10–11 exons) makes them difficult to identify during an annotation procedure.
Expansion of grape MYB genes affecting wine quality parameters
After the sequencing of the grapevine genome, the expansion of large gene families involved in wine characteristics was observed . Such was the case of the Stilbene Synthase Family, responsible for the synthesis of the phytoalexin resveratrol, and the Terpene Synthase Family, involved in aroma production. Other studies have also revealed that genes belonging to terpene and flavonoid metabolic pathways have suffered selective amplifications in the grape genome, in contrast to other plant genomes . If this situation is evident for some biosynthetic genes of the grape phenylpropanoid pathway (PAL, CHS, F3H, F3'5'H, FLS and LAR), it is reasonable to expect that MYB genes controlling this pathway could also be expanded. Nevertheless, it is also possible to suggest that these genes were maintained in the grapevine genome by means of human domestication, and lost in other plant genomes. For example, many genes were lost in Arabidopsis after its two wide genome duplication (WGD) events .
Two different scenarios have been described regarding possible WGD episodes in grape: the first states there is no evidence for such event , while in the second, a recent WGD event might have occurred in close proximity to the Vitis speciation event . In our analysis we observed that MYB genes from the 'epidermal cell fate' clade involved in the synthesis of molecules which confer quality to grapes and wines are largely expanded in the Vitis genome.
Grape and wine colour is defined by glycosylated anthocyanins. Any perturbation in their cytoplasmic synthesis, conjugation or vacuole transport thus directly affects berry colour. Several MYB genes have been characterised in Arabidopsis, maize, petunia, antirrhinum, and apple, which regulate some of the processes concerning anthocyanin production. It was discovered that white grapes arose from multiallelic mutations of the VvMYBA1/A2 genes [57–60], which controls the last biosynthetic step of anthocyanin synthesis, a glycosylation reaction mediated by the UDP-GLUCOSE FLAVONOID 3-O-GLUCOSYLTRANSFERASE (UFGT) enzyme .
The Vitis-Arabidopsis phylogenetic tree shows that there are nine anthocyanin-related MYB gene models distributed on chromosomes 2 and 14 from grape, the most abundant MYB-containing chromosomes (Figure 5A). An additional gene model (GSVIVT 9898001) has yet to be linked to a specific chromosome. Recently, a grape colour locus was described , which we now localised on chromosome 2. In the previous study, four R2R3 MYB genes were found at this locus, but only two (VvMYBA1 and VvMYBA2) had a regulatory effect on fruit colour. These proteins differ in the lengths of their C-regions: MYBA2 has a duplicated C-terminal repeat when compared to MYBA1, while MYBA3 almost completely lacks this region. These differences could cause novel or cooperative functions, and it has been suggested that VvMYBA3 could control anthocyanin synthesis in grape tissues other than the fruits . In addition to the VvMYBA1-3 genes, three other R2R3 MYB gene models were detected in the grape colour locus located on chromosome 2 (Figure 5B, between 12 and 13 Mb. Note that MYBA1 was annotated as two contiguous gene models: GSVIVT38763001 and 38762001).
Chromosome 14 contains three anthocyanin related MYBs, some of them being predicted as possible orthologues of the Arabidopsis anthocyanin-related AtMYB75 and AtMYB90 proteins (Figure 5A). Further experiments are required to determine whether these genes are functional, in order to verify if this new locus influences the colouring of fruits or other plant organs.
In concordance with the Arabidopsis chromosomal duplications classified by Bowers et al , MYB gene evolution is though to involve segmental duplications . The Arabidopsis 'anthocyanin-related' subclade, which includes AtMYB75, 90 , 113 and 114  (all positioned on chromosome 1), was described as a cluster which may have arise by tandem duplications (see Supplementary material in ). The anthocyanin and trichome-related subclades within the Arabidopsis 'epidermal cell fate' clade may have experienced a major rapid expansion after diverging from monocots but before separating from other dicots . Supporting this hypothesis is the finding that these two clusters are absent in the rice subfamily (Additional file 5). Considering the homology and orientation of the nine anthocyanin related MYBs found in the two grape loci, it is possible that some of the original MYBs could have undergone new multiple tandem duplication events somewhere before or after the separation of the Arabidopsis-Vitis lineages. However, if this event occurred before lineage separation, this second expansion has clearly been lost during the subsequent evolution of the Arabidopsis genome.
Seed pigmentation is defined by the accumulation and oxidation of flavan-3-ols, which are transported in the same way as anthocyanins, but are additionally polymerised into proanthocyanidins (PA) and condensed tannins inside vacuoles (reviewed in ). In wines, these molecules define astringency and bitterness. AtMYB123 (also known as TT2) controls the final steps of PA synthesis . Our analysis suggests that grapevine possesses several putative TT2 orthologues (Figure 5A), many of which are more closely related to TT2 than the previously characterised MYB5a, MYB5b and MYBPA1. These genes control PA synthesis and also regulate other biosynthetic genes of the flavonoid pathway [24–26]. Chromosomes 4, 9 and 11 possess four, one and two grape TT2 homologue genes, respectively. MYBPA1/GSVIVT00006679001 and MYB5a/b are shown as dichotomous pairs. Taken together, MYBs related to PA synthesis also appear expanded. The fact that in Lotus japonicus the TT2 lineage compromises a three-member subclade , addresses the question whether this expansion is exclusive to grape. When including the rice MYB subfamily into this analysis, only two models appeared as putative TT2 orthologues (Additional file 5), indicating that grape has substantially more TT2 loci than other species characterised to date.
The observed new expansion events are specific to the anthocyanin related subclade and the TT2 related lineage in grapevine, but not to the trichome subclade. Both the grapevine adaptive responses (because of the adaptative functions of flavonoids) and human breeding and domestication practices could have helped in the selection of grape genotypes with expanded MYB genes involved in colour and astringency.
Selection of grape gene models for expression analysis and modelling of their DNA binding domains
In order to test the annotation procedure and validate gene models, six annotated genes were randomly selected, isolated and their expression patterns were screened in different grapevine organs. VvMYB4 (Genbank accession EF113078), a AtMYB4 repressor orthologue recently characterised in our group was also included. Complete coding regions of VvMYB12-like, 14-like, 24-like, 30-like, 60-like and a C2motif containing repressor-like gene were isolated from different organs of cv. Cabernet-Sauvignon by RT-PCR (Genbank accessions in Additional file 6). The cDNA sequences were 95–100% similar to their respective gene models (annotated from cv. Pinot noir).
Although some Arabidopsis MYBs have yet to be functionally characterised, global gene expression data throughout plant development is available on the Arabidopsis Electronic Fluorescent Pictograph (eFP) Browser [66, 67]. All selected grape MYB genes are highly expressed in reproductive organs and revealed organ-expression profiles remarkably similar to their putative Arabidopsis orthologues (Figure 6 and Additional file 7). In addition, VvMYB24 and 14-like genes are expressed in a tissue-specific manner in ripped berries/inflorescences and seeds, respectively, in a similar way as in Arabidopsis. Subsequent experiments are needed to examine the precise function of these putative orthologies.
A comparative modelling approach (SWISSMODEL) using crystallised mammalian c-MYB (1GV2.PDB) as the structural template was conducted in order to study DNA binding domains from the previously selected MYBs. The in silico 3D model of these R2R3 domains suggest significant structural homology with rat R1R2R3 c-MYB (Additional file 8). Structural studies of the mammalian c-MYB factor show that its DNA binding domain requires several trypthophan residues which participate in the hydrophobic grouping within the DNA-protein interaction interface . These residues are maintained in the grape homologues (as displayed in asterisks in Additional file 8A).
Some subtle differences between the animal and plant MYBs are also present. Animal MYB domains have a single cysteine residue that needs to be reduced for DNA binding and transcriptional activity , contrary to most plant R2R3 MYB domains (including grapes) which contain two cysteines (Additional file 8A,C). Under non-reducing conditions, these cysteines form a disulfide bond preventing DNA binding . The second cysteine residue, conserved in both animals and plants, is essential for regulatory activity in animals. Comparatively, the first cysteine, present only in plants, is crucial for MYB function, demonstrating structural and functional differences between plants and animals that could influence the dynamics of the protein-DNA interaction. The residues surrounding these cysteines are highly similar in the grape MYB domains analysed, further suggesting the importance of their conservation.
Turns between the α2 and α3 helices in each repeat (Additional file 8B, middle and right panels) are essential for conferring globular-like structures (Additional file 8C). Grape R2R3 MYB domains are mostly conserved in the DNA-interacting R2α3, R3α1 and R3α2 helices and turns, as well as in their charge. The most distinct domains in terms of general architecture and positive net charge (yellow regions in each of the structures) were those from VvMYB4 and VvMYBC2 repressor-like proteins, possibly due to their mechanism of action; it has been postulated that MYB repressors affect other MYBs in their binding by either protein-protein interactions or as dominant negative factors .
The sequencing of the grape genome constituted a powerful tool for gene search and evolution studies. Many protein families have been studied in terms of genome conservation across species, in order to establish a genome evolution scenario regarding duplication events and functional diversification. As an example, exon-intron boundaries, expression patterns and regulatory elements in plant D-type cyclin promoters were found to be strongly conserved between Arabidopsis and poplar genomes . The identification of grape MYB genes and their comparative analysis with the Arabidopsis genome suggests strong conservation but also expansions of particular functional clades. Despite recent events of diversification, the exon/intron organisation of MYB genes is highly conserved, suggesting restrictions in R2R3 domain variations while allowing C-terminal modifications for gaining new or cooperative functions.
The identification of MYB genes in grapevine opens the possibility of modulating MYB gene expression in order to control specific aspects of grape physiology and development. In particular, genes involved in flower and seed development or in flavonoid synthesis might be interesting targets for functional characterisation in light of expression similarities in Arabidopsis and their possible duplication histories.
Search for MYB homologues in the Grape Genome and Isolation of Candidate genes
The consensus MYB R2R3 DNA binding domain sequence, obtained from the alignment of previously-isolated and characterised genes from Vitis vinifera and Arabidopsis (Additional file 1), was used in a BLATsearch to identify homologous gene models in the Grape Genome Browser, designed by Genoscope (the French National Sequencing Center ). The results obtained were compared to tentative consensus (TC) sequences available at the Grape TIGR Gene Index EST-database.
Phylogeny reconstruction and bootstrap analysis
To place the obtained sequences in a phylogenetic context, a phylogeny reconstruction was performed with the single DNA-binding domain and full predicted amino acid sequences of grape MYB genes and then integrated in the Arabidopsis phylogenetic tree. Alignments were performed using the BLOSUM matrix (Gap opening and extension penalties of 25 and 1, respectively) using the ClustalW algorithm-based AlignX module from Mega4 software . The phylogenetic tree was constructed using the Neighbour Joining Tree Method in Mega4 and confirmed with Mega3. Tree nodes were evaluated by bootstrap analysis for 1000 or 2000 replicates. All positions containing gaps and missing data were eliminated from the dataset (Complete deletion option in Mega4 software). Evolutionary distances were computed using the Poisson correction method and are represented as amino acid substitutions per site. Sequences of all 126 Arabidopsis R2R3 MYB proteins were downloaded from the TAIR Arabidopsis genome annotation version 7.0, released on April 2007.
Plant materials, RNA isolation and cDNA synthesis
Seven MYB candidate genes were isolated from different cDNAs using specific primers (Additional file 7). Different vegetative and reproductive tissues were collected from grapevine plants (Vitis vinifera L. cv. Cabernet-Sauvignon) growing in commercial fields in central Chile. Total RNA was isolated according to the procedure of Reid et al , using a CTAB-spermidine extraction buffer. For cDNA synthesis, one μg of total RNA was reverse transcribed with random hexamer primers in an 18 μl reaction mixture using the StrataScript® reverse transcriptase (Statagene) according to the manufacturer's instructions. Standard PCR profiles (35 cycles of 30 sec. at 95°C for denaturing, 30 sec. at 57°C for annealing and 1 min at 72°C for extension) were used to amplify MYB genes. Amplified fragments were cloned in TOPO-SD cloning vector (Invitrogen, USA), checked for insertion by PCR amplification and sequenced.
Expression of grape MYB genes in grape organs
Relative transcript quantification of MYB genes was achieved by real time RT-PCR, using the Brilliant® SYBR® Green QPCR Master Reagent Kit (Stratagene) and the Mx3000P detection system (Stratagene) as described in the manufacturer's manual. PCR conditions, standard quantification curves for each gene and relative gene expression calculations were conducted according to . Amplification of the UBIQUITIN1 gene (99 bp; TC53702, TIGR database, VvGi5) was used for normalisation. The ratio between the gene of interest (GOI) and UBIQUITIN expression was calculated using the equation :
where E corresponds to each primer amplification efficiency value. Gene expression levels were normalised to that of the organ sample with the lowest expression in order to obtain a calibrated ΔCt for each gene.
Comparative modelling of grape R2R3 MYB DNA binding domain structures
Using the Ratus sp. c-MYB crystal structure (1GV2.pdb) as a template, grape MYB structures were produced by SWISSMODEL  and viewed using 3D Molecule Viewer. Charge distributions were analysed in each structure generated.
Yang T, Perasso R, Baroin-Tourancheau A: MYB genes in ciliates: a common origin with the MYB proto-oncogene?. Protist. 2003, 154: 229-238.
Jiang C, Gu J, Chopra S, Gu X, Peterson T: Ordered origin of the typical two and three repeat MYB genes. Gene. 2004, 326: 13-22.
Weston K: Myb proteins in life, death and differentiation. Curr Opin Genet Dev. 1998, 8: 76-81.
Martin C, Paz-Ares J: MYB transcription factors in plants. Trends Genet. 1997, 13: 67-73.
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-24.
Ramsay N, Glover B: MYB-bHLH-WD40 protein complex and the evolution of cellular diversity. Trends Plant Sci. 2005, 10 (2): 63-70.
Nesi N, Jond C, Debeaujon I, Caboche M, Lepiniec L: The Arabidopsis TT2 gene encodes an R2R3 MYB domain protein that acts as a key determinant for proanthocyanidin accumulation in developing seed. Plant Cell. 2001, 13: 2099-2114.
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.
Penfield S, Meissner RC, Shoue DA, Carpita NC, Bevan MW: MYB61 is required for mucilage deposition and extrusion in the Arabidopsis seed coat. Plant Cell. 2001, 13: 2777-2791.
Abe H, Yamaguchi-Shinozaki K, Urao T, Iwasaki T, Hosokawa D, Shinozaki K: Role of Arabidopsis MYC and MYB homologs in drought and abscisic acid-regulated gene expression. Plant Cell. 1997, 9: 1859-1868.
Agarwal M, Hao Y, Kapoor A, Dong CH, Fujii H, Zheng X, Zhu JK: A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance. J Biol Chem. 2006, 281 (49): 37636-45.
Zhu J, Verslues PE, Zheng X, Lee BH, Zhan X, Manabe Y, Sokolchik I, Zhu Y, Dong CH, Zhu JK, Hasegawa PM, Bressan RA: HOS10 encodes an R2R3-type MYB transcription factor essential for cold acclimation in plants. Proc Natl Acad Sci USA. 2005, 102 (28): 9966-71.
Lee MW, Qi M, Yang Y: A novel jasmonic acid-inducible rice MYB gene associates with fungal infection and host cell death. Mol Plant Microbe Interact. 2001, 14: 527-35.
Vailleau F, Daniel X, Tronchet M, Montillet JL, Triantaphylidès C, Roby D: A R2R3-MYB gene, AtMYB30, acts as a positive regulator of the hypersensitive cell death program in plants in response to pathogen attack. Proc Natl Acad Sci USA. 2002, 99 (15): 10179-84.
Mengiste T, Chen X, Salmeron J, Dietrich R: The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell. 2003, 15 (11): 2551-65.
Cominelli E, Galbiati M, Vavasseur A, Conti L, Sala T, Vuylsteke M, Leonhardt N, Dellaporta SL, Tonelli C: A guard-cell-specific MYB transcription factor regulates stomatal movements and plant drought tolerance. Curr Biol. 2005, 15: 1196-1200.
Liang YK, Dubos C, Dodd IC, Holroyd GH, Hetherington AM, Campbell MM: AtMYB61, an R2R3-MYB transcription factor controlling stomatal aperture in Arabidopsis thaliana. Curr Biol. 2005, 15 (13): 1201-6.
Ballesteros ML, Bolle C, Lois LM, Moore JM, Vielle-Calzada JP, Grossniklaus U, Chua NH: LAF1, a MYB transcription activator for phytochrome A signaling. Genes Dev. 2001, 15 (19): 2613-25.
Teng S, Keurentjes J, Bentsink L, Koornneef M, Smeekens : Sucrose-specific induction of anthocyanin biosynthesis in Arabidopsis requires the MYB75/PAP1 gene. Plant Physiol. 2005, 139 (4): 1840-52.
Romero I, Fuertes A, Benito MJ, Malpica JM, Leyva A, Paz-Ares J: More than 80 R2R3-MYB regulatory genes in the genome of Arabidopsis thaliana. Plant J. 1998, 14: 273-284.
Jin H, Martin C: Multifunctionality and diversity within the plant MYB-gene family. Plant Mol Biol. 1999, 41: 577-585.
Stracke R, Werber M, Weisshaar B: The R2R3-MYB gene family in Arabidopsis thaliana. Curr Opin Plant Biol. 2001, 4: 447-456.
Kobayashi S, Ishimaru M, Hiraoka K, Honda C: Myb-related genes of the Kyoho grape (Vitis labruscana) regulate anthocyanin biosynthesis. Planta. 2002, 215: 924-933.
Deluc L, Barrieu F, Marchive C, Lauvergeat V, Decendit A, Richard T, Carde J-P, Mérillon JM, Hamdi S: Characterization of a grapevine R2R3-MYB transcription factor that regulates the phenylpropanoid pathway. Plant Physiol. 2006, 140: 499-511.
Deluc L, Bogs J, Walker A, Ferrier T, Decendit A, Merillon JM, Robinson S, Barrieu F: The transcription factor VvMYB5b contributes to the regulation of anthocyanin and proanthocyanidin biosynthesis in developing grape berries. Plant Physiol. 2008
Bogs J, Jaffe FW, Takos AM, Walker AR, Robinson SP: The grapevine transcription factor VvMYBPA1 regulates proanthocyanidin synthesis during fruit development. Plant Physiol. 2007, 143 (3): 1347-61.
Payne CT, Zhang F, Lloyd AM: GL3 encodes a bHLH protein that regulates trichome development in Arabidopsis through interaction with GL1 and TTG1. Genetics. 2000, 156: 1349-1362.
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-80.
Grape Genome Browser. [http://www.genoscope.cns.fr]
Jaillon , et al: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449 (7161): 463-7.
Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, Grando MS, Velasco R: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring Pinot Noir bacterial artificial chromosome contigs. Genetics. 2007, 176 (4): 2637-50.
Velasco R, et al: A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS ONE. 2007, 2 (12): e1326-
Tuskan G, et al: The Genome of Black Cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313 (5793): 1596-1604.
Riechmann JL, Ratcliffe OJ: A genomic perspective on plant transcription factors. Curr Opin Plant Biol. 2000, 3: 423-434.
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-20.
Kranz H, Denekamp M, Greco R, Jin H, Leyva A, Meissner R, 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 (2): 263-276.
Celenza JL, Quiel JA, Smolen GA, Merrikh H, Silvestro A, Normanly J, Bender J: The Arabidopsis ATR1 Myb transcription factor controls indolic glucosinolate homeostasis. Plant Physiol. 2005, 137: 253-262.
Gigolashvili T, Berger B, Mock HP, Müller C, Weisshaar B, Flügge UI: The transcription factor HIG1/MYB51 regulates indolic glucosinolate biosynthesis in Arabidopsis thaliana. Plant J. 2007, 50 (5): 886-901.
Gigolashvili T, Yatusevich R, Berger B, Müller C, Flügge UI: The R2R3-MYB transcription factor HAG1/MYB28 is a regulator of methionine-derived glucosinolate biosynthesis in Arabidopsis thaliana. Plant J. 2007, 51 (2): 247-61.
Gigolashvili T, Engqvist M, Yatusevich R, Müller C, Flügge 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-42.
Sønderby 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-
Higginson T, Li SF, Parish RW: AtMYB103 regulates tapetum and trichome development in Arabidopsis thaliana. Plant J. 2003, 35 (2): 177-92.
Zhang ZB, Zhu J, Gao JF, Wang C, Li H, Li H, Zhang HQ, Zhang S, Wang DM, Wang QX, Huang H, Xia HJ, Yang ZN: Transcription factor AtMYB103 is required for anther development by regulating tapetum development, callose dissolution and exine formation in Arabidopsis. Plant J. 2007, 52 (3): 528-38.
Yang C, Xu Z, Song J, Conner K, Vizcay Barrena G, Wilson ZA: Arabidopsis MYB26/MALE STERILE35 regulates secondary thickening in the endothecium and is essential for anther dehiscence. Plant Cell. 2007, 19 (2): 534-48.
Yang XY, Li JG, Pei M, Gu H, Chen ZL, Qu LJ: Over-expression of a flower-specific transcription factor gene AtMYB24 causes aberrant anther development. Plant Cell Rep. 2007, 26 (2): 219-28.
Mandaokar A, Thines B, Shin B, Lange BM, Choi G, Koo YJ, Yoo YJ, Choi YD, Choi G, Browse J: Transcriptional regulators of stamen development in Arabidopsis identified by transcriptional profiling. Plant J. 2006, 46 (6): 984-1008.
Preston J, Wheeler J, Heazlewood J, Li SF, Parish RW: AtMYB32 is required for normal pollen development in Arabidopsis thaliana. Plant J. 2004, 40 (6): 979-95.
Millar AA, Gubler F: The Arabidopsis GAMYB-like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell. 2005, 17 (3): 705-21.
Reyes JL, Chua NH: ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 2007, 49 (4): 592-606.
Li J, Li X, Guo L, Lu F, Feng X, He K, Wei L, Chen Z, Qu LJ, Gu H: A subgroup of MYB transcription factor genes undergoes highly conserved alternative splicing in Arabidopsis and rice. J Exp Bot. 2006, 57 (6): 1263-73.
Grotewold E, Sainz MB, Tagliani L, Hernandez JM, Bowen B, Chandler VL: Identification of the residues in the Myb domain of maize C1 that specify the interaction with the bHLH cofactor R. Proc Natl Acad Sci USA. 2000, 97 (25): 13579-84.
Jin H, Cominelli E, Bailey P, Parr A, Mehrtens F, Jones J, Tonelli C, Weisshaar B, Martin C: Transcriptional repression by AtMYB4 controls production of UV-protecting sunscreens in Arabidopsis. EMBO J. 2000, 19 (22): 6150-61.
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: 660-677.
Grubb CD, Abel S: Glucosinolate metabolism and its control. Trends Plant Sci. 2006, 11 (2): 89-100.
Halkier BA, Gershenzon J: Biology and biochemistry of glucosinolates. Annu Rev Plant Biol. 2006, 57: 303-33.
Ku H-M, Vision T, Liu J, Tanksley SD: Comparing sequenced segments of the tomato and Arabidopsis genomes: Large-scale duplication followed by selective gene loss creates a network of synteny. Proc Natl Acad Sci USA. 2000, 97: 9121-9126.
Kobayashi S, Goto-Yamamoto N, Hirochika H: Retrotransposon-Induced Mutations in Grape Skin Color. Science. 2004, 304 (5673): 982-
Lijavetzky D, Ruiz-García L, Cabezas JA, De Andrés MT, Bravo G, Ibáñez A, Carreño J, Cabello F, Ibáñez J, Martínez-Zapater JM: Molecular genetics of berry color variation in table grape. Mol Genet Genomics. 2006, 276 (5): 427-35.
This P, Lacombe T, Cadle-Davidson M, Owens CL: Wine grape (Vitis vinifera L.) color associates with allelic variation in the domestication gene VvmybA1. Theor Appl Genet. 2007, 114 (4): 723-30.
Walker AR, Lee E, Bogs J, McDavid DA, Thomas MR, Robinson SP: White grapes arose through the mutation of two similar and adjacent regulatory genes. Plant J. 2007, 49 (5): 772-85.
Bowers JE, Chapman BA, Rong J, Paterson AH: Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events. Nature. 2003, 422: 433-438.
Borevitz JO, Xia Y, Blount J, Dixon RA, Lamb C: Activation tagging identifies a conserved MYB regulator of phenylpropanoid biosynthesis. Plant Cell. 2000, 12 (12): 2383-2394.
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. 2007
Lepiniec L, Debeaujon I, Routaboul JM, Baudry A, Pourcel L, Nesi N, Caboche M: Genetics and biochemistry of seed flavonoids. Annu Rev Plant Biol. 2006, 57: 405-30.
Yoshida K, Iwasaka R, Kaneko T, Sato S, Tabata S, Sakuta M: Functional Differentiation of Lotus japonicus TT2s, R2R3-MYB Transcription Factors Comprising a Multigene Family. Plant and Cell Physiology. 2008, 49: 157-169.
eFP Browser. [http://www.bar.utoronto.ca/efp]
Schmid M, Davison TS, Henz SR, Pape UJ, Demar M, Vingron M, Schölkopf B, Weigel D, Lohmann JU: A gene expression map of Arabidopsis thaliana development. Nat Genet. 2005, 37 (5): 501-6.
Ogata K, Morikawa S, Nakamura H, Hojo H, Yoshimura S, Zhang R, Aimoto S, Ametani Y, Hirata Z, Sarai A: Comparison of the free and DNA-complexed forms of the DNA-binding domain from c-Myb. Nat Struct Biol. 1995, 2: 309-320.
Heine GF, Hernandez JM, Grotewold E: Two cysteines in plant R2R3 MYB domains participate in REDOX-dependent DNA binding. J Biol Chem. 2004, 279 (36): 37878-85.
Quattrocchio F, Baudry A, Lepiniec L, Grotewold E: The regulation of flavonoid biosynthesis. The Science of flavonoids. Edited by: Grotewold E. 2006, New York. Springer Science+Business Media, 97-122.
Menges M, Pavesi G, Morandini P, Bögre L, Murray JA: Genomic organization and evolutionary conservation of plant d-type cyclins. Plant Physiol. 2007, 145 (4): 1558-76.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007
Reid K, Olsson N, Schlosser J, Peng F, Lund S: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6: 27-
Poupin MJ, Federici F, Medina C, Matus JT, Timmermann T, Arce-Johnson P: Isolation of the three grape sub-lineages of B-class MADS-box TM6, PISTILLATA and APETALA3 genes which are differentially expressed during flower and fruit development. Gene. 2007, 404 (1–2): 10-24.
Schwede T, Kopp J, Guex N, Peitsch MC: SWISSMODEL: an automated protein homology-modelling server. Nucleic Acids Res. 2003, 31: 3381-3385.
Thanks to Diego Gutiérrez for his contribution in gene isolation. This work was supported by the Chilean Wine Consortium 05CTE01-03, the Fruit Consortium, 07Genoma01, Millennium Nucleus for Plant Functional Genomics (P06-009-F) and by Fellowships awarded to JTM (MECESUP and CONICYT AT24060171).
JTM carried out the genome search, sequence, phylogenetic, expression and structure analyses and drafted the manuscript. FA contributed to phylogenetic analysis, the comparative modelling and together with JTM contributed to the conception of the study. PA–J was involved in revising the manuscript critically for important intellectual content and gave final approval of the version to be published. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Multiple alignment of 20 representative R2R3 MYB domains from Arabidopsis and six MYB domains from characterised grape MYB genes. Identical amino acid residues are shaded in yellow and the blue and white boxes indicate the extent of the R2 and R3 repeats. The consensus sequence shown under the alignment was used to search for MYB homologues in the Grape Genome. (TIFF 1 MB)
Additional file 4: Complete phylogenetic tree of Vitis and Arabidopsis MYB proteins, using the DNA-binding domain (A) or full protein sequences (B). The Neighbour Joining (NJ) tree method was used. Numbers above nodes represent bootstrap values for 2000 replicates. Asterisks indicate Vitis-specific gene pairs and red and blue letters next to gene model identifiers refer to Genoscope orthologue and Blastp homologue predictions, respectively. (TIFF 1 MB)
Additional file 5: Phylogenetic tree of the Vitis , Arabidopsis and Rice R2R3 MYB Subfamily. The R2R3 DNA-binding domain sequences were used for the construction of a parsimony phylogeny tree with Mega4 software using the Neighbour Joining (NJ) tree method. Numbers above nodes represent bootstrap values for 2000 replicates. (TIFF 493 KB)
Additional file 7: Affymetrix (ATH1 Gene Chip) MYB expression data during Arabidopsis organ development from the putative homologues of the selected grape genes. Data was collected using the Arabidopsis Electronic Fluorescent Pictograph (eFP) Browser. Colours refer to an absolute expression unit, calculated independently for each gene and normalised by the RMA or GCOS method (eFP Browser by Vinegar, drawn by J.Alls and N. Provart. Data from Gene Expression Map of Arabidopsis Development: [66, 67], and the Nambara lab for the imbibed and dry seed stages). (PNG )
Additional file 8: Three dimensional modelling of grape MYB R2R3 domains. A) Vitis and Arabidopsis genes were aligned with the Ratus c-MYB R2R3 domain using the ClustalW module (VECTOR NTI, Invitrogen). Identical amino acid residues are shaded in yellow, and bars indicate the positions of the three α-helices from each repeat. The six asterisks indicate the constantly spaced tryptophan residues. The arrowhead indicates residues previously described as necessary or dependent for MYB-bHLH interaction. B) Ribbon diagram of the Ratus c-MYB R2R3 domain. The DNA-interaction interface is shown (left panel), in which basic residues (coloured in yellow) are exposed outside the R2-α3, R3-α2 and R3-α3 helices (middle panel). Left and middle panels show R2R3 domain in the backwise orientation, as indicated by a black arrow showing the initial R2 residue. C) Charge distributional ribbon diagrams of the R2R3 domain from grape MYBCS-1, MYBC2 repressor like and MYB14 like proteins, obtained by comparative modelling. Residues in volume correspond to cys (green), arg (orange), leu (brown), lys (yellow), gly (blue), thr (cyan), glu (magenta) and ala (grey). Non-volume residues in yellow correspond to positively charged amino acids (arg, lys, his). (TIFF 384 KB)
Authors’ original submitted files for images
About this article
Cite this article
Matus, J.T., 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 8, 83 (2008). https://0-doi-org.brum.beds.ac.uk/10.1186/1471-2229-8-83
- Grape Genome
- Wide Genome Duplication
- R2R3 Domain
- Vitis Genome
- Epidermal Cell Fate