Next Article in Journal
Primary Root and Mesocotyl Elongation in Maize Seedlings: Two Organs with Antagonistic Growth below the Soil Surface
Previous Article in Journal
Validation of Reference Genes for Studying Different Abiotic Stresses in Oat (Avena sativa L.) by RT-qPCR
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Variation in Pollination Deficits in an Insect-Pollinated Dioecious Crop

Centre for Functional Ecology, Department of Life Sciences, University of Coimbra, Calçada Martim de Freitas, 3000-456 Coimbra, Portugal
*
Author to whom correspondence should be addressed.
Plants 2021, 10(7), 1273; https://doi.org/10.3390/plants10071273
Submission received: 27 May 2021 / Revised: 9 June 2021 / Accepted: 18 June 2021 / Published: 22 June 2021

Abstract

:
Inadequate quantity and quality of pollen reaching the stigmas decreases the sexual reproductive output of plants, compromising yield. Still, the current extent of pollen limitation affecting yield (i.e., pollination deficits) is poorly quantified. This study is aimed at quantifying pollination deficits in kiwifruit orchards, a dioecious plant with a fruit caliber and market value largely dependent on pollination services. For that, we set up a pollination experiment and quantified services and yield provided by current pollination vectors, and under optimal pollination, over two years in a total of twenty-three orchards covering the kiwifruit production range in Portugal. We characterized nine fruit traits and used: (1) fruit weight to calculate pollination deficits and relate them with pollinator diversity and abundance, and environmental variables; and (2) production values, fruit caliber, and market values to calculate economic impact of pollination deficits. Results showed that pollination deficits were variable in time and space and were significantly and negatively correlated with pollinator abundance, while the opposite pattern was obtained for production, supporting the notion that a higher pollinator’s abundance is related to lower pollination deficits and higher yields. Understanding the factors affecting pollination deficits is crucial to depict the need for nature-based solutions promoting pollinators and to resort to management practices assisting pollination.

1. Introduction

The yield and quality of over 75% of crops worldwide is directly affected by animal pollination [1], and the area occupied by pollinator dependent crops has increased over the last decades [2]. Pollination is therefore an important biodiversity-dependent service supporting food provisioning. Pollination is the simple transfer of pollen from the anthers to the stigmas, culminating in fertilization. However, it is far from being simple, as plants rely on mutualistic interactions with animals to carry their pollen, frequently suffering pollen limitation [3]. Pollination deficit can be defined as the inadequate quantity and quality of pollen reaching the stigmas, which decreases the sexual reproductive output of plants [4]. It can result from factors such as insufficient or inefficient pollen transfer due to limited pollinator abundance and diversity, low pollinator activity or efficiency, and/or flowering asynchrony [4]. In this context, pollinators’ decline as a result of major global changes such as land use changes (e.g., fragmentation and agriculture intensification), pesticide use, biological invasions, and eutrophication constitutes a threat to pollination services supporting food production [4,5]. Pollination deficits in crops have been assessed indirectly though the quantification of pollinator populations [4,6,7], but direct estimates of pollination deficits and their impact on crop productivity are scarcer (see, however, [8,9,10]). Because pollination deficits are the direct result of pollinator activity and available pollinator populations, quantifying pollinator diversity and abundance and relating it with crop yield is crucial to understanding the observed pollination deficits. Additionally, the direct estimates of pollination deficits measure the effect of limited pollination services on reproductive output from a plant perspective (fitness) and, if linked with productivity, it may also consider the farmer’s perspective, which is more focused on agronomic and economic yields [6].
Kiwifruit (Actinidia spp., Actinidiaceae) is a dioecious crop mainly pollinated by insects, with buzz pollinators being pointed out as more efficient [11] and, to a lesser extent, by wind [12,13]. Kiwifruit caliber and market value are largely influenced by the number of seeds formed, which in turn depend on the number and quality of pollen grains reaching the stigmas [14]. Therefore, efficient pollination is a key aspect in kiwifruit production and its economic viability. High quality market fruits (>100 g) have over 1000 seeds, and to ensure good pollination each flower should receive an estimated number of 2000–3000 viable pollen grains [12,15,16]. In addition, kiwifruit pollination is highly dependent on weather conditions during flowering time, which often limit insect pollination (e.g., [12,17,18]). Finally, the short flowering period, the lack of synchrony between male and female plants, and the low number and/or inadequate distribution of males within the orchard further contribute to an inefficient pollination [19]. While wind pollination can contribute to pollination, enabling the attainment of a high fruit set, its effectiveness is considered highly inconsistent and insufficient [12,13]. As a result of these hurdles, kiwifruit producers frequently include management practices such as the installation of honeybee colonies or the application of artificial pollen [16,20]. However, these practices have two drawbacks. First, honeybees are not the most efficient pollinator and are more sensitive to poor weather conditions when compared to wild pollinators [11,21]. Second, purchased pollen has a high cost, and contributes to the dissemination of diseases, such as the Pseudomonas syringae pv. actinidiae, during pollen application [20,22].
In key kiwifruit production areas worldwide (e.g., New Zealand and Italy) there are well-established practices and guidelines for producers (e.g., [23,24,25]). In these areas, artificial pollen application is a common and cost-effective practice, as pollination resources are scarce due to the vast plantation areas [25]. However, the reality is different in countries with smaller but still significant contribution to the kiwifruit worldwide market, such as Portugal and Spain. In Portugal, kiwifruit was introduced in 1973, the production started in the 1980’s [26], and the production area has grown steadily over the last decades (currently representing around 2736 ha). This region is characterized by small-scale orchards in a mosaic landscape [27]. Landscapes composed of different mosaics of land-uses may provide multiple resources for pollinator populations and promote wild pollinator communities [28,29]. Consequently, pollination requirements may vary significantly from those in massive plantations, setting different challenges for management practices. However, and despite the current knowledge of the limitations to an efficient pollination, the magnitude of pollination deficits in kiwifruit orchards (and in many other crops), its relation with available pollinator communities, and subsequent economic impact are completely unknown. Such information is important to assist decision making and optimize crop production.
Our study aims to quantify pollination services in kiwifruit orchards representing the entire production area of Portugal in order to provide real estimates on pollination deficits, relate them with available pollinator communities, and assess potential impacts on kiwifruit production and market value. We hypothesize that: (1) wind has a significant contribution to fruit set, but it is insufficient to attain high quality market fruits; (2) if current pollination services are inadequate, fruit set and fruit caliber from supplemented flowers will be higher than those from open pollinated flowers, and this will be reflected on overall orchard production and monetary gain; and (3) pollination services, and hence the levels of pollination deficits, are influenced by pollinator diversity and abundance and environmental variables which may explain temporal and spatial variability. To test our hypotheses, we set up a pollination experiment over two years in a total of twenty-three orchards and quantified nine fruit traits and pollination deficits. Environmental and pollinator variables were used to explain pollination deficits and production. Finally, production values, fruit caliber, and market values were used to calculate economic impact of pollination deficits.

2. Results

2.1. Contribution of Wind to Kiwifruit Pollination

Wind pollination produced variable results in fruit set with overall non-significant differences between open and wind pollination and a significantly lower fruit set when compared to supplementary pollination (Table 1). Within an orchard, wind pollination resulted in significantly lower fruit set when compared to supplementary pollination or to both open and supplementary pollination in three of the orchards, while no differences were detected in the remaining orchards and varieties (Table S1).
Overall, wind pollination produced small, nonmarketable fruits with weights significantly lower when compared to open and supplementary pollination (Table 1; Figure 1; Table S2). A significant effect of wind pollination was also obtained on dry matter, flesh firmness, acidity, S:L, and TSS (Table 1). Wind pollinated fruits showed, on average, lower dry matter content when compared to supplementary pollination, lower acidity and TSS when compared to open pollination, lower flesh firmness when compared to both open and supplementary pollination, and higher S:L when compared to supplementary pollination (Table 1). Within orchards, non-significant differences among treatments were observed for most fruit variables (p > 0.05; Table S2).

2.2. Pollination Deficits

Fruit set differed significantly between pollination treatments within year, while no differences were observed between years (Table 2). Fruit set was higher in supplemented flowers compared to open pollinated ones (Table 2). Within orchards and variety, fruit set of open pollinated flowers was significantly higher than that of supplementary pollinated flowers in five orchards (one in 2018 and four in 2019; Table S1).
Fruit weight varied significantly between pollination treatments within year, while no differences were observed between years (Table 2). Overall, fruits resulting from supplementary pollination were heavier than those resulting from open pollination (Table 2). Within orchards and varieties, a similar pattern was observed. Supplemented flowers produced heavier fruits than the ones produced by open pollinated flowers, although this was significant only in eight orchards/varieties (two in 2018 and six in 2019; Tables S2 and S3). For the remaining fruit variables, dry matter, flesh firmness, and S:L differed significantly between years, and none of these fruit variables showed variation between pollination treatment within a year (Table 2). Within the orchards, we found significant differences only for flesh firmness in two orchards and for S:L in one orchard in 2019 (Tables S2 and S3).
Overall, there was a significant decrease in pollination deficits in 2019 when compared to 2018 (Table 2, Figure 2). Within year, we found significant pollination deficits in six orchards out of nine in 2018 (Figure 3A). In 2019, we found pollination deficits in ten orchards out of twenty-two (Figure 3B). In both years, pollination deficits varied among orchards, from −0.03 to 0.43 in 2018, and from −0.05 to 0.27 in 2019.

2.3. Linking Pollination Deficits and Production with Environmental and Pollinator’s Related Variables

Regression analysis showed a significant negative relation between honeybee abundance and pollination deficit (Table 3), indicating that higher abundances of honeybees were associated with lower pollination deficit. Additionally, a significant and positive relation was obtained between the abundance of honeybees and Bombus species and productivity, indicating that higher abundances of these pollinators were associated with higher productivity levels. In contrast, a significant negative relation was found between pollinator richness and productivity, indicating an opposite pattern (Table 3). The two environmental variables included in the regression analyses showed non-significant relations with both pollination deficit and production (Table 3).

2.4. Linking Pollination Deficit with Production and Market Value

Data show high variability in fruit weight (Tables S2 and S3) within orchards. Consequently, fruits were distributed over several caliber categories (Figure S1). A significant effect of pollination treatment within years on the proportion of small and large fruits was observed, while no differences were observed between years (Table 2). Supplementary pollination resulted in a lower proportion of small fruits and a higher proportion of large fruits in both years when compared to open pollination (Table 2). When analyzed individually, orchards that showed pollination deficit also showed clear and similar patterns as the ones described above for the percentage of fruit per caliber category (Table S3). Overall, supplementary pollinated flowers produced a significantly (p < 0.05) or marginally significant (p < 0.10) lower proportion of small fruits compared to open pollination in three out of nine orchards in 2018, and in six out of twenty-two in 2019 (Tables S2 and S3). Supplementary pollination produced significantly (p < 0.05) or marginally significant (p < 0.10) higher proportion of large fruits compared to open pollination in four orchards out of nine in 2018, in five orchards out of twenty-two in 2019 (Tables S2 and S3).
Overall, kiwifruit production was not significantly affected by year, although there was a tendency toward an increase in production in 2019 when compared to 2018 (Table 2). Significant differences were observed between pollination treatments within years (Table 2), with supplementary pollination contributing to higher production values than open pollination in 2018 and 2019 (Table 2). Data further support that production can be increased by adequate pollen supply (Figure 4), particularly in orchards showing pollination deficit as was the case of orchard L in 2018, where a clear increase in production resulting from supplementary pollination can be observed (Figure 4).
Our calculations revealed that differences in the percentage of small and large fruits, resulting from adequate pollen supply (Tables S2 and S3), can translate into monetary gains (Figure S2). The monetary gain was not significantly affected by year, although there is a tendency toward an increase in monetary gain in 2019 when compared to 2018 (Table 2). Significant differences were observed between pollination treatments within years (Table 2), with an increase in the amount of euros under supplementary when compared to open pollination. When analyzed individually, the pattern is clear for several orchards in both years (Figure S2).

3. Discussion

Like many worldwide crops, efficient pollination is a key aspect in kiwifruit production and its economic viability [1,16,18]. In particular, kiwifruit is a dioecious plant that needs high pollen loads (which are transported by biotic and abiotic vectors) to attain marketable fruits [14,16]. This study provides the first spatio-temporal quantification of pollination deficits in kiwifruit orchards across a wide production area. The study region is characterized by productive systems with features that are very different from those of the high kiwifruit yield systems in leading countries, where most studies were performed and where most methodologies have been optimized (e.g., [23,24,25]; but see [18,19] in Spain). Overall, our results show an inefficient pollination by wind and variation in both space and time in pollination services, although significant levels of pollination deficits were observed only in some orchards and/or years. The increased abundance of main pollinator groups (i.e., honeybee and Bombus) were significantly related to a decreasing pollination deficit and increasing production levels. Moreover, pollen supply significantly improved the production in several orchards and, consequently, may result in increased economic revenue to kiwifruit producers. The study, while being mostly descriptive and focused on kiwifruit, is of global interest, provides information that is applicable to other crops, and opens way for further works exploring orchard and landscape factors affecting pollination services to crops. Below, we discuss the results in detail, linking with orchard management practices.

3.1. Contribution of Wind to Kiwifruit Pollination

Kiwifruit flowers have two main pollination syndromes: (1) wind pollination marked by pendulous flowers with large and fleshy stigmas, production of high quantities of pollen, and synchronous mass flower production lasting a short period of time [24,30]; and (2) insect-pollination, specifically buzz pollination, marked by the production of attractive flowers, floral scent, high amounts of pollen as reward, gradually maturing anthers, and a large number of ovules [11,30,31,32]. Although representing one of the two pollination syndromes in kiwifruit, our results show that wind pollination is not sufficient to guarantee adequate pollination for commercial production. The results are consistent with other studies, which reported that, while wind contributes to kiwifruit pollination, it is highly inconsistent and insufficient to attain high quality market fruits (e.g., [12,13,18]). Kiwifruit pollen is dry and carried in the air, and the use of large fans to blow pollen from male flowers and maximize wind pollination is a common practice used by producers in leading kiwifruit producing countries such as Italy [24]. In Portugal, this practice is not common and fans were used in five orchards only. In one of the orchards that used fans, wind pollination produced only 22.7% of fruits in the large and medium fruit categories against 83.9% resulting from open pollination, suggesting that this management practice is still insufficient to attain successful pollination levels. Altogether, the results obtained here highlight the importance of entomophilous pollination to attain marketable kiwifruits.

3.2. Spatio-Temporal Variation of Pollination Deficit

Sampling kiwifruit orchards in two consecutive years and doing so over a wide area of the production range allowed us to observe the existence of both spatial and temporal variation in pollination deficits. Pollination deficits can result from the lack of pollen and/or inefficient pollen transport, both of which can be affected by the available pollinator communities, management factors, and/or weather conditions [18,19,33,34]. While precipitation and temperature did not have a significant effect, the abundance of the main pollinator’s groups was negatively correlated with pollination deficits and positively correlated with orchard production, supporting that variation in pollinator communities is one of the factors driving productivity. In particular, increased abundances of honeybees led to lower pollination deficits, while both honeybee and Bombus species abundances contributed to higher productivity in kiwifruit orchards. Honeybees and Bombus species are indeed the main pollinators in kiwifruit orchards, with honeybees usually attaining high densities which may compensate for their lower kiwifruit pollination efficiencies in comparison to Bombus species [11,18]. Miñarro and Twizell [18] observed that honeybees were more abundant and visited more kiwifruit flowers per unit of time, but Bombus species were more efficient on a per visit basis, which resulted in an overall non-significant effect of pollinator species on fruit weight. Contrary to other studies that show increased pollination services with increased pollinator richness (e.g., [35]), in our case pollinator richness was negatively related to kiwifruit production. This may be explained by: (1) the efficiency of individual pollinator species, as demonstrated by MacInnis and Forrest [36] for strawberries, who found that bee identity had the strongest effect on yield; and (2) by the low abundance of other pollinators [27], which results in these having lower contribution to kiwifruit pollination and production regardless of their contribution to overall pollinator’s richness. Thus, in the case of kiwifruit, management practices and landscape contexts that promote efficient pollinator groups will likely promote pollination effectiveness and productivity.
Interestingly, for a significant number of orchards in both years, we did not find significant differences in fruit weight between open and supplementary pollination, indicating that the available pollinator communities were enough to ensure adequate pollination for commercial production. When the orchards have a good provision of pollination services (either through the wild pollinator communities alone or by its complementation with Apis mellifera colonies) and of pollen (through an adequate male to female ratio, proper distribution of males and flowering synchronization), natural pollination levels can be sufficient to produce high quality marketable fruits. Several studies have shown that pollination services are affected by the interplay between landscape composition and heterogeneity and farm management, and that landscape simplification and intensive farm management have a synergistically negative effect on pollination services [37,38,39,40,41]. In Portugal, kiwifruit orchards are of a small size, with low pesticide input and low insect non-friendly management practices, and are imbedded in a diversified landscape [27]. These factors may contribute to the maintenance of wild pollinator communities in the studied orchards, contributing to the successful pollination of kiwifruit. Future studies should explore the in-field and landscape factors determining pollination services in kiwifruit orchards in heterogenous landscapes such as the ones observed in our study region.
As mentioned above, regularly maintained management practices that impact insects, such as pesticide use, vegetation management, or hail net coverage steadily affect pollinator communities and their behavior [19,42,43]. Four orchards showed pollination deficits in both years and were signaled by the producers as regularly having problems with pollination (author’s field observations). These orchards deserve a thorough evaluation of the management practices that may be impacting pollination services, as it can be one of the main factors affecting productivity. For instance, two of these orchards are covered by a hail net, which is an important factor affecting pollination, as these structures reduce ventilation and pollen movement and restrict movements from pollinators and other beneficial insects and can change or reduce the visual cues necessary for their orientation [20,43]. Also, in New Zealand, a high number of bees were observed foraging uncovered orchards, with flowers in these orchards receiving more visits when compared to the visitation rate observed in covered kiwifruits orchards [43]. Additionally, orchard A3, was very affected by the bacterial disease Pseudomonas syringae pv. Actinidiae, which is an important disease in kiwifruit and contributes to pollination deficit through significant impacts in male flowering and pollen availability, with negative consequences in kiwifruit production [20].

3.3. Pollination Services, Production and Market Value

In our study, pollen supplementation tended to increase productivity and the economic value of kiwifruits, either by increasing total fruit production (evident in several orchards in 2018) or by increasing fruit caliber (evident in several orchards in 2019), both affecting the total economic revenue for the producer. This was more evident in orchards that showed pollination deficits, further highlighting the importance of adequate pollination for kiwifruit production. The importance of pollination for kiwifruit weight and size, which determine economic values, is well documented and has been demonstrated in several studies (e.g., [13,18,44,45]), although none was made at such a large scale. Also, and despite its importance, pollination is still considered the least understood management factor in kiwifruit orchards [24]. A better understanding of the pollination needs at a local scale and across large geographic areas, as we did here, is thus fundamental to improve production and to a more efficient use of resources, which might ultimately result in higher profit for kiwifruit producers.
The results presented in this study reflect a maximum potential gain following optimal pollination. Future studies should incorporate a cost-benefit analysis that considers not only the cost of the various management practices but also potential benefits to pollination services. As observed in our study, it is also important to consider that pollination deficits may vary in space and time due to variable pollinator communities, which may also be a determinant of the management practices needed to improve pollination in some orchards (e.g., the need for artificial pollination). Nevertheless, it is clear that kiwifruit production, and hence producers’ monetary gain, can be increased by adequate pollination. Natural pollination could potentially be enhanced by the use of insect friendly practices such as increasing wild bee nesting sites, the maintenance of diverse hedgerow, or the implementation of green infrastructures as food supplies before and after the crop flowering period [40,46,47]. Additionally, producers may consider complementing wild pollinator’s community with managed bee colonies or, ultimately, and in very particular cases (due to the costs associated and other problems), complementing natural pollination with artificial pollen application [16,20].

4. Materials and Methods

4.1. Studied Orchards

The pollination experiments were conducted in 9 orchards in the spring of 2018 and in 22 orchards in the spring of 2019, 8 of which were surveyed in both years (Table S4). The orchards are located in the North and Centre of Portugal (Figure S3). All but one orchard (orchard G with organic farming) were conventionally managed. Nineteen orchards were planted with Actinidia deliciosa and three orchards were planted with A. chinensis (both dioecious species). Details of kiwifruit orchards and kiwifruit male and female varieties used in this study are given in Table S4. For each orchard locality, environmental variables were extracted from climate data obtained from IPMA (Portuguese Institute for Sea and Atmosphere). Environmental variables extracted include precipitation during the flowering period, because of its effect on pollen movement and insect activity, and average mean temperature in the 1.5 months before flowering peak because of its effect on flowering synchrony and duration [33].
In the study region, the flowering period generally starts in mid-May, but in 2018 the flowering was delayed and occurred in late May–early June. In 2019, flowering occurred in the second half of May but showed high heterogeneity among orchards. The most abundant pollinators found in kiwifruit orchards are Apis mellifera, followed by Bombus spp., and more rarely other wild bees (e.g., Halictus spp.) and Syrphidae (e.g., Eristalis, Melanostoma and Sphaerophoria species; [18,27]). Pollinator richness and abundance was shown to vary among the 22 sampled orchards [27]. A parallel study focused on pollinator community was undertaken by [27] during 2019 at the same 22 orchards samples here. Thus, data on overall pollinator richness and abundance, and on the abundance of the two main pollinator groups, honeybees and Bombus, were taken from [27] and covers the 22 orchards sampled in 2019.

4.2. Pollination Experiments

To quantify pollination deficits in kiwifruit orchards we set up a pollination experiment in the flowering seasons between 2018 and 2019, with the following pollination treatments: (1) open pollination (O), i.e., unmanipulated flowers left open for natural levels of pollination, thereby enabling to quantify the services provided by the available pollinator community; and (2) supplementary pollination (S), i.e., flowers left open to natural levels of pollination and pollinated with fresh pollen from up to eight male flowers of the varieties available in the orchard, thereby enabling to quantify yield under optimal pollination services. In 2018, the following additional treatment was conducted to quantify the contribution of wind: (3) emasculated open pollination (W), i.e., flowers in which petals and anthers were removed at the bud stage and left open, thereby enabling to quantify the services provided by wind. Preliminary pollinator visitation observations revealed that removing anthers and petals (insect attractants) prevented insect visitation to emasculated flowers (author’s field observations). Based on the consistent results obtained in 2018, on the reduced contribution of wind reported in the literature (e.g., [12,13]) and on the increased number of fields surveyed in 2019, in the second year we decided to exclude the evaluation of wind contribution. In detail, for each orchard and variety, at peak of orchard flowering, we selected 30 female plants separated by at least 3 m, in a line at the middle of the orchard, in most cases covering the entire length of the orchard. Per plant, 2–3 flowers, depending on whether two or three pollination treatment were applied, located in the same branch, to guarantee similar resource availability, were marked and assigned to a given treatment. Receptive flowers with petals starting to fall and sticky stigmatic surfaces [16] were selected for both open and supplementary pollinations. Supplementary pollinations were performed by gently rubbing anthers of up to eight freshly collected male flowers with dehiscent pollen (visible by whitish anther pores and visible clouds of pollen when rubbing the flowers). Male flowers were collected from male varieties present in the orchard.

4.3. Fruit Processing

The collection of kiwifruits was coordinated with each orchard producer so that fruits could be collected as close as possible to the harvest date. Thus, fruits were collected, on average two days before the harvesting of the orchard. At harvest time, the fruits from marked flowers were collected and the percentage of flowers that set fruit was calculated (fruit set). All fruits resulting from the flowers marked for pollination treatments were weighed, and measured for short (at the equatorial zone) and long (apical-basal axis) diameters with a digital caliper; all deformations and skin damages were also registered. Additionally, we determined flesh firmness, °Brix, acidity (n = 10 and all fruits per treatment in 2018 and 2019, respectively), and dry matter content (n = 10 and 15 fruits per treatment in 2018 and 2019, respectively). Flesh firmness was measured by puncture with a handheld digital penetrometer (Fruit hardness tester, STEPS Systems) fitted with a flat 8 mm diameter plunger after removing 1 mm skin at the fruit equator. °Brix was measured with a hand-held digital refractometer (Kern Optics), and acidity was measured with an ATAGO PAL-Easy ACID F5 by taking a juice sample from longitudinally cut fruit half of each fruit. For acidity, the juice sample was diluted in distilled water (1:50) according to the manufactures’ instructions. The second half of each fruit was used to determine dry matter content, recording the fresh weight and the dry weight after drying at 60 °C until achieving a constant weight. Dry matter content was given by: (dry weight/fresh weight) * 100 and is expressed in percentage.
In 2018, the fruit half used for °Brix and acidity was also used to determine total soluble sugars (TSS), starch and total antioxidant capacity (TAA). The samples were triturated with a blender and three subsamples of 10 mg were weighed, placed in extraction buffer, and stored at −20 °C until analysis. Total soluble sugars were quantified following Irigoyen et al. [48]. For this, one subsample per fruit was homogenized with ethanol (80%, v/v) and incubated at 80 °C for 1 h; after centrifugation at 5000× g for 10 min at 4 °C, one aliquot of the supernatant was mixed with an anthrone solution (40 mg anthrone + 1 mL dH2O + 20 mL H2SO4) and incubated for 10 min at 100 °C; samples were cooled and centrifuged as described above. Starch was quantified according to Osaki et al. [49]. For this, another subsample was homogenized with perchloric acid (30%, v/v) and incubated at 60 °C for 1 h. After centrifugation at 10,000× g for 10 min at 4 °C, one aliquot of the supernatant was incubated with the anthrone solution (as described for TSS), followed by another centrifugation (5000× g for 10 min and 4 °C). For both TSS and starch, the absorbance of the supernatant was read at 625 nm using a Multilabel Reader (EnSpire™, PerkinElmer, Norwalk, USA) and the contents calculated using the standard curve of glucose (y = 5.1061x + 0.2279, r2 = 0.9995, for TSS and y = 1.0514x + 0.2719, r2 = 0.9919 for starch). For TAA, the last subsample was homogenized with 1 mL of methanol and incubated at 40 °C for 30 min. The mixture was centrifuged at 15,000× g for 15 min at 4 °C and one aliquot of the supernatant was mixed with ABTS [2,20-azino-bis (3-ethylbenzothiazoline-6-sulphonic acid)]. The absorbance was read at 734 nm (Re et al., 1999) and TAA was calculated using a gallic acid standard curve (y = 1094.5x − 38.977, r2 = 0.9855).

4.4. Pollination Deficit

Pollination deficit was calculated using the commonly used formula to calculate pollen limitation, and was obtained for each plant based on fruit weight, according to Larson and Barrett [50]: PL = 1 − O/S, where O is the fruit weight of the open pollination treatment, and S is the fruit weight of the supplementary pollination treatment. Positive values resulting from higher fruit weight in supplementary versus open pollination treatments indicate pollination deficit and, consequently, the occurrence of pollination deficits, while zero or negative values indicate no pollination deficit. In cases where we obtained negative values for pollination deficit, we considered that these indicate the absence of pollination deficit and were treated as zero in subsequent calculations in this study. Despite our attempts to minimize them, negative values may result from within individual differences (e.g., ovary size, resource allocation, shadow) known to affect fruit weight within plants [51,52,53].

4.5. Commercial Grading

Commercial grading of every sampled kiwifruit was assigned following standard grading tables provided by the Portuguese Association of Kiwifruit producers (APK—Associação Portuguesa de Kiwicultores). The tables assign fruits into classes and calibers. First, kiwifruits were assigned to class I or class II based on their short to long diameter ratio (S:L), deformation, and skin damages. Fruits with an S:L below 0.75, deformed or bearing skin damage, were assigned to class II, while fruits with an S:L equal or above 0.75, well-formed and without skin damages, were assigned to class I. Secondly, within classes, fruits were assigned to one of 11 caliber categories (namely, 18, 20, 23, 25, 27, 30, 33, 36, 39, 42, or 46 category) based on their weight (each category representing the average number of fruits of that size to attain 1 kg of kiwifruit). Fruits with weights lower than 65 g were considered with no market value.

4.6. Economic Analysis

We estimated effects of pollination deficit on kiwifruit production by comparing yields under natural (open pollination) and optimal pollination (supplementary pollination). Orchard production values of 2018 and 2019 were provided by each kiwifruit producer and taken as the values corresponding to production under natural pollination. Production under optimal pollination was estimated, taking into account the gain in fruit set in supplementary pollinated fruits compared to natural pollination and the pollination deficit values obtained for each orchard according to the following equation: Popt = Po + (Po * PL) + (Po * Fs), where Popt is the production under optimal pollination, Po is the production under open pollination, PL is the pollination deficit, and Fs is the difference in fruit set between supplementary and open pollination treatments.
We estimated the economic effects of pollination deficit by comparing the monetary gain under open and optimal pollination. Monetary gain was obtained by calculating the amount, in Euros, corresponding to each orchard production, both under open and supplementary pollination, considering the percentage of fruits in each class and caliber, and the corresponding monetary value. Average prices paid to the producer according to class and caliber assignments were provided by APK. Calculations were done according to the following equation: Monetary gain (euros) = ∑ P * Ci * Ei, where P is the production (t/ha) of a given field under natural or optimal pollination (calculated as described above), Ci is the proportion of fruits in each class and caliber combination, and Ei is the price, in Euros, payed per kg of fruit for each class and caliber combination in a given year. These were then used to calculate the difference in monetary gain resulting from supplementary pollination.

4.7. Statistical Analyses

Generalized linear models (GLM) were used to explore differences between pollination treatments within orchards on fruit parameters (namely, fruit set, weight, dry matter, flesh firmness, °Brix, acidity, S:L, TSS, starch, TAA, and fruit caliber). Fruit weight and quality parameters were analyzed using a Gaussian distribution with the identity link function, and fruit set and fruit distribution by caliber categories were analyzed using a binomial distribution with the logit link function. For the GLM analysis, fruit caliber was grouped into three categories, adapted from NZKGI [25]: (1) large fruits, which are fruits with high market quality and includes calibers 18 to 30; (2) medium fruits, which includes calibers 33 to 36; and (3) small fruits, which includes fruits in calibers 39, 42, 46, fruits in class II, and non-marketable fruits. The effect of pollination treatment on fruit set and fruit parameters across orchards was analyzed using generalized linear mix models (GLMM) with pollination treatment as a fixed factor and orchard and kiwifruit variety as random factors. The S:L was transformed using arcsin (sqrt) and the TAA was transformed using sqrt. The effect of year and pollination treatment within year on fruit set, fruit parameters, and pollination deficit, production and monetary gain across orchards was analyzed using GLMM with year and pollination treatment nested within year as fixed factors, and orchard, kiwifruit variety, precipitation, and average mean temperature as random factors. Model validation was performed on the residuals by checking heteroscedasticity and normality [54] and response variables were transformed when necessary to achieve such requirements. One-sample t-test was used to test if pollination deficits within orchard differed significantly from zero. Regression models using pollinator richness, honeybee abundance, Bombus spp. abundance, precipitation, and temperature as predictors were used to explain patterns in pollination deficits and production across orchards. Only orchards for which pollinator variables and production were obtained in the same year were included in this analysis. Other pollinators’ related variables (e.g., overall pollinator abundance, wild pollinator richness and abundance) showed significant correlations and, thus, were excluded from the regression analyses.
All analyses were done using R version 3.3.2 [55] using the package “car” for Type-III analysis of variance [56] and “lme4” for generalized linear models and generalized linear mixed models [57], and “multcomp” for multiple comparisons after Type-III analysis of variance in GLMM [58] and TukeyHSD for multiple comparisons after analysis of variance in GLM.

5. Conclusions

The quantification of pollination deficits across the production area of Portugal enabled us to conclude that pollination deficits are variable both in time and space, likely resulting from a combination of multiple factors that determine pollen availability and pollinator’s communities. Although pollen supplementation increased productivity and the economic values of kiwifruits, positively affecting the producer’s economic revenue, in general, a representative number of kiwifruit orchards showed low pollination deficits, suggesting sufficient pollination services. Still, consistent pollination deficits in some orchards suggest the existence of factors that consistently affect pollination success. Consequently, future studies should address how orchard management practices and landscape factors may affect kiwifruit production, contributing with relevant data to design guidelines towards productive and sustainable agroecosystems.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10071273/s1, Table S1: Fruit set; Tables S2 and S3: Fruit weight and quality parameters in 2018 and 2019, respectively; Figure S1: Distribution of kiwifruits by caliber category; Figure S2: Change in monetary gain; Table S4: Details of kiwifruit orchards; Figure S3: Distribution of sampling orchards.

Author Contributions

Experiment performance and sample collection, H.C., C.S., V.C., H.G., M.C., S.C., J.L.; chemical analyses, H.C., M.C.D., C.S.; statistical analysis of data, H.C.; writing—original draft preparation, H.C.; writing—review and editing, H.C., C.S., V.C., H.G., M.C., S.C., J.L., M.C.D. Funding acquisition, S.C., J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Operational group I9Kiwi—Developing strategies for the sustainability of kiwifruit production through creation of an added value product, funded by PDR2020. Project RENATURE—“Programa Operacional Regional do Centro 2014–2020 (Centro2020)—CENTRO-01-0145-FEDER-000007” financed the work of HC, SC and MC; SC was also financed by CULTIVAR project (CENTRO-01-0145-FEDER-000020), co-financed by the Regional Operational Programme Centro 2020, Portugal 2020 and European Union, through European Fund for Regional Development (ERDF). MCD (SFRH/BPD/100865/2014) was funded by national funds (OE), through FCT, I.P., in the scope of the framework contract foreseen in the numbers 4–6 of the article 23, of the Decree-Law 57/2016, 29 August, changed by Law 57/2017, 19 July. HC as also financed by the later. This work was carried out at the R&D Unit Center for Functional Ecology—Science for People and the Planet (CFE), with reference UIDB/04004/2020, financed by FCT/MCTES through national funds (PIDDAC).

Data Availability Statement

Data are available from the authors upon request.

Acknowledgments

The authors are thankful to APK (Associação Portuguesa de Kiwicultores) and to all kiwifruit producers and field technicians involved in this study for collaboration and allowing development of the experiments in their orchards and providing the necessary information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Klein, A.-M.; Vaissière, B.E.; Cane, J.H.; Steffan-Dewenter, I.; Cunningham, S.A.; Kremen, C.; Tscharntke, T. Importance of pollinators in changing landscapes for world crops. Proc. R. Soc. B Biol. Sci. 2007, 274, 303–313. [Google Scholar] [CrossRef] [Green Version]
  2. Aizen, M.A.; Garibaldi, L.A.; Cunningham, S.A.; Klein, A.M. How much does agriculture depend on pollinators? Lessons from long-term trends in crop production. Ann. Bot. 2009, 103, 1579–1588. [Google Scholar] [CrossRef] [PubMed]
  3. Knight, T.M.; Steets, J.A.; Vamosi, J.C.; Mazer, S.J.; Burd, M.; Campbell, D.R.; Dudash, M.R.; Johnston, M.O.; Mitchell, R.J.; Ashman, T.-L. Pollen Limitation of Plant Reproduction: Pattern and Process. Annu. Rev. Ecol. Evol. Syst. 2005, 36, 467–497. [Google Scholar] [CrossRef] [Green Version]
  4. Vaissière, B.; Freitas, B.; Gemmill-Herren, B. Protocol to Detect and Assess Pollination Deficits in Crops: A Handbook for Its Use; Food and Agiculture Organization of the United Nations: Rome, Italy, 2011. [Google Scholar]
  5. Potts, S.G.; Biesmeijer, J.C.; Kremen, C.; Neumann, P.; Schweiger, O.; Kunin, W.E. Global pollinator declines: Trends, impacts and drivers. Trends Ecol. Evol. 2010, 25, 345–353. [Google Scholar] [CrossRef]
  6. Freitas, B.M.; Vaissiére, B.E.; Saraiva, A.M.; Carvalheiro, L.G.; Garibaldi, L.A.; Ngo, H. Identifying and assessing pollina-tion deficits in crops. In Pollination Services to Agriculture: Sustaining and Enhancing a Key Ecosystem Service; Taylor & Francis: Oxfordshire, UK, 2016; pp. 17–32. ISBN 9781317445678. [Google Scholar]
  7. Lebuhn, G.; Droege, S.; Connor, E.; Gemmill-Herren, B.; Azzu, N. Protocol to Detect and Monitor Pollinator Communities–Guidance for Practitioners; Food and Agiculture Organization of the United Nations: Rome, Italy, 2016; ISBN 9789251089781. [Google Scholar]
  8. Garratt, M.P.D.; Evans, D.J.; Moss, R.L.; Dodson, E.D.; Jenner, C.; Potts, J.C. Pollination deficits in UK apple orchards. J. Pollinat. Ecol. 2014, 12, 9–14. [Google Scholar]
  9. Knapp, J.L.; Osborne, J.L. Courgette Production: Pollination Demand, Supply, and Value. J. Econ. Èntomol. 2017, 110, 1973–1979. [Google Scholar] [CrossRef] [Green Version]
  10. Webber, S.M.; Garratt, M.P.D.; Lukac, M.; Bailey, A.P.; Huxley, T.; Potts, S.G. Quantifying crop pollinator-dependence and pollination deficits: The effects of experimental scale on yield and quality assessments. Agric. Ecosyst. Environ. 2020, 304, 107106. [Google Scholar] [CrossRef]
  11. Pomeroy, N.; Fisher, R.M. Pollination of kiwifruit (Actinidia deliciosa) by bumble bees (Bombus terrestris): Effects of bee density and patterns of flower visitation. N. Z. Èntomol. 2002, 25, 41–49. [Google Scholar] [CrossRef]
  12. Testolin, R.; Vizzotio, G.; Costa, G. Kiwifruit pollination by wind and insects in Italy. N. Z. J. Crop. Hortic. Sci. 1991, 19, 381–384. [Google Scholar] [CrossRef]
  13. Costa, G.; Testolin, R.; Vizzotto, G. Kiwifruit pollination: An unbiased estimate of wind and bee contribution. N. Z. J. Crop. Hortic. Sci. 1993, 21, 189–195. [Google Scholar] [CrossRef]
  14. Ferguson, A.R. Kiwifruit: The Wild and the Cultivated Plants. Adv. Food Nutr. Res. 2013, 68, 15–32. [Google Scholar] [CrossRef]
  15. Ferguson, A.R. Kiwifruit: A Botanical Review. In Horticultural Reviews; Janick, J., Ed.; AVI Publishing Company, INC: Westport, CT, USA, 1984; pp. 1–64. [Google Scholar]
  16. Tacconi, G.; Michelotti, V.; Cacioppo, O.; Vittone, G. Kiwifruit pollination: The interaction between pollen quality, pollination systems and flowering stage. J. Berry Res. 2016, 6, 417–426. [Google Scholar] [CrossRef] [Green Version]
  17. Oliveira, M.; Veloso, A.; Teixeira, M.; Veloso, F.; Antunes, D. Avaliação do efeito da polinização complementar em actinídea deliciosa CV. “Hayward” na Região de Entre—Douro e Minho. In Proceedings of the I Simpósium Nacional de Fruticultura, Alcobaça, Portugal, 12–13 October 2006; Ramos, A., Martins, J.M.S., Eds.; Associação Portuguesa de Horticultura: Lisbon, Portugal, 2009; pp. 97–101. [Google Scholar]
  18. Miñarro, M.; Twizell, K.W. Pollination services provided by wild insects to kiwifruit (Actinidia deliciosa). Apidologie 2015, 46, 276–285. [Google Scholar] [CrossRef] [Green Version]
  19. González, M.V.; Coque, M.; Herrero, M. Influence of pollination systems on fruit set and fruit quality in kiwifruit (Actinidia deliciosa). Ann. Appl. Biol. 1998, 132, 349–355. [Google Scholar] [CrossRef]
  20. Tacconi, G.; Michelotti, V. Artificial Pollination in Kiwifruit and Olive Trees. In Pollination in Plants; Segneanu, A.E., Orbeci, C., Lazau, C., Sfirloaga, P., Vlazan, P., Bandas, C., Grozescu, I., Eds.; InTech: London, UK, 2018; p. 29. ISBN 978-1-78923-237-0. [Google Scholar]
  21. Vicens, N.; Bosch, J. Weather-Dependent Pollinator Activity in an Apple Orchard, with Special Reference to Osmia cornuta and Apis mellifera (Hymenoptera: Megachilidae and Apidae). Environ. Èntomol. 2000, 29, 413–420. [Google Scholar] [CrossRef]
  22. Donati, I.; Cellini, A.; Buriani, G.; Mauri, S.; Kay, C.; Tacconi, G.; Spinelli, F. Pathways of flower infection and pollen-mediated dispersion of Pseudomonas syringae pv. actinidiae, the causal agent of kiwifruit bacterial canker. Hortic. Res. 2018, 5, 56. [Google Scholar] [CrossRef] [Green Version]
  23. Testolin, R.; Ferguson, A.R. Kiwifruit (Actinidia spp.) Production and Marketing in Italy. N. Z. J. Crop Hortic. Sci. 2009, 37, 1–32. [Google Scholar] [CrossRef]
  24. Goodwin, M. Pollination of Crops in Australia and New Zealand; Rural Industries Research and Development Corporation: Wagga Wagga, Australia, 2012; ISBN 9781742544021. [Google Scholar]
  25. NZKGI. Kiwifruit Book; NZKGI: Mount Maunganui, New Zealand, 2016. [Google Scholar]
  26. Antunes, M.D.; Terra, C. Da Produção à Comercialização; Universidade do Algarve: Algarve, Portugal, 2008. [Google Scholar]
  27. Gaspar, H. Impact of Insect Diversity and Abundance in Kiwi Orchard’s Production; University of Coimbra: Coimbra, Portugal, 2020. [Google Scholar]
  28. Martins, K.T.; Gonzalez, A.; Lechowicz, M.J. Pollination services are mediated by bee functional diversity and landscape context. Agric. Ecosyst. Environ. 2015, 200, 12–20. [Google Scholar] [CrossRef] [Green Version]
  29. Fijen, T.P.M.; Scheper, J.A.; Boekelo, B.; Raemakers, I.; Kleijn, D. Effects of landscape complexity on pollinators are moderated by pollinators’ association with mass-flowering crops. Proc. R. Soc. B Boil. Sci. 2019, 286, 20190387. [Google Scholar] [CrossRef] [Green Version]
  30. Craig, J.L.; Stewart, A.M.; Pomeroy, N.; Heath, A.C.G.; Goodwin, R.M. A review of kiwifruit pollination: Where to next? N. Z. J. Exp. Agric. 1988, 16, 385–399. [Google Scholar] [CrossRef]
  31. Corbet, S.A.; Chapman, H.; Saville, N. Vibratory Pollen Collection and Flower Form: Bumble-Bees on Actinidia, Symphytum, Borago and Polygonatum. Funct. Ecol. 1988, 2, 147. [Google Scholar] [CrossRef]
  32. Goodwin, R.M.; Congdon, N.M. Recognition and attractiveness of staminate and pistillate kiwifruit flowers (Actinidia deliciosa var. deliciosa) by honey bees (Apis mellifera L.). N. Z. J. Crop. Hortic. Sci. 2017, 46, 72–80. [Google Scholar] [CrossRef]
  33. McPherson, H.G.; Hall, A.J.; Stanley, C.J. Seasonal and regional variation in budbreak and flowering of kiwifruit vines in New Zealand. N. Z. J. Crop. Hortic. Sci. 1994, 22, 263–276. [Google Scholar] [CrossRef]
  34. Tuell, J.K.; Isaacs, R. Weather during bloom affects pollination and yield of highbush blueberry. J. Econ. Èntomol. 2010, 103, 557–562. [Google Scholar] [CrossRef]
  35. Blitzer, E.J.; Gibbs, J.; Park, M.G.; Danforth, B. Pollination services for apple are dependent on diverse wild bee communities. Agric. Ecosyst. Environ. 2016, 221, 1–7. [Google Scholar] [CrossRef] [Green Version]
  36. MacInnis, G.; Forrest, J.R.K. Pollination by wild bees yields larger strawberries than pollination by honey bees. J. Appl. Ecol. 2019, 56, 824–832. [Google Scholar] [CrossRef]
  37. Connelly, H.; Poveda, K.; Loeb, G. Landscape simplification decreases wild bee pollination services to strawberry. Agric. Ecosyst. Environ. 2015, 211, 51–56. [Google Scholar] [CrossRef]
  38. Nicholson, C.C.; Koh, I.; Richardson, L.L.; Beauchemin, A.; Ricketts, T.H. Farm and landscape factors interact to affect the supply of pollination services. Agric. Ecosyst. Environ. 2017, 250, 113–122. [Google Scholar] [CrossRef]
  39. Bartual, A.M.; Sutter, L.; Bocci, G.; Moonen, A.-C.; Cresswell, J.; Entling, M.; Giffard, B.; Jacot, K.; Jeanneret, P.; Holland, J.; et al. The potential of different semi-natural habitats to sustain pollinators and natural enemies in European agricultural landscapes. Agric. Ecosyst. Environ. 2019, 279, 43–52. [Google Scholar] [CrossRef]
  40. Eeraerts, M.; Smagghe, G.; Meeus, I. Pollinator diversity, floral resources and semi-natural habitat, instead of honey bees and intensive agriculture, enhance pollination service to sweet cherry. Agric. Ecosyst. Environ. 2019, 284, 106586. [Google Scholar] [CrossRef]
  41. Lázaro, A.; Alomar, D. Landscape heterogeneity increases the spatial stability of pollination services to almond trees through the stability of pollinator visits. Agric. Ecosyst. Environ. 2019, 279, 149–155. [Google Scholar] [CrossRef]
  42. Bartholomée, O.; Aullo, A.; Becquet, J.; Vannier, C.; Lavorel, S. Pollinator presence in orchards depends on landscape-scale habitats more than in-field flower resources. Agric. Ecosyst. Environ. 2020, 293, 106806. [Google Scholar] [CrossRef]
  43. Evans, L.J.; Cutting, B.T.; Jochym, M.; Janke, M.A.; Felman, C.; Cross, S.; Jacob, M.; Goodwin, M. Netted crop covers reduce honeybee foraging activity and colony strength in a mass flowering crop. Ecol. Evol. 2019, 9, 5708–5719. [Google Scholar] [CrossRef]
  44. Antunes, M.; Oliveira, M.; Teixeira, M.; Veloso, A.; Veloso, F.; Panagopoulos, T. Evaluation of the effect of complementary pollination on Actinidia deliciosa ’hayward’ in northwest portugal. Acta Hortic. 2007, 753, 347–352. [Google Scholar] [CrossRef]
  45. Goodwin, R.; McBrydie, H.; Taylor, M. Wind and honey bee pollination of kiwifruit (Actinidia chinensis‘HORT16A’). N. Z. J. Bot. 2013, 51, 229–240. [Google Scholar] [CrossRef] [Green Version]
  46. Castle, D.; Grass, I.; Westphal, C. Fruit quantity and quality of strawberries benefit from enhanced pollinator abundance at hedgerows in agricultural landscapes. Agric. Ecosyst. Environ. 2019, 275, 14–22. [Google Scholar] [CrossRef]
  47. Karamaouna, F.; Kati, V.; Volakakis, N.; Varikou, K.; Garantonakis, N.; Economou, L.; Birouraki, A.; Markellou, E.; Liberopoulou, S.; Edwards, M. Ground cover management with mixtures of flowering plants to enhance insect pollinators and natural enemies of pests in olive groves. Agric. Ecosyst. Environ. 2019, 274, 76–89. [Google Scholar] [CrossRef]
  48. Irigoyen, J.J.; Einerich, D.W.; Sanchez-Diaz, M. Water stress induced changes in concentrations of proline and total soluble sugars in nodulated alfalfa (Medicago sativd) plants. Physiol. Plant. 1992, 84, 55–60. [Google Scholar] [CrossRef]
  49. Osaki, M.; Shinano, T.; Tadano, T. Redistribution of carbon and nitrogen compounds from the shoot to the harvesting organs during maturation in field crops. Soil Sci. Plant Nutr. 1991, 37, 117–128. [Google Scholar] [CrossRef]
  50. Larson, B.M.H.; Barrett, S.C.H. A comparative analysis of pollen limitation in flowering plants. Biol. J. Linn. Soc. Linn. Soc. 2000, 69, 503–520. [Google Scholar] [CrossRef]
  51. Lawes, G.S.; Woolley, D.J.; Lai, R. Seeds and Other Factors Affecting Fruit Size in Kiwifruit. Acta Hortic. 1990, 257–264. [Google Scholar] [CrossRef]
  52. Vasilakakis, M.; Papadopoulos, K.; Papageorgiou, E. Factors affecting the fruit size of “hayward” kiwifruit. Acta Hortic. 1997, 419–424. [Google Scholar] [CrossRef]
  53. McPherson, H.G.; Richardson, A.C.; Snelgar, W.P.; Patterson, K.J.; Currie, M.B. Flower quality and fruit size in kiwifruit (Actinidia deliciosa). N. Z. J. Crop. Hortic. Sci. 2001, 29, 93–101. [Google Scholar] [CrossRef] [Green Version]
  54. Zuur, A.; Ieno, E.N.; Walker, N.; Saveliev, A.A.; Smith, G.M. Mixed Effects Models and Extensions in Ecology with R; Statistics for Biology and Health; Springer: New York, NY, USA, 2009; ISBN 978-0-387-87457-9. [Google Scholar]
  55. R Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  56. Fox, J.; Weisberg, S. An R Companion to Applied Regression, 3rd ed.; Sage: Thousand Oaks, CA, USA, 2019. [Google Scholar]
  57. Bates, D.; Mächler, M. Package ‘lme4’: Linear mixed-effects models using eigen and S4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  58. Hothorn, T.; Bretz, F.; Westfall, P. Simultaneous Inference in General Parametric Models. Biom. J. 2008, 50, 346–363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Mean fruit weight (g) ± SE of kiwifruits resulting from wind (W), open (O) and supplementary (S) pollination treatments applied in the nine orchards surveyed in 2018.
Figure 1. Mean fruit weight (g) ± SE of kiwifruits resulting from wind (W), open (O) and supplementary (S) pollination treatments applied in the nine orchards surveyed in 2018.
Plants 10 01273 g001
Figure 2. Yearly variation in pollination deficit in the orchards surveyed both in 2018 and 2019.
Figure 2. Yearly variation in pollination deficit in the orchards surveyed both in 2018 and 2019.
Plants 10 01273 g002
Figure 3. Pollination deficit (mean ± SE) in the nine orchards surveyed in 2018 (A) and in the twenty-two orchards surveyed in 2019 (B). Statistically significant differences (p < 0.05) in pollination deficit resulting from one-sample t-test analyses are indicated by an asterisk.
Figure 3. Pollination deficit (mean ± SE) in the nine orchards surveyed in 2018 (A) and in the twenty-two orchards surveyed in 2019 (B). Statistically significant differences (p < 0.05) in pollination deficit resulting from one-sample t-test analyses are indicated by an asterisk.
Plants 10 01273 g003
Figure 4. Kiwifruit production (t/ha) corresponding to current orchard pollination services (i.e., open pollination, O) and estimated production under optimal pollination services (i.e., supplementary pollination, S) for the nine orchards surveyed in 2018 (left) and the twenty-two orchards surveyed in 2019 (right).
Figure 4. Kiwifruit production (t/ha) corresponding to current orchard pollination services (i.e., open pollination, O) and estimated production under optimal pollination services (i.e., supplementary pollination, S) for the nine orchards surveyed in 2018 (left) and the twenty-two orchards surveyed in 2019 (right).
Plants 10 01273 g004
Table 1. Generalized mixed-effect model analyses of the effect of pollination treatments (O, Open; S, Supplementary; W, Wind) on fruit set and fruit parameters (S:L, fruit short to long diameter ratio; TSS, total soluble sugars; TAA, total antioxidant capacity and starch). Statistically significant differences are indicated by different letters and highlighted in bold.
Table 1. Generalized mixed-effect model analyses of the effect of pollination treatments (O, Open; S, Supplementary; W, Wind) on fruit set and fruit parameters (S:L, fruit short to long diameter ratio; TSS, total soluble sugars; TAA, total antioxidant capacity and starch). Statistically significant differences are indicated by different letters and highlighted in bold.
VariableMeans ± SEX2 or Fp Values
WindOpenSupplement
Fruit set (%)66.52 ± 3.18 (b)73.20 ± 2.81 (b)85.60 ± 2.22 (a)X22 = 24.72<0.001
Weight (g)61.82 ± 2.02 (c)89.76 ± 2.01 (b)107.17 ± 1.87 (a)F2, 527.6 = 164.51<0.001
Dry matter (%)14.75 ± 0.28 (b)15.14 ± 0.21 (a, b)15.66 ± 0.21 (a)F2, 214.0 = 164.51<0.001
Flesh firmness (N)40.59 ± 2.08 (b)42.74 ± 1.85 (a)43.96 ± 1.75 (a)F2, 215.03 = 6.120.003
°Brix (%)6.03 ± 0.236.12 ± 0.206.21 ± 0.20F2, 215.0 = 0.2460.783
Acidity (%)1.46 ± 0.03 (b)1.55 ± 0.03 (a)1.50 ± 0.03 (a, b)F2, 208.25 = 3.1750.044
S:L0.88 ± 0.00 (a)0.87 ± 0.00 (a, b)0.86 ± 0.00 (b)F2, 515.9 = 3.080.047
TAA (µmol eq. galic acid/g)3.08 ± 0.123.07 ± 0.102.98 ± 0.10F2, 211.0 = 0.3540.702
TSS (mg/g)4.81 ± 0.18 (b)5.33 ± 0.18 (a)5.31 ± 0.19 (a, b)F2, 201.2 = 3.610.029
Starch (mg/g)2.89 ± 0.062.93 ± 0.042.90 ± 0.05F2, 204.3 = 0.1180.889
Table 2. Generalized mixed-effect model analyses of the effect of year or the effects of year and pollination treatment (Open vs. Supplementary pollination) nested within year on pollination deficit, fruit set and fruit parameters (S:L, fruit short to long diameter ratio). Values for year and pollination treatments are given as means and SE of the mean. Statistically significant differences are highlighted in bold.
Table 2. Generalized mixed-effect model analyses of the effect of year or the effects of year and pollination treatment (Open vs. Supplementary pollination) nested within year on pollination deficit, fruit set and fruit parameters (S:L, fruit short to long diameter ratio). Values for year and pollination treatments are given as means and SE of the mean. Statistically significant differences are highlighted in bold.
FactorsYearYear: Pollination Treatment
Variables20182019X2 or Fp Values20182019X2 or Fp Values
OpenSupplementOpenSupplement
Fruit set (%)78.68 ± 1.5785.94 ± 0.88X21 = 1.510.21973.24 ± 2.4084.12 ± 1.9982.99 ± 1.3488.92 ± 1.13X22 = 26.27<0.001
Weight (g)97.95 ± 1.22103.65 ± 0.75F1, 20.26 = 0.2560.61991.89 ± 1.61104.56 ± 1.53100.57 ± 1.08106.67 ± 1.02F2, 1865 =44.30<0.001
Dry matter (%)15.56 ± 0.1414.03 ± 0.07F1, 20.0 = 0.063<0.00115.45 ± 0.1915.88 ± 0.1914.00 ± 0.1014.07 ± 0.10F2, 1140 = 2.720.066
Flesh firmness (N)47.24 ± 0.9070.38 ± 0.35F1, 7.9 =65.65<0.00146.89 ± 1.3047.57 ± 1.2470.14 ± 0.7870.62 ± 1.24F2, 1519 = 1.260.284
°Brix (%)6.19 ± 0.126.56± 0.05F1, 19.0 = 1.000.3316.34 ± 0.166.26 ± 0.166.42 ± 0.076.57 ± 0.07F2, 1479 = 0.210.814
Acidity (%)1.52 ± 0.021.54 ± 0.01F1, 12.4 = 0.130.7281.54 ± 0.031.51 ± 0.021.54 ± 0.011.55 ±0.01F2, 1511= 0.690.501
S:L0.88 ± 0.000.90 ± 0.00F1, 10.6 = 11.520.0060.88 ± 0.000.90 ± 0.000.88 ± 0.000.90 ± 0.00F2, 1789= 1.290.274
Pollination deficit0.11 ± 0.010.06 ± 0.01F1, 8.4 = 5.470.046------
Large fruits (%)50.98 ± 2.1158.35 ± 1.34X21 = 0.350.55442.86 ± 2.8959.93 ± 3.0054.45 ± 1.9262.19 ± 1.85X22 = 29.33<0.001
Medium fruits (%)24.06 ± 1.8120.53 ± 1.10X21 = 1.170.28224.49 ± 2.5123.60 ± 2.6020.18 ± 1.5520.88 ± 1.55X22 = 0.210.901
Small fruits (%)24.96 ± 1.8321.12 ± 1.11X21 = 0.670.41232.65 ± 2.7416.48 ± 2.2725.37 ± 1.6816.93 ± 1.43X22 = 38.40<0.001
Production (t/ha)15.97 ± 1.7618.07 ± 1.63F1, 7.3 = 0.070.80114.56 ± 2.4617.38 ± 2.5717.18 ± 2.2718.97 ± 2.38F1, 28 = 11.71<0.001
Monetary gain (€/ha)14,842.25± 1820.7318,479.98± 2105.59F1, 7.3 = 0.010.94812,385.04
± 2912.32
17,299.47
± 2664.70
16,935.04
± 2340.69
20,024.43
± 3075.26
F1, 28 = 7.350.003
Table 3. Regression analyses of the effect of environmental (precipitation and temperature) and pollinator’s related variables (pollinator richness, honeybee abundance, and Bombus spp. abundance) on pollination deficit and production. Statistically significant differences are highlighted in bold.
Table 3. Regression analyses of the effect of environmental (precipitation and temperature) and pollinator’s related variables (pollinator richness, honeybee abundance, and Bombus spp. abundance) on pollination deficit and production. Statistically significant differences are highlighted in bold.
VariablesPollination DeficitProduction
Estimatet ValuepEstimatet Valuep
Pollinator richness0.00070.2670.792−0.603−2.6000.021
Honeybee abundance−0.0004−2.2860.0340.0392.3990.031
Bombus spp. abundance−0.0003−0.8220.4210.0772.2060.045
Precipitation0.00050.5020.6220.1061.2740.224
Temperature0.01860.6980.494−2.086−0.7260.480
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Castro, H.; Siopa, C.; Casais, V.; Castro, M.; Loureiro, J.; Gaspar, H.; Dias, M.C.; Castro, S. Spatiotemporal Variation in Pollination Deficits in an Insect-Pollinated Dioecious Crop. Plants 2021, 10, 1273. https://doi.org/10.3390/plants10071273

AMA Style

Castro H, Siopa C, Casais V, Castro M, Loureiro J, Gaspar H, Dias MC, Castro S. Spatiotemporal Variation in Pollination Deficits in an Insect-Pollinated Dioecious Crop. Plants. 2021; 10(7):1273. https://doi.org/10.3390/plants10071273

Chicago/Turabian Style

Castro, Helena, Catarina Siopa, Vinícius Casais, Mariana Castro, João Loureiro, Hugo Gaspar, Maria Celeste Dias, and Sílvia Castro. 2021. "Spatiotemporal Variation in Pollination Deficits in an Insect-Pollinated Dioecious Crop" Plants 10, no. 7: 1273. https://doi.org/10.3390/plants10071273

APA Style

Castro, H., Siopa, C., Casais, V., Castro, M., Loureiro, J., Gaspar, H., Dias, M. C., & Castro, S. (2021). Spatiotemporal Variation in Pollination Deficits in an Insect-Pollinated Dioecious Crop. Plants, 10(7), 1273. https://doi.org/10.3390/plants10071273

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop