Skip to main content

A new image-based tool for the high throughput phenotyping of pollen viability: evaluation of inter- and intra-cultivar diversity in grapevine

Abstract

Background

Low pollen viability may limit grapevine yield under certain conditions, causing relevant economic losses to grape-growers. It is usually evaluated by the quantification of the number of viable and non-viable pollen grains that are present in a sample after an adequate pollen grain staining procedure. Although the manual counting of both types of grains is the simplest and most sensitive approach, it is a laborious and time-demanding process. In this regard, novel image-based approaches can assist in the objective, accurate and cost-effective phenotyping of this trait.

Results

Here, we introduce PollenCounter, an open-source macro implemented as a customizable Fiji tool for the high-throughput phenotyping of pollen viability. This tool splits RGB images of stained pollen grains into its primary channels, retaining red and green color fractionated images (which contain information on total and only viable pollen grains, respectively) for the subsequent isolation and counting of the regions of interest (pollen grains). This framework was successfully used for the analysis of pollen viability of a high number of samples collected in a large collection of grapevine cultivars. Results revealed a great genetic variability, from cultivars having very low pollen viability (like Corinto Bianco; viability: 14.1 ± 1.3%) to others with a very low presence of sterile pollen grains (Cuelga; viability: 98.2 ± 0.5%). A wide range of variability was also observed among several clones of cv. Tempranillo Tinto (from 97.9 ± 0.9 to 60.6 ± 5.9%, in the first season). Interestingly, the evaluation of this trait in a second season revealed differential genotype-specific sensitivity to environment.

Conclusions

The use of PollenCounter is expected to aid in different areas, including genetics research studies, crop improvement and breeding strategies that need of fast, precise and accurate results. Considering its flexibility, it can be used not only in grapevine, but also in other species showing a differential staining of viable and non-viable pollen grains. The wide phenotypic diversity observed at a species level, together with the identification of specific cultivars and clones largely differing in this trait, pave the way of further analyses aimed to understand the physiological and genetic causes driving to male sterility in grapevine.

Background

The commercial production of wine and table grapes (Vitis vinifera L.) is mostly based on the cultivation of plants with hermaphroditic flowers, in which self-fertilization is the major route for pollination. Fruit-set rate (i.e. proportion of ovaries in an inflorescence that become fruits) can be highly variable in grapevine [1], and it is generally considered as normal if greater than 50%, and insufficient if less than 30% [2]. In fact, poor fruit-set is said to limit crop yield under certain conditions, causing relevant economic losses to grape-growers [1]. Fruit-set is affected by many factors, including climate (e.g. light, temperature, rain), vineyard management (e.g. pruning system, row configuration, training system) and genetic factors [3,4,5,6,7]. Regarding the latter, low pollen viability [5, 8], low ovule fertility [6], low germination rate [9], or the development of anomalous flowers [10] are some of the reasons explaining poor fruit-set and the excessive abscission of flowers (termed shatter, shedding or by the French word “coulure”) [11].

Grapevine pollen formation and development follow a time-course process [12]. Under optimal climatic conditions, sporogenous cells can be detected in the anthers as soon as inflorescences are clearly visible in the growing shoot. Meiosis takes place after some days, giving rise to tetrads which rapidly release microspores in the loculus. Microspores vacuolate thereafter, and pollen mitosis take place just before anthesis, when pollen grains are mature enough to be released from the anthers [13, 14]. In this regard, genotype-specific differences in pollen starch reserves affect this developmental process, having an effect in the successful obtaining of fully mature pollen grains at anthesis [13]. Additionally, pollen production is highly variable between grapevine cultivars. In a work performed in a set of table grape cultivars, Kelen et al. [7] report a variation ranging between ca. 3000 and 9000 pollen grains per flower. These figures suggest that pollination and fruit set might be ensured even at very low rates of pollen viability. Nevertheless, pollen sterility may be accompanied by a certain degree of ovule sterility [15], becoming a more limiting factor for successful berry set and subsequent grape production, since there are only four ovules per flower [11], and at least one fertilized ovule is needed for a normal berry set.

Pollen viability indicates the ability of the pollen grain to deliver sperm cells to the embryo sac following compatible pollination [16]. Consequently, it can be considered as a measurement of pollen grain quality and vigor [17]. In this regard, a positive relationship between pollen viability and pollen germination capability has been suggested [7], and Royo et al. [18] have recently reported the incapability of the non-viable Corinto Bianco grapevine pollen grains to germinate in vitro, indicating that viable pollen is needed for a successful germination and so for seed and fruit development and growth. Pollen viability has been studied by different methods, including the evaluation of seed- and fruit-set, diverse in vitro and in vivo germination tests, and by different staining techniques. Alexander’s staining [19] is one of the most common methods, since it allows the distinction between viable and non-viable (sterile, aborted) pollen grains based on the differential staining of the pollen protoplasm and the cellulose contained in pollen walls. The method has been suggested to be appropriate for large scale breeding programs, as it correlates with pollen germination rates and it is simple and easy to use [20]. After pollen grains staining, viability can be rated by the manual (eye) counting of non-viable (stained in pale turquoise blue) and viable (stained dark blue or purple) pollen grains. Although this manual approach is the simplest and most sensitive one to evaluate pollen viability, it is a slow and laborious process that might result in diminished accuracy and precision (reproducibility) rates. In an effort to automate this task, Kelly et al. [21] suggested the use of pollen grain size variation as an indirect evaluation of pollen viability, considering that viable grains are considerably bigger than non-viable grains. Authors found a significant positive correlation between the proportion of viable and non-viable grains (after differential staining), and the so-called pollen size index (PSI), which relativizes the number of pollen grains above a certain size by the total number of grains. This approach has been used for the automatic estimation of pollen viability in species like Mimulus guttatus, Collinsia verna and Beta vulgaris [21, 22].

Image-based phenotyping is revolutionizing many areas of plant research [23]. The use of high-throughput systems is being used to extract biological information from image data for subsequent decision support. The analysis of microscope preparations can be sped-up by the collection of color digital microscope images (RGB images) for their further automate processing with specific image analysis systems [24]. As an example, an image processing workflow for the automated analysis of images of stained sections of maize internodes has been recently proposed to discriminate and quantify lignified and non-lignified tissues [25]. In this regard, the development of objective, accurate and rapid cost-effective methods for the high-throughput evaluation of pollen viability will be helpful in different palynological areas, including genetics research, crop improvement and breeding studies [17]. In addition, it will be useful to evaluate the effect of storage on pollen viability and longevity in genebanks [26], and to assess the effect of different environmental factors on pollen performance [27]. Although image based technologies have been already applied in palynological analyses [28,29,30], their use has been mainly restricted to the automatic counting of pollen grains of several plant species, mostly for their morphological classification according to predefined categories.

Here, we introduce a customizable macro developed for the open-source Fiji image-analysis platform [31] for the automatic analysis of digital microscope RGB images of pollen grains after Alexander’s modified staining [32], as a solution to enable an automated and fast pollen viability phenotyping. This macro, PollenCounter, splits digital images into its primary channels, and isolates the regions of interest (pollen grains) from the background in the red and green color fractionated images to address the independent counting of total and viable pollen grains. The PollenCounter macro has been designed and validated in a sample of grapevine cultivars that show a high degree of pollen viability variability, providing accurate and precise information in all the tested scenarios. In addition, it has been applied to evaluate if there is some difference in pollen viability within different inflorescences of the same clone, as well as to show the range of variability among different clones of cv. Tempranillo Tinto and a high number of table- and wine-grape cultivars.

Methods

Plant material

Grapevine cultivars and flowers sampling

In this study, 120 grapevine accessions have been analyzed (Additional file 1). These accessions belong to four different entities, and they were studied on their own experimental site: (I) the Grapevine Germplasm Collection of the Instituto de Ciencias de la Vid y del Vino (ICVV; FAO Institute Code ESP-217), in Logroño (La Rioja, Spain); (II) the Experimental Grapevine Collection of the Universität für Bodenkultur Wien (BOKU), in Tulln an der Donau (Niederösterreich, Austria); (III) the cv. Tempranillo Tinto clones collection of a commercial grapevine nursery (VP, Viveros Provedo S.A.), also in Logroño; and (IV) a private vineyard of cv. Tempranillo Blanco (FZ), in Fonzaleche (La Rioja, Spain). Accessions from BOKU were used to evaluate sampling effect, accessions from ICVV were used to evaluate the inter-cultivar genetic diversity, and Tempranillo Tinto clones from VP were used to evaluate the intra-cultivar genetic variability. In these four plots, plants are maintained following standard conditions. Cultivar identity was confirmed by SSRs and/or SNPs analyses [33]. If available, cultivar prime name according to the VIVC database (http://www.vivc.de) has been used in the manuscript.

BOKU sampling

To evaluate within-clone differences in pollen viability, flowers from the basal and second-order inflorescences from the proximal and distal shoots were collected early in the morning in three plants of the same clone. The analysis was done one season (2017) in four wine cultivars (Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc) (Additional file 1). In BOKU plot, plants are trained as a Guyot system with a single cane of 8 buds and one spur of 2 retained buds.

ICVV, FZ and VP sampling

For the inter- and intra-cultivar analyses of pollen viability variability, the first inflorescences from a random shoot from two plants were sampled early in the morning at ICVV, VP and FZ plots in 2015. This study involved 116 grapevine accessions (111 in ICVV, 4 in VP and 1 in FZ). This analysis was repeated in 2017 on 24 selected varieties (ICVV) and on Tempranillo Tinto clones (VP) (Additional file 1). In ICVV, VP and FZ plots, plants are trained on a Royat system, composed of a double spur-pruned permanent cordon with 6 spurs of 2 buds per vine.

Pollen sample preparation and imaging

To determine pollen viability, inflorescences were collected at full bloom (modified E-L 23 stage [34]), when 50% of flower caps have fallen. After collection, they were transported to the laboratory and maintained at room temperature until their analysis, which was finished within the same day. Recently opened flowers (with erect stamen filaments and light yellowish anthers) from the central section of the inflorescence were then collected and stained following Alexander’s modified staining [32], which can differentiate between aborted (non-viable) and non-aborted (viable) grains (Fig. 1a). Briefly, three-four flowers were selected, soaked in 40 μL of staining solution and shaken vigorously for 15 s to facilitate pollen release. 20 μL of the solution were then transferred onto a pre-heated microscope slide and observed to take digital images. Samples from ICVV, FZ and VP samples were observed using a Zeiss SteREO Discovery V20 stereo-microscope, and images were taken using a Zeiss AxioCam camera and AxioVision software (v. 4.8, Zeiss). Samples from BOKU were visualized with an Olympus CX41 microscope equipped with an Olympus SC50 digital camera using CellSens Entry software (v. 1.14, Olympus). Image contrast and saturation was adapted to maximize the differentiation between aborted and non-aborted pollen grains. In total, 4126 RGB images saved as JPEG files with an average number of 130 pollen grains per image were taken for further analysis.

Fig. 1
figure 1

PollenCounter pipeline for the automated detection and counting of total and viable pollen grains. In a, an input RGB image of grapevine pollen grains after Alexander’s modified staining is shown (dark pollen grains are viable, whilst pale grains are sterile). This image is firstly segmented into its primary colors, and the color fractionated grayscale pictures for the red (b) and green (c) channels are taken for their parallel analysis (red and green image frame lines, respectively). Images are then segmented into black regions of interest (pollen grains) and a homogenous white background (d, e). Regions of interest are closed and filled with black pixels (f, g) before the use of a watershed algorithm for the separation of touching regions of interest (h, i). Identified pollen grains are then counted, and output images (j, k) are automatically saved as.jpg files. Note that in these images, counted pollen grains are automatically colored in cyan. In l, an example of the output pivot table with the counting results in both channels is shown. The main commands implemented in PollenCounter are indicated in yellowish boxes (grey arrow). For simplification, only a representative area of one of the processed images from cv. Vermentino is presented

Automated analysis of pollen images for pollen grains counting

The workflow designed for the automatic counting of viable and non-viable pollen grains is shown in Fig. 1. The open-source Fiji platform [31] was used for macro development and subsequent image analysis. It is based on the popular open-source software focused on biological-image processing ImageJ v. 1.51 (US National Institutes of Health). This sequential analysis was executed in batch mode on the entire set of images stored in a predefined directory. Briefly, each input RGB image is split into its three primary color channels (red, green and blue) as color-fractionated grayscale pictures. Considering that they hold independent information regarding total and viable pollen grains (Fig. 1b, c), only red and green channels are retained for their analysis. The threshold of both color-channel images is automatically set to segment images into background and individual regions of interest (in our case, pollen grains), and every image is black and white-converted. As a result, images with pollen grains transformed to black regions of interest and a homogenous white background are obtained (Fig. 1d, e). Next, the Dilate and Fill holes operations are applied to generate solid particles, by adding black pixels to the edges and inner sections of the black regions of interest. Additional pixels from the edges are then removed using the Erode operation (Fig. 1f, g). The automatic separation of touching regions is performed with a watershed segmentation algorithm. It calculates the Euclidean distance map and looks for the ultimate eroded points, which are then dilated as far as possible (either until the edge of the particle is reached, or it touches another region of interest) (Fig. 1h, i). Regions of interest were then measured (number, surface areas and diameters) using the Analyze Particles command (Fig. 1j,k). In this regard, only particles with a surface area ranging between 60 and 800 pixels2 were counted in the red-channel images (Pollen grainsR), whereas only those with a surface area between 100 and 800 pixels2 were considered in the green-channel images (Pollen grainsG). For both channels, circularity (calculated from diameter data) range was set to 0.40–1.00. After processing, automatic counting results are shown in a pivot table that can be saved by user as a .txt or .xls file (Fig. 1l). Additionally, two.JPEG files corresponding to the analysis of red and green channels are created in the directory. In these images, the detected and counted pollen grains are shown in cyan and the discarded structures in black.

The custom Fiji macro for the automatic counting of total and viable pollen grains, called PollenCounter, and a set of images are available as Additional files 2 and 3.

Evaluation and validation of the automatic approach

For the validation of the novel approach, a subset of 392 RGB images from 19 cultivars grown in ICVV and FZ plots (Albillo Real, Bouschet Petit, Cuelga, Dominga, Gamay noir, Gewuerztraminer, Mencía, Molar, Pinot Noir, Planta Mula, Ruby Cabernet, Rubired, Sauvignon Blanc, Schiava Grossa, Silvaner Gruen, Syrah, Tempranillo Blanco, Trousseau Noir and Vijiriega Común) were chosen. These images were specifically selected for their variable number of pollen grains (from 20 to 1042 pollen grains) and viability rate (from ca. 20 to 100%).

Manual analysis of pollen images for viability estimation

The number of viable and non-viable pollen grains was manually counted in each image considering their differential staining (pale blue = non-viable/dark blue = viable). Manual pollen viability per image [MPVimg (%)] was calculated as follows:

$$MPV_{img}\, \left( \% \right) = \frac{Dark\,blue\,pollen\,grains}{All\,pollen\,grains} \times 100$$

Automatic analysis of pollen images for viability estimation

Automatic pollen viability per image [APVimg (%)] was calculated using the information obtained with the automatic approach as follows:

$$APV_{img}\, \left( \% \right) = \frac{{Pollen\,grains_{G} }}{{Pollen\,grains_{R} }} \times 100$$

Pollen size index (PSI) determination

As suggested by Kelly et al. [21], we calculated the pollen size index per image (PSIimg) in a subsample of 38 images from the 19 previously listed cultivars (2 images per cultivar). To this aim, the surface area of all pollen grains in each image was individually obtained (in pixels2) using the Analyze Particles command of Fiji, and viable and non-viable pollen grain surface areas were compared to determine the optimal surface area threshold capable to differentiate both types (Additional file 4A). Although the overlapping observed between the distribution of the area of viable and non-viable pollen grains suggested that this method would not be useful in grapevine, a general threshold of 118 pixels2 was established for PSI (%) calculation. This value corresponds to the average value of the area of the smallest viable pollen grain detected in each of the 38 images. PSI (%) was calculated as follows:

$$PSI_{img}\, \left( \% \right) = \frac{{Pollen\,grains\,over\,118\,pixels^{2} }}{Total\,pollen\,grains} \times 100$$

Grapevine pollen viability assessment

To estimate pollen viability in a sample (inflorescence), automatic pollen viability [APV (%)] was calculated as follows:

$$APV\,\left( \% \right) = \frac{{\sum Pollen\,grains_{G} }}{{\sum Pollen\,grains_{R} }} \times 100,$$

where ∑Pollen grainsG and ∑Pollen grainsR indicate the sum of automatically-counted pollen grains in all the green- and red-channel images taken from a sample, respectively. In average, more than 1000 pollen grains per sample were analyzed.

Statistical analysis

The performance of the automatic method was evaluated by comparing the automatic results (total counting and viable counting) with those manually obtained (n = 392 images) through linear regression models and associated coefficients of determination (R2). In addition, Bland–Altman plots were calculated as previously indicated [35] to estimate the accuracy of both approaches. Box-plots combined with independent t-tests were used to detect significant differences (p ≤ 0.05) between pollen viability (%) rates obtained manually [MPVimg (%)] and by the automatic method [APVimg (%)]. The performance of the PSI method was evaluated by comparing the automatic result [PSIimg (%)] with the manual value [MPVimg (%)], through a linear regression model and its associated coefficient of determination (R2) (n = 38 images).

To analyze within-clone differences in APV (%), BOKU samples were used (n = 384 images). First, a one-way ANOVA was individually calculated for each cultivar to evaluate if APV (%) varied significantly (p ≤ 0.05) between plants of the same clone (cultivar). Then, a three-way ANOVA with Fisher’s LSD post hoc tests was calculated in the whole set of samples to check if any of the three factors considered (genotype, shoot position and inflorescence order) played a significant effect (p ≤ 0.05) on APV (%). For this analysis, four levels were considered for genotype (Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc), two for shoot position (proximal and distal) and two for inflorescence position (first and second order) factors.

A Fisher’s LSD-test was performed to assess differences in APV (%) between Tempranillo Tinto clones, being considered significant at p ≤ 0.05.

The normality of the data was graphically assessed using histograms. All calculations were done using SPSS v. 21.0 (IBM, Chicago, IL, USA).

Results

Automatic pollen counting tool designing and validation

Two typical RGB images of two cultivars with low (cv. Albillo Real) and high (cv. Forcallat Tinta) pollen viability rates after Alexander’s modified staining approach [32] are shown in Fig. 2b, e. As previously reported [21, 22], a difference in size between viable and non-viable pollen grains was observed for the different cultivars evaluated (see Figs. 1a, 2b and Additional file 5 for some examples), with non-viable pollen grains slightly smaller. Nonetheless, the calculation of PSIimg yield a non-significant correlation with the manual evaluation of pollen viability (MPVimg; R2 = 0.15; p > 0.01, Additional file 4B) in a subsample of 38 images of 19 grapevine cultivars, and so its usefulness is limited when considering multi-cultivar genetic frameworks.

Fig. 2
figure 2

Original and processed images obtained for two typical pollen samples with low (ac, cv. Albillo Real) and high (df, cv. Forcallat Tinta) viability. Original images are shown in b and e. Images obtained in the red (a, d) and green (c, f) channels after PollenCounter use are shown with the number of pollen grains automatically identified (lower right corners). For simplification, only representative areas of the processed images are shown

On the other hand, viable and non-viable pollen grains can be easily identified according to their differential staining in samples with both types of pollen grains (Figs. 1a, 2b and Additional file 5). The three dyes used in the Alexander’s modified staining solution (malachite green, acid fuchsin and orange G) split into very different red, green and blue image profiles (see Additional file 6). In this regard, it is easy to identify malachite green dye (which stains cellulose in pollen walls, present in all pollen grains) in the red channel, as happens with acid fuchsin dye (which stains pollen protoplasm, only present in mature viable pollen grains) in the green one (Additional file 5). Orange G (clearly identified in the blue channel) reinforces the differentiation of the other two dyes, aiding in the distinction of viable and non-viable pollen grains. These differences pave the way to test the use of the red and green channels for the automatic differentiation and counting of total and viable pollen grains through an image-based approach. Consequently, we developed a processing tool, called PollenCounter, to be used in the open-source Fiji platform for the segmentation of the three channels in our set of images, retaining green and red layers (Fig. 2) for their independent analysis, as detailed in the pipeline indicated in “Methods” (Fig. 1).

The use of PollenCounter allowed the differential counting of total and viable pollen grains with absolute reproducibility while substantially reducing the time needed for the analysis. Considering the subset of 392 images, an average of 240 s of dedicated work were needed per image for the manual counting of viable and non-viable pollen grains, whereas this time was reduced to 3.5 s when PollenCounter was run in an IntelCore i3 laptop.

PollenCounter accuracy was evaluated by comparing the number of total and viable pollen grains counted manually and automatically in the individual images taken from 19 cultivars (n = 392, Fig. 3). A significant correlation between manual and automatic values (R2 = 0.98, p ≤ 0.01) was obtained both for total and viable pollen grains, and regression lines nearly match identity lines (x = y) (Fig. 3a, c). Similarly, the Bland and Altman approach showed a high agreement between the automatic and manual number of total and viable pollen grains, with a mean difference of 2.9 and 1.4 pollen grains between both systems, respectively (Fig. 3b, d). The individual differences between these two counting systems were well distributed within the interval limits and no bias was observed, especially in the range between 30 and 284 pollen grains, which constitute percentile 10 and 90 of the whole set of images analyzed, respectively. Moreover, 95% confidence intervals were small enough to sustain that the automatic method could substitute the manual one. Our approach also yield high correlation results (R2 = 0.98; p ≤ 0.01) between APVimg (%) estimations and MPVimg (%) values in the same subsample of images (n = 38) used for testing the usefulness of the pollen size index (Additional file 4C).

Fig. 3
figure 3

Correlation analyses (a, c) and Bland–Altman plots (b, d) comparing manual (MC) and automatic (AC) counting of total (a, b) and viable (c, d) pollen grains. In a and c, red solid lines shows linear regression line, whereas the identity line (x = y) is shown as a blue solid line. In b and d, red solid lines indicate the mean difference between both approaches, while upper and lower solid black lines represent limits of agreement

On a next step, the validity of PollenCounter to estimate pollen viability was checked in the 19 cultivars to evaluate potential genotype-specific inconsistencies. In this regard, APVimg (%) and MPVimg (%) values were compared for each cultivar independently, using the information obtained from each processed image. No significant differences were found on pollen viability rates between the manual and the automatic approach in any of the grape cultivars (Fig. 4), even when analyzing genotypes with very low (e.g. Tempranillo Blanco) or very high (e.g. Cuelga) pollen viability.

Fig. 4
figure 4

Viability rates (%) obtained using automatic (APVimg, blue) and manual (MPVimg, yellow) pollen grain counting in 19 grapevine cultivars (AR: Albillo Real; BP: Bouschet Petit; CU: Cuelga; DM: Dominga; GN: Gamay noir; GW: Gewuerztraminer; MN: Mencia; MO: Molar; PM: Planta Mula; PN: Pinot Noir; RC: Ruby Cabernet; RU: Rubired; SB: Sauvignon Blanc; SG: Schiava Grossa; SI: Silvaner Gruen; SY: Syrah; TB: Tempranillo Blanco; TN: Trousseau Noir; VC: Vijiriega Común). No significant differences between manual and automatic approaches were observed for any cultivar (t test, p > 0.05). Whiskers indicate standard deviations between images

As a whole, PollenCounter provided precise, accurate and valid counts of total and viable grains from pollen samples stained using Alexander’s protocol, which could be used to estimate pollen viability rate in diverse grapevine cultivars. Consequently, pollen viability results presented in the next sections have been obtained using PollenCounter.

Evaluation of between- and within-vines pollen viability differences

One way ANOVA revealed no significant differences in APV (%) between the three plants sampled for any of the cultivars assayed: Cabernet Franc (p = 0.10), Chardonnay (p = 0.16), Gewuerztraminer (p = 0.09) and Sauvignon Blanc (p = 0.81). In a more detailed evaluation, we tested the effect of shoot position (proximal/distal) and inflorescence order (first/second) on APV (%). The three-way ANOVA indicated that genotype is the only factor contributing significantly to variability (F = 45.16; p ≤ 0.05), and no significant effect was released for shoot position (F = 0.06; p = 0.80) nor inflorescence order (F = 0.25; p = 0.62). Similarly, no significant two-way or three-way interactions for the factors analyzed were obtained. Additional Fisher’s LSD post hoc tests revealed a significant (p ≤ 0.05) lower APV (%) values for Chardonnay inflorescences (in average, 88.4 ± 3.9%) compared to the other three cultivars (96.7 ± 1.9, 98.0 ± 1.5 and 98.5 ± 1.2% for Cabernet Franc, Gewuerztraminer and Sauvignon Blanc inflorescences, respectively), but no significant differences between any of the inflorescences of the same genotype were found (p > 0.05) (Fig. 5).

Fig. 5
figure 5

Automatic pollen viability rates [APV (%)] obtained in different inflorescences of four different grapevine cultivars with PollenCounter. Different lowercase letters indicate significant differences between cultivars (p ≤ 0.05). P1: First inflorescence of the proximal shoot (yellow); P2: Second inflorescence of the proximal shoot (orange); D1: First inflorescence of the distal shoot (light blue); D2: second inflorescence of the distal shoot (dark blue). Whiskers indicate standard deviations between inflorescences

Pollen viability in grapevine clones and cultivars

On a next level of comparison, APV (%) values from 4 different clones of the same cultivar (Tempranillo Tinto) sampled in the same experimental plot (VP) were compared (Fig. 6). We found two Tempranillo clones with very high APV (%) values (RJ-51 and VP-2, > 90%), one clone with high APV (%) (VP-25, 75–90%), and one clone with medium APV (%) (VP-11, 50–75%). These APV (%) values were rather stable considering 2015 and 2017 seasons.

Fig. 6
figure 6

Automatic pollen viability rates [APV (%)] obtained for four Tempranillo clones (VP-11, VP-25, RJ-51, VP-2) in two different seasons—2015 (blue bars) and 2017 (yellow bars)—with PollenCounter. Different lowercase letters indicate significant differences between samples (p ≤ 0.05). Whiskers indicate standard deviations between inflorescences

Next, we compared pollen viability between different cultivars. Table 1 indicates the average and standard deviation of APV (%) obtained for the 111 accessions (107 cultivars) analyzed from the ICVV Grapevine Germplasm Collection in 2015 using PollenCounter. In general, we found high average APV (%) values for the cultivars analyzed in this work, with 81.0% of individuals (90 accessions) over 75%, and 49.5% (55 accessions) over 90%. In contrast, the lowest pollen viability value was found for the parthenocarpic cultivar Corinto Bianco, with an APV (%) of 14.1 ± 1.3%. In addition, we found other five accessions with average pollen viability lower than 50%: the interspecific hybrid Flot Rouge (35.1 ± 11.1), and the V. vinifera cultivars Planta Mula (44.9 ± 6.5), Molar (45.4 ± 9.8), Tempranillo Tinto (49.9 ± 0.8) and Vermentino (49.9 ± 3.2). Considering that they were grown under similar conditions in the same location and plot, the wide range of diversity observed between cultivars evidences a clear genetic effect on pollen formation and/or development, driving to different pollen viability rates. Additionally, no effect of the pruning system (Guyot vs. Royat) on flower vigor and pollen viability was observed when both systems were used for a same cultivar (Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc), although they were tested in two different locations and two different clones (Fig. 5 and Table 1).

Table 1 List of grapevine accessions (n = 111), corresponding to 107 different cultivars, sampled in the ICVV Grapevine Germplasm Collection for this study

Lastly, 24 cultivars were selected to phenotype this trait in a second season (2017). They included cultivars which in 2015 showed either low/very low pollen viability (like Corinto Bianco, Flot Rouge and Planta Mula) or high pollen viability difference between inflorescences (like Dominga, Maturana Blanca, Ondarrabi Beltza), as well as some randomly chosen cultivars with very high pollen viability (like Bobal, Castelao, Cabernet Sauvignon). On average, an increase in pollen viability (9.05%) was observed in this set of cultivars in 2017 compared to 2015, although important differences between cultivars were obtained (Fig. 7). They can be split into two groups according to their variation between seasons. One group includes 8 cultivars (Flot Rouge, Muscat a Petits Grains Blancs, Ondarrabi Beltza, Planta Mula, Tempranillo Tinto, Trousseau Noir, Verdil and Vijiriega Común) with an APV (%) difference between seasons larger than 10%. The rest of cultivars (Airén, Albillo Real, Bobal, Cabernet Sauvignon, Castelao, Colombard, Corinto Bianco, Dominga, Maturana Blanca, Mollar Cano, Naparo, Palomino de Jerez, Palomino Fino, Planta Nova, Tempranillo Blanco and Vermentino) show a lower difference between seasons. In the second group, we found diverse cultivars in terms of their APV (%) value, including very low: Corinto Bianco; low: Vermentino; medium: Albillo Real, Dominga, Mollar Cano and Tempranillo Blanco; high: Colombard, Naparo and Palomino de Jerez; and very high: Airén, Bobal, Cabernet Sauvignon, Castelao and Planta Nova.

Fig. 7
figure 7

Automatic pollen viability rates [APV (%)] for 24 grapevine cultivars in 2015 (blue) and 2017 (red) obtained with PollenCounter. Dots indicate average values, whereas whiskers show standard deviations between inflorescences

Discussion

PollenCounter: a new tool for pollen viability phenotyping

Phenotyping is a critical component of plant genetic research and crop improvement. Accurate data collection, if combined with genetic information, can greatly accelerate progress in breeding for yield and quality traits of new and better adapted cultivars [23]. Recent works have reported the use of high-throughput technologies for the phenotyping of diverse grapevine traits, including plant phenology [36], crop yield components [36,37,38], grape quality [35, 39, 40] and fungal disease resistance [41,42,43]. In this context, the use of efficient and objective image-based systems for pollen viability estimation may provide an alternative solution for the time-consuming phenotyping of this trait. Noticing the variable size of viable and non-viable pollen grains in M. guttatus and C. verna, Kelly et al. [21] proposed the use of the pollen size index (PSI) as an indicator of pollen viability. In grapevine, we also observed a smaller size in non-viable pollen grains compared to viable ones, but the use of the PSI indicator yielded a non-significant correlation with the manual measurement of pollen viability (Additional file 4B), probably due to the substantial genetic variation in pollen grains size. As an example, the area of the smallest viable pollen grain ranged from 100 to 181 pixels2 (for cvs. Syrah and Vijiriega Común, respectively) in a subsample of 38 images from 19 cultivars. Although rather small, the detailed analysis of this subsample was enough to confirm that the establishment of a global surface area threshold for the automatic classification of pollen grains as viable or non-viable is not satisfactory enough to automate PSI calculation (and its subsequent pollen viability estimation) in grapevine wide genetic frameworks. Consequently, a genotype-specific evaluation of viable and non-viable pollen grains surface area distribution will be needed to set case-to-case size thresholds, hindering the automation of the process and its application to inter-genotypic studies.

In contrast to the great multi-cultivar variability detected for the pollen grain size, we observed that pollen grains looked similarly in all cultivars after Alexander’s staining. In all cultivars, viable and non-viable pollen grains were easily visually differentiated (dark blue = viable/pale blue = non-viable), suggesting the possibility of automating their differential counting through the use of novel image-based technologies on color digital microscope images. Any color digital image is the result of the combination of three primary color-channels [red (R), green (G) and blue (B)], which imitate the way humans perceive color [44]. These three channels store complementary information, and their analysis might yield different biological information for subsequent decision support. In our case, R, G and B channels contain differential information on viable and non-viable pollen grains (which correspond to the differential R, G and B profiles of the stains included in Alexander’s modified staining) which can be used for its separation and further independent analysis. In this regard, the separation algorithm of the RGB images into independent colors channels is a critical step towards automatic image techniques, and so different separation processes have been proposed [44]. Here, we used the separation of the red, green and blue primary channels of the original RGB color digital microscope images to represent the intensity of the stains. Although color separation was not perfect (the color spectrum of both malachite green and acid fuchsin dyes spread over other color channels), we obtained enough selective contrast between these two dyes in the red and green layers to provide an appropriate difference between light blue- and dark purple-stained structures. In this regard, diverse processes (like stain deconvolution) have been proposed as alternative approaches to improve the separation of stains [45], but prior information is needed about the red, green and blue components of the pure stains to separate, which is not always feasible [46]. The simple separation approach used in our work and implemented in PollenCounter proved to be sufficient to make an appropriate counting of viable and total pollen grains in a diverse genetic framework, yielding highly accurate viability measurements and replacing labor-intensive and time-demanding manual processes.

PollenCounter can be easily run by users by opening the macro (Additional file 2) with Fiji and by choosing the source directory where the pollen digital images are stored. Then, the macro automatically analyses the whole set of images, releasing the number of pollen grains counted in the red and green color fractionated images in a pivot table for pollen viability calculation. Obviously, an adequate contrast between viable and non-viable pollen grains is needed when acquiring the RGB images, and settings should be optimized by the user to obtain the best differentiation of both type of pollen grains. The macro can be adapted to user needs by modifying different parameters related to pollen grain morphology, like their size and circularity. In fact, this flexibility opens its use to the evaluation of pollen viability in other crops with different pollen grain size and shape. As default settings, we set a surface area range of 60–800 pixels2 for the total counting of grapevine pollen grains (either viable or non-viable) in the red channel. In contrast, and aware of the general higher size of viable grains compared to non-viable, this range was slightly shortened to 100–800 pixels2 for the counting of viable pollen grains in the green channel, since no viable pollen grains were generally observed under 100 pixels2 and some small dark-stained particles (other than pollen grains) can appear as the result of the staining procedure. Thus, the presence of an adjacent small acid fuchsin-stained balloon-like structure next to some viable pollen grains was found to be common in some cultivars like Doña Blanca, Muscat Ottonel or Tinto Velasco (Additional file 7A). This structure could correspond to some protoplasm content released after intine rupture due to pollen grain hydration during the staining procedure, as previously suggested [47]. Although these structures are discarded during the processing of red-channel derived images due to their color features (Additional files 7B and 7D), they are detected in the green channel (Additional file 7C), so potentially counted as viable pollen grains if no area filter is included. The area filter included in our tool allows discarding these particles (Additional file 7E), decreasing the number of false positives and improving the accuracy of our tool.

Additionally, circularity factor can be modified in PollenCounter to adapt to pollen shape features. In general, we observed that non-viable pollen grains tend to be less spherical than viable grains (see Additional file 5 for an example). Nonetheless, such difference was not wide enough to establish it as an additional criterion for pollen grains differentiation, so we opted to maintain a fixed wide range of 0.40–1.00 to ensure a wide-spectrum counting of pollen grains.

Pollen viability shows a high uniformity between different inflorescences of the same clone

Nowadays, certified grapevine clones are asexually multiplied in nurseries by vegetative propagation to obtain genetically identical copies. During this process, the frequency and rate of spontaneous mutations driving to divergent phenotypes are suggested to be low [48], and so the vast majority of the new plant material obtained shows identical phenotypic features to the original plant. Accordingly, we observed similar pollen viability values for the different plants sampled for Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc genotypes, which respectively derived from a single clone. Besides, diverse works indicate that shoot and inflorescence position can influence grapevine reproductive performance. As some examples, basal clusters on Vidal blanc and Concord hybrids have been reported to be significantly heavier than second order clusters [49, 50], and similar differences have been found in Cabernet Sauvignon, Merlot and Sauvignon blanc V. vinifera cultivars [51, 52]. Phenology is another trait suggested to be affected by cluster and shoot positions, with lower inflorescences emerging before those in an upper position, and distal shoots flowering before proximal ones [52, 53]. Here, we aimed to evaluate the contribution of shoot and inflorescence position to pollen viability variability in four different genetic backgrounds. Our results indicate that pollen viability values of Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc are stable between different inflorescences, and we did not find any evidence of differences between basal and second clusters, nor between clusters from proximal and distal shoots. It is generally accepted that carbohydrate reserves are needed for pollen grain formation and development [13, 54], which could be compromised in some flowers due to the known preferential sugar transport to specific inflorescences during flowering time [52]. Unfortunately, carbohydrate reserves were not evaluated in our work, so it was not possible to establish a correlation between sugar level and pollen viability in inflorescences of different order/position. Anyway, it could be expected that carbohydrate reserves and/or supply during pollen formation were sufficient to ensure an adequate pollen development in these cases.

To our knowledge, the current work is the first one aimed to determine sources of variability within-plant (cluster-to-cluster differences on a shoot and shoot-to-shoot differences on a vine) in grapevine pollen viability through its direct quantification. The observed uniformity of pollen viability in the different inflorescences along the shoot and the cordon has relevant implications in sampling strategies. Our findings suggest that under optimal growing conditions, a random sampling of flowers from inflorescences from selected shoots is enough to obtain unbiased results to explore the genetic variability in pollen viability in diverse grapevine genotypes.

Tempranillo Tinto clones show high pollen viability variability

Somatic variation can occasionally affect major phenotypic traits, and somatic variants can be exploited for clonal selection programs for the improvement of elite cultivars [48]. As seen for other traits like grape aroma [55], berry skin color [56], cluster compactness [57] or water use efficiency [58], we observed a wide intra-cultivar genetic variability for pollen viability. In this regard, clones of Tempranillo Tinto with very high (RJ-51 and VP-2) and much lower (VP-11 and VP-25) pollen viability values were identified within the same plot (VP). These four Tempranillo Tinto clones were intentionally selected for their different viticulture performances [33]: whereas RJ-51 and VP-2 are two high-yielding clones with compact clusters and a high number of berries per cluster, VP-11 and VP-25 are two less productive and less compact clones with a minor number of berries per cluster. As previously suggested for a set of Cabernet Sauvignon clones [8], our results indicate a direct link between pollen viability and crop yield. In this regard, high pollen sterility (and the probably linked ovule sterility [15]) may lower fruit set rates and berry number per cluster, reducing cluster weight and consequently, crop yield [1, 5]. Similarly, the reduction of the berry number will reduce cluster compactness, which in turn might improve grape and wine quality via cluster microclimate improvement [59].

Specific comparative cytological and transcriptome analyses between sterile mutants and fertile wild-types have been useful to uncover the role of different genetic networks in the pollen development of the model plant A. thaliana [60, 61]. In grapevine, such analyses are hindered by the lack of mutant collections [62]. An alternative approach is the comparative analysis of clones differing in the trait of interest, like recently done for the in-depth analysis of cluster compactness [57]. In this regard, the identification and further characterization of the Tempranillo Tinto clones described here may be very useful to understand the genetic mechanisms involved in male sterility determination in grapevine.

Grapevine cultivars show large variability for pollen viability and different susceptibility to environment

The assessment of genetic diversity by exploring natural variation is the basis for the improvement of any crop. Over the last years, an increasing number of studies have dealt with the evaluation of grapevine diversity for many relevant traits, such as berry size [63], cluster structure [33] and water use efficiency [58]. The quality of the pollen produced by grapevine flowers is an important component of its reproductive behavior [5], but the information available is limited to a low number of genotypes [5, 7, 18, 64,65,66]. In this regard, our work aimed to explore the genetic diversity of this trait among a wide set of table- and wine-grape cultivars of diverse origin, showing for the first time the wide variability present in the cultivated grapevine.

In contrast to the uniformity observed for the pollen viability between plants in Cabernet Franc, Chardonnay, Gewuerztraminer and Sauvignon Blanc, we observed some cultivars with a high level of variation between the two plants (inflorescences) sampled in the ICVV plot. These plants were generated in 2009 from the duplication of an older grapevine collection grown in another plot [33]. During this process, different plants of the same cultivar were used to obtain the scions, which could explain part of the differences observed. Nevertheless, other additional genetic factors reducing pollen viability uniformity between inflorescences cannot be discarded.

As expected, most of the cultivars analyzed in this work showed high or very high pollen viability values, including very relevant wine grape varieties like Airén, Cabernet Sauvignon and Pinot Noir. During grapevine domestication and selection, humans focused on maintaining those genotypes capable to yield regular grape production, and so they unconsciously selected beneficial alleles for fertility, flower drop and productivity [67]. In this regard, self-fruitful hermaphroditic plants were preferentially selected over the dioecious wild form, since they allowed the obtaining of much more fruit on a predictable basis [68]. During this selection process, it is likely that grapevine genotypes with low pollen viability could have been progressively unwittingly discarded, maintaining those with better reproductive performance. Exceptionally, individuals with a low capability to generate viable pollen grains could have been selected for grape production despite this disadvantage. This can be the case of Corinto Bianco, a somatic variant of the seeded cultivar Pedro Ximenes [69] with very low pollen viability, but able to produce seedless berries without fertilization, which are highly appreciated for dried fruit consumption (raisins). Recently, an array of causes has been described to understand the genetic origin of parthenocarpy of Corinto Bianco [18], which include the absence of viable pollen. In contrast, we observed a low presence of viable pollen grains, which would explain the occasional formation of seeded berries after fertilization [18, 69]. Genetic and ploidy analysis of regenerated plantlets from Corinto Bianco seeds revealed their polyploidy (3n or 4n), probably derived by the fusion of unreduced gametes [18]. Interestingly, we observed an abnormal higher size for the viable pollen grains of this cultivar compared to non-viable pollen grains, especially in 2015 (Additional file 8). Considering the generally accepted correlation between pollen grain size and ploidy level [70], these stained and bigger pollen grains may correspond to viable diploid (2n) pollen grains of Corinto Bianco.

Pollen development and maturation is known to be sensitive to environmental conditions, and both abiotic (heat, cold, drought) and biotic factors may affect male reproductive performance [54]. In our work, 2015 and 2017 environmental conditions might have played specific effects in the different grapevine cultivars analyzed, identifying more and less-sensitive cultivars. Regarding the sensitive ones (those with a difference in APV (%) between seasons above 10%), we found higher pollen viability in 2017 than in 2015, except for the Spanish cultivar Verdil. Pre-flowering climate conditions were slightly different in La Rioja during these two seasons (2015 was warmer and drier than 2017, Additional file 9) and it could have been one of the factors affecting pollen development, as generally accepted for other crops [54, 71, 72]. In this light, previous reports indicate that the effect of temperature on pollen viability is genotype-specific [9, 27, 73], and very high temperatures (42 °C for 4 h) can jeopardize the regular formation of the outer wall of pollen grains in cv. Touriga Nacional, driving to an overall reduced viability [27]. Nevertheless, specific works are needed to understand how environmental factors individually affect grapevine pollen formation and/or maturation. Besides, and although many factors are involved, the sensitivity of these cultivars to temperature could explain part of the seasonal variance observed for some yield components for these genotypes, like previously reported for cv. Trousseau Noir [74].

Conclusions

Pollen viability assessment is critical for understanding grapevine reproductive development and for grapevine genetic improvement. Here, we introduce PollenCounter, a new processing tool capable to do an automatic counting of pollen grains after Alexander’s modified staining. Its high precision and accuracy make it suitable to estimate pollen viability in a diverse genetic framework, substituting time-demanding manual processes. The automatic analysis of pollen viability in a wide collection of cultivars revealed a great genetic variability for this trait, as well as differences in their susceptibility to the environment. A large variation was also identified among Tempranillo Tinto clones, but very small within single plants of the same clone of diverse cultivars. Some of the insights given in this work will guide future analyses aimed to unravel the physiological and genetic mechanisms affecting pollen viability in grapevine.

References

  1. Dry PR, Longbottom ML, McLoughlin S, Johnson TE, Collins C. Classification of reproductive performance of ten winegrape varieties. Aust J Grape Wine Res. 2010;16:47–55.

    Article  Google Scholar 

  2. Bessis R. La maîtrise des rendements. Revue Des Oenologues. 1993;19:7–10.

    Google Scholar 

  3. Collins C, Dry PR. Response of fruitset and other yield components to shoot topping and 2-chlorethyltrimethyl-ammonium chloride application. Aust J Grape Wine Res. 2009;15:256–67.

    Article  CAS  Google Scholar 

  4. Acimovic D, Tozzini L, Green A, Sivilotti P, Sabbatini P. Identification of a defoliation severity threshold for changing fruitset, bunch morphology and fruit composition in Pinot Noir. Aust J Grape Wine Res. 2016;22(3):399–408.

    Article  CAS  Google Scholar 

  5. Baby T, Gilliham M, Tyerman SD, Collins C. Differential fruitset between grapevine cultivars is related to differences in pollen viability and amine concentration in flowers. Aust J Grape Wine Res. 2016;22(1):149–58.

    Article  CAS  Google Scholar 

  6. Ewart A, Kliewer WM. Effects of controlled day and night temperatures and nitrogen on fruit-set, ovule fertility, and fruit composition of several wine grape cultivars. Am J Enol Vitic. 1977;28(2):88–95.

    CAS  Google Scholar 

  7. Kelen M, Demirtas I. Pollen viability, germination capability and pollen production level of some grape varieties (Vitis vinifera L.). Acta Physiol Plant. 2003;25(3):229–33.

    Article  Google Scholar 

  8. Bowed PA, Kliewer WM. Influence of clonal variation, pruning severity, and cane structure on yield component development in “Cabernet Sauvignon” grapevines. J Am Soc Hort Sci. 1990;115(4):530–4.

    Google Scholar 

  9. Ebadi A, May P, Sedgley M, Coombe BG. Effect of low temperature near flowering time on ovule development and pollen tube growth in the grapevine (Vitis vinifera L.), cvs Chardonnay and Shiraz. Aust J Grape Wine Res. 1995;1(1):11–8.

    Article  Google Scholar 

  10. Longbottom ML, Dry PR, Sedgley M. A research note on the occurrence of ‘star’ flowers in grapevines: observations during the 2003–2004 growing season. Aust J Grape Wine Res. 2004;10(3):199–202.

    Article  Google Scholar 

  11. Keller M. Botany and anatomy. In: Keller M, editor. The science of grapevines: anatomy and physiology. Amsterdam: Elsevier; 2010.

    Google Scholar 

  12. Cresti M, Ciampolini F. Ultrastructural characteristics of pollen development in Vitis vinifera L. (cv. Sangiovese). Vitis. 1999;38(4):141–4.

    Google Scholar 

  13. Lebon G, Duchene E, Brun O, Magne C, Clement C. Flower abscission and inflorescence carbohydrates in sensitive and non-sensitive cultivars of grapevine. Sex Plant Reprod. 2004;17(2):71–9.

    Article  CAS  Google Scholar 

  14. Lebon G, Duchene E, Brun O, Clement C. Phenology of flowering and starch accumulation in grape (Vitis vinifera L.) cuttings and vines. Ann Bot. 2005;95(6):943–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Iyer CPA, Randhawa GS. Induction of pollen sterility in grapes (Vitis vinifera). Vitis. 1966;5:433–45.

    Google Scholar 

  16. Shivanna KR, Linskens HF, Cresti M. Pollen viability and pollen vigor. Theor Appl Genet. 1991;81:34–42.

    Article  Google Scholar 

  17. Dafni A, Firmage D. Pollen viability and longevity: practical, ecological and evolutionary implications. Plant Syst Evol. 2000;222:113–32.

    Article  Google Scholar 

  18. Royo C, Carbonell-Bejerano P, Torres-Pérez R, Nebish A, Martinez O, Rey M, Aroutiounian R, Ibanez J, Martinez-Zapater JM. Developmental, transcriptome, and genetic alterations associated with parthenocarpy in the grapevine seedless somatic variant Corinto bianco. J Exp Biol. 2016;67(1):259–73.

    CAS  Google Scholar 

  19. Alexander MP. Differential staining of aborted and nonaborted pollen. Stain Technol. 1969;44(3):117–22.

    Article  CAS  PubMed  Google Scholar 

  20. Atlagic J, Terzic S, Marjanovic-Jeromela A. Staining and fluorescent microscopy methods for pollen viability determination in sunflower and other plant species. Ind Crops Prod. 2012;35:88–91.

    Article  Google Scholar 

  21. Kelly JK, Rasch A, Kalisz S. A method to estimate pollen viability from pollen size variation. Am J Bot. 2002;89(6):1021–3.

    Article  PubMed  Google Scholar 

  22. Dufay M, Vaudey V, De CauwerI, Touzet P, Cuguen J, Arnaud J-F. Variation in pollen production and pollen viability in natural populations of gynodioecious Beta vulgaris ssp. maritima: evidence for a cost of restoration of male function? J Evol Biol. 2008;21(1):202–12.

    Article  CAS  PubMed  Google Scholar 

  23. Walter A, Liebisch F, Hund A. Plant phenotyping: from bean weighing to image analysis. Plant Methods. 2015;11:14.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Lamprecht MR, Sabatini DM, Carpenter AE. Cell Profiler™: free, versatile software for automated biological image analysis. Biotechniques. 2007;42(1):71–5.

    Article  CAS  PubMed  Google Scholar 

  25. Legland D, El-Hage F, Méchin V, Reymond M. Histological quantification of maize stem sections from FASGA-stained images. Plant Methods. 2017;13:84.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Volk G. Collecting pollen for genetic resources conservation. In: Guarino L, Ramanatha VR, Goldberg E, editors. Collecting plant genetic diversity: technical guidelines 2011 update. Rome: CAB International; 2011.

    Google Scholar 

  27. Pereira HS, Delgado M, Avó AP, Barao A, Serrano I, Viegas W. Pollen grain development is highly sensitive to temperature stress in Vitis vinifera. Aust J Grape Wine Res. 2014;20(3):474–84.

    Article  CAS  Google Scholar 

  28. Costa CM, Yang S. Counting pollen grains using readily available, free image processing and analysis software. Ann Bot. 2009;104:1005–10.

    Article  PubMed  PubMed Central  Google Scholar 

  29. France I, Duller AWG, Duller GAT, Lamb HF. A new approach to automated pollen analysis. Quat Sci Rev. 2000;19:537–46.

    Article  Google Scholar 

  30. Holt K, Allen G, Hodgson R, Marsland S, Flenley J. Progress towards an automated trainable pollen location and classifier system for use in the palynology laboratory. Rev Palaeobot Palynol. 2011;167:175–83.

    Article  Google Scholar 

  31. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, et al. Fiji: an open-source platform for biological-image analysis‎. Nat Methods. 2012;9(7):676–82.

    Article  CAS  PubMed  Google Scholar 

  32. Peterson R, Slovin JP, Chen C. A simplified method for differential staining of aborted and non-aborted pollen grains. Int J Plant Biol. 2010;1:e13.

    Article  Google Scholar 

  33. Tello J, Aguirrezábal R, Hernaiz S, Larreina B, Montemayor MI, Vaquero E, Ibáñez J. Multicultivar and multivariate study of the natural variation for grapevine bunch compactness. Aust J Grape Wine Res. 2015;21(2):277–89.

    Article  Google Scholar 

  34. Coombe BG. Adoption of a system for identifying grapevine growth stages. Aust J Grape Wine Res. 1995;1(2):104–10.

    Article  Google Scholar 

  35. Tello J, Cubero S, Blasco J, Tardaguila J, Aleixos N, Ibáñez J. Application of 2D and 3D image technologies to characterize morphological attributes of grapevine clusters. J Sci Food Agric. 2016;96(13):4575–83.

    Article  CAS  PubMed  Google Scholar 

  36. Klodt M, Herzog K, Töpfer R, Cremers D. Field phenotyping of grapevine growth using dense stereo reconstruction. BMC Bioinform. 2015;16:143.

    Article  Google Scholar 

  37. Diago MP, Tardaguila J, Aleixos N, Millán B, Prats-Montalbán JM, Cubero S, Blasco J. Assessment of cluster yield components by image analysis. J Sci Food Agric. 2015;95(6):1274–82.

    Article  CAS  PubMed  Google Scholar 

  38. Herrero-Huerta M, González-Aguilera D, Rodríguez-Gonzalvez P, Hernández-López D. Vineyard yield estimation by automatic 3D bunch modelling in field conditions. Comput Electron Agric. 2015;110:17–26.

    Article  Google Scholar 

  39. Cubero S, Diago MP, Blasco J, Tardaguila J, Prats-Montalbán JM, Ibáñez J, Tello J, Aleixos N. A new method for assessment of bunch compactness using automated image analysis. Aust J Grape Wine Res. 2015;21(1):101–9.

    Article  Google Scholar 

  40. Portalés C, Ribes-Gómez E. An image-based system to preliminary assess the quality of grape harvest batches on arrival at the winery. Comput Ind. 2015;68:105–15.

    Article  Google Scholar 

  41. Hill GN, Evans KJ, Beresford RM, Dambergs RG. Comparison of methods for the quantification of botrytis bunch rot in white wine grapes. Aust J Grape Wine Res. 2014;20:432–41.

    Article  CAS  Google Scholar 

  42. Oerke E-C, Herzog K, Töpfer R. Hyperspectral phenotyping of the reaction of grapevine genotypes to Plasmopara viticola. J Exp Biol. 2016;67(18):5529–43.

    CAS  Google Scholar 

  43. Knauer U, Matros A, Petrovic T, Zanker T, Scott ES, Seiffert U. Improved classification accuracy of powdery mildew infection levels of wine grapes by spatial-spectral analysis of hyperspectral images. Plant Methods. 2017;13:47.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Cheng HD, Jiang XH, Sun Y, Wang J. Color image segmentation: advances and prospects. Pattern Recognit. 2001;34:2259–81.

    Article  Google Scholar 

  45. Alsubiae N, Trahearn N, Raza SEA, Snead D, Rajpoot NM. Stain deconvolution using statistical analysis of multi-resolution stain colour representation. PLoS ONE. 2017;12(1):e0169875.

    Article  Google Scholar 

  46. Tadrous PJ. Digital stain separation for histological images. J Microsc. 2010;240(2):164–72.

    Article  CAS  PubMed  Google Scholar 

  47. Chichiriccò G, Pacini E. Cupressus arizonica pollen wall zonation and in vitro hydration. Pl Syst Evol. 2008;270:231–42.

    Article  Google Scholar 

  48. Pelsy F. Molecular and cellular mechanisms of diversity within grapevine varieties. Heredity. 2010;104(4):331–40.

    Article  CAS  PubMed  Google Scholar 

  49. Wolpert JA, Howell GS, Cress CE. Sampling strategies for estimates of cluster weight, soluble solids and acidity of ‘Concord’ grapes. J Am Soc Hort Sci. 1981;105:434–8.

    Google Scholar 

  50. Wolpert JA, Howell GS, Mansfield TK. Sampling Vidal blanc grapes. I. Effects of training system, pruning severity, shoot exposure, shoot origin, and cluster thinning on cluster weight and fruit quality. Am J Enol Vitic. 1983;34(2):72–6.

    Google Scholar 

  51. Tarter ME, Keuter SE. Shoot-based sampling of Vitis vinifera cultivars. Am J Enol Vitic. 2008;59(1):55–60.

    Google Scholar 

  52. Trought MCT, Naylor AP, Frampton C. Effect of row orientation, trellis type, shoot and bunch position on the variability of Sauvignon Blanc (Vitis vinifera L.) juice composition. Aust J Grape Wine Res. 2017;23:240–50.

    Article  CAS  Google Scholar 

  53. Eltom M, Trought MCT, Agnew R, Parker A, Winefield CS. Pre-budburst temperature influences the inner and outer arm morphology, phenology, flower number, fruitset, TSS accumulation and variability of Vitis vinifera L. Sauvignon Blanc bunches. Aust J Grape Wine Res. 2017;23:280–6.

    Article  CAS  Google Scholar 

  54. De Storme N, Geelen D. The impact of environmental stress on male reproductive development in plants: biological processes and molecular mechanisms. Plant Cell Environ. 2014;37:1–18.

    Article  PubMed  Google Scholar 

  55. Duchene E, Legras JL, Karst F, Merdinoglu D, Claudel P, Jaegli N, Pelsy F. Variation of linalool and geraniol content within two pairs of aromatic and non-aromatic grapevine clones. Aust J Grape Wine Res. 2009;15(2):120–30.

    Article  CAS  Google Scholar 

  56. Kobayashi S, Yamamoto NG, Hirochika H. Association of VvmybA1 gene expression with anthocyanin production in grape (Vitis vinifera) skin–color mutants. J Jpn Soc Hort Sci. 2005;74(3):196–203.

    Article  CAS  Google Scholar 

  57. Grimplet J, Tello J, Laguna N, Ibáñez J. Differences in flower transcriptome between grapevine clones are related to their cluster compactness, fruitfulness and berry size. Front Plant Sci. 2017;8:632.

    Article  PubMed  PubMed Central  Google Scholar 

  58. Tortosa I, Escalona JM, Bota J, Tomás M, Hernández E, García-Escudero E, Medrano H. Exploring the genetic variability in water use efficiency: evaluation of inter and intra cultivar genetic diversity in grapevines. Plant Sci. 2016;251:35–43.

    Article  CAS  PubMed  Google Scholar 

  59. Tello J, Ibáñez J. What do we know about grapevine bunch compactness? A state-of-the-art review. Aust J Grape Wine Res. 2017. https://0-doi-org.brum.beds.ac.uk/10.1111/ajgw.12310.

    Google Scholar 

  60. Jakobsen MK, Poulsen LR, Schulz A, Fleurat-Lessard P, Moller A, Husted S, Schiott M, Amtmann A, Palmgren MG. Pollen development and fertilization in Arabidopsis is dependent on the MALE GAMETOGENESIS IMPAIRED ANTHERS gene encoding a Type V P-type ATPase. Genes Dev. 2005;19:2757–69.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  61. Feng B, Lu D, Ma X, Peng Y, Sun Y, Ning G, Ma H. Regulation of the Arabidopsis anther transcriptome by DYT1 for pollen development. Plant J. 2012;72:612–24.

    Article  CAS  PubMed  Google Scholar 

  62. Santos-Rosa M, Poutaraus A, Merdinoglu D, Mestre P. Development of a transient expression system in grapevine via agro-infiltration. Plant Cell Rep. 2008;27(6):1053–63.

    Article  CAS  PubMed  Google Scholar 

  63. Houel C, Martin-Magniette ML, Nicolas SD, Lacombe T, Le Cunff L, Franck D, Torregrosa L, Conejero G, Lalet S, This P, et al. Genetic variability of berry size in the grapevine (Vitis vinifera L.). Aust J Grape Wine Res. 2013;19(2):208–20.

    Article  Google Scholar 

  64. Abreu I, Costa I, Oliveira M, Cunha M, de Castro R. Ultrastructure and germination of Vitis vinifera cv. Loureiro pollen. Protoplasma. 2006;228(1–3):131–5.

    Article  CAS  PubMed  Google Scholar 

  65. Carbonell-Bejerano P, Royo C, Torres-Pérez R, Grimplet J, Fernandez L, Franco-Zorrilla JM, Lijavetzky D, Baroja E, Martinez J, García-Escudero E, et al. Catastrophic unbalanced genome rearrangements cause somatic loss of berry color in grapevine. Plant Physiol. 2017;175(2):786–802.

    PubMed  Google Scholar 

  66. Liu SM, Sykes SR, Clingeleffer PR. Pollen fertility and berry setting behaviour of the grape variety Carina. Aust J Exp Agric. 2007;47(7):877–82.

    Article  Google Scholar 

  67. Bacilieri R, Lacombe T, Le Cunff L, Di Vecchi-Staraz M, Laucou V, Genna B, Péros JP, This P, Boursiquot J-M. Genetic structure in cultivated grapevine is linked to geography and human selection. BMC Plant Biol. 2013;13:25.

    Article  PubMed  PubMed Central  Google Scholar 

  68. This P, Lacombe T, Thomas MR. Historical origins and genetic diversity of wine grapes. Trends Genet. 2006;22(9):511–9.

    Article  CAS  PubMed  Google Scholar 

  69. Vargas AM, Velez MD, de Andres MT, Laucou V, Lacombe T, Boursiquot JM, Borrego J, Ibáñez J. Corinto Bianco: a seedless mutant of Pedro Ximenes. Am J Enol Vitic. 2007;58(4):540–3.

    Google Scholar 

  70. De Storme N, Zamariola L, Mau M, Sharbel TF, Geelen D. Volume-based pollen size analysis: an advanced method to assess somatic and gametophytic ploidy in flowering plants. Plant Reprod. 2013;26(2):65–84.

    Article  CAS  PubMed  Google Scholar 

  71. Sheoran IS, Saini HS. Drought-indiced male sterility in rice: changes in carbohydrate levels and enzyme activities associated with the inhibition of starch accumulation in pollen. Sex Plant Reprod. 1996;9(3):161–9.

    Article  Google Scholar 

  72. Zinn KE, Tunc-Ozdemir M, Harper JF. Temperature stress and plant sexual reproduction: uncovering the weakest links. J Exp Biol. 2010;61(7):1959–68.

    CAS  Google Scholar 

  73. Malaescu M, Blidariu C, Dobrei A, Dobrei A, Nistor E, Dragunescu A, Olaru D, Velicevici G, Dobromir D. Research concerning the effect of temperature on pollen germination capacity from several grape cultivars grown in Timisoara Didactic Station area. J Hortic For Biotechnol. 2015;19(4):25–35.

    Google Scholar 

  74. Shellie K, Cragin J, Serpe M. Performance of alternative European wine grape cultivars in Southwestern Idaho: cold hardiness, berry maturity, and yield. Hortechnology. 2014;24(1):138–47.

    Google Scholar 

Download references

Authors’ contributions

Conceived and designed this work: JI and JT. Designed phenotyping tool: JI and JT. Conducted the experiments: JT and MIM. Analyzed the data: JT. Drafted the manuscript: JT. Contributed reagents/materials: JI and AF. Critically revised the manuscript: JT, AF and JI. All authors read and approved the final manuscript.

Acknowledgements

We sincerely acknowledge Viveros Provedo for providing part of the plant material used in this work, Dr. P. Carbonell-Bejerano for Tempranillo blanco (FZ) pollen images, and the Servicio de Investigación Vitivinícola (SIVV, Gobierno de La Rioja), J.L. Pérez Sotés, A. Nebish, K. Frank, V. Dockner and G. Wiedeschitz for their technical assistance and/or plant maintenance.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

Consent for publication

The authors consent to the publication of the manuscript.

Ethics approval and consent to participate

All authors have read and accept the content of this manuscript. Field studies were conducted in accordance with local legislations.

Funding

This work was co-funded by the Spanish Ministerio de Economía y Competitividad (MINECO) through Project AGL2014-59171-R (co-funded by FEDER). Open access funding provided by University of Natural Resources and Life Sciences Vienna (BOKU).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Javier Tello.

Additional files

Additional file 1.

List of the 120 grapevine accessions selected for this work. Accession code, cultivar name, and season sampled are indicated.

Additional file 2.

PollenCounter Fiji-compatible macro.

Additional file 3.

Pollen images.

Additional file 4.

Analysis of two different automatic methods to estimate pollen viability. In A, histograms are shown for non-viable (NV, blue) and viable (V, green) pollen grains area (pixels2). Pollen grains (n = 6082) contained in 38 images from 19 different cultivars are considered. Red broken line indicates the average value of the smallest viable pollen grain detected in each image (118 pixels2). In B, pollen size index (PSIimg) is compared to the manual pollen viability (MPVimg) assessment. In B, the manual approach is compared to the automatic system designed in this work (APVimg). For B and C, 38 images are considered.

Additional file 5.

Separation of the red (R), green (G) and blue (B) layers of a RGB image of six grapevine pollen grains (three viable, three non-viable). To obtain red, green and blue color fractionated grayscale pictures, the original RGB image was processed using the “split channels” tool in Fiji. Pollen grains correspond to cv. Chardonnay.

Additional file 6.

Separation of RGB images of the three dyes used in the Alexander’s modified staining solution (AF: Acid fuchsin; MG: Malachite green; OG: Orange G) on their red (R), green (G) and blue (B) basic layers. To obtain R, G and B color fractionated grayscale pictures, RGB images of pure dyes were processed using the “split channels” tool in Fiji.

Additional file 7.

Original (A), red- and green-color fractionated grayscale pictures (B, C), and PollenCounter output images (D, E) obtained for a pollen sample of cv. Muscat Ottonel. An example of protoplasm content released from a viable pollen grain is indicated with a red arrow. Note that it is not detected in the red-channel derived images (B, D). Although detected in the green channel (C), the structure is not considered as a valid pollen grain to be counted (black region, E). For simplification, only a representative area of the processed image is shown. B and C images were obtained using the “split channels” tool in Fiji. The number of pollen grains automatically counted by PollenCounter (in cyan) is indicated in the lower right corners (D, E).

Additional file 8.

Corinto Bianco pollen grains after Alexander’s modified staining obtained in 2015 (A) and 2017 (B). Red arrows indicate viable pollen grains. For simplification, only representative areas are shown.

Additional file 9.

Pre- and post-flowering climate conditions in the Grapevine Germplasm Collection of the Instituto de Ciencias de la Vid y del Vino in 2015 and 2017. Days are shown according to the date when the first sample was collected in 2015 (2-June) and 2017 (23-May), indicated as a “0”. Sample collection extended for 14 and 10 days in 2015 and 2017, respectively. Light and dark red lines indicate mean temperatures for 2015 and 2017, respectively. Light and dark blue columns indicate accumulate daily rainfalls for 2015 and 2017, respectively. Data were obtained from La Rioja Government website (http://www.larioja.org/siar).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Tello, J., Montemayor, M.I., Forneck, A. et al. A new image-based tool for the high throughput phenotyping of pollen viability: evaluation of inter- and intra-cultivar diversity in grapevine. Plant Methods 14, 3 (2018). https://0-doi-org.brum.beds.ac.uk/10.1186/s13007-017-0267-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13007-017-0267-2

Keywords