Impact Factor 4.106 | CiteScore 4.47
More on impact ›

Original Research ARTICLE

Front. Plant Sci., 08 October 2019 |

Genome-Wide Analysis of the YABBY Gene Family in Grapevine and Functional Characterization of VvYABBY4

Songlin Zhang1,2†, Li Wang1,3†, Xiaomeng Sun1,2, Yunduan Li1,2, Jin Yao1,2, Steve van Nocker4 and Xiping Wang1,2*
  • 1State Key Laboratory of Crop Stress Biology in Arid Areas, College of Horticulture, Northwest A&F University, Yangling, China
  • 2Key Laboratory of Horticultural Plant Biology and Germplasm Innovation in Northwest China, Ministry of Agriculture, Northwest A&F University, Yangling, China
  • 3College of Horticulture, Agricultural University of Hebei, Baoding, China
  • 4Department of Horticulture, Michigan State University, East Lansing, MI, United States

Genes of the plant-specific YABBY transcription factor family have various roles, including lateral organ development, establishment of dorsoventral polarity, and response to abiotic stress. In this study, we carried out a genomic census of YABBY genes in grapevine (Vitis vinifera) and characterized their expression pattern during ovule development. We identified seven YABBY genes and classified them into five subfamilies, based on peptide sequence, similarity of exon–intron structure and composition of peptide sequence motifs. Analysis of YABBY gene expression in various grapevine structures and organs suggested that these genes function in diverse aspects of development and physiology. Analysis of expression during ovule development in four cultivars showed that one gene, VvYABBY4, was preferentially expressed during the period of ovule abortion in seedless cultivars. Transgenic expression of VvYABBY4 in tomato conferred reduced plant stature, dark green leaves, elongated pistil, and reduced size of fruit and seeds. Reduced seed size was associated with smaller endosperm cells. Expression of VvYABBY4 also affected expression of numerous tomato genes with presumed roles in seed development. These data suggest the potential for VvYABBY4 to influence seed development in grapevine, which may impact seedless grape breeding.


The plant-specific YABBY transcription factors play important roles in the establishment of adaxial–abaxial polarity (Kumaran et al., 2002), development of lateral organs (Bowman et al., 1989), and stress response (Zhao et al., 2017). YABBY proteins comprise two highly conserved DNA-binding domains: a zinc finger-like domain (C2C2) and helix–loop–helix domain (YABBY) (Jang et al., 2004). The genome of Arabidopsis thaliana (Arabidopsis) encodes for six YABBY members, which have been classified into five subfamilies, designated FIL/YAB3, CRC, INO, YAB2, and YAB5. FILAMENTOUS FLOWER (FIL), YAB3, YAB2, and YAB5 act redundantly to promote lateral organ development. CRABS CLAW (CRC) participates in establishing polarity of developing carpels and nectaries (Siegfried et al., 1999). INNER NO OUTER (INO) plays a seemingly distinct role to promote development of the ovule outer integument, the layer of cells surrounding the nucleus, into the seed coat (Villanueva et al., 1999; Bowman, 2000).

Based on their critical functions in development, the genomic organization and expression of YABBY genes have been studied widely in plants. In maize, there are 13 YABBY genes, and at least a subset of these is strongly expressed in the shoot meristem and leaf primordia, suggesting a function in organ formation and polarity (Sawa et al., 1999; Siegfried et al., 1999). Rice contains eight YABBY genes, including the well-studied DROOPING LEAF (DL) (Tanaka et al., 2017). In a loss of function dl mutant, carpels are transformed to stamens, and formation of leaf veins is defective (Yamada et al., 2004; Tanaka and Hirano, 2012). DL is orthologous with CRC, and other members of the CRC subfamily have been reported to be associated with polarity (Nagasawa et al., 2003; Ohmori et al., 2008). Tomato contains nine YABBY genes, and Cong found that FASCIATED (SlYABBY2b) mainly controls the formation of carpel during the flowering and fruit development of tomato (Cong et al., 2008). SlYABBY1a was highly expressed in the epidermis, and it is speculated that the SlYABBY1a gene may affect the development of distal pericarp cells and play an important role in the regulation of fruit ripening (Bartley and Ishida, 2003; Bartley and Ishida, 2007). Previously, Xiang et al. cloned two genes from the wild Chinese grapevine species, Vitis pseudoreticulata, and studied their function by heterologous expression in Arabidopsis. They found that the expression of VpYABBY1 caused the loss of dorso-ventral polarity in the leaf blade, similar to previous observations for FIL and YABBY3 (Siegfried et al., 1999; Xiang et al., 2013).

Most reports on YABBY genes have examined potential functions in symmetry of shoot apical and floral meristems in Arabidopsis (Siegfried et al., 1999; Kumaran et al., 2002). By contrast, there have been few previous investigations into the YABBY genes in common grapevine (Vitis vinifera L.). In this study, we investigated the genomic organization and expression of the YABBY gene family in grapevine, focusing on those genes that may have a role in seed development or abortion. Our results have potential significance for molecular breeding programs of seedless grapes.

Materials and Methods

Plant Material

Four grapevine cultivars, including two seedless cultivars, Flame Seedless (V. vinifera) and Thompson Seedless (V. vinifera), and two seeded cultivars, Red Globe (V. vinifera) and Kyoho (V. vinifera × Vitis labrusca), were used in this study. Plants were maintained in the grape germplasm resource collection of Northwest A&F University, Yangling, China (34°20′N108°24′E). Roots, stems, leaves, flowers (at full bloom stage), tendrils, and fruit (42 days postanthesis) were collected from 4-year-old “Red Globe” and “Thompson Seedless” plants. As noted previously (Li et al., 2017b), fertilized ovules in seedless cultivars typically abort at 30 to 40 days after full bloom (DAF). Consequently, ovules and seeds were dissected from developing fruit at 27 DAF, 30 DAF, 33 DAF, 36 DAF, 39 DAF, and 42 DAF. Samples were frozen immediately in liquid nitrogen and stored at −80°C for analysis of gene expression.

Genetically modified tomato (genotype micro-Tom) and the harvested seeds (T1) were cultivated in a greenhouse under 14-h light/10-h dark photoperiods, 25°C day/20°C night temperatures, and relative air humidity of 60%. Plant nutrient solution was applied once a week.

Identification of YABBY Genes in the Grapevine Genome

The YABBY signature domain (PF04690) was identified with SMART ( (Lu et al., 2017). The Hidden Markov Model (HMM) profile of the YABBY DNAbinding domain (PF04690) was obtained from the Pfam protein families database ( (Punta et al., 2012). Candidate genes were obtained from the 12 coverage assembly of the V. vinifera PN40024 genome (12×; by using HMM online software (Eddy, 1998; Finn et al., 2015). Open-reading frames (ORFs) were obtained by ORF finder ( Detailed information, including gene ID, protein sequence, and coding sequence, was retrieved from the Grape Genome Database and NCBI GenBank. The number of amino acids and isoelectric point were obtained by ExPASy (

Multiple-Sequence Alignment and Phylogenetic Analyses

Multiple alignment of YABBY protein sequences from grapevine, Arabidopsis (Bowman, 2000), tomato (Huang et al., 2013), rice (Toriba et al., 2007), and maize (Juarez et al., 2004) were performed using ClustalX 2.1 (Thompson et al., 2002). Phylogenetic tree was constructed with the neighbor-joining (NJ) method using MEGA 5.0 (Stracke et al., 2001). Bootstrap analysis was performed using 1000 replicates with the following parameters: “p distance,” “complete deletion,” and “gap setting” (Sun et al., 2017).

Distribution of Conserved Domains and Exon–Intron Structure Analysis of VvYABBY Genes

Conserved sequence motifs in YABBY genes were identified using MEME 4.11.2 ( et al., 2006); protein domains were identified using SMART ( et al., 2017a; Zhang et al., 2017). Exon–intron structures were determined according to the alignments of their transcribed sequences with corresponding genomic sequences, and the diagram was generated with the online Gene Structure Display Server ( et al., 2015).

Determination of Chromosomal Localization and Synteny Analysis

Genomic positions of the YABBY genes were as previously annotated in the Grape Genome Database (12×; Tandemly duplicated genes were defined as adjacent to homologous YABBY genes, with no more than one intervening gene (Holub, 2001). Syntenic blocks within the grape genome, as well as between grape and Arabidopsis, were identified from the Plant Genome Duplication Database ( (Zhang et al., 2012b).

Gene Expression Analysis

Total RNA was isolated from 100 to 200 mg of all the samples using an EZNA Plant RNA Kit (R6827-01; Omega Bio-Tek, Norcross, GA, USA), according to the manufacturer’s instructions. Agarose gel electrophoresis was used to check the integrity of RNA. Moreover, total RNA was quantified with Nanodrop8000 Spectrophotometer (Thermo Scientific, Waltham, MA, USA). One microgram of total RNA was used to synthesize first-strand cDNA using the Prime Script RT reagent Kit (TaKaRa Biotechnology, Dalian, China) in a final volume of 20 μl. The products were diluted sixfold for further experiments.

Quantitative real-time PCR analysis was performed with an IQ5 real-time PCR machine (Bio-Rad, Hercules, CA, USA). Each reaction was conducted in triplicate with a final volume of 20 μl, including 1.6-μl primer (1.0 μM), 1.0-μl cDNA template, 10.0-μl SYBR green (TaKaRa, Bio Inc.), and 7.4-μl sterile distilled H2O. The PCR amplification profile used was following a previous study (Wang et al., 2015). Relative expression levels were analyzed using IQ5 software and the normalized expression method (Guo et al., 2014). The sequence of the gene-specific primers used for each of the YABBY genes in grapevine, as well as two house-keeping VvACTIN (GenBank Accession number AY680701) and VvEF1-α (GenBank Accession number EC931777) gene, is listed in Table S1. The CT values of the VvACTIN and VvEF1-α were geometrically averaged as the new CT (internal control gene). The relative transcript level of a gene of interest is calculated as 2−ΔΔCT [ΔΔCT = ΔCT (candidated sample) − ΔCT (control sample); ΔCT = CT (interesting gene) − CT (internal control gene)] (Guo et al., 2018).

Subcellular Localization and Transactivation Assay of VvYABBY4

The coding region of VvYABBY4 gene without terminator code was amplified and connected with pBI221–GFP vector. The recombinant plasmid pBI221–VvYABBY4–GFP and empty vector were transformed into Agrobacterium tumefaciens strain GV3101. The leaves of 5-week-old tobacco (Nicotiana benthamiana) were injected with bacteria solution. After 72-h treatment, leaves were observed under a Nikon A1R/A1 confocal microscope (Nikon, Tokyo, Japan). GFP was excited at 514 nm. DAPI was excited at 405 nm (Fang et al., 2016).

Yeast GAL4 system was used for transactivation assay. Since p53 and large T-antigen are known to interact in a yeast two-hybrid assay (Zhu, 1997; Li and Fields, 1993), it can be used as a positive control. A negative control should also be performed using pGBKT7-Lam (which encodes the Gal4 BD fused with lamin) and pGADT7-T. Then, select the colonies as the positive and negative control, respectively. These diploids are useful as reference strains for checking new batches of growth media (pGBKT7-53–positive control bait plasmid, pGBKT7-Lam–negative control bait plasmid, and pGADT7-T–positive control prey plasmid purchased from TAKARA company, 630489).The full-length ORF of VvYABBY4 was cloned into the pGBKT7 vector with Sma I and Sal I cleavage sites. The empty pGBKT7 vector and recombinant vector pGBKT7–VvYABBY4 were transformed into the Y2H yeast strain for the transcriptional activation analysis. These constructs were painted on the same SD/Trp−, SD/Trp−/X-α-gal or SD/Trp−/X-α-gal/AbA medium plates at 30°C for 3 to 4 days. (Reference Matchmaker™ Gold Yeast Two-Hybrid System User Manual; clontech, TAKARA).

Molecular Cloning of the VvYBBY4 Gene and Expression in Transgenic Tomato

The VvYABBY4 cDNA was engineered into the pCAMBIA2300 vector under control of the CaMV 35S promoter, making use of the primers shown in Table S2. The resulting pCAMBIA2300–VvYABBY4 plasmid was introduced into Agrobacterium strain GV3101. Transformation was carried out using the leaf disc method (Mccormick et al., 1986). Primers were designed according to the sequence of pCAMBIA2300 vectors and used for PCR identification of transgenic lines.

Analysis of Transgenic Plants

Pollen viability and germination were evaluated in accordance with the method used in previous study (Zhang et al., 2012a; Pham et al., 2015). Pollen morphology was observed under a scanning electron microscope (SEM). Sample preparation for SEM was carried out according to previous study (Xiang et al., 2013). Cell measurements were performed using Image-pro software.


Identification of YABBY Gene Family Members in V. Vinifera

Based on the presence of an ORF encoding a single, conserved YABBY domain, we identified seven YABBY genes in a draft genome sequence of grapevine. We designated these genes V. vinifera YABBY1YABBY7 according to their consecutive chromosomal positions. The parameters used to describe the VvYABBY proteins are listed in Table 1 and included gene locus identifiers, gene accession numbers, chromosomal position, pI, and length of the coding sequence. The deduced length of VvYABBY proteins ranged from a minimum of 160 amino acids (VvYABBY1) to a maximum of 211 residues (VvYABBY3), whereas the isoelectric point (pI) values ranged from 7.03 (VvYABBY2) to 8.98 (VvYABBY1). This variability suggests that different VvYABBY proteins might play a part in different microenvironments.


Table 1 Information of grape YABBY genes/proteins.

Multiple Sequence Alignment and Structural Analysis of YABBY Genes in Grape

A sequence alignment of the YABBY proteins is shown in Figure 1A. The ∼60-amino acid YABBY domain is situated near the carboxyl terminus, similar to YABBY proteins in other plants where characterized. All grape YABBYs also contain a highly conserved C2C2 zinc finger domain near the amino terminus. Within this domain, the spatial configuration of cysteine (C) and histidine (H) residues directly involved with Zn2+ binding are conserved (C-X2-C-X20-C-X1-HC). Analysis of additional conserved segments within the VvYABBY1–VvYABBY7 protein sequences revealed a short domain bordering the YABBY domain in VvYABBY3, 4, 5, 6, and 7, and a carboxyl-terminal domain in VvYABBY3, 4, 5, and 7 (Figure 1B). Consistent with the similarity in domain composition, a phylogeny based on peptide sequence revealed a close relationship between VvYABBY4 and VvYABBY5, and VvYABBY3 and VvYABBY7 (Figures 1B, C). This was further supported by analysis of the exon–intron structures of the seven genes, which identified similar exon partitions between VvYABBY4 and VvYABBY5, and VvYABBY3 and VvYABBY7 (Figure 1D).


Figure 1 YABBY genes in grapevine. (A) Alignment of open-reading frame translations from the Vitis vinifera YABBY genes. The conserved C2C2 and YABBY domains are marked. (B) Phylogenetic relationships among the ORF translations of the seven grapevine YABBY genes. Numbers near the nodes represent bootstrap values. (C) Motif analysis of grapevine YABBY ORF translations. Motifs were designated as 1–4 and distinguished by different colors. (D) Exon–intron structure of grapevine YABBY genes. Exons are indicated as yellow boxes. Black lines linking two exons represent introns.

Phylogenetic and Synteny Analyses of YABBY Genes

To better understand the evolutionary relationships among YABBY genes, a phylogenetic tree was constructed with a total of 43 YABBY genes from different species (Figure 2A). As important model plants, the functions of Arabidopsis and tomato YABBY genes have been well discussed, allowing inference of function for phylogenetically related proteins from grape. This analysis divided the grapevine YABBY genes into five well-defined clades or subfamilies, consistent with previous analysis in Arabidopsis and tomato: 1) FIL/YAB3, 2) YAB2, 3) YAB5, 4) INO, and 5) CRC. It is evident that YABBY genes of the monocotyledon (rice and maize) clustered together. As expected, grape YABBY genes showed a closer relationship to dicotyledons (Arabidopsis and tomato). In clade YAB5, no genes of monocotyledon species was included. In addition, as a distinct branch, ZmYABBY10 and OsYABBY1 are highly homologous but not homologous with any dicotyledonous plants. It shows that the YABBY gene is not conservative in the evolution of species, and they have diversified during their evolution. The clear orthology observed within each of these clades among the five species suggests that these subfamilies were distinguished prior to the evolutionary separation of these plants.


Figure 2 Phylogenetic relationships and synteny analysis of grapevine YABBY genes. (A) Phylogenetic analysis of YABBY genes from grapevine and other species. Yellow circles represent Arabidopsis thaliana, red circles represent Solanum lycopersicum, green circles represent Vitis vinifera. Blue and purple circles denote Zea mays and Oryza sativa, respectively. Numbers at the nodes indicate bootstrap values; values lower than 30% are not shown. (B) Synteny analysis of YABBY genes between grape and Arabidopsis. Colored bars linking two chromosomal parts represent syntenic regions. Chr, chromosomes.

As shown in Figure 2B, two pairs of grape YABBY genes (VvYABBY3/VvYABBY7 and VvYABBY4/VvYABBY5) are associated with segmental duplication (Table S3). Moreover, studies of syntenic relationships often allow transferring available genomic information from a model species to a less-studied species (Lyons et al., 2008). Some of the predicted YABBY domains from grapevine and Arabidopsis showed pairwise relationships, such as those between VvYABBY1AtCRC and VvYABBY2AtINO (Table S4), and this is consistent with the grouping in Figure 2A. This suggested that they came from a common ancestor and duplicated prior to the divergence of grape and Arabidopsis.

Expression Profiles of YABBY Genes in Diverse Structures

As the first step to define function of the grapevine YABBY genes, we evaluated their expression levels in stems, young leaves, roots, flowers, tendrils, and fruit. In the seeded cultivar “Red Globe” and seedless cultivar “Thompson Seedless”,: all of the genes showed relatively high expression in the flower (Figure 3). VvYABBY3, 4, 5, 6, and 7 were also relatively highly expressed in the leaf. In addition, VvYABBY4, 5, and 7 showed relatively high expression in fruit. VvYABBY1, 3, and 7 were expressed relatively strongly in root. Comparing these results with our phylogenetic and synteny analyses, it is interesting that some homologous genes (VvYABBY3 and VvYABBY7, for example) displayed a similar expression pattern in all structures.


Figure 3 Analysis of grape YABBY genes expression in various grapevine structures by quantitative PCR. “T.S.” stand for “Thompson Seedless” and “R.G.” represent “Red Globe.” Using VvACTIN and VvEF1-α genes followed by geometric averaging for determining a normalization factor. The means ± SD of three biological replicates are presented. Asterisks represent statistical significance (*P < 0.05, one-way ANOVA).

Expression Patterns of Grape YABBY Genes During Ovule Development Among Seeded and Seedless Cultivars

The expression of YABBY genes in four grapevine cultivars, including two seedless cultivars, Flame Seedless (V. vinifera) and Thompson Seedless (V. vinifera), and two seeded grapevines, Red Globe (V. vinifera) and Kyoho (V. vinifera × V. labrusca), are shown in Figure 4. Some of the genes, like VvYABBY1, 2, and 6, exhibited no striking differences in expression pattern during ovule development between seedless and seeded grape cultivars. On closer inspection, however, it becomes clear that VvYABBY3/VvYABBY7 and VvYABBY4/VvYABBY5 also had the same pattern of expression during the development of ovule. In addition, we can conclude that YABBY genes participate in the development of ovules, because all of them were expressed during ovule growth.


Figure 4 Expression pattern of grape YABBY genes in different stages of ovule development in four grapevine cultivars. Seedless cultivars: Thompson Seedless (T.S) and Flame Seedless (F). Seeded cultivars: Red Globe (R.G.) and Kyoho (K). The value of the x - axis represents the days after full bloom. Using VvACTIN and VvEF1-α genes followed by geometric averaging for determining a normalization factor. The means ± SD of three biological replicates are presented. Values with different lowercase superscripts mean significant difference (P < 0.05).

It is noteworthy that VvYABBY4 and VvYABBY5 were expressed to relatively high levels in the two seedless cultivars. According to previous studies, ovule abortion usually occurs between 30 and 40 DAF (Li et al., 2017b). During this period, the level of VvYABBY4 gene expression increased with the beginning of ovule abortion and then decreased. However, the expression level of VvYABBY5 was more stable compared with VvYABBY4. It suggests that ovule abortion begins with the up-regulation of VvYABBY4 gene. This hypothesis, though, needs further proof.

VvYABBY4 Is Positioned in the Nucleus and Functions as a Transcriptional Activator

Recombinant plasmid pBI221–VvYABBY4–GFP was transfected into tobacco leaf. In the case of DAPI, the fluorescence signal was alone detected in the nucleus; therefore, the pBI221–VvYABBY4–GFP fusion protein was localized to the nucleus (Figure 5A). To analyze the transcriptional activation of VvYABBY4, full-length VvYABBY4 was inserted into the pGBKT7 vector. Results showed that yeast transformed with pGBKT7-VvYABBY4 had activation ability and presented β-galactosidase activity (Figure 5B). These results confirm that VvYABBY4 is an active transcription activator which localized to the nucleus.


Figure 5 Subcellular localization and transactivation analysis of VvYABBY4. (A) Subcellular localization of the 35S::VvYABBY4–eGFP fusion protein in Nicotiana benthamiana leaves. Scale bar: 20 μm. (B). Transactivation activity analysis of full length VvYABBY4 in yeast.

Characterization of the Tomato Expressing VvYABBY4

To further define a functional role for the grapevine YABBY genes in development, we engineered and characterized transgenic tomato constitutively expressing VvYABBY4. PCR identification showed that the VvYABBY4 gene has been successfully transformed in the tomato (Figure S1). Compared with non-transgenic plants, transgenic lines were shorter and displayed reduced apical dominance (Figure 6A). The leaves of transgenic lines were also darker green compared with the control (Figures 6B, F). Microscopic observation revealed that leaf mesophyll cells of transgenic plants were present at a higher density, with the increased concentration of chloroplasts potentially contributing to the dark green color. Moreover, mesophyll cells of transgenic plants displayed significantly enhanced polarity (Figures 6C, G). In the adaxial region, long square columnar palisade tissue cells were closely arranged, with narrow intercellular space. In contrast, in the abaxial region, the spongy tissue cells displayed a spherical shape. In addition, a closer view of the transgenic leaves revealed obvious long hair on the surface, which rarely appeared on leaves of non-transgenic control plants.


Figure 6 Phenotypic analysis of effects of VvYABBY4 expression in transgenic tomato. (A) Observation of non-transgenic control plant (NT; left) and representatives of two VvYABBY4 expressing transgenic lines (OE-1 and OE-2). (B) Leaf from non-transgenic plant. (C). Longitudinal sections of non-transgenic plant leaf. (D) Flower from non-transgenic plant. (E). High magnification of the non-transgenic plant pistil. (FI) The observation of VvYABBY4 transgenic plant corresponding to (B, C, D, E), respectively. The red arrow indicates the protruding part of the pistil. Bars:(A), 1 cm (B, D, F, H), 10 mm (C, G), 200 μm (E, I), 1 mm.

Finally, 35S::VvYABBY4 transgenic plants had larger flowers (Figures 6D, H), and the pistil was longer than non-transgenic control plants at anthesis (Figures 6E, I). Besides that, there were no obvious difference between transgenic and control plants in other respects, such as bud, sepal, or stamen development.

Pollen Viability, Germination Percentage, and Pollen Morphology

Pollen viability in transgenic plants was similar (∼85%) to non-transgenic plants, as evaluated using TTC staining (Figures 7A, D). In addition, transgenic and non-transgenic plants exhibited similar pollen size and shape. However, we found that the mature pollen grains from transgenic plants exhibited irregularly shaped protuberances attached to the exine (Figure 7F). In contrast, the pollen of non-transgenic plants appeared glabrous (Figure 7C). We collected pollen at the same time from transgenic and control plants and cultured the pollen under identical conditions. As shown in Figures 7B, E, the germination of pollen from 35S::VvYABBY4 plants was slower than the control.


Figure 7 Experimentation of pollen energy, pollen germination and pollen grains. (A, D). Pollen viability of non-transgenic control and 35S::VvYABBY4 transgenic plant using TTC staining 2 days before full-bloom respectively. (B, E). Pollen germination of non-transgenic control and 35S::VvYABBY4 transgenic plant using aniline blue staining respectively. (C, F). Observation of non-transgenic control and 35S::VvYABBY4 transgenic plant pollen under scanning electron microscope (SEM). Bars: (A, B, D, E), 100 μm. (C, F), 30 μm.

Effect of VvYABBY4 Expression on Transgenic Tomato Fruit and Seed Development and Morphology

Compared with non-transgenic control plants, transgenic plants, constitutively expressing VvYABBY4, exhibited slightly reduced fruit size, approximately 90.3% that of the control (Figures 8A,E). Relative to the seeds of non-transgenic plants, which were plump and wide, the seeds of transgenic lines were significantly smaller (Figure 8B). The hundred grain weight of seed from transgenic plants was reduced to 78.2% that of the controls (Figure 8G). However, the number of seeds per fruit was similar to controls (Figure 8F).


Figure 8 Effect of VvYABBY4 expression on tomato fruit and seed development and morphology. (A) Fruit from non-transgenic control and transgenic plant (from top to bottom), “d” represents diameter. (B) Seeds from non-transgenic control and transgenic plant (from top to bottom).(C) and (D). Histological observation of non-transgenic control and over-expression of VvYABBY4 plants seeds. (E) Average fruit diameters of non-transgenic control (NT) and over-expression of VvYABBY4 lines (OE-1 and OE-2). (F) Number of seeds per fruit in non-transgenic control (NT) and over-expression of VvYABBY4 lines (OE-1 and OE-2). (G) Hundred grain weight of non-transgenic control (NT) and over-expression of VvYABBY4 lines (OE-1 and OE-2). Values in (E-G) are given as means ± SEs, *P < 0.05, compared with the non-transgenic control. Three independent repeats were set with five fruits in each group and calculated the average of the three sets of data, total 15 fruits were used in one line test. The diameter of each single fruit was measured three times. Bars: (A, B), 1 mm (C, D), 20 μm.

To obtain clues about the potential influence of VvYABBY4 on seed development, we observed the seeds during period of mature green using microscopy (Figures 8C, D). Histological longitudinal section of transgenic plants seeds showed that endosperm cells were more closely arranged, and were narrower than those of nontransgenic plants. A potential explanation is that the development of the endosperm cells limits the seed size.

Effect of VvYABBY4 Expression in Transgenic Tomato on Genes Related to Seed and Fruit Development

To index growth and development of the transgenic tomato seed and fruit at the molecular level, we analyzed the expression of four previously characterized tomato genes involved in seed or fruit growth and development. Tomato FW2.2 encodes an unknown functional transmembrane protein and interacts with CKII kinase to form a cell cycle regulatory complex that inhibits cell division and participates in cell cycle signaling pathways (Mei and Simmons, 2011; Chakrabarti et al., 2013). Tomato WEE1 participates in the control of cell size during fruit development through the negative regulation of CDK activity, resulted in a changes of plant and fruit size (Nathalie et al., 2004; Nathalie et al., 2007). As shown in Figure 9, the expression patterns of SlWEE1 in transgenic plants are different from SlFW2.2. At all the stages of fruit development, the expression of SlFW2.2 continues to be highly expressed compared with the control. However, the expression of SlWEE1 has no obvious difference compared with non-transgenic plant during fruit development. This result indicated that the reduced fruit size of transgenic plant was not caused by the endoreduplication which was regulated by SlWEE1, but by the high expression of SlFW2.2 caused by VvYABBY4, thus negatively regulating the fruit size.


Figure 9 Expression of tomato genes linked with seed and fruit growth in VvYABBY4 transgenic plants. Two independent lines (line 1 and line 2) were evaluated. Relative transcript accumulation of related genes in tomato fruits and seeds at different fruit development stages. Transcripts were normalized to the expression of the SlACTIN1 gene. dpa, days postanthesis; F, fruit; S, seed; The means ± SD of three biological replicates are presented. Asterisks represent statistical significance (*P < 0.05, one-way ANOVA).

As for seed size, we selected two genes that are involved in the development of tomato seeds. Tomato TM29 functions in floral meristem, flower, and seed development (Sutherland et al., 2002). Tomato FUL1 (SlFUL1, formerly known as TDR4) is expressed during ripening (Busi et al., 2003), where it interacts with RIN to mediate ethylene-independent aspects of fruit and seed ripening (Bemer et al., 2012). The expression of these two genes was assessed at three different seed developmental stages. As shown in Figure 9, at each stage, the expression of the SlFUL1 in transgenic plants was higher than the control. SlTM29 is highly expressed in the seeds of 30 days postanthesis compared with the controls, and the amount of expression dropped with the maturity of seeds. It has been reported that high expression of SlFUL1 produced a smaller and reduced number of tomato seeds (Lozano et al., 1998; Czerednik et al., 2015) and decreased the expression of SlTM29 conferred seedlessness in tomato (Sutherland et al., 2002). Our results are consistent with these observations and suggest that VvYABBY4 gene may affect the expression of genes associated with seed growth to co-regulate seed development.


The YABBY family of plant-specific transcription factor plays an important role in abaxial–adaxial polarity and development of lateral organs. In addition, YABBY genes have been implicated as important components of plant response to abiotic stress. In this study, we identified seven YABBY genes in grapevine. Genomic studies have documented the existence of six YABBY members in Arabidopsis, nine in tomato, eight in rice, 17 in soybean, and 13 in maize. The relatively few YABBY genes found in grapevine could reflect the lack of expansion of this family during evolution.

To understand grape YABBY evolutionary relationships, we constructed a phylogenetic tree, including YABBY proteins from five species, including monocotyledons and dicotyledons. In all these plants, YABBY members were separated into five subfamilies. Genes in the same subfamily, such as VvYABBY3 and VvYABBY7, or VvYABBY4 and VvYABBY5, show similar motif and exon–intron structure. This apparent recent evolution implies that these gene pairs may perform a similar function. On the other hand, evolutionary differentiation of gene structure implies distinction of gene function(Wan et al., 2018). For example, in Arabidopsis, CRC is involved in the formation of nectaries and carpels, but INO is related to the outer integument development (Bowman, 2000).

Gene duplication and divergence events have been suggested to be the main contributors to evolutionary momentum (Ohno et al., 2010). We did not find any tandemly repeated gene events for the YABBY family. However, we identified two pairs of segmental duplications. Due to the genome of grapevine that originated from a hexaploidization (Jaillon et al., 2007), several synteny blocks would be expected to exist in this genome. Segmental duplicates may be more often retained due to subfunctionalization, without increasing the likelihood of gene rearrangement (Lynch and Conery, 2000; Lynch and Conery, 2003). Meanwhile, a total of two pairs of homologous genes were identified between grape and Arabidopsis, suggesting that they were derived from the same ancestor gene and became distinguished prior to species differentiation.

The function of Arabidopsis YABBY genes is related to development of abaxial tissues in lateral organs. For example, FIL, YAB2, and YAB3 are expressed in cotyledons, leaves, flowers, and other lateral organs, whereas CRC is expressed in carpels. Similarly, YABBY genes in grapevine are widely expressed in many structures. Specifically, the relative expression of grape YABBY genes in floral organs is more obvious, suggesting that these YABBY genes are involved in the development of the flower. The observation that YABBY genes in the same subfamily showed similar patterns of expression, such as VvYABBY3 and VvYABBY7, suggests that they perform similar functions. In our study, to complement previous research, we focused on identifying genes involved in seed development, as well as ovule abortion. Thus, we documented the expression of the seven genes in different stages of ovule development, in two seeded (Red Globe, Kyoho) and two seedless (Thompson Seedless, Flame Seedless) cultivars.

Quantitative results indicate that the expression of VvYABBY4 and VvYABBY5 showed significant differences between seeded and seedless grape varieties, this is also consistent with the published transcriptome data (Table S5) (Wang et al., 2016). And it is certain that due to different grape materials and sampling periods, there are some differences between us, just like in his results, there was no expression of VvYABBY1 or 2 during the three stages. However, VvYABBY3, 4, 5, and 7 are significantly differentially expressed in seedless and seedless grapes. And these results are consistent with each other. Recent reports indicated that the key period of grape ovule abortion falls between 30 and 40 DAF (Wang et al., 2016; Li et al., 2017b). During this time, the expression level of VvYABBY5 was not influenced by ovule abortion, although its expression in seedless grape varieties is excellent. By contrast, the level of VvYABBY4 expression increased with the beginning of ovule abortion and then decreased with the end of procedure, and so we speculate that both VvYABBY4 and VvYABBY5 play a role in seed development, but that VvYABBY4 plays an additional role in seed abortion.

It is worth noting that considering the sequence divergencies in four varieties (Genova et al., 2014), we cloned these seven YABBY genes, and the results showed that there was no difference in VvYABBY1, 2, 3, 4, 6 sequences among the four cultivars (Figure S2). In addition, two points were found to be different from the reference genome in VvYABBY5, but they were exactly the same among the four species, which we doubt the accuracy of the reference genome at this site. The sequence alignment of YABBY7 in four varieties shows that YABBY7 has the same sequence in Kyoho and Pinot noir (reference genome), but two points are different from the other three varieties. So we designed the quantitative primer based on the sequencing results.

As for housekeeping genes in grape, grape berries or seeds undergo significant metabolic changes throughout their development, orchestrated in part by the up- and down-regulation of transcripts (Coombe, 1992). The ability to identify transcripts that are resistant to growth fluctuations or stresses is challenging, so to date, no single candidate housekeeping genes has been shown to be universally acceptable. In 2006, Reid demonstrated that the most stable reference genes: actin, GAPDH, EF1-α, and SAND ranked among the top genes when data from three independent statistical approaches were evaluated during grape development (Gonçalves et al., 2005; Reid et al., 2006). For more accurate normalization, use of at least two reference genes followed by geometric averaging is recommended for determining a normalization factor (Vandesompele et al., 2002; Reid et al., 2006). Therefore, two reference genes, VvACTIN and VvEF1-α, were used in our study.

Moreover, we carried out a preliminary exploration of the function of the VvYABBY4 gene. The reduced apical dominance of the transgenic tomato plants is difficult to explain; however, it has been shown that overexpression of FW2.2 will dwarf corn plants (Mei et al., 2010), so we speculated that the overexpression of VvYABBY4 led to the increase of FW2.2 expression, which inhibited the growth of the plant, given the presumed function of the YABBY genes. The more tightly packed arrangement of palisade cells in the 35S:VvYABBY4 plant may lead to the apparent deeper green color. The surface of pollen grains showed irregular bulges, but this does not seem to affect normal pollination and fertilization. In addition, as shown in Figure 2A, VvYABBY4 belong to the YAB2 subfamily, which is known to contribute to the development of lateral organs. Interestingly, we found that heterologous expression of VvYABBY4 in tomato negatively impacts seed size, suggesting it may participate in seed development in grapevine. Similarly, the function of the rice YAB2-clade gene DROOPING LEAF is unrelated to lateral organ polarity (Jang et al., 2004; Yamada et al., 2004). Besides that, ectopic expression of VpYABBY2 in Arabidopsis did not lead to any alteration in polarity. Rather it may play a specific role in carpel development and grape berry morphogenesis (Xiang et al., 2013). From these perspectives, we can conclude that during evolution of angiosperms, the functions of YABBY genes became increasingly diversified.

In addition, it was reported that seed size in higher plants can be affected by the development of the maternal tissue (Li et al., 2011) or promoting cell proliferation in the integuments of ovules (Adamski et al., 2009). Our finding that transgenic tomato plants constitutively expressing VvYABBY4 formed smaller seeds compared with wild-type suggests that VvYABBY4 may mainly function to restrict endosperm cell expansion during seed growth.

In conclusion, seven grape YABBY genes were analyzed systematically in terms of evolution and expression pattern during ovule development. Functional characterization suggested that VvYABBY4 may play a role in growth of endosperm cells during seed development, which may provide new insights on genetic pathways involving YABBY genes as potential candidates for seedless grape breeding.

Data Availability Statement

All datasets for this study are included in the manuscript/Supplementary Files.

Author Contributions

XW and SZ designed the study. SZ and XS contributed to most experiments. YL and JY assisted with the analysis of the results. XW and LW provided guidance on the study. SN modified the language of the manuscript. SZ and XW wrote the manuscript. All of the authors approved the final version of this manuscript.


This work was supported by the Joint Funds of National Natural Science Foundation of China (U1603234) and the Program for Innovative Research Team of Grape Germplasm Resources and Breeding (2013KCT-25).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at:


Adamski, N. M., Anastasiou, E., Eriksson, S., O’Neill, C. M., Lenhard, M. (2009). Local maternal control of seed size by KLUH/CYP78A5-dependent growth signaling. Proc. Natl. Acad. Sci. U. S. A. 106 (47), 20115–20120. doi: 10.1073/pnas.0907024106

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartley, G. E., Ishida, B. K. (2003). Developmental gene regulation during tomato fruit ripening and in-vitro sepal morphogenesis. BMC Plant Biol. 3 (1), 4. doi: 10.1186/1471-2229-3-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Bartley, G. E., Ishida, B. K. (2007). Ethylene-sensitive and insensitive regulation of transcription factor expression during in vitro tomato sepal ripening. J. Exp. Bot. 58 (8), 2043–2051. doi: 10.1093/jxb/erm075

PubMed Abstract | CrossRef Full Text | Google Scholar

Bemer, M., Karlova, R., Ballester, A. R., Tikunov, Y. M., Bovy, A. G., Wolters-Arts, M., et al. (2012). The tomato FRUITFULL homologs TDR4/FUL1 and MBP7/FUL2 regulate ethylene-independent aspects of fruit ripening. Plant Cell 24 (11), 4437–4451. doi: 10.1105/tpc.112.103283

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowman, J. L. (2000). The YABBY gene family and abaxial cell fate. Curr. Opin. Plant Biol. 3 (1), 17–22. doi: 10.1016/S1369-5266(99)00035-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowman, J. L., Smyth, D. R., Meyerowitz, E. M. (1989). Genes directing flower development in Arabidopsis. Plant Cell 1 (1), 37–52. doi: 10.1105/tpc.1.1.37

PubMed Abstract | CrossRef Full Text | Google Scholar

Busi, M. V., Bustamante, C., D’Angelo, C., Hidalgocuevas, M., Boggio, S. B., Valle, E. M., et al. (2003). MADS-box genes expressed during tomato seed and fruit development. Plant Mol. Biol. 52 (4), 801–815. doi: 10.1023/a:1025001402838

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakrabarti, M., Zhang, N., Sauvage, C., Muños, S., Blanca, J., Cañizares, J., et al. (2013). A cytochrome P450 regulates a domestication trait in cultivated tomato. Proc. Natl. Acad. Sci. U. S. A. 110 (42), 17125–17130. doi: 10.1073/pnas.1307313110

PubMed Abstract | CrossRef Full Text | Google Scholar

Cong, B., Barrero, L. S., Tanksley, S. D. (2008). Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication. Nat. Genet. 40 (6), 800. doi: 10.1038/ng.144

PubMed Abstract | CrossRef Full Text | Google Scholar

Coombe, B. G. (1992). Research on development and ripening of the grape berry. Am. J. Enol. Vitic. 43 (1), 101–110. doi: 10.1016/S0065-2164(08)70256-9

CrossRef Full Text | Google Scholar

Czerednik, A., Busscher, M., Angenent, G. C., de Maagd, R. A. (2015). The cell size distribution of tomato fruit can be changed by overexpression of CDKA1. Plant Biotechnol. J. 13 (2), 259–268. doi: 10.1111/pbi.12268

PubMed Abstract | CrossRef Full Text | Google Scholar

Eddy, S. R. (1998). Profile Hidden Markov Models;review. Bioinformatics 14 (9), 755–763. doi: 10.1093/bioinformatics/14.9.755

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, L., Su, L., Sun, X., Li, X., Sun, M., Karanja, K. S., et al. (2016). Expression ofVitis amurensis NAC26in Arabidopsis enhances drought tolerance by modulating jasmonic acid synthesis. J. Exp. Bot. 67 (9), 2829–2845. doi: 10.1093/jxb/erw122

PubMed Abstract | CrossRef Full Text | Google Scholar

Finn, R. D., Jody, C., William, A., Miller, B. L., Wheeler, T. J., Fabian, S., et al. (2015). HMMER web server: 2015 update. Nucleic Acids Res. 43 (1), 30–38. doi: 10.1093/nar/gkv397

CrossRef Full Text | Google Scholar

Genova, A. D., Almeida, A. M., Muñoz-Espinoza, C., Vizoso, P., Travisany, D., Moraga, C., et al. (2014). Whole genome comparison between table and wine grapes reveals a comprehensive catalog of structural variants. BMC Plant Biol. 14 (1), 7–7. doi: 10.1186/1471-2229-14-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonçalves, S., Cairney, J., Maroco, J., Oliveira, M. M., Miguel, C. (2005). Evaluation of control transcripts in real-time RT-PCR expression analysis during maritime pine embryogenesis. Planta 222 (3), 556–563. doi: 10.1007/s00425-005-1562-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, C., Guo, R., Xu, X., Gao, M., Li, X., Song, J., et al. (2014). Evolution and expression analysis of the grape (Vitis vinifera L.) WRKY gene family. J. Exp. Bot. 65 (6), 1513. doi: 10.1093/jxb/eru007

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, S., Dai, S., Singh, P. K., Wang, H., Wang, Y., Tan, J. L. H., et al. (2018). A membrane-bound NAC-like transcription factor OsNTL5 represses the glowering in Oryza sativa. Front. Plant Sci. 9, 555–. doi: 10.3389/fpls.2018.00555

CrossRef Full Text | Google Scholar

Holub, E. B. (2001). The arms race is ancient history in Arabidopsis, the wildflower. Nat. Rev. Genet. 2 (7), 516–527. doi: 10.1038/35080508

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Houten, J. V., Gonzalez, G., Han, X., Knaap, E. V. D. (2013). Genome-wide identification, phylogeny and expression analysis of SUN, OFP and YABBY gene family in tomato. Mol. Genet. Genomics 288 (3–4), 111–129. doi: 10.1007/s00438-013-0733-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaillon, O., Aury, J. M., Noel, B., Policriti, A., Clepet, C., Casagrande, A., et al. (2007). The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature

Google Scholar

Jang, S., Hur, J., Kim, S. J., Han, M. J., Kim, S. R., An, G. (2004). Ectopic expression of OsYAB1 causes extra stamens and carpels in rice. Plant Mol. Biol. 56 (1), 133–143. doi: 10.1007/s11103-004-2648-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Juarez, M. T., Twigg, R. W., Timmermans, M. C. (2004). Specification of adaxial cell fate during maize leaf development. Development 131 (18), 4533–4544. doi: 10.1242/dev.01328

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumaran, M. K., Bowman, J. L., Sundaresan, V. (2002). YABBY polarity genes mediate the repression of KNOX homeobox genes in Arabidopsis. Plant Cell 14 (11), 2761–2770. doi: 10.1105/tpc.004911

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B., Fields, S. (1993). Identification of mutations in p53 that affect its binding to SV40 large T antigen by using the yeast two-hybrid system. FASEB J. 7 (10), 957–963. doi: 10.1096/fasebj.7.10.8344494

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Zheng, L., Corke, F., Smith, C., Bevan, M. W. (2011). Control of final seed and organ size by the DA1 gene family in Arabidopsis thaliana. Development 22 (20), 4545–4554. doi: 10.1101/gad.463608

CrossRef Full Text | Google Scholar

Li, Y., Zhu, Y., Yao, J., Zhang, S., Wang, L., Guo, C., et al. (2017a). Genome-wide identification and expression analyses of the homeobox transcription factor family during ovule development in seedless and seeded grapes. Sci. Rep. 7 (1), 12638. doi: 10.1038/s41598-017-16207-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Zhang, C., Guo, Y., Niu, W., Wang, Y., Xu, Y. (2017b). Evolution and expression analysis reveal the potential role of the HD-Zip gene family in regulation of embryo abortion in grapes (Vitis vinifera L.). BMC Genomics 18 (1), 744. doi: 10.1186/s12864-017-4110-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Lozano, R., Angosto, T., Gomez, P., Payan, C., Capel, J., Huijser, P., et al. (1998). Tomato flower abnormalities induced by low temperatures are associated with changes of expression of MADS-Box genes. Plant Physiol. 117 (1), 91. doi: 10.1104/pp.117.1.91

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, T., Zhang, G., Sun, L. (2017). Genome-wide identification of CBL family and expression analysis ofCBLsin response to potassium deficiency in cotton. PeerJ 5 (8), e3653. doi: 10.7717/peerj.3653

CrossRef Full Text | Google Scholar

Lynch, M., Conery, J. S. (2000). The evolutionary fate and consequences of duplicate genes. Science 290 (5494), 1151–1155. doi: 10.1126/science.290.5494.1151

PubMed Abstract | CrossRef Full Text | Google Scholar

Lynch, M., Conery, J. S. (2003). The evolutionary demography of duplicate genes. J. Struct. Funct. Genomics 3 (1–4), 35–44. doi: 10.1023/A:1022696612931

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyons, E., Pedersen, B., Kane, J., Alam, M., Ming, R., Tang, H., et al. (2008). Finding and comparing syntenic regions among Arabidopsis and the outgroups papaya, poplar, and grape: CoGe with Rosids. Plant Physiol. 148 (4), 1772–1781. doi: 10.1104/pp.108.124867

PubMed Abstract | CrossRef Full Text | Google Scholar

Mccormick, S., Niedermeyer, J., Fry, J., Barnason, A., Horsch, R., Fraley, R. (1986). Leaf disc transformation of cultivated tomato (L. esculentum) using Agrobacterium tumefaciens. Plant Cell Rep. 5 (2), 81–84. doi: 10.1007/BF00269239

PubMed Abstract | CrossRef Full Text | Google Scholar

Mei, G., Rupe, M. A., Jo Ann, D., Jijun, Z., Daniel, S., Duncan, K. E., et al. (2010). Cell number regulator1 affects plant and organ size in maize: implications for crop yield enhancement and heterosis. Plant Cell 22 (4), 1057–1073. doi: 10.1105/tpc.109.073676

PubMed Abstract | CrossRef Full Text | Google Scholar

Mei, G., Simmons, C. R. (2011). Cell number counts—the fw2.2 and CNR genes and implications for controlling plant fruit and organ size. Plant Sci. 181 (1), 1–7. doi: 10.1016/j.plantsci.2011.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Nagasawa, N., Miyoshi, M., Sano, Y., Satoh, H., Hirano, H., Sakai, H., et al. (2003). SUPERWOMAN1 and DROOPING LEAF genes control floral organ identity in rice. Development 130 (4), 705–718. doi: 10.1242/dev.00294

PubMed Abstract | CrossRef Full Text | Google Scholar

Nakano, T., Suzuki, K., Fujimura, T., Shinshi, H. (2006). Genome-wide analysis of the ERF gene family in Arabidopsis and rice. Plant Physiol. 140 (2), 411. doi: 10.1104/pp.105.073783

PubMed Abstract | CrossRef Full Text | Google Scholar

Nathalie, G., Frédéric, G., Michel, H., Christian, C., Armand, M. (2007). The cell cycle-associated protein kinase WEE1 regulates cell size in relation to endoreduplication in developing tomato fruit. Plant J. 51 (4), 642–655. doi: 10.1111/j.1365-313X.2007.03167.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Nathalie, G., Michel, H., Frédéric, D., Frédéric, G., Philippe, D., Mathilde, C., et al. (2004). Molecular characterization of a WEE1 gene homologue in tomato (Lycopersicon esculentum Mill.). Plant Mol. Biol. 56 (6), 849–861. doi: 10.1007/s11103-004-5110-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohmori, Y., Abiko, M., Horibata, A., Hirano, H. Y. (2008). A transposon, Ping, is integrated into intron 4 of the DROOPING LEAF gene of rice, weakly reducing its expression and causing a mild drooping leaf phenotype. Plant Cell Physiol. 49 (8), 1176–1184. doi: 10.1093/pcp/pcn093

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohno, S., Wolf, U., Atkin, N. B. (2010). Evolution from fish to mammals by gene duplication. Hereditas 59 (1), 169–187. doi: 10.1111/j.1601-5223.1968.tb02169.x

CrossRef Full Text | Google Scholar

Pham, V. T., Herrero, M., Hormaza, J. I. (2015). Effect of temperature on pollen germination and pollen tube growth in longan (Dimocarpus longan Lour.). Sci. Hortic. 197, 470–475. doi: 10.1016/j.scienta.2015.10.007

CrossRef Full Text | Google Scholar

Punta, M., Coggill, P. C., Eberhardt, R. Y., Mistry, J., Tate, J., Boursnell, C., et al. (2012). The Pfam protein families database. Nucleic Acids Res. 40 (database issue), D290. doi: 10.1093/nar/gkr1065

PubMed Abstract | CrossRef Full Text | Google Scholar

Reid, K. E., Olsson, N., Schlosser, J., Peng, F., Lund, S. T. (2006). An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 6 (1), 27–27. doi: 10.1186/1471-2229-6-27

PubMed Abstract | CrossRef Full Text | Google Scholar

Sawa, S., Watanabe, K., Goto, K., Kanaya, E., Morita, E. H., Okada, K. (1999). FILAMENTOUS FLOWER, a meristem and organ identity gene of Arabidopsis, encodes a protein with a zinc finger and HMG-related domains. Genes Dev. 13 (9), 1079–1088. doi: 10.1101/gad.13.9.1079

PubMed Abstract | CrossRef Full Text | Google Scholar

Siegfried, K. R., Eshed, Y., Baum, S. F., Otsuga, D., Drews, G. N., Bowman, J. L. (1999). Members of the YABBY gene family specify abaxial cell fate in Arabidopsis. Development 126 (18), 4117–4128. doi: 10.1021/ie020097t

PubMed Abstract | CrossRef Full Text | Google Scholar

Stracke, R., Werber, M., Weisshaar, B. (2001). The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 4 (5), 447–456. doi: 10.1016/S1369-5266(00)00199-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Q., Wang, G., Zhang, X., Zhang, X., Qiao, P., Long, L., et al. (2017). Genome-wide identification of the TIFY gene family in three cultivated Gossypium species and the expression of JAZ genes. Sci. Rep. 7, 42418. doi: 10.1038/srep42418

PubMed Abstract | CrossRef Full Text | Google Scholar

Sutherland, P., Veit, B., Yao, J. L. (2002). Down-regulation of TM29, a tomato SEPALLATA homolog, causes parthenocarpic fruit development and floral reversion. Plant Physiol. 130 (2), 605–617. doi: 10.1104/pp.005223

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, W., Hirano, H. Y. (2012). The YABBY gene TONGARI-BOUSHI1 is involved in lateral organ development and maintenance of meristem organization in the rice spikelet. Plant Cell 24 (1), 80–95. doi: 10.1105/tpc.111.094797

PubMed Abstract | CrossRef Full Text | Google Scholar

Tanaka, W., Toriba, T., Hirano, H. Y. (2017). Three TOB1-related YABBY genes are required to maintain proper function of the spikelet and branch meristems in rice. New Phytol. 215 (2), 825–839. doi: 10.1111/nph.14617

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, J. D., Gibson, T. J., Higgins, D. G.. (2002). Multiple Sequence Alignment Using ClustalW and ClustalX. Curr Protoc Bioinformatics Chapter 2 (Unit 2):Unit 2.3. doi: 10.1002/0471250953.bi0203s00

CrossRef Full Text | Google Scholar

Toriba, T., Harada, K., Takamura, A., Nakamura, H., Ichikawa, H., Suzaki, T., et al. (2007). Molecular characterization the YABBY gene family in Oryza sativa and expression analysis of OsYABBY1. Mol. Genet. Genomics 277 (5), 457–468. doi: 10.1007/s00438-006-0202-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3. doi: 10.1186/gb-2002-3-7-research0034

CrossRef Full Text | Google Scholar

Villanueva, J. M., Broadhvest, J., Hauser, B. A., Meister, R. J., Schneitz, K., Gasser, C. S. (1999). INNER NO OUTER regulates abaxial–adaxial patterning in Arabidopsis ovules. Genes Dev. 13 (23), 3160. doi: 10.1101/gad.13.23.3160

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, T., Liu, Z. M., Li, L. F., Leitch, A. R., Leitch, I. J., Lohaus, R., et al. (2018). A genome for gnetophytes and early evolution of seed plants. Nat. Plants 4 (2), 82. doi: 10.1038/s41477-017-0097-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Hu, X., Jiao, C., Li, Z., Fei, Z., Yan, X., et al. (2016). Transcriptome analyses of seed development in grape hybrids reveals a possible mechanism influencing seed size. BMC Genomics 17 (1), 898. doi: 10.1186/s12864-016-3193-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Yin, X., Cheng, C., Wang, H., Guo, R., Xu, X., et al. (2015). Evolutionary and expression analysis of a MADS-box gene superfamily involved in ovule development of seeded and seedless grapevines. Mol. Genet. Genomics 290 (3), 825–846. doi: 10.1007/s00438-014-0961-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiang, J., Liu, R. Q., Li, T. M., Han, L. J., Zou, Y., Xu, T. F., et al. (2013). Isolation and characterization of two VpYABBY genes from wild Chinese Vitis pseudoreticulata. Protoplasma 250 (6), 1315–1325. doi: 10.1007/s00709-013-0514-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamada, T., Ito, M., Kato, M. (2004). YABBY2-homologue expression in lateral organs of Amborella trichopoda(Amborellaceae). Int. J. Plant Sci. 165 (6), 917–924. doi: 10.1086/423793

CrossRef Full Text | Google Scholar

Zhang, G., Lu, T., Miao, W., Sun, L., Tian, M., Wang, J., et al. (2017). Genome-wide identification of ABA receptor PYL family and expression analysis ofPYLsin response to ABA and osmotic stress inGossypium. PeerJ 5 (Web Server issue), e4126. doi: 10.7717/peerj.4126

CrossRef Full Text | Google Scholar

Zhang, S. W., Huang, G. X., Ding, F., He, X. H., Pan, J. C. (2012a). Mechanism of seedlessness in a new lemon cultivar ‘Xiangshui’ [Citrus limon (L.) Burm. F]. Sex. Plant Reprod. 25 (4), 337–345. doi: 10.1007/s00497-012-0201-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Gao, M., Singer, S. D., Fei, Z., Wang, H., Wang, X. (2012b). Correction: genome-wide identification and analysis of the TIFY gene family in grape. PLoS One 7 (9), e44465. doi: 10.1371/journal.pone.0044465

CrossRef Full Text | Google Scholar

Zhao, S. P., Lu, D., Yu, T. F., Ji, Y. J., Zheng, W. J., Zhang, S. X., et al. (2017). Genome-wide analysis of the YABBY family in soybean and functional identification of GmYABBY10 involvement in high salt and drought stresses. Plant Physiol. Biochem. 119, 132. doi: 10.1016/j.plaphy.2017.08.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, L. (1997). Yeast GAL4 Two-Hybrid System. Recombinant Protein Protocols. Humana Press. 63, 173–196. doi: 10.1385/0-89603-481-X:173

CrossRef Full Text | Google Scholar

Keywords: grapevine, seedlessness, transcription factor, YABBY, gene function

Citation: Zhang S, Wang L, Sun X, Li Y, Yao J, Nocker Sv and Wang X (2019) Genome-Wide Analysis of the YABBY Gene Family in Grapevine and Functional Characterization of VvYABBY4. Front. Plant Sci. 10:1207. doi: 10.3389/fpls.2019.01207

Received: 28 March 2019; Accepted: 03 September 2019;
Published: 08 October 2019.

Edited by:

Takaki Yamauchi, PRESTO, Japan Science and Technology Agency, Japan

Reviewed by:

Patricio Hinrichsen, INIA Chile, Chile
Liangsheng Zhang, Fujian Agriculture and Forestry University, China

Copyright © 2019 Zhang, Wang, Sun, Li, Yao, Nocker and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiping Wang,

These authors have contributed equally to this work