Next Article in Journal
Closed-Loop Optimal Control of Greenhouse Cultivation Based on Two-Time-Scale Decomposition: A Simulation Study in Lhasa
Next Article in Special Issue
Plant Composition and Feed Value of First Cut Permanent Meadows
Previous Article in Journal
Comparative Assessment of Genetic Variability Realised in Doubled Haploids Induced from F1 and F2 Plants for Response to Fusarium Stalk Rot and Yield Traits in Maize (Zea mays L.)
Previous Article in Special Issue
Evaluating Growth, Biomass and Cannabinoid Profiles of Floral Hemp Varieties under Different Planting Dates in Organic Soils of Florida
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dissecting the Genotype × Environment Interaction for Potato Tuber Yield and Components

1
Department of Agriculture, Food and Environment (Di3A), University of Catania, 95123 Catania, Italy
2
Institute of BioEconomy, National Research Council (CNR-IBE), 95126 Catania, Italy
*
Author to whom correspondence should be addressed.
Agronomy 2023, 13(1), 101; https://doi.org/10.3390/agronomy13010101
Submission received: 30 November 2022 / Revised: 20 December 2022 / Accepted: 26 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Sustainable Management of Herbaceous Field Crops)

Abstract

:
In the Mediterranean Basin, potato (Solanum tuberosum L.) is a major staple crop, the yield performance of which can vary dramatically based on pedo-climatic conditions and genotype. Hence, dissecting the genotype × environment interaction (GEI) in this region is mandatory for the setup of high-yielding and stable potato genotypes, also considering its importance for local potato breeding and the development of protected geographical identifications. Therefore, this research evaluated the marketable tuber yield (MY) and several yield components of five potato genotypes (Arizona, Generosa, Levante, Paradiso, and Vogue) over 4 years (2019, 2020, 2021, and 2022) and three locations in Sicily (Southern Italy) by additive main effects and multiplicative interactions (AMMI) and genotype main effects + genotype × environment interaction (GGE) biplot analyses. From combined ANOVA emerged a high significance of GEI, with the environment that provided the most considerable extent of variation for the most of the productive traits. The AMMI and GGE analyses identified Arizona as the best leading genotype in the studied area by virtue of its high productivity (44.5 t ha−1 of mean MY) coupled with stability, followed by Generosa (46.5 t ha−1). Ideal environments (location × year) were highly dependent on the productive trait, but most of them belonged to Acireale, characterized by fertile soils. According to our results, this approach could be recommended for breeding programs and commercial cultivation in the studied regions, along with the setup of potato protected geographical identifications.

1. Introduction

The potato (Solanum tuberosum L.) is recognized as the fourth-highest crop in terms of global production quantity, with about 359 million tons of tubers deriving from 16.5 million hectares [1]. Its present use is mainly focused on human nutrition in the form of both fresh produce consumption and various industrially processed versions. Potatoes, in fact, provide a lot of energy and are an excellent source of carbohydrates (starch and free sugars), good-quality proteins, and numerous minerals and antioxidants [2,3]. Thanks to the potato’s remarkable ability to adapt to the elements of the climate, it is grown on all continents, in developed and developing countries, at different altitudes (from the plains to the mountains), and at different times of the year. It is a staple food in the Mediterranean Basin, where the surface area dedicated to this crop is ~3 million hectares with an overall production of 28 Mt of tubers [1]. In this area, it is cultivated not only in the usual cycle (spring–summer) but also in a winter–spring cycle (from November–January to March–early June) for early production, especially in several coastal countries such as southern Italy, Tunisia, Egypt, Cyprus, Israel, Lebanon and Turkey [4,5,6].
The large variability in tuber yield is a consequence of a multitude of factors, namely genotype, soil properties, climate, and agronomic management [4,5,7,8,9,10,11]. In a given environment, a cultivar provides a maximum yield when it is well-adapted to the prevailing climate conditions and belongs to a maturity class that matches the available number of growing days. Stable and high-yielding cultivars are recommended for cultivation, however tuber yield traits of cultivars can vary from environment to environment and within a location for different growing seasons. The basic cause of low stability is a wide occurrence of genotype × environment interaction (GEI). GEI is a complex phenomenon involving agroecological, climatic, and agronomic conditions, as well as all the physiological and genetic factors that determine plant growth and development [12]. It results in the inconsistent performance of genotypes across environments, caused either by differential responses of the same set of genes to changes in the environment or by the expression of different sets of genes in different environments [13]. In order to prepare high-yielding and stable genotypes, genotype × environment interaction (GEI) is one of the most important and complicated issues for breeding programs [14]. Therefore, stability and adaptation studies across multi-environmental trials (METs) are imperative, also considering that most of the traits of agronomic importance are hereditary quantitative. Tuber yield, for instance, which is a character that presents continuous distribution, has polygenic inheritance and is highly influenced by variations in the environment [15]. The farmers or growers need cultivars that show high performance for yield and other essential agronomic traits and are stable, i.e., perform well over both different environmental conditions and years. Usefulness in dissecting the GEI is not only related to the development of widely adapted and stable genotypes across environments, but also to the improvement of yield traits of existing cultivars under conditions that positively influence such traits [16].
Several statistical techniques have been developed to study the GEI effects, measure stability, and facilitate varietal recommendations in METs. In addition to linear mixed models (LMMs) [17], two actual methods are widely used to dissect the GEI: the additive main effects and multiplicative interaction (AMMI) analysis and the genotype main effect plus genotype × environment interaction (GGE) biplot analysis [14,18]. AMMI and GGE analyses for yield traits stability have been studied in several crops including maize [19], rice [20], sugarcane [21], sorghum [22], Bambara groundnut [23], etc. In potatoes, AMMI analysis was used for cultivar recommendation in Peru on purple-pulp colored genotypes [24] and in Poland on post-registration series cultivars [25]; the GGE biplot analysis was performed on Russet Burbank, Umatilla Russet, and Cal White cultivars in Canada [26] and eight potato genotypes in Ethiopia [27]. The only study in which AMMI and GGE analysis has been conducted simultaneously was realized in Kenya to estimate the magnitude of GEI for potato tuber yield and bacterial wilt resistance in 42 breeding families [28]. However, this research was mainly focused on bacterial wilt resistance. To the best of our knowledge, although the GEI on potato tuber yield was partially investigated in previous research [26,27,29,30], a systematic study on potato yield by AMMI and GGE biplot analyses is lacking. Furthermore, there is a need to incorporate the GEI into local potato breeding programs for the setup of Protected Geographical Identifications. In a given geographical area, this could be useful to identify well-adapted and stable potato genotypes with productive characteristics that comply with quality standards, with the goal of better matching the needs of farmers and distinct market demands. For this reason, the present study aimed to investigate the GEI and determine the most stable and productive potato genotypes (five) over four seasons and three locations in Sicily based on AMMI and GGE analysis.

2. Materials and Methods

2.1. Study Sites

The present experiment was performed throughout three locations sited in Sicily (Southern Italy) on the southeast foothills of Mount Etna, an area traditionally devoted to potato cultivation: Acireale (37°60′ N 15°16′ E, 80 m a.s.l.), Giarre (37°74′ N 15°18′ E, 50 m a.s.l.), and Mascali (37°77′ N 15°11′ E, 60 m a.s.l.). According to the USDA soil classification [31], the soils over the three locations are Typic Xerochrepts type, with sandy texture and low water field capacity. Under a chemical point of view, in Acireale the soil had higher organic matter and macro (P, N, K, Ca, and Mg) and micro (Zn, Cu, and Mn) element contents, as well as a higher cation exchange capacity than Giarre and Mascali; in Giarre, the soil stood out for the high amount of both Fe and exchangeable Na. Table 1 shows the detailed physical and chemical characteristics of the −30 cm soil layer of the three locations under study, carried out in accordance with the procedures approved by the Italian Society of Soil Science [32]. The experimental trials last four growing seasons from 2019 to 2022, hereafter referred to as Season I, Season II, Season III, and Season IV. The zone has a typical Mediterranean climate, with mild–wet winters, springs with low rainfall and raising temperatures, and with most rainfall falling during the autumn–early winter period, which allows the growth of the winter–spring cycle potato.

2.2. Genotypes

Five potato genotypes—namely Arizona, Generosa, Levante, Paradiso, and Vogue (hereafter referred to as G1, G2, G3, G4, and G5, respectively)—were evaluated in this study. Except for Arizona, developed in 2009, the other genotypes are new cultivars recently introduced (between 2015 and 2018), never evaluated in the studied area, and whose pedigrees have been taken from van Berloo et al. [33]. Arizona (UK 150-19D22 × Mascotte) and Generosa are AB cooking type (based on the European Association for Potato Research cooking classification) cultivars with a medium-to-late cycle, high vigor, and low resistance to common scab. Levante (AR 01-3218 × Almera) is a B cooking type featuring excellent organoleptic characteristics with a medium-to-late cycle and high vigor, and high tolerance to late blight and potato cyst nematodes (PCN) 1 and 4. Paradiso (Rodeo × Valor) is a B cooking type characterized by very large tubers and high tolerance to late blight. Vogue (Marabel × Amorosa) is an AB cooking type characterized by high tuberification and tuber-enlargement speed, with medium vigor and high tolerance to common scab and PCNs 1 and 4. All cultivars are Dutch and high yielding, skin and flesh yellow colored, oval and long-oval shape tubers, and are potentially suitable for the early potato market.

2.3. Design and Agronomic Management

The experimental design used at each location was a 3.75 × 4.5 m randomized block with three replications. Each block consisted of five genotypes (Arizona, Generosa, Levante, Paradiso, and Vogue), i.e., six rows each consisting of 90 plants at 5.3 plants m−2 (0.25 m × 0.75 m). The external rows and two plants of each row-end were considered as the border. Planting was done by hand in the first week of December throughout the four growing seasons, using disease-free and certified A-class sprouted seed tubers. Agronomic management (tillage, fertilization, irrigation, and weed control) was carried out as per recommendations for the local farming practices and with uniformity criteria across the three locations. In detail, a 30 cm depth plowing followed by harrowing in October was performed, together with a fertilization program consisting of a pre-planting application of an NPK synthetic fertilizer at the dose of 50-80-120 kg ha−1 of N (as ammonium nitrate), P2O5 (as double perfosphate), and K2O (as potassium sulphate), respectively, plus a further 70 kg ha−1 application of N (as ammonium nitrate) as top dressing. About 70 mm, 100 mm, 50 mm, and 120 mm of irrigation water were provided by furrow irrigation to fully satisfy crop water requirements in seasons I II, III, and IV, respectively.

2.4. Yield and Yield Components

Potato yield and yield components were determined for each location and all growing seasons by manually harvesting tubers at about 120 days after planting (DAP), i.e., when the leaves and stems were bleached and dried, from 52 plants from the central rows of each plot. Initially, tubers were classified as marketable (unitary weight > 20 g) and unmarketable (unitary weight < 20 g, or deformed, green, and diseased), with each category counted and weighed separately. Only marketable tubers were used to determine marketable yield (MY), the number of marketable tubers per plant (NMTP), and the average marketable tuber weight (AMTW). Based on the measurement of the transverse diameter, marketable tubers were also distributed into three dimensional classes: <40 mm, 40–60 mm, and >60 mm. Then, the percentage incidence of each size class on the overall marketable production was calculated as caliber tubers (CT). Discarded production (unmarketable tuber, UT) was determined as the % weight of unmarketable tubers of the total tuber yield (marketable + unmarketable). A representative sample of marketable tubers—at least fifty tubers per replicate for each genotype and location—was washed and dried at 65 °C up to a constant weight to evaluate dry matter content. Tuber specific gravity (TSG) was determined from a sample weight of 5 kg by calculating the quotient weight in air divided by the weight in air minus the weight in water [34].

2.5. Weather Description

In each location, meteorological measurements (rainfall and minimum and maximum air temperatures) were recorded daily using a data logger (CR21, Campbell Scientific, Logan, UT, USA) positioned about 20 m from the experimental fields. The growing degree days (GDD) were also calculated with the formula GDD = [(Tmax − Tmin)/2] − Tbase, where Tmax and Tmin are the daily maximum and minimum air temperatures, respectively, and Tbase is the temperature below which the crop growth does not progress, which is 7 °C. The three locations under study showed different climatic characteristics. Considering the months in which the growing season (December–May) takes place during the 10-year period 2013–2022, Acireale was characterized by the lowest average minimum and maximum temperatures (10.1 and 16.5 °C, respectively) and accumulated GDDs (1152), whereas Mascali showed the highest values (11.0 and 17.5 °C, and 1322, respectively). The temperature differences between Acireale and Mascali were more marked in March, with monthly accumulated GDD values of 166 and 199, respectively. Regarding rainfall, the lowest fell in Acireale (831 mm) and the highest in Giarre (1040 mm). Over the three locations, about 80% of rainfall was concentrated during the December–March period.
With regard to the weather conditions during the experimental field trials across the three locations (Figure 1), compared to the 10-year period, there emerged:
  • higher minimum and maximum temperatures in March (+1.4 and +0.9 °C in Acireale, +1.1 and +0.9 °C in Giarre, +1.1 and +0.9 °C in Mascali) in Season I;
  • lower rainfall throughout the growing season (612 vs. 831 mm in Acireale, 732 vs. 1040 in Giarre, and 702 vs. 960 in Mascali) in Season II, but a high concentration in March when rainfall was 41%, 37%, and 34% of the total for the growing season in Acireale, Giarre, and Mascali, respectively;
  • higher rainfall (1150 vs. 831 mm in Acireale, 1176 vs. 1040 mm in Giarre, and 1206 vs. 960 mm in Mascali) in Season III, mostly concentrated in December and March, as well as higher minimum and maximum temperatures in April and May, which led to a strong increase of accumulated GDD (+171 in Acireale, +193 in Giarre, and +194 in Mascali, respectively);
  • lower minimum and maximum temperatures in March and lower rainfall (by about half) throughout the growing season across the three locations in Season IV.
However, the variations in weather conditions between the seasons were more pronounced in Acireale, where there was a greater variation in the accumulated GDD throughout the growing seasons, essentially due to the greater variation in the average minimum temperatures.

2.6. Statistical Analysis

The quantitative traits under study were statistically analyzed following a three-stage procedure consisting of an ANOVA followed by two multivariate models, i.e., AMMI and GGE biplot analysis, both based on the singular value decomposition (SVD). The ANOVA was adopted to estimate the variations among genotypes, environments, and GEI using CoStat® version 6.003 (CoHort Software, Monterey, CA, USA). For the purpose of this study and to simplify the estimation of variance components, location × year contributions were considered to be environments, consistent with Affleck et al. [26] and Gedif and Yigzaw [27]. The ANOVA was firstly computed for individual environments to check the genotype effect in each environment and then combined to check the GEI, considering both the environments and genotypes as fixed factors. One-way ANOVAs were also applied to pooled data of both environments and genotypes to test the significance of each variable. Before ANOVA and in accordance with Lombardo et al. [35], percentage data (UT and CT) and NMTP data were Bliss- and square-root transformed, respectively, to obtain normalized distribution. Homogeneity of variance and normality were respectively verified with the Bartlett and Shapiro–Wilk tests. Differences between means were assessed using Fisher’s protected least significant difference (LSD) test at the 95% probability level.
Once the significance of the GEI was assessed, it was further elucidated through stability analysis to graphically generate biplots. In the AMMI model, combining additive and multiplicative parameters into a single analysis, the GEI is divided into interaction principal component (IPC) scores and residual values according to the equation described by Zobel et al. [36]. In detail, the AMMI firstly involves an ANOVA to analyze the genotype and environment main effects without the interaction and, then, a principal component analysis to standardized residuals that include the experimental error and the effect of the GEI. The GGE model absorbs the main effects of genotype plus the GEI into the multiplicative component and uses the SVD of the first two principal components (PCs) [18,37]. The GGE biplot is based on the sites regression (SREG) linear-bilinear model indicated by Crossa et al. [38]. Model equations of AMMI and GGE are reported in Scavo et al. [39]. Moreover, following Scavo et al. [39], the results of the AMMI and GGE analyses were graphically displayed via five biplot patterns with scaling = 0, centering = 2, and environment-focused partitioning (SVP = 2): AMMI1 (IPC1 scores vs. trait main effect), AMMI2 (IPC1 vs. IPC2), “which-won-where” (mega-environment analysis), “ranking genotypes” (genotype evaluation), and “ranking environments” (environment evaluation). The PBTools version 1.4 (IRRI, Los Baños, Laguna, Philippines) software was used to perform AMMI and GGE analysis.

3. Results and Discussion

With the aim of identifying productive and stable potato genotypes under diverse environmental conditions in Sicily, in this work, a multi-environmental trial was performed considering five potato genotypes, three locations, and four growing seasons. Overall, 12 environments resulting from the combination of three locations × four growing seasons were considered. In particular, a three-stage approach was employed. It started with a combined ANOVA on pooled data to estimate the variance contribution of genotype and environment and to test the significance of GEI. Then, stability analysis was conducted by AMMI and GGE biplot analysis to elucidate and partition the GEI, graphically observe the interrelationships among tested genotypes and environments, recognize mega-environments (if any) within the studied area, and rank the genotypes and environments for each yield trait.

3.1. Combined ANOVA

From the ANOVA emerged that the main effects of genotype and environment were a significant source of variation (p ≤ 0.001) for all productive traits (Table 2). MY, NMTP, AMTW, and CT 40–60 mm showed the most considerable extent of variation due to environment, which respectively accounted for 68.8%, 85.9%, 54.9%, and 70.2% of the total variation (G + E + GEI + blocks). This is not surprising since quantitative traits are well known to be closely dependent on environmental conditions [20,21,40]. On the contrary, genotype was the largest source of variation for TSG (77.9%), CT < 40 mm (47.3%), and CT > 60 mm (60.8%) (Table 2). Potato’s TSG, in fact, is closely correlated to dry matter content [34], which in turn is highly genotype-dependent. Partitioning the environment into predictable (location) and unpredictable (year) variance components (Table S1), the three-way ANOVA (L × Y × G) showed that all interaction effects were significant at the 99% probability level, in accordance with Movahedi et al. [41], and that year was the most important source of variation for productive traits, except for UT. This indicated the higher influence of weather conditions than soil characteristics on potato yield. An opposite trend was found by Gedif and Yigzaw [27] in mid-high-altitude areas of Ethiopia. Specifically, in our study, the differences in weather conditions recorded between seasons were higher than those between locations in the same season (Figure 1). Weather influences potato yield and productive characteristics throughout all the phenological stages. For instance, tuber yield integrates development, total growth, and partitioning processes, and each of these processes is differently affected by temperature [42].
Regardless of environments, one-way ANOVAs indicated that Generosa showed the highest MY (46.5 t ha−1), followed by Arizona (44.5 t ha−1) (Table 2). For Arizona, MY can be attributed to NMTP (on average 7.4), while for Generosa it was mainly due to AMTW (137.9 g). This finding agreed with Lombardo et al. [34], who also reported high Arizona yields due to NMTP in a similar environment. On the contrary, Levante showed the lowest MY (35.7 t ha−1), due to the lowest AMTW (99.5 g) mainly caused by the highest percentage of small caliber tubers (18.2% of <40 mm tubers) and by the greatest UT (2.3%).
Pooling over genotypes (Table 2), E4 had the greatest MY (54.7 t ha−1), followed by E5 (49.2 t ha−1), which showed also the highest TSG (1059.8 kg m−3) and CT > 60 mm (36.5%). E10, despite the greatest AMTW (185.2 g) and a high percentage of 40–60 mm and >60 mm CTs, showed a very low MY (34.0 t ha−1) due to the lowest NMTP (3.6), like E2 and E6 which both had low MYs. The highest NMTP was found in E1 (11.0), which at the same time showed the lowest AMTW (77.8 g) and the greatest CT < 40 mm percentage (24.8%). The poor NMTP of E2, E6, and E10 (i.e., Season II) may be attributable to the very low rainfall in January and February 2020 (also considering that irrigation was performed from March), which coincides with the phase of tuber differentiation. Indeed, NMTP is highly dependent on environmental conditions and water availability, among others, which play a key role during the critical period around tuberization [43]. According to the previous study on tuber quality characteristics [39], the variance magnitudes of the environment explained the differentiated performances of genotypes between environments, indicating a higher influence of seasonal climatic conditions than location. Although it had a low contribution to total variation, as also reported by Khan et al. [23], the GEI was highly significant (p ≤ 0.001) for MY and yield components, so further analysis was justified.

3.2. Stability Analysis

AMMI and GGE analyses divided the GEI sum of squares into IPCs (Table S2) and PCs (Table S3) to unfold the variance components. Both multivariate models were significant for potato yield and yield-related traits, except for CT < 40 mm and CT > 60 mm (Table 3), which therefore were not further analyzed. About AMMI, the cumulative contribution of IPC1 and IPC2 captured 70.2% of the GEI sum of squares for MY, 82.4% for UT, 80.4% for NMTP, 72.2% for AMTW, 75.1% for TSG, and 73.7% for CT 40–60 mm. Similarly, the first two PCs explained 67.5% of the GEI sum of squares for MY, 76.8% for UT, 77.7% for NMTP, 78.6% for AMTW, 87.2% for TSG, and 69.9% for CT 40–60 mm in the GGE model. These percentages demonstrated the validity of the AMMI and GGE models in representing the GEI in the present study. In fact, Thangavel et al. [44] asserted that an IPCA1 + IPCA2 cumulative variation > 66% is a reliable interpretation of the GEI for the AMMI models, while according to Yang et al. [45] a PC1 + PC2 variation > 60% is desirable for GGE. GEI interactions were graphically represented through five biplot patterns.

3.2.1. AMMI Analysis

According to the FR test [46], the model better describing the GEI in this study is AMMI3, except for UT, since the IPC3 was the last component with significant residue (Table 3). However, the AMMI2 model was selected to explain the GEI effect because it always provided a total variation greater than 70%. Through the AMMI1 biplot shown in Figure 2, it is possible to identify superior genotypes (i.e., those with IPC1 scores nearing zero and above-average mean values) across the environments. At the same time, environments with IPC1 scores close to zero have a low interaction effect, thus ensuring good performances for all genotypes (Table S2). Overall, results obtained in the AMMI1 biplots are in line with the combined ANOVA, and highlighted G1 (Arizona) and G2 (Generosa) as the superior genotypes for potato production throughout the studied environments. In particular, G2 showed the highest MY and AMTW, although it was very unstable, as well as high CT 40–60 mm and TSG, and low UT. G1 was preferred for MY because it showed above-average values and an IPC1 score near zero, in addition to high and stable AMTW and low UT. G3 (Levante) was the superior genotype for TSG and NMTP. It is a new cultivar recently introduced on the market, so no scientific literature is available. However, the Romagnoli F.lli Spa (https://www.romagnolipatate.it; Bologna, Italy; accessed on 29 October 2021) reports similar values of dry matter. Except for AMTW, G5 (Vogue) had below-average main effects and was the most unstable; hence, it was not suitable in the studied area. In addition, G5 and G3 showed the highest UT, even though they were highly influenced by the environment (large IPC1 scores).
Regarding environments, E4 and E5 were the most favorable ones for MY, allowing both high performances and low interactive effects, while the lowest MY occurred in Season II (E2, E6, and E10). Seasons I and II (E1, E2, E5, E6, and E9) allowed also the highest TSG values, with E5 showing the greatest stability. Except for 2022 (E4), Acireale had above-average UT, whereas E9 and E10 had the lowest UT for all genotypes. E10 and E1 showed, respectively, the highest and the lowest performance for AMTW. Moreover, except for Season IV (E8), Giarre allowed a high AMTW, likely due to a twofold reason: one is that, compared to Acireale, during April–May (the period of tuber growth) Giarre experienced higher minima and maxima air temperatures which are known to favor (to some extent) the development of tubers; the other reason may be related to the higher soil total N, available P2O5 and exchangeable K2O of Giarre than Mascali.
The AMMI2 biplot (Figure 3) was more accurate to extract the GEI variation since it captured a higher portion of variance than AMMI1. This statement is corroborated by the reports of Shahriari et al. [47] and Tena et al. [21]. The vertical and horizontal lines that passed through the origin (0, 0) divide the biplot into four sectors. Similar to the “which-won-where” GGE biplot, genotypes far from the origin are considered winning genotypes in the environments that fall in that sector [23]. In this study, G1 wins in E1, E2, and E9 for MY, mainly due to AMTW; G2 wins in E12 for MY, AMTW, and TSG; G3 wins in E2, E4, E9, and E11 for TSG; the same genotype wins in E2, E5, and E6 for CT 40–60 mm; G5 wins in E6, E7, E8, and E10 for MY. AMMI analysis is increasingly used for the evaluation of genotypes in multi-environmental trials, often reporting specific adaptations as in the present study [19,23]. The distance from the biplot origin controls the extent of GEI. Except for G3 and G4 in MY, G4 and G5 in AMTW, and G3 in TSG, all genotypes were far from the origin, hence they are more sensitive to the GEI and therefore specifically adapted to environments [48]. Considering that the three locations under study are close together, the fact that G3 and G4 show lower productive performances than G1 and G2 but high stability for all yield-related traits means that they could be introduced in the studied area. Neither for MY nor yield-related traits, genotypes clustered together, meaning that they did not behave similarly across the environments. Furthermore, except for UT and NMTP, E12 had a long vector; this indicated that it elicited strong interactive forces. On the contrary, environments with IPC1 scores near zero had small interaction effects, which indicated that all genotypes performed well in them, thus indicating they are favorable environments for all genotypes under study [20]. From the angle between the vectors of environment and genotype is possibly highlighting the type and degree of correlation, keeping in mind that it always points in the direction of a positive interaction. In this study, G2 positively correlated to E12 for all traits unless UT and CT 40–60 mm, G5 positively correlated to E8 unless UT, TSG, and CT 40–60 mm, and G4 negatively correlated to E9 unless NMTP and CT 40–60 mm (Figure 3). Regarding correlations between environments, E12 positively correlated to E4 and E11 for MY, UT, AMTW, and only E4 for NMTP, since they formed an acute angle. In the cases of obtuse and right angles, there is a negative correlation and no correlation, respectively [23,40]. Here, no clear correlations between environments and no environment clusters were observed.

3.2.2. GGE Analysis

The results from GGE were similar to those obtained from AMMI and combined ANOVA. A “which-won-where” biplot was constructed for each trait to explore the possible existence of mega-environments within the studied region and to identify winning genotypes in each mega-environment (Figure 4). In this biplot, a polygon is drawn joining the genotypes that are located farthest from the origin, containing all other genotypes within the polygon area. Each side of the polygon is crossed by a perpendicular line that divides the biplot into sectors. Vertex genotypes are the best performers in the environments encompassed within the sector (mega-environment) and, oppositely, vertex genotypes placed in a sector where there are no environments are treated as poor genotypes for all tested environments [21,40]. Here, in line with the AMMI2 biplot, no clear differentiation of mega-environments with similar characteristics was found, but results were closely dependent on the specific yield trait. For MY, three mega-environments have been identified: one encompassing E4, E5, E11, and E12 with G2 as the winning genotype, another encompassing E1, E2, E3, E9, and E10 with G1 as the winning genotype, and the last mega-environment comprising E7 and E8 with G5 at the vertex. G1 was also the winning genotype for NMTP in a mega-environment encompassing E1, E3, E9, and E11, as well as for AMTW in a mega-environment including E1, E2, and E9. The fact that G1 was the winning genotype for yield and other yield-related traits in Acireale in three out of four seasons (E1, E2, E3) may indicate that Arizona is able to optimally use soil fertility, given the higher soil nutrient element and organic matter levels of Acireale, thus suggesting its introduction in zones with high soil inputs. About AMTW, most environments (E4, E5, E6, E7, E8, E10, and E12) were encompassed in a mega-environment with G2 as the winner. At the same time, G1 and G2 were at the vertex in a sector with only E8 for UT; therefore, they are genotypes with low UT in the studied region. For TSG, excluding E7, all environments formed a mega-environment with G3 positioned at the vertex, in accordance with what was observed in AMMI2.
Figure 5 shows the “ranking genotypes” biplots, which compare the genotypes to an ideal genotype located at the center of the concentric circles. An ideal genotype has high performance coupled with high stability across environments [18,37]. In accordance with Rakshit et al. [22], the biplot was created by environment-focused partitioning of data (SPV = 2) without scaling to study the interrelationships among test environments. The average environment coordination (AEC) method was employed, in which an average environment axis (AEA) passes through the biplot origin and the average environment, represented on the AEA at the center of concentric circles [37]. Genotypes on the right side of AEC abscissa are above-average performers, whereas those on the left side underperform the overall mean. Moreover, the lower the length of the projection of a genotype, the more stable it is. From this study, it can be stated that G2 and G1 were close to the ideal genotype for productive components, despite a low stability of G2 for MY (Figure 5). They are also the worst genotypes for UT, which is an undesirable trait of potatoes. Cadersa et al. [29], Gedif and Yigzaw [27], and Muthoni et al. [28] reported similar findings for potato yield in Mauritius, Ethiopia, and Kenya, respectively. Summarizing, G1 (Arizona) is a variety specially bred for the South European and North African markets, and it is no coincidence that it can be noted as the best leading genotype in the studied area by virtue of its high productivity, as a consequence of both yield components and high stability. G2 (Generosa), despite the highest MY due to AMTW, can be considered less suitable than G1 because of the lower stability than G1. In addition, G3 was confirmed as the superior genotype for TSG; this, coupled with an MY stability (despite not being very high), suggests its introduction into the studied area for protected geographical identifications. Indeed, the choice of superior genotypes for protected geographical identifications needs the genotypes’ assessment over diverse environments in a specific region.
Like in Figure 5, Figure 6 shows the “ranking environments” biplots, which help in evaluating the environments with respect to an ideal environment. Similarly to the ideal genotype, the ideal environment, a crucial element in multi-environmental trials, has the highest discriminativeness (power of genotype discrimination, visualized through the length of environment vectors, i.e., high PC1 scores) and the lowest representativeness (representativeness of all other environments, visualized through a small angle with the AEA, i.e., zero PC2 scores) [18]. Hence, a desirable environment should be near the ideal environment, thus having both high PC1 and low PC2 scores (Table S3), graphically represented by long vectors that form a short angle with the AEA (Figure 6). In this study, E4 for MY, E5 for UT and TSG, E1 for NMTP, E10 (followed by E4) for AMTW, and E3 for CT 40–60 mm were recorded as the most suitable environments for the studied potato genotypes in Sicily. It could be noted that the majority of ideal environments in terms of yield and yield components belong to Acireale, likely due to higher soil fertility (Table 1). A similar approach was used by Rakshit et al. [22] in sorghum by Khan et al. [23] in Bambara groundnut and by Gedif and Yigzaw [27] for potatoes in Ethiopia.

4. Conclusions

Once a new genotype is proposed for cultivation, its GEI with a specific area should be assessed. For this reason, the present multi-environmental trial was carried out to evaluate the productive performances of five recently introduced potato genotypes over four years and three locations in a Mediterranean climate. The GEI of all productive traits was significant, with the environment (and in particular the year) explaining the highest amount of variation, except for TSG and CT < 40 mm and >60 mm. AMMI and GGE biplot analyses recognized Arizona and Generosa as the superior genotypes in the studied region by virtue of their above-average mean performances and stability throughout the environments. In particular, Arizona and Generosa were the winning genotypes in most environments for MY, NMTP, and AMTW, whereas Levante was recognized as the leading genotype for TSG. The environments belonging to Acireale were the closest to the ideal environment, likely due to its higher soil fertility, while Season II (E2, E6, and E10) showed unfavorable climatic conditions, as demonstrated by the low NMTP and high UT. In conclusion, both AMMI and GGE proved to be effective tools for local potato breeding, the selection of protected geographical identifications, and the commercial introduction of high-yielding and stable potato cultivars across the Mediterranean climate.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/agronomy13010101/s1: Table S1: F values of main effects and their interactions resulting from three-way analysis of variance (ANOVA) of potato yield and yield-related traits; Table S2: IPC1 and IPC2 scores derived from AMMI analysis of the five potato genotypes and 12 environments under study; Table S3: PC1 and PC2 scores derived from GGE analysis of the five potato genotypes and 12 environments under study.

Author Contributions

Conceptualization, G.M. and A.I.; methodology, A.S. and A.I.; software, data curation and formal analysis, A.S.; investigation, A.S., G.M. and A.I.; writing—original draft preparation, A.S. and A.I.; writing—review and editing, A.S., G.M. and A.I.; supervision, G.M. and A.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data are available via email request to the corresponding author.

Acknowledgments

The authors thank S. La Rosa and A. Pellegrino (CNR, IBE) for their excellent technical assistance. We are also grateful to Giovanni Continella for logistically organizing field experiments and Romagnoli F.lli Spa for providing the tuber ‘seed’ material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO Statistical Database. Production. Crops. Rome, Italy. 2020. Available online: http://faostat.fao.org (accessed on 30 October 2021).
  2. Ierna, A.; Mauromicale, G. How irrigation water saving strategy can affect tuber growth and nutritional composition of potato. Sci. Hortic. 2022, 299, 111034. [Google Scholar] [CrossRef]
  3. Lombardo, S.; Pandino, G.; Mauromicale, G. The influence of growing environment on the antioxidant and mineral content of “early” crop potato. J. Food Compos. Anal. 2013, 32, 28–35. [Google Scholar] [CrossRef]
  4. Lombardo, S.; Scavo, A.; Abbate, C.; Pandino, G.; Parisi, B.; Mauromicale, G. Mycorrhizal inoculation improves mineral content of organic potatoes grown under calcareous soil. Agriculture 2021, 11, 333. [Google Scholar] [CrossRef]
  5. Ierna, A.; Parisi, B. Crop growth and tuber yield of “early” potato crop under organic and conventional farming. Sci. Hortic. 2014, 165, 260–265. [Google Scholar] [CrossRef]
  6. Foti, S.; Mauromicale, G.; Ierna, A. Influence of irrigation regimes on growth and yield of potato cv. Spunta. Potato Res. 1995, 38, 307–317. [Google Scholar] [CrossRef]
  7. Ierna, A.; Mauromicale, G. Sustainable and profitable nitrogen fertilization management of potato. Agronomy 2019, 9, 582. [Google Scholar] [CrossRef] [Green Version]
  8. Ierna, A.; Mauromicale, G. How moderate water stress can affect water use efficiency indices in potato. Agronomy 2020, 10, 1034. [Google Scholar] [CrossRef]
  9. Ierna, A.; Parisi, B.; Melilli, M.G. Overall quality of “early” potato tubers as affected by organic cultivation. Agronomy 2022, 12, 296. [Google Scholar] [CrossRef]
  10. Lombardo, S.; Pandino, G.; Mauromicae, G. The mineral profile in organically and conventionally grown “early” crop potato tubers. Sci. Hortic. 2014, 167, 169–173. [Google Scholar] [CrossRef]
  11. Samaniego, I.; Espin, S.; Cuesta, X.; Arias, V.; Rubio, A.; Llerena, W.; Angós, I.; Carrillo, W. Analysis of environmental conditions effect in the phytochemical composition of potato (Solanum tuberosum) cultivars. Plants 2020, 9, 815. [Google Scholar] [CrossRef]
  12. Pour-Aboughadareh, A.; Khalili, M.; Poczai, P.; Olivoto, T. Stability Indices to Deciphering the Genotype-by-Environment Interaction (GEI) Effect: An Applicable Review for Use in Plant Breeding Programs. Plants 2022, 11, 414. [Google Scholar] [CrossRef] [PubMed]
  13. Suzuki, D.T.; Griffiths, A.J. An Introduction to Genetic Analysis, 7th ed.; WH Freeman and Company: New York, NY, USA, 2000. [Google Scholar]
  14. Gauch, H.G. Statistical analysis of yield trials by AMMI and GGE. Crop Sci. 2006, 46, 1488–1500. [Google Scholar] [CrossRef]
  15. Bradshaw, J.E. Genetics of Agrihorticultural Traits. In Handbook of Potato Production, Improvement and Postharvest Management; Gopal, J., Khurana, S.M.P., Eds.; Food Products Press: Binghampton, NY, USA, 2006; pp. 41–75. [Google Scholar]
  16. Yan, W.; Kang, M.S.; Ma, B.; Woods, S.; Cornelius, P.L. GGE biplot vs. AMMI analysis of genotype-by-genotype environment data. Crop Sci. 2007, 47, 643–655. [Google Scholar] [CrossRef] [Green Version]
  17. Flis, B.; Domański, L.; Zimnoch-Guzowska, E.; Polgar, Z.; Pousa, S.Á.; Pawlak, A. Stability Analysis of Agronomic Traits in Potato Cultivars of Different Origin. Am. J. Potato Res. 2014, 91, 404–413. [Google Scholar] [CrossRef] [Green Version]
  18. Yan, W.; Hunt, L.A.; Sheng, Q.; Szlavnics, Z. Cultivar evaluation and mega-environment investigation based on the GGE biplot. Crop Sci. 2000, 40, 597–605. [Google Scholar] [CrossRef]
  19. Katsenios, N.; Sparangis, P.; Leonidakis, D.; Katsaros, G.; Kakabouki, I.; Vlachakis, D.; Efthimiadou, A. Effect of Genotype × Environment Interaction on Yield of Maize Hybrids in Greece Using AMMI Analysis. Agronomy 2021, 11, 479. [Google Scholar] [CrossRef]
  20. Senguttuvel, P.; Sravanraju, N.; Jaldhani, V.; Divya, B.; Beulah, P.; Nagaraju, P.; Manasa, Y.; Prasad, A.S.; Brajendra, P.; Gireesh, C.; et al. Evaluation of genotype by environment interaction and adaptability in lowland irrigated rice hybrids for grain yield under high temperature. Sci. Rep. 2021, 11, 15825. [Google Scholar] [CrossRef]
  21. Tena, E.; Goshu, F.; Mohamad, H.; Tesfa, M.; Tesfaye, D.; Seife, A. Genotype × environment interaction by AMMI and GGE-biplot analysis for sugar yield in three crop cycles of sugarcane (Saccharum officinirum L.) clones in Ethiopia. Cogent. Food. Agric. 2019, 5, 1651925. [Google Scholar] [CrossRef]
  22. Rakshit, S.; Ganapathy, K.N.; Gomashe, S.S.; Rathore, A.; Ghorade, R.B.; Kumar, M.V.; Ganesmurthy, K.; Jain, S.K.; Kamtar, M.Y.; Sachan, J.S.; et al. GGE biplot analysis to evaluate genotype, environment and their interactions in sorghum multi-location data. Euphytica 2012, 185, 465–479. [Google Scholar] [CrossRef]
  23. Khan, M.M.H.; Rafii, M.Y.; Ramlee, S.I.; Jusoh, M.; Al Mamun, M. AMMI and GGE biplot analysis for yield performance and stability assessment of selected Bambara groundnut (Vigna subterranea L. Verdc.) genotypes under the multi-environmental trials (METs). Sci. Rep. 2021, 11, 22791. [Google Scholar] [CrossRef]
  24. Tirado, R.M.; Tirado, R.L.; Mendoza, J.C. Genotype × environment interaction in potato (Solanum tuberosum L.) yield with pigmented flesh in Cutervo, Peru. Chilean J. Agric. Anim. Sci. 2018, 34, 191–198. [Google Scholar]
  25. Lenartowicz, T.; Piepho, H.P.; Przystalski, M. Stability analysis of tuber yield and starch yield in mid-late and late maturing starch cultivars of potato (Solanum tuberosum). Potato Res. 2020, 63, 179–197. [Google Scholar] [CrossRef]
  26. Affleck, I.; Sullivan, J.A.; Tarn, R.; Falk, D.E. Genotype by environment interaction effect on yield and quality of potatoes. Can. J. Plant Sci. 2008, 88, 1099–1107. [Google Scholar] [CrossRef]
  27. Gedif, M.; Yigzaw, D. Genotype by environment interaction analysis for tuber yield of potato (solanum tuberosum L.) using a GGE biplot method in Amhara Region, Ethiopia. Agric. Sci. 2014, 5, 239–249. [Google Scholar] [CrossRef] [Green Version]
  28. Muthoni, J.; Shimelis, H.; Melis, R. Genotype x Environment Interaction and Stability of Potato Tuber Yield and Bacterial Wilt Resistance in Kenya. Am. J. Potato Res. 2015, 92, 367–378. [Google Scholar] [CrossRef]
  29. Cadersa, Y.; Santchurn, D.; Govinden Soulange, J.; Saumtally, S.; Parmessur, Y. Genotype- by-environment interaction for marketable tuber yield in advanced potato clones using AMMI and GGE methods. Afr. Crop Sci. J. 2022, 30, 331–346. [Google Scholar] [CrossRef]
  30. Sood, S.; Bhardwaj, V.; Kumar, V. BLUP and stability analysis of multi-environment trials of potato varieties in sub-tropical Indian conditions. Heliyon 2020, 6, e05525. [Google Scholar] [CrossRef]
  31. Soil Survey Staff. Soil Taxonomy: A Basic System of Soil Classification for Making and Interpreting Soil Surveys, 2nd ed.; U.S. Government Publishing Office: Washington, DC, USA, 1999.
  32. Italian Society of Soil Science. Metodi Normalizzati Di Analisi Del Suolo; Edagricole: Bologna, Italy, 1985. [Google Scholar]
  33. van Berloo, R.; Hutten, R.C.B.; van Eck, H.J.; Visser, R.G.F. An Online Potato Pedigree Database Resource. Potato Res. 2007, 50, 45–57. [Google Scholar] [CrossRef] [Green Version]
  34. Schippers, P.A. The relationship between specific gravity and percentage dry matter in potato tubers. Am. Potato J. 1976, 53, 111–122. [Google Scholar] [CrossRef]
  35. Lombardo, S.; Abbate, C.; Pandino, G.; Parisi, B.; Scavo, A.; Mauromicale, G. Productive and physiological response of organic potato grown under highly calcareous soils to fertilization and mycorrhization management. Agronomy 2020, 10, 1200. [Google Scholar] [CrossRef]
  36. Zobel, R.W.; Wright, M.J.; Gauch Jr, H.G. Statistical analysis of a yield trial. Agron. J. 1988, 80, 388–393. [Google Scholar] [CrossRef]
  37. Yan, W.; Kang, M.S. GGE Biplot Analysis: A Graphical Tool for Breeders, Geneticists, and Agronomists; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  38. Crossa, J.; Cornelius, P.L.; Yan, W. Biplots of linear-bilinear models for studying crossover genotype × environment interaction. Crop Sci. 2002, 42, 619–633. [Google Scholar] [CrossRef]
  39. Scavo, A.; Mauromicale, G.; Ierna, A. Genotype × environment interactions of potato tuber quality characteristics by AMMI and GGE biplot analysis. Sci. Hortic. 2023, 310, 111750. [Google Scholar] [CrossRef]
  40. Balakrishnan, D.; Subrahmanyam, D.; Badri, J.; Raju, A.K.; Rao, Y.V.; Beerelli, K.; Mesapogu, S.; Surapaneni, M.; Ponnuswamy, R.; Padmavathi, G.; et al. Genotype × environment interactions of yield traits in backcross introgression lines derived from Oryza sativa cv. Swarna/Oryza nivara. Front. Plant Sci. 2016, 7, 1530. [Google Scholar] [CrossRef] [Green Version]
  41. Movahedi, H.; Mostafavi, K.; Shams, M.; Golparvar, A.R. AMMI analysis of genotype × environment interaction on grain yield of sesame (Sesamum indicum L.) genotypes in Iran. Biotechnol. Biotechnol. Equip. 2020, 34, 1013–1018. [Google Scholar] [CrossRef]
  42. Lizana, C.X.; Sandaña, P.; Behn, A.; Ávila-Valdés, A.; Ramírez, D.A.; Soratto, R.P.; Campos, H. Potato. In Crop Physiology: Case Histories for Major Crops; Sadras, V.O., Calderini, D.F., Eds.; Academic Press: Cambridge, MA, USA, 2021; Volume 18, pp. 551–587. [Google Scholar] [CrossRef]
  43. Celis-Gamboa, C.; Struik, P.C.; Jacobsen, E.; Visser, R.G.F. Temporal dynamics of tuber formation and related processes in a crossing population of potato (Solanum tuberosum). Ann. Appl. Biol. 2003, 143, 175–186. [Google Scholar] [CrossRef]
  44. Thangavel, P.; Anandan, A.; Eswaran, R. AMMI analysis to comprehend genotype-by-environment (G × E) interactions in rainfed grown mungbean (Vigna radiata L.). Aust. J. Crop Sci. 2011, 5, 1767–1775. [Google Scholar]
  45. Yang, R.C.; Crossa, J.; Cornelius, P.L.; Burgueño, J. Biplot analysis of genotype × environment interaction: Proceed with caution. Crop Sci. 2009, 49, 1564–1576. [Google Scholar] [CrossRef] [Green Version]
  46. Piepho, H.P. Robustness of statistical tests for multiplicative terms in the additive main effects and multiplicative interaction model for cultivar trials. Theor. Appl. Genet. 1995, 90, 438–443. [Google Scholar] [CrossRef] [PubMed]
  47. Shahriari, Z.; Heidari, B.; Dadkhodaie, A. Dissection of genotype × environment interactions for mucilage and seed yield in Plantago species: Application of AMMI and GGE biplot analyses. PLoS ONE 2018, 13, e0196095. [Google Scholar] [CrossRef] [Green Version]
  48. Koundinya, A.V.V.; Pandit, M.K.; Ramesh, D.; Mishra, P. Phenotypic stability of eggplant for yield and quality through AMMI, GGE and cluster analyses. Sci. Hortic. 2019, 247, 216–223. [Google Scholar] [CrossRef]
Figure 1. Monthly rainfall, average minima and maxima monthly air temperatures, and 10-year (2013–2022) climatic trend in Acireale (a), Giarre (b), and Mascali (c) throughout Season I (2018–2019), Season II (2019–2020), Season III (2020–2021), and Season IV (2021–2022).
Figure 1. Monthly rainfall, average minima and maxima monthly air temperatures, and 10-year (2013–2022) climatic trend in Acireale (a), Giarre (b), and Mascali (c) throughout Season I (2018–2019), Season II (2019–2020), Season III (2020–2021), and Season IV (2021–2022).
Agronomy 13 00101 g001aAgronomy 13 00101 g001b
Figure 2. ‘AMMI 1’ biplot showing the mean trait effect and the first interaction principal component (IPC1) of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Figure 2. ‘AMMI 1’ biplot showing the mean trait effect and the first interaction principal component (IPC1) of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Agronomy 13 00101 g002
Figure 3. ‘AMMI 2’ biplot showing the first two interaction principal components (IPC1 and IPC2) effects of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Figure 3. ‘AMMI 2’ biplot showing the first two interaction principal components (IPC1 and IPC2) effects of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Agronomy 13 00101 g003
Figure 4. Polygon view of the GGE biplot illustrating the “which-won-where” analysis of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Figure 4. Polygon view of the GGE biplot illustrating the “which-won-where” analysis of five potato genotypes (G) over 12 environments (E). The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Agronomy 13 00101 g004
Figure 5. ‘Ranking genotypes’ GGE biplots illustrating the performance of five potato genotypes (G) over 12 environments (E) in comparison to an ideal genotype. The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Figure 5. ‘Ranking genotypes’ GGE biplots illustrating the performance of five potato genotypes (G) over 12 environments (E) in comparison to an ideal genotype. The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Agronomy 13 00101 g005
Figure 6. ‘Ranking environments’ GGE biplots showing the performance of 12 environments (E) and five potato genotypes (G) in comparison to an ideal environment. The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Figure 6. ‘Ranking environments’ GGE biplots showing the performance of 12 environments (E) and five potato genotypes (G) in comparison to an ideal environment. The biplots were constructed with scaling = 0, centering = 2, and SVP = 2. G1: Arizona; G2: Generosa; G3: Levante; G4: Paradiso; G5: Vogue; E1–E4: Acireale (seasons I to IV); E5–E8: Giarre (seasons I to IV); E9–E12: Mascali (seasons I to IV).
Agronomy 13 00101 g006
Table 1. Description of soil characteristics (−30 cm depth) of the three locations under study.
Table 1. Description of soil characteristics (−30 cm depth) of the three locations under study.
Soil CharacteristicUnit of MeasurementAcirealeGiarreMascali
Sand%828590
Silt%14117
Clay%443
Organic matter%4.92.52.4
C/N/9.58.910.7
Total N%0.300.160.13
P2O5 availableppm428236149
K2O exchangeableppm1193723672
pH/7.37.47.4
Electrical conductivitymS cm−11.962.260.60
Cation exchange capacitymeq 100 g−120.212.511.9
Ca exchangeablemeq 100 g−113.07.28.0
Mg exchangeablemeq 100 g−13.81.92.2
Na exchangeablemeq 100 g−10.81.90.3
K exchangeablemeq 100 g−12.51.51.4
Znppm5786
Cuppm30819
Feppm10311078
Mnppm301110
Table 2. Joint analysis of variance (ANOVA, F-values) of five potato genotypes (G) and 12 environments (E) for marketable yield (MY, t ha−1), unmarketable tubers (UT, %), number of marketable tubers plant−1 (NMTP), average marketable tuber weight (AMTW, g), tuber specific gravity (TSG, kg m−3), and caliber tubers (CT, %).
Table 2. Joint analysis of variance (ANOVA, F-values) of five potato genotypes (G) and 12 environments (E) for marketable yield (MY, t ha−1), unmarketable tubers (UT, %), number of marketable tubers plant−1 (NMTP), average marketable tuber weight (AMTW, g), tuber specific gravity (TSG, kg m−3), and caliber tubers (CT, %).
Source of VariationMYNMTPAMTWUTTSGCT < 40 mmCT 40–60 mmCT > 60 mm
GArizona44.5 (2.0) b7.4 (0.5) a122.5 (5.1) c1.3 (0.3) b1048.2 (1.0) c9.3 (1.5) bc70.9 (2.2) bc19.8 (1.9) bc
Generosa46.5 (2.4) a6.9 (0.4) bc137.9 (6.7) a1.034 (0.2) b1055.7 (1.1) b9.8 (1.5) b72.6 (2.4) b17.6 (2.5) c
Levante35.7 (1.6) d7.3 (0.5) ab99.5 (3.9) e2.3 (0.84) a1062.2 (1.1) a18.2 (1.7) a76.1 (1.7) a5.7 (1.3) d
Paradiso41.5 (1.8) c6.5 (0.4) c128.4 (5.5) b1.5 (0.3) b1049.3 (1.2) c6.3 (0.9) d69.8 (1.9) c23.9 (1.9) a
Vogue40.0 (2.1) c6.8 (0.4) bc115.9 (5.7) d2.9 (0.7) a1049.4 (1.1) c8.5 (1.1) cd71.1 (2.6) bc20.4 (2.8) b
EE143.4 (1.7) c11.0 (0.4) a77.8 (4.0) f2.2 (0.6) bc1057.6 (1.5) ab24.8 (3.2) a57.6 (3.1) f17.5 (2.5) cd
E223.7 (1.3) g3.9 (0.2) hi123.6 (7.4) d2.3 (0.6) ab1054.4 (2.1) cde5.2 (1.5) e72.4 (2.7) cd22.4 (3.7) bc
E339.8 (2.9) de7.3 (0.6) ef101.5 (3.3) e4.6 (1.6) a1050.4 (0.9) f9.6 (1.1) d75.2 (1.1) bc15.1 (1.5) d
E454.7 (2.8) a8.1 (0.2) cd127.5 (5.9) cd1.1 (0.3) ef1051.7 (1.8) ef13.3 (1.6) b61.8 (3.5) f24.9 (3.5) b
E549.2 (2.7) b6.9 (0.4) ef134.3 (3.1) bc0.6 (0.21) fg1059.8 (1.3) a3.6 (0.9) e59.7 (2.9) f36.6 (3.1) a
E633.3 (1.7) f4.5 (0.2) h139.8 (6.7) b2.5 (0.6) ab1056.1 (1.9) bc2.6 (1.1) f87.3 (1.5) a10.1 (2.1) e
E738.8 (1.5) e6.2 (0.5) g123.1 (5.8) d1.4 (0.3) cde1050.5 (2.3) f11.9 (1.7) bcd84.8 (1.3) a3.2 (1.2) f
E847.1 (2.9) b8.8 (0.4) bc101.6 (43) e2.0 (0.5) bcd1052.8 (2.0) def11.3 (1.1) cd66.6 (3.5) e22.1 (3.9) bc
E947.1 (2.2) b9.1 (0.3) b98.0 (3.5) e0.6 (0.2) fg1054.7 (2.1) cd14.6 (2.5) b73.5 (1.6) cd11.9 (2.8) e
E1034.0 (1.4) f3.6 (0.2) i185.2 (8.0) a0.5 (0.2) g1052.9 (2.8) def3.3 (1.6) f77.9 (2.5) b18.8 (2.7) bc
E1147.1 (4.1) b7.5 (0.4) de118.6 (8.0) d1.2 (0.3) def1043.3 (2.1) g13.0 (1.7) bc78.1 (1.4) b8.9 (2.0) e
E1242.0 (3.8) cd6.7 (0.3) fg119.1 (10.2) d1.7 (0.5) cde1051.2 (2.5) f11.9 (2.1) bcd69.8 (3.5) de18.2 (5.1) cd
ANOVAG a41.1 ***5.1 ***49.4 ***10.0 ***85.1 ***57.5 ***5.9 ***70.3 ***
E b70.0 ***71.87 ***70.6 ***10.2 ***17.7 ***56.4 ***35.5 ***34.3 ***
G × E c18.2 ***4.1 ***7.5 ***5.5 ***4.7 ***6.3 ***6.5 ***8.4 ***
Blocks d1.0 ns2.6 ns1.1 ns1.0 ns1.8 ns1.1 ns2.7 ns2.7 ns
Values are means with SEM (standard error of the mean, n = 36 for G, n = 15 for E) in brackets. Values within a column with different letters indicate significant differences (p ≤ 0.05) using the LSD test. a 4 d.f.; b 11 d.f.; c 44 d.f.; d 2 d.f.; *** and ns indicate significant effect at p ≤ 0.001 and not significant, respectively.
Table 3. Additive main effects and multiplicative interaction (AMMI) and G + GE (GGE) analyses in potato yield and yield components.
Table 3. Additive main effects and multiplicative interaction (AMMI) and G + GE (GGE) analyses in potato yield and yield components.
Marketable YieldN. of Marketable Tubers Plant−1Average Marketable Tuber WeightUnmarketable Tubers
SSMSSS% of GEISSMSSS% of GEISSMSSS% of GEISSMSSS% of GEI
AMMI analysis
IPC1 a3414.4243.9 ***42.069.14.9 ***57.814471.71033.6 **43.4769.154.9 **62.0
IPC2 b2294.3191.1 ***28.227.12.1 ***22.69612.1801.0 **28.8252.721.1 *20.4
IPC3 c1806.4180.6 ***22.218.11.7 ***15.17052.6705.3 *21.2135.613.6 ns10.9
GGE analysis
PC1 a4177.4298.4 ***42.770.95.1 ***53.127903.21993.1 **52.3798.757.0 **55.3
PC2 b2429.3202.5 ***24.832.82.6 ***24.614039.11169.9 **26.3310.625.9 *21.5
PC3 c2186.0218.5 ***22.324.12.3 ***18.08337.3833.6 *15.6215.521.5 ns14.9
Tuber specific gravityCaliber tubers < 40 mmCaliber tubers 40–60 mmCaliber tubers > 60 mm
SSMSSS% of GEISSMSSS% of GEISSMSSS% of GEISSMSSS% of GEI
AMMI analysis
IPC1 a973.469.5 ***47.01284.291.6 ns51.11435.0102.4 ***47.63606.6257.5 ns52.4
IPC2 b581.248.4 ***28.1623.952.0 ns24.8785.765.5 ***26.11416.4118.0 ns20.6
IPC3 c337.433.6 ***16.3481.948.2 ns19.2620.662.1 ***20.6971.097.0 ns14.1
GGE analysis
PC1 a3820.4272.9 ***69.92948.4210.6 ns63.91490.2106.4 ***45.66388.5456.2 ns53.0
PC2 b947.078.8 ***17.3975.681.3 ns21.2794.266.2 ***24.33353.7279.5 ns27.8
PC3 c443.544.4 ***8.1541.154.0 ns11.7661.166.1 ***20.21415.5141.5 ns11.7
a 14 d.f.; b 12 d.f.; c 10 d.f.; Notes: IPC: interaction principal component; SS: sum of squares; MS: mean square; SS% of GEI: variance proportion of genotype × environment interaction; ***, **, * and ns indicate significant effect at p ≤ 0.001, p ≤ 0.01, p ≤ 0.05 and not significant, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Scavo, A.; Mauromicale, G.; Ierna, A. Dissecting the Genotype × Environment Interaction for Potato Tuber Yield and Components. Agronomy 2023, 13, 101. https://doi.org/10.3390/agronomy13010101

AMA Style

Scavo A, Mauromicale G, Ierna A. Dissecting the Genotype × Environment Interaction for Potato Tuber Yield and Components. Agronomy. 2023; 13(1):101. https://doi.org/10.3390/agronomy13010101

Chicago/Turabian Style

Scavo, Aurelio, Giovanni Mauromicale, and Anita Ierna. 2023. "Dissecting the Genotype × Environment Interaction for Potato Tuber Yield and Components" Agronomy 13, no. 1: 101. https://doi.org/10.3390/agronomy13010101

APA Style

Scavo, A., Mauromicale, G., & Ierna, A. (2023). Dissecting the Genotype × Environment Interaction for Potato Tuber Yield and Components. Agronomy, 13(1), 101. https://doi.org/10.3390/agronomy13010101

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