Abstract
In the recent past, genetic analyses of grapevine focused mainly on the identification of resistance loci for major diseases such as powdery and downy mildew. Currently, breeding programs make intensive use of these results by applying molecular markers linked to the resistance traits. However, modern genetics also allows to address additional agronomic traits that have considerable impact on the selection of grapevine cultivars. In this study, we have used linkage mapping for the identification and characterization of flowering time and ripening traits in a mapping population from a cross of V3125 (‘Schiava Grossa’ × ‘Riesling’) and the interspecific rootstock cultivar ‘Börner’ (Vitis riparia × Vitis cinerea). Comparison of the flowering time QTL mapping with data derived from a second independent segregating population identified several common QTLs. Especially a large region on linkage group 1 proved to be of special interest given the genetic divergence of the parents of the two populations. The proximity of the QTL region contains two CONSTANS-like genes. In accordance with data from other plants such as Arabidopsis thaliana and Oryza sativa, we hypothesize that these genes are major contributors to control the time of flowering in Vitis.
Similar content being viewed by others
Introduction
Grapevine is a crop used worldwide for fermentation into wine, for fresh and dried fruit supply as well as for juice production. Most grapevine cultivars suffer seriously from various pests and diseases. Thus, a major aim in grapevine breeding is the introgression of traits from Vitis species that carry resistances found in species native to Eastern Asia and North America. However, traits such as time of flowering and veraison (beginning of berry ripening) have recently also gained increasing interest due to shifts of plant phenology observed in the context of global warming. The development of molecular markers linked to loci or genes affecting these traits has thus become an important research topic. Such markers will facilitate the early selection of plants on the genetic level, simplifying the introgression of desirable and excluding undesired traits for optimal adaptation of plant development to current climatic conditions.
Several agronomic traits important for grapevine cultivation have been genetically dissected during the past few years. The traits addressed include berry color, berry size and seedlessness (Mejía et al. 2011; Cabezas et al. 2006; Fischer et al. 2004; Doligez et al. 2002), phenology-related traits such as veraison and flowering time (Costantini et al. 2008; Mejía et al. 2007), inflorescence number per shoot (Marguerit et al. 2009) or flower sex (Fechter et al. 2012; Marguerit et al. 2009; Lowe and Walker 2006; Riaz et al. 2004; Dalbó et al. 2000). However, the results usually vary greatly between the genotypes analyzed. Consequently, the identification of loci and genes determining these traits at species level has so far proven difficult. Therefore, the development of robust trait-linked genetic markers remains challenging.
Complex traits that are controlled by many loci such as regulation of flowering time (Costantini et al. 2008) are particularly difficult to analyze, although the published Vitis reference genome sequences (Adam-Blondon et al. 2011; Jaillon et al. 2007; Velasco et al. 2007) allow rather easy access to grapevine genes. The genome annotation also enables to compare grapevine genes with genetic data and gene functions assigned in other plants.
Few findings have so far been reported for flowering time control in grapevine. By studying one of the rare mutants known in grapevine, Boss and Thomas (2002) suggested that gibberellic acid inhibits the development of inflorescences. As the authors point out, this is in contradiction with results obtained from model plants such as Arabidopsis thaliana where gibberellic acid promotes flowering (Song et al. 2013; Blazquez et al. 1998, Blazquez and Weigel 2000; Langridge 1957). This finding may indicate that greater differences in phenological control mechanisms between annual plants and woody perennial plants like grapevine are to be expected. However, since it seems likely that similar genes are involved, research in Vitis has focused on the identification and characterization of candidate genes and gene families for flowering time homologous to those found in A. thaliana and other well-studied plants. Using this strategy, the VvFT/TFL1 gene family was already examined in detail and the VvFT gene was shown to promote flowering in grapevine as expected from experiments run with model plants (Joly et al. 2004; Boss et al. 2006; Sreekantan and Thomas 2006; Carmona et al. 2007). In addition, the MADS-box gene VvMADS8 has been characterized in detail (Sreekantan and Thomas 2006). Three members of this subfamily of MADS transcription factor genes are present in the grapevine genome (Carmona et al. 2008). Other work focused on NAP (no-apical meristem) homologs from grapevine, being targets of the floral homeotic genes PISTILLATA/APETALA3 in A. thaliana (Fernandez et al. 2006). However, compared to the detailed knowledge of the pathways controlling flowering in A. thaliana (e.g. Ehrenreich et al. 2009; Andres and Coupland 2012), information about the genes involved in grapevine flowering is still very limited.
The construction of detailed genetic maps in combination with intensive phenotyping and subsequent QTL analysis provides an excellent way to identify genomic regions and define loci involved in expression of the traits investigated. We have chosen this strategy, extended by candidate gene prediction in the target regions, to localize loci controlling flowering and veraison time. The basis for mapping is a cross between the high-quality Vitis vinifera breeding line V3125 (‘Schiava Grossa’ × ‘Riesling’) and the interspecific cultivar ‘Börner’ (Vitis riparia × Vitis cinerea). We searched for flowering time QTLs and closely linked markers using a significantly extended marker set. By analyzing the genotypic and phenotypic segregation in the progeny of V3125 × ‘Börner’, we detected highly significant QTLs on linkage groups 1 and 14. The same QTLs became evident in a second independent mapping population, suggesting an overall role of these regions in the control of flowering time. Genes annotated in the delimited QTL regions include several CONSTANS-like genes as well as other genes that are candidates for a contribution to the control of flowering time.
Materials and methods
Plant material
Mapping population 1: V3125 × ‘Börner’
The mapping population of 202 individuals was obtained by crossing the cultivars V3125 (‘Schiava Grossa’ × ‘Riesling’) and ‘Börner’ (V. riparia Gm183 × V. cinerea Arnold) in the years 1998 and 2001 and was planted in the vineyard at Geilweilerhof. V3125 is an elite F1 individual derived from two traditional Vitis vinifera cultivars and represents the highly pathogen-susceptible, high-quality V. vinifera type. ‘Börner’ is a hybrid of the two American Vitis species V. riparia and V. cinerea and was bred for use as a rootstock cultivar. This cultivar was shown to be highly resistant to different pathogens, including grape phylloxera, downy and powdery mildew and black rot. The F1 progeny of the cross V3125 × ‘Börner’ segregates for several agronomic, morphological and resistance-related traits. Phenotypic and genotypic data of the mapping population were used for the calculation of genetic maps and QTL analysis.
Mapping population 2: GF.GA-47-42 × ‘Villard blanc’
This second mapping population was used to confirm selected QTL loci for flowering time found in the V3125 × ‘Börner’ population. It is derived from the cross of the breeding line GF.GA-47-42 (‘Bacchus weiss’ x ‘Seyval’) × ‘Villard blanc’ (Seibel 6468 × ‘Subereux’), comprises 151 F1 individuals and segregates considerably for flowering time. The cross was performed in 1989 and the F1 progeny was planted in the vineyards at Geilweilerhof in 1996. The maternal breeding line GF.GA-47-42 is early flowering. In contrast, ‘Villard blanc’ is an offspring of two late-flowering hybrids and flowers very late. A first genetic map for this mapping population was published in 2006 (Zyprian et al. 2006). This map was improved with additional markers (E. Zyprian, in preparation) and used for QTL analysis.
Phenotyping of flowering time and veraison
The progeny of the mapping population V3125 × ‘Börner’ was evaluated in 2010–2012 for flowering time and veraison time (2010–2011) (Fig. 1). Flowering time was scored according to the BBCH scale (Lorenz et al. 1995) as follows: beginning of flowering (stage 61: 10 % of flowerhoods fallen) and full bloom (stage 65: 50 % of flowerhoods fallen). The date at which each individual reached stages 61 and 65 was recorded. The date of veraison time (begin of ripening) was recorded according to the BBCH scale stage 81 (beginning of ripening: berries begin to develop variety-specific color).
The GF.GA-47-42 × ‘Villard blanc’ population was phenotyped for flowering time (full bloom only) in 5 years (1999, 2009–2012). Since these are both white varieties, veraison time was determined as beginning of fruit softening (10–15 % of fruits soft).
For subsequent QTL analysis, the phenotypic data were classified according to five stages for flowering time following OIV descriptor 302 (OIV 2009) (1 = early flowering; 2 = medium early flowering; 3 = medium flowering; 4 = medium late flowering; 5 = late flowering). Assigning the data into these classes proved to be advantageous in comparison to the use of direct dates for QTL analysis as unavoidable differences in flowering dates and periods in the respective years were masked. Using classes thus led to more defined QTL regions with higher LOD scores and smaller confidence intervals. The whole period between the first and last flowering individual was counted in days—starting with the first flowering individual—and divided equally into five periods. Depending on the respective year, each class thus represented 2–4 days of the flowering period. Individuals were then assigned to the according class.
For veraison time, the phenotypic data was classified using a similar system as described above basically following OIV descriptor 303 (OIV 2009). However, it showed that only the application of a 6-class system led to a normal distribution of the data, which—in contrast to the 5-class OIV system—was thus classified to six stages with 1 = very early veraison, 2 = early veraison, 3 = medium veraison, 4 = medium late veraison, 5 = late veraison and 6 = very late veraison.
DNA extraction
Young leaves were collected and immediately frozen at −70 °C. After lyophilisation of the material, DNA extraction was carried out using the method of Lemke et al. (2011). The DNA was quantified with a spectral NanoPhotometer (Implen GmbH, München, Germany). Samples were adjusted to a DNA concentration of 1 ng/µl for subsequent PCR analyses.
SSR marker development and analysis
Simple sequence repeat (SSR) markers were designed based on anchored scaffolds from the 8x and 12x version of the published grapevine genome sequence PN40024 internet release (Jaillon et al. 2007; http://www.genoscope.cns.fr/cgi-bin/ggb/vitis/12X/gbrowse/vitis/; Adam-Blondon et al. 2011), respectively. Special focus was laid on regions which had not been covered sufficiently with markers in the previous map elaborated by Zhang et al. (2009). Particular emphasis was laid on the terminal regions of the chromosomes as well as on newly assembled genome regions in the 12x version for which no or only few markers were available before. Supplementary Table S1 gives an overview of all newly designed markers, their source of reference genome (8x or 12x version) and the amplimers. Sequences in regions of interest were scanned for SSRs using the program WebSat (http://wsmartins.net/websat/). Primers flanking the SSRs were developed using the Primer3 software included in the WebSat program. For multiplexing purposes, the length of the amplimers was varied from 80 to 400 bp. Annealing temperatures were set to 60 °C. Forward primers were labeled at the 5′ end with either one of four different fluorescent dyes (HEX, FAM, TAMRA, ROX).
SSR markers were screened for informative segregation using a small subset of genotypes including the parents, grandparents and ten F1 progeny. Additional markers were chosen from literature and also analyzed as described above: VChr- (Cipriani et al. 2008), VMC- (Vitis Microsatellite consortium; http://www.agrogene.com), VVMD- (Bowers and Meredith 1996; Bowers et al. 1999), VrZAG- (Sefc et al. 1999) and UDV-(Di Gaspero et al. 2005). All informative markers were combined in multiplexes with 6–12 markers, taking into account the different fragment sizes and fluorescence labels. PCR amplifications were carried out using the Qiagen Multiplex kit (Qiagen GmbH, Hilden, Germany) following the supplier’s instructions. The lengths of the resulting amplicons were analyzed on an ABI 3110xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Marker data and allele calls were collected in flat file format for subsequent analyses.
Map construction and QTL analysis
For construction of genetic maps, the double pseudo-testcross strategy (Grattapaglia and Sederoff 1994) was applied, resulting in two separate maps for the parents and a consensus map, respectively. The genotypic information was evaluated and the maps were calculated using the program JoinMap 4.0 (Van Ooijen 2006). Markers were tested for their goodness of fit segregation ratio with a Chi-square test. Samples with very high Chi-square values >50 or more than 20 % of missing marker data were excluded from the calculation. Map distances were estimated by applying the Kosambi function (Kosambi 1944). LOD score thresholds of 5 or higher were used for linkage group determination. To avoid map inflation, the reliability of the terminal marker positions was checked with special care by evaluating their data and deduced locations. Suspicious markers that were placed on the map by JoinMap at the outermost positions of the chromosomes with an unreasonable long distance (e.g. 30 or more cM) to all other markers were excluded from the map. In addition, a few markers that were placed at wrong positions with regard to marker phases were also removed from the final dataset, leaving 374 integrated markers.
The data points that had been collected and assigned to the respective phenotypic classes in each season were used for QTL detection. In addition, the median values of the phenotypic classes were also calculated over all years and used as data sets. In our analyses, medians were used instead of means as for some individuals the scored date of flowering time differed greatly in one of the years compared to the other years. Statistically speaking, the median is more robust against these outliers than the mean and therefore represents data that is less skewed. QTL analyses were carried out using the program MapQTL 5.0 (Van Ooijen 2004). Putative QTLs were identified via interval mapping (IM) and subsequent multiple QTL mapping (MQM), using QTL-flanking markers as co-factors. For non-normally distributed traits the non-parametric Kruskal–Wallis rank sum test was applied. Both linkage group specific and genome-wide significance thresholds were determined by permutation tests (1,000 permutations) and the thresholds were set at a value of p < 0.05. QTLs were considered reliable within the population if they were identified in at least two growing seasons. Markers associated with peak LOD score values were additionally evaluated for their significance with the Kruskal–Wallis rank sum test.
Results
Marker development based on the reference sequence
In total, 215 SSR markers were newly designed based on the reference genome sequence PN40024 (Jaillon et al. 2007; Adam-Blondon et al. 2011) and tested for their segregation type in V3125 × ‘Börner’ (Table 1; supplementary Table S1). Roughly 2/3 (=67 %) of all tested markers could be assigned to one of the segregation types shown in Table 1. Overall, fifteen percent of the markers proved to be monomorphic, while 18 % gave no or unclear signals and were thus excluded from further analysis. 47 % of the usable markers were fully informative segregating with four or three alleles from both parents (ab × cd, ef × eg) and 53 % segregated from only one heterozygous parent (maternal: lm × ll, paternal: nn × np).
Map construction and extension: the genetic maps
One hundred forty-one newly designed markers and 50 additional markers transferred from other public sources (see “Materials and methods”) could be integrated into the previously developed consensus map of V3125 × ‘Börner’ (Zhang et al. 2009). In total, 374 markers could be assigned to 19 linkage groups. These were numbered according to the International Grape Genome Program (IGGP) rules of nomenclature (Adam-Blondon et al. 2004; Riaz et al. 2004). The extended consensus map covered 1365 cM with an average distance of 3.9 cM between the markers (Table 2; Fig. 2).
A high degree of co-linearity could be observed between the marker order given in the PN40024 reference genome sequence and the new map presented here. Almost all of the newly designed markers could be assigned to the targeted chromosome region. Only three markers did not localize in the region they were designed from, namely marker GF10-06 (expected on LG10, localized on LG 06), GF02-02 (expected on LG02, localized on LG 10) and GF02-24 (expected on LG02, mapped on LGs 14 and 16).
Using marker development based on the newly assembled and other chromosome parts in the 12x version of the reference genome sequence from PN40024, the map could be extended by 209 cM. Extension of map coverage was achieved particularly at chromosome ends. It should be noted that the percentage of newly identified informative SSR markers varied strongly for the different chromosomes. It proved especially difficult to develop markers for chromosome 16, on which only 20 % of all tested markers showed informative segregation. In addition, only one of them, namely GF16-01, was fully informative while the other markers on this linkage group were segregating only from one (heterozygous) parent yielding limited information.
Phenotypic data
Flowering time variation in the mapping population V3125 × ‘Börner’ was recorded twice per season, firstly at the onset of opening of flowers and a second time at full bloom. Begin of flowering and full bloom started at different time points between May and June in the years 2010–2012, with an interval of 2–4 days between begin of flowering and full bloom (data not shown). The time interval until all genotypes had reached the full bloom stage also varied between 13 days in the years 2010/2011 and 19 days in 2012 (Fig. 3). The respective data correlated with r = 0.4–0.5 (Pearson’s correlation coefficient) between all the years. The frequency distribution of full bloom time over all years is presented in Fig. 4a.
Veraison time (begin of berry ripening) was recorded 39–79 days after full bloom for the mapping population V3125 × ‘Börner’. The data correlated with r = 0.3 between the years 2010–2011 and 2011–2012 and r = 0.5 between years 2010 and 2012. The frequency distribution over all years is presented in Fig. 4b.
The progeny of the mapping population GF.GA-47-42 × ‘Villard Blanc’ was phenotyped for full bloom only. In this population, flowering started 3–8 days later than in the first mapping population V3125 × ‘Börner’ in the respective years (Fig. 3).
QTL mapping for flowering time and ripening traits
Using the extended genetic map and the phenotypic data captured from the 202 genotypes in the V3125 × ‘Börner’ population, a detailed QTL analysis was performed. Several QTLs were detected for the agronomic traits considered. Flanking markers of primarily identified major QTLs were used as co-factors for subsequent MQM mapping. Few major QTLs that explain more than 20 % of the phenotypic variance were identified. In addition, a number of weaker QTLs were detected that varied in their position and relevance throughout the years. These “smaller” QTLs were only considered reliable if they were detected in at least two or more different seasons. The phenotypic data did not always exhibit a normal distribution. Additional Kruskal–Wallis tests were thus carried out for evaluation of the QTL results. An overview of the results of all QTL analyses based on the median values of the data sets is given in Table 3. The segregation pattern of the respective trait is indicated in regard to its parental origin.
For the start of flowering, genomic regions showing significant effects were detected on chromosomes 1, 10, 11, 14 and 17 (see Table 3; marked fl-s). The same positions were determined for factors affecting the time of full bloom (fl-bl). Additional reliable QTLs for time of full bloom were detected on linkage groups 16 and 19. However, the QTL on linkage group 1 showed the most consistent effects over all years and explained 12.7–12.8 % of phenotypic variance. The 1-LOD interval for the time of full bloom spanned 7.2 cM around the markers GF01-19 and GF01-20. A second LOD peak with a confidence interval of 3.3 and 4.8 cM is found for time of full bloom and the beginning of flowering, respectively, with the same corresponding marker GF01-22. The large QTL region on chromosome 1 covers an area of 20.8 cM for full bloom. Referring to the reference genome sequence PN40024, the physical location for the QTL is between 2.8 and 8.2 Mb. An additional two-peak QTL was detected for the beginning of flowering on LG 14 with the associated markers VMC2a5 and UDV_123. Searching genetic factors affecting full bloom, two adjacent loci on the same chromosome were detected. They stretch over a 1-LOD confidence interval of 9.8 respective 18.5 cM with the corresponding markers VMC2H12 and UDV_057a.
Several regions relevant to determine veraison time were detected on chromosomes 1, 11 and 13 (Table 3; marked ver). The QTLs on chromosomes 1 and 11 covered the same genomic regions as those for full bloom. The QTL on chromosome 1 explained 17.7 % resp. 19.7 % of the phenotypic variance. While the QTL region exceeding the significance threshold was extending over 31.3 cM and showed three distinct peaks in 2010, the 1-LOD interval was drastically reduced to 3.4 cM in the following year with marker GF01-19 lying closest to this QTL. Altogether, the same QTL was identified for veraison time and for flowering time. For the flowering–veraison interval, no significant effects could be determined. All results concerning flowering time and ripening-related traits in the mapping population of V3125 × ‘Börner’ are summarized in Fig. 5.
For the second mapping population GF.GA-47-42 × ‘Villard Blanc’, QTLs for flowering time (full bloom) were detected on the same chromosomes 1, 14, 17 and 19. In addition, further QTLs were located on chromosomes 4 and 8 (Table 4).
For chromosome 1, two QTL peaks could be determined with the corresponding markers GF01-16 and VrZAG29. They explain 17.0 respective 14.6 % of the phenotypic variance. The physical positions of the flanking markers of their 1-LOD confidence intervals indicate that these regions are located between 3.46 and 5.98 Mb with the associated marker VrZAG29 at a position of 5.28 Mb. The two regions overlap to a great extend with the large QTL region on the same chromosome found in the mapping population V3125 × ‘Börner’ (2.8–8.2 Mb).
On LG 14, the QTL around marker VMC2A5 explains 21.1 % of the phenotypic variance for the time of full bloom. The same marker was found to be closest to the QTL peak for the start of flowering in the mapping population V3125 × ‘Börner’ but remained undetected for the time of full bloom. An overview of all QTLs identified in the mapping population GF.GA-47-42 × ‘Villard blanc’ is given in Table 4.
Discussion
In the recent past, genetic analyses have led to the identification of numerous resistance loci for the most important grapevine diseases powdery and downy mildew (for review see Töpfer et al. 2011). These loci have been characterized, and some of them are used via marker-assisted selection (MAS) in breeding programs. Despite the fact that there is further work required to optimize markers for durable resistance, breeders also demand for molecular markers addressing other agronomic and resistance-related traits to support breeding. Up to now, there are only few suitable markers available for genetically well-described quantitative or polygenic agronomic traits. QTL mapping is one of the most successful approaches for trait-linked marker development. We have thus evaluated mapping populations for phenological traits (flowering time, veraison) using this method. The results were compared with the outcomes of other research groups trying to identify QTLs which can reliably be characterized in different genetic backgrounds.
Map construction
The published genome sequence of the nearly homozygous inbred line PN40024 (http://www.genoscope.cns.fr/externe/GenomeBrowser/Vitis/) was used for SSR marker development to increase saturation in the genetic map of the interspecific cultivar ‘Börner’, a hybrid of V. riparia × V. cinerea.
Testing the newly deduced markers showed that an unexpectedly high amount of 82 % was amplification-positive in pretests and two-thirds were informative in segregation analysis. Moreover, most of the newly designed markers could be assigned to the targeted chromosome regions, which confirms the high synteny within the genus Vitis and the usefulness and reliability of the reference genome sequence for marker design. Genomes of different Vitis species have been shown to have a high degree of variability with an average of 3.4 million SNPs (1 every 140 bp) per genotype species (Le Paslier et al. 2013). We therefore expected a much higher rate of amplification failure than observed, especially as this work was carried out with a progeny of an interspecific cross of the type V. vinifera × (V. riparia × V. cinerea). Some of the linkage groups, however, proved difficult to be covered with new markers. Especially linkage group 16, on which only 20 % of all tested markers showed segregation, could not be covered as satisfactorily as the other linkage groups. Segregation was only detectable from the male parent, giving rise to only partially informative markers. As the mapping population genetically represents three different Vitis species, these results could be due to altered levels of heterozygosity in the corresponding chromosome or chromosomal segments. Chromosome 16 seems to be highly homozygous in V. vinifera and thus only shows a few detectable length polymorphisms for the two American species of the ‘Börner’ genome. In addition, chromosome 16 seems to contain large stretches of repetitive sequences, making it difficult to design unique markers. However, the overall results and statistics prove that the strategy of designing SSR marker assays (primer sequences) from an already existing reference sequence is a powerful and easy-to-use tool for future marker design.
Flowering and ripening-related parameters
The most consistent QTL in this study was a region on LG 1 affecting both flowering and veraison time, with two distinct peaks which together cover a broad confidence interval of 20.8 cM respective 5.4 Mb of the chromosome in the first mapping population V3125 × ‘Börner’ (Fig. 6). QTL analyses in the second segregating population of GF.GA-47-42 × ‘Villard blanc’ indicated that especially this region on chromosome 1 might be of interest regarding different phenology traits including the time of full bloom. For this developmental stage, the same QTL regions on chromosome 1 were found independently in both mapping populations, using different sets of markers and phenotypic data.
QTLs for flowering and veraison time and the flowering–veraison interval were also previously reported for chromosomes 1 and 2 by Costantini et al. (2008). Their closest marker linked to flowering time is located approximately 6 Mb away from the QTL detected in this study. In the mapping population V3125 × ‘Börner’, we could not identify reproducible significant QTLs for the flowering–veraison interval. The reason for this might be due to the high consistency between flowering and veraison time QTLs in this mapping population as early-flowering plants also reached veraison time earlier than the late-flowering plants. Just one additional chromosome, namely chromosome 13, was found to be only related with veraison time in comparison to the time of flowering, meaning that there are almost no differences in the intervals between the two traits which could be mapped as a QTL.
Additional QTLs for flowering time were observed on chromosomes 10, 11, 14, 16, 17 and 19 in the mapping population V3125 × ‘Börner’ (Fig. 5), underlining the hypothesis of a complex trait with numerous genes involved in flower formation and the onset of flowering. The effects on chromosome 14 and 17 were also evident in the second mapping population. On both of these LGs, loci have already been proposed for the flowering time trait (Duchene et al. 2012; Mejía et al. 2007). The QTL on LG14 identified by Duchene et al. (2012) is located in close proximity to the one we found in this study. The authors propose the gene VvCOL2 as a candidate gene for flowering time. All results indicate that a combination of several genomic regions influences flowering parameters, with each of them contributing only a more or less small effect to the overall performance of the plants.
The difference in flowering time (early vs. late) is inherited from different genetic backgrounds in the two mapping populations, leading to the shared QTLs shown in Tables 3 and 4. For the progeny of V3125 × ‘Börner’, genetic factors inherited from V. cinerea contributed most to the late-flowering phenotypes on LG 1, 10, 11 and 19 while the alleles from V. riparia led to early flowering. For the QTLs on LG 14 and 17, effects were observed for ‘Schiava Grossa’ (early flowering) and ‘Riesling’ (late flowering) (Table 3). In the pedigree of the second mapping population GF.GA-47-42 × ‘Villard blanc’, however, early flowering was exclusively transmitted by the ancestors of GF.GA-47-42, suggesting that early flowering is a dominant effect. Interspecific French hybrids are the origin of the late-flowering ‘Villard blanc’ (Seibel 6468 and ‘Subereux’). These include V. rupestris in their pedigrees. Although genetically unrelated, the alleles inherited from the parents of the two mapping populations led to the same effect of early respective late flowering (cf. “Materials and methods” section “Plant material”).
In regard to the control of flowering time, complex mechanisms have been revealed in model plants such as A. thaliana (e.g. Ehrenreich et al. 2009), giving a hint at genes that might also contribute to flowering and ripening time control in Vitis species. With the availability of a reference genome sequence for grapevine, it has become possible to postulate candidate genes homologous to A. thaliana or other plant flowering time control genes in the grapevine genome. These include flowering integrator genes such as members of the FT/TFL1 family (Joly et al. 2004; Boss et al. 2006; Carmona et al. 2007; Sreekantan and Thomas 2006) and flower meristem identity genes, e.g. transcription factor LEAFY (VFY in grapevine) (Carmona et al. 2002; Joly et al. 2004; Boss et al. 2006) or the MADS-box genes APELATA1 (VAP1) and FRUITFULL (VFUL-L) (Calonje et al. 2004). Based on the reference sequence annotation, several genes homologous to known plant genes—especially from A. thaliana and O. sativa—are located within the 1-LOD confidence interval of the QTL on chromosome 1 that might contribute to flowering time (Tables 5, 6; Fig. 5). Directly at the QTL peak for the start of flowering around marker VChr01a, the phytochrome and flowering time regulatory protein 1 (PFT1) is annotated in the reference sequence (Fig. 6). It therefore represents an important candidate gene for flowering time control. As a subunit of the plant mediator complex, the most important targets of PFT1 concerning flowering time control have been shown to be the important flowering genes CONSTANS (CO) and Flowering Locus T (FT). PFT1 serves as an independent transcriptional activator of both these genes (Iñigo et al. 2012). The QTL LOD value is highest around the region of VvPFT1 for the onset of flowering. VvPFT1 may initiate a pathway leading to the activation of CO and FT, subsequently followed by the action of other genes.
During further development of the flowers, the highest LOD value of the QTL switches slightly to a region further downstream around markers GF01-19/GF01-20 (Fig. 6). This region can also be found in the second mapping population with a QTL for the same trait. Within the associated 1-LOD confidence interval, several other genes are annotated that might also have an impact on flowering time and flower development. However, the most interesting candidates are probably two CONSTANS-like genes, situated within the LOD ± 1 confidence interval near markers VMC4f9.2 and GF01-19, respectively. In addition to several other zinc finger transcription factors, CONSTANS-like genes are known to play a key role in the control of flowering via the photoperiodic pathway in A. thaliana (Putterill et al. 1995). Two CONSTANS homologous genes have yet been identified in grapevine (VvCO and VvCOL1) and their expression patterns characterized (Almada et al. 2009). A spatial and temporal relationship in the expression of VvCO, VvFY and VvMADS8 (equivalent to the Arabidopsis genes LEAFY and SOC1) was observed during bud development. Taking into account that the respective QTL region is found in two genetically highly diverse mapping populations, an influence of the two CONSTANS-like genes on flowering time near the QTL peaks identified in this study seems very likely.
The gene VvGAI1 (gibberellic acid-insensitive mutant protein 1) is located in the close proximity of the QTL peak, showing 63.5 % identity (76.5 % positives) to the RGL2 gene from A. thaliana. Grapevine plants with a mutation in this gene have an altered response to gibberellic acid causing a dwarf habitus and the production of numerous inflorescences while normal tendril development is suppressed (Boss and Thomas 2002). Reversely speaking, one of the functions of gibberellic acid signals in grapevine seems to be the suppression of floral meristem production. This gene therefore represents another candidate gene for the control of the onset of flowering in grapevine.
Near the second QTL peak around marker GF01-22 (Fig. 6), several other annotated genes are present in the reference genome sequence that might function in flower development or flowering time control. These include an auxin-independent growth promoter, the gene for protein SPIRAL1-like 2, RPT2 and the gene encoding the homeobox protein BEL1 homolog KNOX (Fig. 6). Several transporter protein genes and genes for non-defined transcription factors can be assigned to this region (Tables 5, 6). However, in this region well-known flowering time control genes like CO or FT do not appear in the reference genome sequence and we could not confirm this locus in the second mapping population. Further studies will thus be necessary to confirm whether these or other yet unknown genes contribute to regulation of flowering and ripening time in Vitis.
This study presents comprehensive data to develop new markers linked to phenological traits in marker-assisted breeding. We demonstrated that due to high synteny within the genus Vitis, it is easily possible to develop markers for fine mapping of specific loci. This relies on the availability of a high-quality reference genome sequence that allows a detailed deduction of candidate genes. Several of the QTLs for flowering time presented in this study were indentified using two independent mapping populations, permitting to identify functional candidate genes to be tested in the near future. A next step will be to unravel in more detail the pathways of phenological appearance in the field, a process that still is a challenge when addressing a highly complex perennial species such as grapevine. Our strategy proved suitable for the development of markers for extended marker-assisted selection, thus opening up new possibilities for fast and efficient modern grapevine breeding.
Author contributions
MD and IF carried out the phenotyping and genotyping steps. LH developed the concept of marker design based on the reference genome sequence. IF designed and carried out the data evaluation, construction of genetic maps for V3125 × ‘Börner’ and subsequent QTL analyses, DH and IF analyzed the sequences regarding candidate gene prediction. EZ provided the genetic map for GF.GA-47-42 × ‘Villard blanc’. LH, BW and RT designed the project and helped discuss and guide through all steps of the experiments. All the authors have read and approved the final manuscript.
References
Adam-Blondon AF, Roux C, Claux D, Butterlin G, Merdinoglu D, This P (2004) Mapping 245 SSR markers on the Vitis vinifera genome: a tool for grape genetics. Theor Appl Genet 109:1017–1027
Adam-Blondon AF, Jaillon O, Vezzulli S, Zharkikh A, Troggio M, Velasco R (2011) Genome sequence initiatives. In: Adam-Blondon AF, Martinez-Zapater JM, Kole C (eds) Genetics, genomics and breeding of grapes. Science Publishers, Enfield, pp 211–234
Almada R, Cabrera N, Casaretto JA, Ruiz-Lara S, González Villanueva E (2009) VvCO and VvCOL1, two CONSTANS homologous genes, are regulated during flower induction and dormancy in grapevine buds. Plant Cell Rep 28:1193–1203
Andres F, Coupland G (2012) The genetic basis of flowering responses to seasonal cues. Nat Rev Genet 13(9):627–639
Blazquez MA, Weigel D (2000) Integration of floral inductive signals in Arabidopsis. Nature 404:889–892
Blazquez MA, Green R, Nilsson O, Sussman MR, Weigel D (1998) Gibberellins promote flowering of Arabidopsis by activating the LEAFY promoter. Plant Cell 10:791–800
Boss PK, Thomas MR (2002) Association of dwarfism and floral induction with a grape “green revolution” mutation. Nature 416:847–850
Boss PK, Sreekantan L, Thomas MR (2006) A grapevine TFL1 homologue can delay flowering and alter floral development when overexpressed in heterologous species. Funct Plant Biol 33:31–41
Bowers JE, Meredith CP (1996) Genetic similarities among wine grape cultivars revealed by restriction fragment length polymorphism (RFLP) analysis. J Am Soc Hortic Sci 121:620–624
Bowers JE, Dangl GS, Meredith CP (1999) Development and characterization of additional microsatellite DNA markers for grape. Am J Enol Vitic 50:243–246
Cabezas JA, Cervera MT, Ruiz-Garcia L, Carreno J, Martínez-Zapater JM (2006) A genetic analysis of seed and berry weight in grapevine. Genome 49:1572–1585
Calonje M, Cubas P, Martinez-Zapater JM, Carmona MJ (2004) Floral meristem identity genes are expressed during tendril development in grapevine. Plant Physiol 135:1491–1501
Carmona MJ, Cubas P, Martinez-Zapater JM (2002) VFL, the grapevine FLORICAULA/LEAFY ortholog, is expressed in meristematic regions independently of their fate. Plant Physiol 130:68–77
Carmona MJ, Calonje M, Martinez-Zapater JM (2007) The FT/TFL1 gene family in grapevine. Plant Mol Biol 3:637–650
Carmona MJ, Chaib J, Martínez-Zapater JM, Thomas MR (2008) A molecular genetic perspective of reproductive development in grapevine. J Exp Bot 59(10):2579–2596
Cipriani G, Marrazzo MT, Di Gaspero G, Pfeiffer A, Morgante M, Testolin R (2008) A set of microsatellite markers with long core repeat optimized for grape (Vitis spp.) genotyping. BMC Plant Biol 8:127
Costantini L, Battilana J, Lamaj F, Fanizza G, Stella Grando M (2008) Berry and phenology-related traits in grapevine (Vitis vinifera L.): from quantitative trait loci to underlying genes. BMC Plant Biol 8:38
Dalbó MA, Ye GN, Weeden NG, Steinkellner H, Sefc KM, Reisch BI (2000) A gene controlling sex in grapevines placed on a molecular-based genetic map. Genome 43:333–340
Di Gaspero G, Cipriani G, Marazzo M, Andreetta D, Prado Castro M, Peterlunger E, Testolin R (2005) Isolation of (AC)n-microsatellites in Vitis vinifera L. and analysis of genetic background in grapevines under marker assisted selection. Mol Breed 15:11–20
Doligez A, Bouquet A, Danglot Y, Lahogue F, Riaz S, Meredith CP, Edwards KJ, This P (2002) Genetic mapping of grapevine (Vitis vinifera L.) applied to the detection of QTLs for seedlessness and berry weight. Theor Appl Genet 105:780–795
Duchene E, Butterlin G, Dumas V, Merdinoglu D (2012) Towards the adaptation of grapevine varieties to climate change: QTLs and candidate genes for developmental stages. Theor Appl Genet 124:623–635
Ehrenreich IM, Hanzawa Y, Chou L, Roe JL, Kover PX, Puruggananet MD (2009) Candidate gene association mapping of Arabidopsis flowering time. Genetics 183:325–335
Fechter I, Hausmann L, Daum M, Rosleff Sörensen T, Viehöver P, Weisshaar B, Töpfer R (2012) Candidate genes within a 143 kb region of the flower sex locus in Vitis. Mol Genet Genomics 287:247–259
Fernandez L, Ageorges A, Torregrosa L (2006) A putative NAP homolog specifically expressed during grapevine flower and berry development. Vitis 45(1):51–52
Fischer BM, Salakhutdinov I, Akkurt M, Eibach R, Edwards KJ, Töpfer R, Zyprian EM (2004) Quantitative trait locus analysis of fungal disease resistance factors on a molecular map of grapevine. Theor Appl Genet 108:501–515
Grattapaglia D, Sederoff R (1994) Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: mapping strategy and RAPD markers. Genetics 137:1121–1137
Iñigo S, Alvarez MJ, Strasser B, Califano A, Cerdán PD (2012) PFT1, the MED25 subunit of the plant mediator complex, promotes flowering through CONSTANS dependent and independent mechanisms in Arabidopsis. Plant J 69(4):601–612
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Le Clainche I, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pè ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P; French-Italian Public Consortium for Grapevine Genome Characterization (2007) The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature 449:463–467
Joly D, Perrin M, Gertz C, Kronenberger J, Demangeat G, Masson JE (2004) Expression analysis of flowering genes from seedling-stage to vineyard life of grapevine cv. Riesling. Plant Sci 166:1427–1436
Kosambi DD (1944) The estimation of map distances from recombination values. Ann Eugenn 12:172–175
Langridge J (1957) Effect of day-length and gibberellic acid on the flowering of Arabidopsis. Nature 180:36–37
Le Paslier M-C, Choisne N, Bacilieri R, Bounon R, Boursiquot J-M, Bras M, Brunel D, Di Gaspero G, Hausmann L, Lacombe T, Laucou V, Launay A, Martinez-Zapater JM, Morgante M, Raj PM, Ponnaiah M, Quesneville H, Scalabrin S, Torres-Perez R, Adam-Blondon A-F (2013) The GrapeReSeq 18k Vitis genotyping chip. In: 9th International symposium grapevine physiology and biotechnology, 21–26 April 2013, La Serena, p 123
Lemke L, Rex M, Zyprian E, Töpfer R (2011) A simple, inexpensive and environmentally friendly method for high throughput DNA extraction from grapevine (Vitis spp.). Vitis 50:7–10
Lorenz DH, Eichhorn KW, Bleiholder H, Klose R, Meier U, Weber E (1995) Phenological growth stages of the grapevine, Vitis vinifera L. ssp. vinifera. Codes and descriptions according to the extended BBCH scale. Aust J Grape Wine Res 1(2):100–103
Lowe KM, Walker MA (2006) Genetic map of the interspecific grape rootstock cross Ramsey (Vitis champinii) × Riparia Gloire (Vitis riparia). Theor Appl Genet 112:1582–1592
Marguerit E, Boury C, Manicki A, Donnart M, Butterlin G, Nemorin A, Wiedemann-Merdinoglu S, Merdinoglu D, Ollat N, Decroocq S (2009) Genetic dissection of sex determinism, inflorescence morphology and downy mildew resistance in grapevine. Theor Appl Genet 118:1261–1278
Mejía N, Gebauer M, Munoz L, Hewstone N, Hinrichsen P (2007) Identification of QTLs for seedlessness, berry size, and ripening date in a seedless × seedless table grape progeny. Am J Enol Vitic 58:499–507
Mejía N, Soto B, Guerrero M, Casanueva X, Houel C, de los Ángeles Miccono M, Ramos R, Le Cunff L, Boursiquot JM, Hinrichsen P, Adam-Blondon AF (2011) Molecular, genetic and transcriptional evidence for a role of VvAGL11 in stenospermocarpic seedlessness in grapevine. BMC Plant Biol 11:57
OIV (2009) OIV descriptor list for grape varieties and Vitis species, 2nd edn. Organisation Internationale de la Vigne et du Vin, Organisation Intergouvernementale, Paris
Putterill J, Robson F, Lee K, Simon R, Coupland G (1995) The CONSTANS gene of Arabidopsis promotes the flowering and encodes a protein showing similarities to zinc finger transcription factors. Cell 80:847–857
Riaz S, Dangl GS, Edwards KJ, Meredith CP (2004) A microsatellite marker based framework linkage map of Vitis vinifera L. Theor Appl Genet 108:864–887
Sefc KM, Regner F, Tureschek E, Glo J, Steinkellner H (1999) Identification of microsatellite sequences in Vitis riparia and their application for genotyping of different Vitis species. Genome 42:367–373
Song YH, Ito S, Imaizumi T (2013) Flowering time regulation: photoperiod- and temperature-sensing in leaves. Trends Plant Sci 18(10):575–583
Sreekantan L, Thomas MR (2006) VvFT and VvMADS8, the grapevine homologues of the floral integrators FT and SOC1, have unique expression patterns in grapevine and hasten flowering in Arabidopsis. Funct Plant Biol 33:1129–1139
Töpfer R, Hausmann L, Harst M, Maul E, Zyprian E, Eibach R (2011) New horizons for grapevine breeding. In: Flachowsky H, Hanke M-V (eds) Methods in temperate fruit breeding. fruit, vegetable and cereal science and biotechnology, vol 5 (special issue 1), pp 79–100
Van Ooijen JW (2004) MapQTL®5, software for the mapping of quantitative trait loci in experimental populations. Kyazma B.V., Wageningen
Van Ooijen JW (2006) JoinMap®4, software for the calculation of genetic linkage maps in experimental populations. Kyazma B.V., Wageningen
Velasco R, Zharkikh A, Troggio M et al (2007) A high quality draft consensus sequence of the genome of a heterozygous grapevine variety. PLoS One 2:e1326
Zhang J, Hausmann L, Eibach R, Welter L, Toepfer R, Zyprian E (2009) A framework map from grapevine V3125 (Vitis vinifera ‘Schiava Grossa’ × ‘Riesling’) × rootstock cultivar ‘Börner’ (Vitis riparia × Vitis cinerea) to localize genetic determinants of phylloxera root resistance. Theor Appl Genet 119:1039–1051
Zyprian E, Eibach R, Töpfer T (2006) Eine neue genetische Karte der Weinrebe aus der Kreuzung “Gf.Ga-47-42” × “Villard blanc”. Deutsch Weinbau Jahrb 2006:151–158
Acknowledgments
We gratefully thank the German Federal Ministry of Education and Research (BMBF) for funding this project (Grant No. 0315460A and 0315460B). Additional support came from funding via the priority program SPP 1530 of the German research foundation (DFG) (Grant No. TO 152/5-1 and WE 1576/16-1).
Our further thanks go to Andreas Preiss for technical assistance, Dr. Friederike Rex for helpful discussions and Dr. Rudolf Eibach and his team for grapevine breeding.
Conflict of interest
The authors declare that they have no conflict of interest.
Ethical standard
The authors declare that the experiments of this study comply with the current laws. We confirm to have the authority to publish this work and that the manuscript has not been published before and is not under consideration for publication elsewhere.
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by A. Charcosset.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.
About this article
Cite this article
Fechter, I., Hausmann, L., Zyprian, E. et al. QTL analysis of flowering time and ripening traits suggests an impact of a genomic region on linkage group 1 in Vitis . Theor Appl Genet 127, 1857–1872 (2014). https://doi.org/10.1007/s00122-014-2310-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00122-014-2310-2