Next Article in Journal
A Multi-Criteria Goal Programming Model to Analyze the Sustainable Goals of India
Next Article in Special Issue
Parameter Uncertainty Analysis of the Life Cycle Inventory Database: Application to Greenhouse Gas Emissions from Brown Rice Production in IDEA
Previous Article in Journal
Disentangling the Complex Effects of Socioeconomic, Climatic, and Urban Form Factors on Air Pollution: A Case Study of China
Previous Article in Special Issue
Intercropping Halophytes to Mitigate Salinity Stress in Watermelon
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Environment Evaluation and Genetic Characterisation of Common Bean Breeding Lines for Organic Farming Systems

1
Dipartimento di Scienze Agrarie, Alimentari e Ambientali (DSA3), Università degli Studi di Perugia, Borgo XX Giugno 74, 06121 Perugia, Italy
2
Enza Zaden Italia Research, S.S. Aurelia Km 96, 710, Tarquinia, 01016 Viterbo, Italy
3
N8 AgriFood, Department of Biology, University of York, Heslington, York YO10 5DD, UK
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(3), 777; https://0-doi-org.brum.beds.ac.uk/10.3390/su10030777
Submission received: 21 January 2018 / Revised: 2 March 2018 / Accepted: 9 March 2018 / Published: 12 March 2018

Abstract

:
It is recognised that one of the main causes for the relative low yields under organic conditions is the use of modern cultivars which are bred for high-input management systems. The work described here aimed to study and test possible breeding strategies to produce cultivars of common bean for organic agriculture. To this purpose, crosses between a traditional Italian landrace named “Gnocchetto” and a cultivar were carried out. The F1 plants obtained were either backcrossed or self-fertilised and the obtained materials subjected to selection for quality traits at different development stages. The resulting lines were tested under four different environmental conditions for three years in order to determine their potential performance. The resulting data were analysed using a Multi-Environment Trial Analysis (MET) approach and different visualisations of the GGE biplot were generated. Furthermore, to assess the level of genetic similarity, the lines were characterised using 25 Simple Sequence Repeat (SSR) molecular markers. Results showed that the breeding approach applied allowed to select lines with the same technological and agronomic characteristics as commercially available cultivars, but with different adaptation abilities that make them suitable for organic agriculture.

1. Introduction

The demand for organic food is currently increasing worldwide and interest in organic agriculture has grown significantly [1]. Behind this trend, there is the consumer perception that organic products are healthier, safer and more beneficial for the environment and biodiversity when compared to those produced under conventional management systems [2].
In the past three decades, the impact of organic management on crop yields has been widely discussed and contrasting results have been highlighted. While many studies showed that organic management resulted in significantly lower yields than conventional ones [3,4,5,6], others reported comparable yields [7,8,9,10,11].
Recently, many reviews have used comprehensive meta-analyses to investigate the contrasting results mentioned above. Considering data from over 350 published studies, De Ponti et al. [12] concluded that organic yields are on average 20% lower than conventional ones. On the other hand, Seufert and colleagues [13] showed that yield differences between organic and conventional agriculture do exist, but that they are highly contextual, depending on the system, site characteristics and crop types. Moreover, Ponisio et al. [14] built a hierarchical meta-analytic model to better dissect and estimate the yield gap between the two management systems. They found that organic yields are, on average, 19.2% (±3.7%) lower than conventional, while, in contrast to the previous studies, their results demonstrated a lack of significant differences in yield gaps for leguminous versus non-leguminous crops and perennials versus annuals; even more interestingly, results showed that agricultural diversification practices such as rotations or multi-cropping can reduce the yield gap between the two management systems.
It has been suggested that one of the main causes for the low yields achievable under organic conditions is the use of modern cultivars which are bred for high-input production systems [15,16]. Currently, most cultivars are specifically bred for conventional agriculture and rely on the use of pesticides and fertilisers for optimal performance [14,17,18]. Nevertheless, there is also evidence that conventionally bred cultivars do not always perform worse than those bred for organic agriculture when tested under organic conditions; for example, genotype by system interactions were not reported in studies carried out on wheat [19] or oats [20]. In comparative trials of 32 dry bean genotypes under organic versus conventional systems, Heiling and Kelly [21] reported that some genotypes were better suited for organic production than others. However, those genotypes poorly performing under organic conditions performed badly also under conventional conditions, confirming the absence of a genotype by system interaction in the studied materials. Studying soft wheat, Annicchiarico et al. [22] highlighted a high correlation between yield in organic and conventional management systems and, in particular, that genotype by location interactions were much more important than genotype by system interactions. A similar result was also obtained by Raggi et al. [23] in barley. Even if the effect of location is more relevant than the management system in determining yield performances, certain traits that are useful for organic farming are often lacking in modern cultivars [24].
Many consumers prefer organic products with certain physico-chemical properties in accordance with their personal values that go beyond the basic ethical criteria established by EU regulation on organic farming (EC 834/2007). In this context, landraces represent a useful genetic resource as starting material for breeding new cultivars for organic agriculture that meet both consumers’ and growers’ needs [25]. Landraces are variable but identifiable populations that lack “formal” improvement and have been cultivated for hundreds of years in restricted areas by small-scale farmers under marginal conditions which are often analogous to those of organic or low-input management [26,27]. For these reasons, landraces are characterised by (i) specific adaptation to certain environmental conditions and (ii) rare or even unique organoleptic properties. However, despite the above-mentioned beneficial qualities, landraces may also have agronomic disadvantages when compared to modern cultivars such as particular plant architecture or inadequate physiological attributes that make cultivation more difficult and expensive. Therefore, obtaining and testing new breeding lines that combine landraces and commercial cultivar qualities are among the main objectives for many plant breeders [15].
Common bean landraces are populations that maintain a relatively high degree of genetic diversity and are able to adapt to different environmental conditions [26]. Genetic diversity itself provides a buffer against environmental stresses caused by a multitude of factors which are generally referred to as “biotic and abiotic stresses” [15]. In a recent study, Klaedtke and colleagues also showed that the genetic diversity present in some historical common bean cultivars was sufficient for short-term local adaptation to occur in two different environmental conditions under organic management [28].
Common bean (Phaseolus vulgaris L.) is one of the most important grain legume species for direct human consumption, being a rich source of proteins, vitamins, minerals and fibre. This crop is also highly valued for the long storage life of its seeds and for its numerous culinary applications [29]. Due to these characteristics, common bean is considered a model food legume [30]. The total world production of dry beans is estimated at 26 million tons/year, harvested from an area covering about 30 million ha [31]. Currently, detailed information about global organic production is lacking for this crop. However, it was recently estimated that organic certified bean production accounts for 2.4% and 1.5% of total bean production in the U.S. and Europe respectively, and these percentages are steadily rising [32,33].
As a leguminous species, common bean contributes a substantial amount of nitrogen to the soil [34]. Nitrogen is one of the most significant nutrients for plant health, and cost-effective nitrogen management practices are necessary for an efficient production. For this reason, common bean is frequently intercropped or included in rotation schemes in organic agriculture and therefore, breeding and trials for high-performing cultivars specifically bred for organic and low-input agriculture are of great importance for optimising production in these systems.
This paper reports on the yield evaluation of 17 common bean lines bred for organic agriculture that were developed from a cross between a landrace and a commercial cultivar. The lines were compared with cultivars of the same commercial class phenotype in multi-environment field trials (MET) under organic management. In order to compare the level of genetic similarity between the breeding and the parental lines, molecular characterization was also carried out using 18 “neutral” and 7 “gene derived” Simple Sequence Repeat (SSR) molecular markers.

2. Materials and Methods

2.1. The Breeding Programme and Lines

The landrace used in the breeding programme, known as “Gnocchetto” or “Fagiolo a pisello” (Gn), is traditionally cultivated in Colle di Tora (Rieti, Italy, 42°12′31.33′′N, 12°56′49.92′′E, 557 m a.s.l.) on the steep hillsides of Lake Turano in the Apennines, mainly under organic conditions [35]. Gn is a landrace characterised by indeterminate growth habit, late flowering, white flowers and small round white seeds, ascribable to the Navy commercial class. Due to its growth habit and climbing ability, cultivation is difficult and expensive: plants can be more than four metres tall, requiring long supports, and only flower in the upper part. The flowering of this landrace continues until late in the cropping season, shifting the harvesting period by approximately a month compared to the most common Italian late-flowering cultivars. Seeds of the Gn landrace have a thin external tegument which means they cook quickly. Furthermore, beyond the technological advantages of this characteristic, it has also been recently demonstrated that fast-cooking common bean cultivars have higher protein and mineral retention compared to slow-cooking ones, suggesting that they have a greater nutritive value [36]. All these characteristics together with its smooth doughy taste [26] make this landrace an interesting genetic resource for a breeding programme.
Coco nano (Cn) is a cultivar characterised by determinate (dwarf) growth habit, early and prolific seed production which was selected under conventional conditions. Since Cn is the same commercial seed type as Gn (i.e., Navy), resultant crosses maintain this market class seed-phenotype.
The final aim of the breeding programme was to develop a panel of dwarf common bean lines characterised by high yield performance under organic management producing fast-cooking, smooth and doughy tasting seeds. Indeed, in P. vulgaris the determinate growth habit has been extensively exploited in breeding programmes to decrease plant biomass and optimise allocation between vegetative and reproductive growth [37]. In fact, the determinate growth habit and the shorter cycle allow many mechanical operations, including the harvest [38].
A schematic representation of the breeding programme carried out is shown in Figure 1. In 1998, 150 manual crosses were performed using the landrace Gn as the recurrent and seed-bearing female parent and the cultivar Cn (male parent) as the donor of the “dwarf” trait. In common bean this trait has been described as controlled by a single locus [39]. The high number of initial crosses performed was necessary to ensure introgression of the existing genetic diversity of Gn in the new breeding lines [35,40]. The F1 plants obtained, and those of the subsequent generations, were separated into families, corresponding to each single cross. With no exceptions, all the F1 plants showed an indeterminate growth habit. To fulfil the breeding objectives, two different approaches were used throughout the following years: (A) 33 F1 plants were backcrossed to the recurrent parent (Gn); (B) 47 F1 plants were self-fertilised (F2).
In the year 2000, lines of the two groups were self-pollinated and those in group B were selected for the dwarf growth habit. In the same year, the seed quality was evaluated by means of a panel test. The principal aim of the panel test was to identify the above-mentioned taste characteristics which are peculiar to Gn together with its fast cooking feature. All the samples tested were weighed and prepared according to the experimental conditions (soaking/cooking) described by Torricelli and colleagues [40]. Only the lines with a score equal to 4 (Gn-like) were selected as detailed in [40]. In 2001, the same procedure was repeated and the panel test was also carried out on the lines in group A (Figure 1). At the end of the same year, the best 18 families were identified according to the results of the panel test and the seed weight (lines that produced seeds with an average weight ≤ 0.38 g were discarded). In the panel test, seeds from the breeding lines were scored by “blind” tasters according to their organoleptic qualities and their response to the soaking/cooking treatment. The panellists were trained according to Garruti and Bourne [41] using seed samples of Gn and Cn; therefore, the breeding lines tested were scored according to the following criteria: 1 = Cn-like; 2 = similar to Cn; 3 = similar to Gn; 4 = Gn-like.
In 2002, 78 lines derived from the 18 previously selected families were subjected to the panel test again and the best 22 lines were selected. In the following breeding step (2003–2004), plants were selected for yield potential (i.e., number of seeds per plant) and 5 low yielding lines were discarded.

2.2. Multi-Environment Field Trials

Seventeen lines derived from the above described breeding programme and two cultivars, Cn and Cb (as commercial controls), were included in multi-environment field trials [42]. Cn is the dwarf cultivar included in the breeding programme, while Cb (Coco Blanc Gautier) is a cultivar, characterised by the same growth habit and belonging to the same commercial class seed phenotype as Cn and Gn.
The Gn landrace was not included in the trials because it shows a vigorous vine growth, with plants that can be as tall as four metres, and consequently are usually spaced 100 × 100 cm and starts flowering late (mid-August). For all these reasons it was impossible to compare Gn with dwarf plants, which require less space and flower much earlier (June). Therefore, all tested materials belonged to the same commercial seed class and plant growth habit according to Voysest and Dessert [43]. The list of assessed materials, with information on their biological status and pedigree notes, is reported in Table 1.
Eight field trials were carried out in four different locations during three consecutive crop seasons between 2010 and 2012. In order to include a wide range of environmental variability, the four experimental sites covered a latitude range of about 10° (from 42°02′35.34′′N to 52°21′14.15′′N) and an altitude from 26 to 310 m a.s.l. Two experimental sites were located in Italy, in San Martino in Campo, Perugia (L1) and Rispescia, Grosseto (L2). The other two were located in Fressingfield, United Kingdom (L3) and Bellegade du Razes, France (L4).
All experiments were carried out following the principles and guidelines for organic management systems established by the International Foundation for Organic Agriculture (IFOAM, 2005). Throughout the trials, average daily temperature (°C), average daily relative humidity (%), and total rainfall precipitations (mm) (from May to August) were recorded using a weather station (Table 2).
The entries were arranged in a completely randomised block design with three replications. For each entry and in each replication, five plants were assessed [44]. The row plots measured 300 cm in length and were spaced 70 cm apart; plants within row plots were spaced 50 cm apart. In all the environments, watering was applied according to Thung [45]. Finally, plants were manually harvested at maturity and the yield of each plant was measured as grams of seeds per plant (g plant−1).

2.3. Statistical Analysis of Yield Data

In the MET, each location-year combination was considered as a different environment (hereafter the environments are designated by a code combining the location ID with the year. e.g., L1-10 = San Martino in Campo 2010). Mean yield was calculated for each environment and for each line as grams of seeds per plant. In each environment, the harvest was carried out manually at maturity. Mean yields of the tested lines in each environment were separated using Duncan’s multiple range test using XLSTAT 2013 (Addinsoft).
In order to quantify the environment (E), the genotype (G) and the genotype by environment interaction (GE) effects on the yield data, a multi-environment analysis of variance (ANOVA) was performed using the same software according to the following mixed model:
yij = µ + αi + βj + φij + εij
where yij is the yield of genotype i in environment j, µ is the grand mean, αi is the genotype effect, βj is the environmental effect, φij is the genotype by environment interaction and εij is the random error.
GGE biplot analyses were carried out using GenStat 12.1 software (VSN International). Different graphical outputs were produced to: (i) identify adaptation patterns by visualising relations between the studied lines and the environments under organic management; (ii) rank the lines on the basis of the mean potential yield performance and stability; and (iii) highlight meaningful mega-environments and attempt to identify the best yielding line in each mega-environment (the so-called “which-won-where” method).

2.4. DNA Isolation, Markers Selection and Genetic Characterization

Young fresh leaves were collected from (i) six individual plants of each breeding line and (ii) a bulk of six plants of each of the parents (Gn and Cn) and of the control cultivar (Cb). Leaf tissues were disrupted using the TissueLyser II (Qiagen, Hilden, Germany) and the Genomic DNA was extracted using the DNeasy 96 Plant Kit (Qiagen, Hilden, Germany) according to the protocol provided by the manufacturer. DNA quality and concentration were assessed by means of UV-Vis spectrophotometry using the NanoDrop 2000 (Thermo Scientific, Waltham, MA, USA) and 1.2% (w/v) agarose gel separation.
A set of 59 Simple Sequence Repeat (SSR) markers (35 genomic and 24 gene-based) covering all 11 linkage groups was initially identified [46,47,48,49,50,51,52] (Table S1). Seven out of the 24 gene-based SSRs were located in genes involved in flowering control in common bean [52]. Flowering related traits are strongly linked to the reproductive success and seed production of a crop in a certain environment. Microsatellite prediction was carried out using WebSat software [53] and primer sequences were designed using Primer3 software [54]. All the other SSRs were retrieved through bibliographic analysis using the National Center for Biotechnology Information (NCBI) databases (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/). The entire set of 59 SSRs was arranged and optimised in multiplex PCRs using MPprimer software [55] and initially tested on eight genotypes (Cn, Cb and six breeding lines randomly chosen from those available).
Amplification of SSR loci was performed by PCR using a 2720 Thermal Cycler (Applied Biosystems, Foster City, CA, USA) in a volume of 12.5 μL, containing the following reagents: 20 ng of template DNA, 1× Type-it® Multiplex PCR Master 1 Mix (Qiagen, Hilden, Germany) and 0.4 µM of each primer. The reaction mixture was initially denatured at 94 °C for 4 min, followed by 12 cycles of amplification composed of: denaturation (94 °C for 30 s), annealing (temperatures are reported in Table S1) for 30 s and extension (72 °C for 30 s). Subsequently, 25 cycles of amplification were carried out as follows: 94 °C for 10 s, annealing temperature (decreased by 1 °C) for 15 s and 72 °C for 20 s. The final extension was carried out for 20 min at 72 °C.
2 μL of the amplified DNA was mixed with 3 μL of GelPilot DNA Loading Dye 5x (Qiagen, Hilden, Germany) and 10 μL of ultra-pure H2O and then loaded into 3% (w/v) agarose gels with 1× TAE buffer. The electrophoresis was performed at 100 V for approximately 1.5 h. All gels were visualised under UV light after staining with ethidium bromide. Taking into account the observed polymorphisms, amplification quality, annealing temperature and amplification fragment sizes, 25 SSR markers were selected for genotyping. The forward primers were 5′-end-labeled with one of the fluorescent tags 6-FAM, HEX, PET and VIC (Applied Biosystems, Foster City, CA, USA). Amplicons were separated and sized according to the internal size standard GeneScan 500 LIZ (Applied Biosystems, Foster City, CA, USA) using the automatic sequencer 3130XL Genetic Analyzer (Applied Biosystem, Foster City, CA, USA). The visualisation and scoring of the amplicons were manually performed using the GENEMAPPER software (Applied Biosystems, Foster City, CA, USA).

2.5. Genetic Diversity Analysis

Number of polymorphic loci (Npol), different alleles (Na), effective alleles (Ne) [56], observed and expected homozygosity (Homo and Home, respectively) and heterozygosity (Ho and He, respectively) of the tested breeding lines were assessed using the GenAlEx [57] and PowerMarker software [58]. The significance of differences among the breeding lines and among individuals within lines was tested using Analysis of Molecular Variance (AMOVA) in Arlequin v 3.5 software [59]. The significance level of the tested sources of variation was assessed using 9999 permutations.
A pairwise accession-by-accession genetic distance matrix was calculated using GenAlEx [60]. In order to obtain information on the genetic similarity/dissimilarity of the breeding lines with parental lines (Gn, Cn) and the control (Cb), the latter were also included in the analysis. The resulting matrix was used to perform a Principal Coordinate Analysis (PCoA) using the covariance matrix with data standardisation and to build a neighbor joining phylogenetic tree, using Mega 6 software [61].

3. Results

3.1. Multi-Environment Field Trials and Statistical Analysis

The mean yields of all the tested lines (17 breeding lines and two controls) in the eight environments (L1-10, L1-11, L1-12, L2-11, L2-12, L3-10, L3-11, L4-10) are reported in Table 3. The highest mean yield across environment was recorded for line 11 (55.66 g plant−1) and the lowest for line 9 (28.71 g plant−1). Among the tested environments, the highest mean yield was obtained in L3-11 (93.60 g plant−1) and the lowest in L1-11 (14.27 g plant−1).
The environments L3-10 and L3-11 were the best yielding ones because they showed the most favourable climatic conditions, characterized by the lowest average maximum and minimum temperatures between flowering and maturity of the pods (Tmin FM = 11.1 and 11.6, Tmax FM = 18.7 and 20.0 °C, respectively, Table 2).
All the sources of variation modelled in the ANOVA: E, G and GE, were significant and explained 55.5%, 4.7% and 7.1% of yield variance, respectively (Table 4). According to ANOVA results, environment was the most important source of yield variation. The relatively high value of GE compared to G suggests the probable existence of different mega-environments (ME) that need to be investigated [42]. Knowledge of MEs, namely, a group of environments defined by similar biotic and abiotic stresses and cropping system requirements, is important to direct more effectively the development of the germplasm involved in a breeding programme.
The GGE biplot explained a high proportion of the total GGE variance (i.e., 83.71%) (Figure 2). In Figure 2a, the cosine of the angle between two vectors approximates the correlation between them: the more acute an angle is, the more positively correlated the vectors are. In the same figure, both genotype (red) and environment (black) vectors are drawn such that the specific relations between a certain genotype and a certain environment can be visualised.
The genotypic vectors of the breeding lines 3, 4, 5, 6, 7 and of the Cn (18) and Cb (19) cultivars, showed obtuse angles (i.e., negative interaction) with vectors of all environments: these genotypes were not able to take advantage of the conditions in any environment. On the contrary, genotypic vectors of lines 8, 9, 13, 14, 15, and 16 formed acute angles with all the environment vectors, showing a positive interaction at all the experimental sites. Breeding lines 1, 2, 10, 11, 12 and 17 showed right angles and/or large acute angles (i.e., poorer interactions) with some environment vectors and positive interaction with others.
In Figure 2b, the GGE biplot reports a ranking of the genotypes based on their mean yield performance and stability. In the figure, the average-environment axis (AEA) abscissa points to higher mean yield across environments and the AEA ordinate points to greater variability (i.e., poorer stability or stronger interaction with a certain environment) in both directions. The average yield of the accessions on the left part of plot was generally lower than the grand mean, whilst that of the accessions on the right was higher.
The cultivars and breeding lines that were characterised by negative interaction with the environments (see above) also showed an average yield lower than the grand mean (Figure 2b). However, the other breeding lines showed mean yields higher than the grand mean. The “which-won-where” GGE biplot allowed us to identify MEs and to distinguish the best yielding line in each of them. In Figure 2c, a polygon is drawn on genotypes that are furthest from the biplot origin and perpendicular lines to each side of the polygon (or to their extension) are drawn, starting from the biplot origin. The latter line divides the biplot into different sectors that are different MEs. The eight studied environments fall into three sectors (ME-1, ME-2 and ME-3), indicating the existence of three meaningful MEs (Figure 2c). In the “which-won-where” biplot, the best yielding lines were 17, 15 and 16, that won in ME-1, ME-2 and ME-3 respectively.

3.2. Genotyping and Genetic Diversity

Among the 17 breeding lines, DNA was successfully extracted from 93 individuals and from the three bulks representing the two parents and the control. Of the 59 SSR markers amplified, 18 neutral and 7 gene-derived were found to be polymorphic and were consequently used for the genotyping (42.4% of the total). The selected SSR markers detected a total of 47 alleles; number of alleles per locus ranged from 1 to 3. Of the 47 identified alleles, 3 (6.38%) were private: two belonged to line 6 and one to line 7. The average number of alleles (Na) ranged from 1.00 to 1.08 and the average number of effective alleles (Ne) from 1.00 to 1.06, confirming that the breeding materials developed are virtually pure lines.
At each studied locus, both Cn and Cb bulks showed the presence of a single allele (Na = 1 and Ne = 1) which agrees with the genetic condition of common bean cultivars (i.e., pure lines). On the contrary, higher values of Na and Ne (both 1.64) were observed for Gn, complying with the known within-population diversity of this landrace (Table S2) [26].
The AMOVA (Table 5) showed that 94.0% of total molecular variance of the breeding lines was due to variation among lines (p < 0.001) and only 3.4% among individual plants within lines; finally, 2.6% of the molecular variance was due to variation among alleles within individuals. The molecular results reported clearly show that the UNIPG lines are genetically differentiated and can be regarded as near-pure lines.
In the PCoA, the first two principal coordinates accounted for 54.86% and 14.79% of the variation, respectively (total 69.65%) (Figure 3). The first axes clearly distinguished three groups of samples: (i) Gn together with most of the breeding lines (i.e., 1, 2, 3, 4, 8, 10, 11, 12, 13, 14, 15, 16, 17); (ii) lines 5, 6, 7, 9 and (iii) Cn and Cb. The second axes mainly distinguished lines 6 and 9 from lines 5 and 7, underlining the existence of a higher level of diversity in the second group in comparison with the other two. Indeed, according to their position on the plot, all the samples in the first group are highly similar. Finally, the PCoA did not highlight any genetic differentiation between the two cultivars (Cn and Cb).
The neighbor joining tree (Figure 4) confirmed that lines of the first group belong to a single cluster that also included the landrace (Gn). All the lines in this cluster derived from breeding strategy A, while lines 5, 6, 7, 9 from breeding strategy B. Cn and Cb are not grouped in any major cluster. Results of the genetic analyses comply with the pedigree of the lines, clearly separating those derived from the backcross (A) from those obtained by the simple cross (B).

4. Discussion

The two principal aims of this study were: (1) to evaluate the performance of 17 breeding lines that share the same organoleptic and fast-cooking features by using a multi-environmental approach under organic conditions and (2) to understand the level of genetic diversity among these lines. The high yields recorded in Fressingfield in 2010 and 2011 were due to favourable weather conditions. Indeed, in comparison to the other environments, the two years (i.e., L3-10 and L3-11) were characterized by the lowest average maximum and minimum temperatures between flowering and maturity of the pods. It is well known that most of the crops are particularly sensitive to the temperature recorded between flowering and maturity and how heat stress can reduce yields [62]. In common bean, it has been reported that high temperature stress causes abortion of buds, flowers and pods, leading to high yield losses [63,64,65]. Moreover, the ANOVA showed that all the sources of variation were significant and, of these, the effect of environment (E) was the greatest. The strong differences among environments (i.e., the combination of location and year) can be considered the main cause of the predominant role of E. Similar results have been reported by studies carried out on different self-pollinated leguminous crops, such as pea [66] and soybean [67], as well as on major cereals [22,68,69]. The effect of genetic variation (G) was also significant, suggesting that the genotypic variability among the tested lines significantly affected their performance. Considering that all the best performing materials were some of the lines obtained by crossing the Gnocchetto landrace and the Coco nano cultivar (a pure line according to our data), it can be inferred that the genomic contribution of the landrace determined their good performance under organic conditions.
To better understand the potential utility of the selected common bean lines under organic and different agro-climatic conditions, genotype main effect (G) plus genotype by environment (GE) interaction biplot analyses were performed. By dissecting the GE interaction, this methodology allows the identification of promising genotypes under different environmental conditions. This method has also been, and is still largely used by many breeding companies and research institutions to evaluate and test the performance of different genotypes in certain areas using conventional agricultural practices. It has already been underlined how GE interaction (positive or negative) is of major importance when a breeding programme is intended for organic and low-input agriculture [70]. Indeed, important traits such as yield and overall quality were reported to be strongly linked to tolerance of various stresses and are complex inherited traits characterised by high GE interactions [71].
The GGE biplot analysis revealed strong differences in terms of potential performance with three of the 17 breeding lines expressing their best performance in different mega-environments. In Figure 2b, both breeding lines 16 and 17, which were shown to be the “winning lines” in ME-1 and ME-3, respectively, are the furthest lines from the AEA abscissa on the vertical axis. This distance is closely related to genotype instability (as it represents the contribution of GE to the genotype score) which, in this case, might be interpreted as a specific interaction with a narrower range of environmental conditions in a certain location. By contrast, line 15, which was the “winning line” in ME-2, together with line 14, showed a broader ability to interact positively with all the study environments since the distance of their genotype score falls very close to the AEA abscissa. It is also noteworthy that in the same figure, line 11, which showed the highest mean yield across all environments (see Table 3), was not ultimately among the “winning lines”. In this case, the effect of the environment, discarded by the analysis, had a strong effect on the mean yield across environments. It is important to stress that Figure 2 explains only a portion of the total variation (i.e., G and GE) and some apparently stable genotypes may not be truly stable as the variation in their performance is not fully explained by the biplot.
The three MEs identified mainly group environments belonging to same geographical location, with very few exceptions. In fact, ME-1 includes two out of the three testing environments of San Martino in Campo, IT (L1-10, L1-12,) and ME-2 includes all the testing environments in the second Italian location (Rispescia) plus the only environment in San Martino in Campo (L1-11) that does not belong to ME-1. ME-1 and ME-2 group all the Italian test environments. In addition, it is noteworthy that in ME-2, the environment L2-12 lies farther from the other two (L1-11, L2-11) and it is characterised by severely dry conditions. Since L2-11 (the same location but in a different year) is characterised by the highest level of total precipitation, ME-2 is definitely characterized by variable rainfall conditions and, consequently, it seems that other environmental factors may have played a stronger role in the definition of this ME. Finally, ME-3 clearly groups all the environments of the two northern European experimental sites.
According to our results, the level of diversity among the breeding lines developed allowed the best performing lines in the different MEs to be identified, which mainly corresponded to the different geographical areas considered. Therefore, the breeding programme was successful in selecting a panel of diverse dwarf common bean breeding lines sharing the organoleptic qualities of the landrace Gn. Moreover, according to our data, the residual genetic diversity among lines sharing the same quality and technological traits (i.e., Gn tasting qualities and fast-cooking features) allowed those with different adaptation ability to be identified. It should be noted that a previous study highlighted the existence of a relatively high level of genetic diversity in Gn [26], which is in line with our results.
As suggested by Fess and colleagues [25], our data confirm that the genetic diversity within a landrace population is likely to be a very important resource for coping with unexpected stresses and/or adapting to different environmental conditions, especially when the population is conserved on-farm for many generations. Indeed, on-farm conservation is a dynamic form of conservation that allows crops to evolve in response to both biotic and abiotic stresses. In this regard, it also important to note that all the lines suitable for organic farming (lines 15, 16, 17) belong to cluster A which also includes Gn (Figure 4). This grouping of the tested materials suggests that the landrace might have contributed with traits for high yields under organic conditions.
Regarding the two different breeding approaches used to develop the trailed lines, apparently, the use of a single backcross (A) did not negatively affect introgression of the quality traits of Gn, indicating that early selection, based on the growth habit and the panel tests, effectively increased the general performance of the breeding programme. This modus operandi allowed the number of backcrosses—generally a substantial investment in terms of time and development costs for a breeding company and/or institution—to be minimised.

5. Conclusions

The breeding strategy adopted in this study allowed the optimisation of time and cost in developing new common bean lines for organic agriculture. Moreover, across a panel of lines that share the same organoleptic and technological qualities, multi-environment trials enabled effective selection of those with the highest yields in the different MEs. As current cultivars for organic and low-input agriculture—characterised by environmental variability—are still limited, these results, together with careful preservation, screening and use of common bean germplasm such as landraces, which are adapted to organic agriculture per se, are an important step forward in fostering and accelerating successful breeding programmes for this growing sector.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2071-1050/10/3/777/s1, Table S1: List of SSR markers used for the genotyping; Table S2: Genetic descriptive statistics of the studied common bean lines.

Acknowledgments

The research leading to these results was funded by the European Community’s Seventh Framework Programme (FP7/2007–2013) under grant agreement No. 245058 SOLIBAM (Strategies for Organic and Low-Input Integrated Breeding and Management). Thanks are due to Frederick Rey for carrying out the experiment in France and for providing data.

Author Contributions

Renzo Torricelli, Lorenzo Raggi and Valeria Negri conceived and designed the experiments; Leonardo Caproni, Lorenzo Raggi, Carlo Tissi and Sally Howlett performed the experiments; Leonardo Caproni, Carlo Tissi and Lorenzo Raggi analysed the data; Valeria Negri provided reagents/materials/analysis tools; all the authors contributed to the paper writing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Willer, H.; Kilcher, L. The World of Organic Agriculture—Statistics and Emerging Trends 2012; Willer, H., Kilcher, L., Eds.; IFOAM: Bonn, Germany; FiBL: Frick, Switzerland, 2012. [Google Scholar]
  2. Shafie, F.A.; Rennie, D. Consumer Perceptions Towards Organic Food. Procedia Soc. Behav. Sci. 2012, 49, 360–367. [Google Scholar] [CrossRef]
  3. Stanhill, G. The comparative productivity of organic agriculture. Agric. Ecosyst. Environ. 1990, 30, 1–26. [Google Scholar] [CrossRef]
  4. Padel, S.; Lampkin, N. Conversion to organic farming: An overview. In The Economics of Organic Farming; Lampkin, N.H., Padel, S., Eds.; CAB: Wallingford, UK, 1994; pp. 295–313. [Google Scholar]
  5. Mäder, P.; Fliessbach, A.; Dubois, D.; Gunst, L.; Fried, P.; Niggli, U. Soil fertility and biodiversity in organic farming. Science 2002, 296, 1694–1697. [Google Scholar] [CrossRef] [PubMed]
  6. Ryan, M.H.; Derrick, J.W.; Dann, P.R. Grain mineral concentrations and yield of wheat grown under organic and conventional management. J. Sci. Food Agric. 2004, 84, 207–216. [Google Scholar] [CrossRef]
  7. Smolik, J.; Dobbs, T. Crop yields and economic returns accompanying the transition to alternative farming systems. J. Prod. Agric. 1991, 4, 153–161. [Google Scholar] [CrossRef]
  8. Clark, M.S.; Horwath, W.R.; Shennan, C.; Scow, K.M.; Lantni, W.T.; Ferris, H. Nitrogen, weeds and water as yield-limiting factors in conventional, low-input, and organic tomato systems. Agric. Ecosyst. Environ. 1999, 73, 257–270. [Google Scholar] [CrossRef]
  9. Reganold, J.P.; Glover, J.D.; Andrews, P.K.; Hinman, H.R. Sustainability of three apple production systems. Nature 2001, 410, 926–930. [Google Scholar] [CrossRef] [PubMed]
  10. Pimentel, D.; Hepperly, P.; Hanson, J.; Douds, D.; Seidel, R. Environmental, Energetic, and Economic Comparisons of Organic and Conventional Farming Systems. Bioscience 2005, 55, 573–582. [Google Scholar] [CrossRef]
  11. Posner, J.L.; Baldock, J.O.; Hedtcke, J.L. Organic and Conventional Production Systems in the Wisconsin Integrated Cropping Systems Trials: I. Productivity 1990–2002. Agron. J. 2008, 100, 253–260. [Google Scholar] [CrossRef]
  12. De Ponti, T.; Rijk, B.; Van Ittersum, M.K. The crop yield gap between organic and conventional agriculture. Agric. Syst. 2012, 108, 1–9. [Google Scholar] [CrossRef]
  13. Seufert, V.; Ramankutty, N.; Foley, J.A. Comparing the yields of organic and conventional agriculture. Nature 2012, 485, 229–232. [Google Scholar] [CrossRef] [PubMed]
  14. Ponisio, L.C.; M’Gonigle, L.K.; Mace, K.C.; Palomino, J.; De Valpine, P.; Kremen, C. Diversification practices reduce organic to conventional yield gap. Proc. R. Soc. B Biol. Sci. 2015, 282, 1–7. [Google Scholar] [CrossRef]
  15. Lammerts Van Bueren, E.T.; Jones, S.S.; Tamm, L.; Murphy, K.M.; Myers, J.R.; Leifert, C.; Messmer, M.M. The need to breed crop varieties suitable for organic farming, using wheat, tomato and broccoli as examples: A review. NJAS Wagening. J. Life Sci. 2011, 58, 193–205. [Google Scholar] [CrossRef]
  16. Murphy, K.M.; Campbell, K.G.; Lyon, S.R.; Jones, S.S. Evidence of varietal adaptation to organic farming systems. Field Crop. Res. 2007, 102, 172–177. [Google Scholar] [CrossRef]
  17. Lammerts van Bueren, E.T.; Struik, P.C.; Jacobsen, E. Ecological concepts in organic farming and their consequences for an organic crop ideotype. NJAS Wagening. J. Life Sci. 2002, 50, 1–26. [Google Scholar] [CrossRef]
  18. Wolfe, M.S.; Baresel, J.P.; Desclaux, D.; Goldringer, I.; Hoad, S.; Kovacs, G.; Löschenberger, F.; Miedaner, T.; Østergård, H.; Lammerts Van Bueren, E.T. Developments in breeding cereals for organic agriculture. Euphytica 2008, 163, 323–346. [Google Scholar] [CrossRef]
  19. Hildermann, I.; Thommen, A.; Dubois, D.; Boller, T.; Wiemken, A.; Mader, P. Yield and baking quality of winter wheat cultivars in different farming systems of the DOK long-term trial. J. Sci. Food Agric. 2009, 2477–2491. [Google Scholar] [CrossRef]
  20. Saastamoinen, M.; Hietaniemi, V.; Pihlava, J.-M.; Merja, E.; Kontturi, M.; Tuuri, H.; Niskanen, M.; Kangas, A. β-Glucan contents of groats of different oat cultivars in official variety, in organic cultivation, and in nitrogen ferilization trials in Finland. Agric. Food Sci. 2004, 13, 68–79. [Google Scholar] [CrossRef]
  21. Heilig, J.A.; Kelly, J.D. Performance of Dry Bean Genotypes Grown under Organic and Conventional Production Systems in Michigan. Agron. J. 2012, 104, 1485. [Google Scholar] [CrossRef]
  22. Annicchiarico, P.; Chiapparino, E.; Perenzin, M. Response of common wheat varieties to organic and conventional production systems across Italian locations, and implications for selection. Field Crop. Res. 2010, 116, 230–238. [Google Scholar] [CrossRef]
  23. Raggi, L.; Ciancaleoni, S.; Torricelli, R.; Terzi, V.; Ceccarelli, S.; Negri, V. Field Crops Research Evolutionary breeding for sustainable agriculture: Selection and multi-environmental evaluation of barley populations and lines. Field Crop. Res. 2017, 204, 76–88. [Google Scholar] [CrossRef]
  24. Lammerts van Bueren, E.T.; Backes, G.; de Vriend, H.; Østergård, H. The role of molecular markers and marker assisted selection in breeding for organic agriculture. Euphytica 2010, 175, 51–64. [Google Scholar] [CrossRef]
  25. Fess, T.L.; Kotcon, J.B.; Benedito, V.A. Crop breeding for low input agriculture: A sustainable response to feed a growing world population. Sustainability 2011, 3, 1742–1772. [Google Scholar] [CrossRef]
  26. Tiranti, B.; Negri, V. Selective microenvironmental effects play a role in shaping genetic diversity and structure in a Phaseolus vulgaris L. landrace: Implications for on-farm conservation. Mol. Ecol. 2007, 16, 4942–4955. [Google Scholar] [CrossRef] [PubMed]
  27. Brush, S. The issues of in situ conservation of crop genetic resources. In Genes in the Field. On Farm Conservation of Crop Diversity; Brush, S., Ed.; IPGRI: Rome, Italy, 2000; pp. 3–28. [Google Scholar]
  28. Klaedtke, S.; Caproni, L.; Klauck, J.; de la Grandville, P.; Dutartre, M.; Stassart, P.; Chable, V.; Negri, V.; Raggi, L. Short-Term Local Adaptation of Historical Common Bean (Phaseolus vulgaris L.) Varieties and Implications for In Situ Management of Bean Diversity. Int. J. Mol. Sci. 2017, 18, 493. [Google Scholar] [CrossRef] [PubMed]
  29. Negri, V.; Tosti, N. Phaseolus genetic diversity maintained on-farm in central Italy. Genet. Resour. Crop Evol. 2002, 49, 511–520. [Google Scholar] [CrossRef]
  30. Broughton, W.; Hernandez, G.; Blair, M. Beans (Phaseolus spp.)—model food legumes. Plant Soil 2003, 252, 55–128. [Google Scholar] [CrossRef]
  31. Food and Agriculture Organization of the United Nations. FAOSTAT Database. Available online: http://www.fao.org/faostat/en/#data/QC (accessed on 20 August 2017).
  32. United States Department of Agriculture (USDA). Economic Research Service. 2013. Available online: https://www.ers.usda.gov/data-products/organic-production.aspx (accessed on 20 August 2017).
  33. European Commission Statistics Database (EUROSTAT). Available online: http://ec.europa.eu/eurostat/data/database (accessed on 20 August 2017).
  34. Graham, P.H.; Ranalli, P. Common bean (Phaseolus vulgaris L.). Field Crop. Res. 1997, 53, 131–146. [Google Scholar] [CrossRef]
  35. Negri, V.; Tiranti, B. Effectiveness of in situ and ex situ conservation of crop diversity. What a Phaseolus vulgaris L. landrace case study can tell us. Genetica 2010, 138, 985–998. [Google Scholar] [CrossRef] [PubMed]
  36. Wiesinger, J.A.; Cichy, K.A.; Glahn, R.P.; Grusak, M.A.; Brick, M.A.; Thompson, H.J.; Tako, E. Demonstrating a Nutritional Advantage to the Fast-Cooking Dry Bean (Phaseolus vulgaris L.). J. Agric. Food Chem. 2016, 64, 8592–8603. [Google Scholar] [CrossRef] [PubMed]
  37. Kwak, M.; Toro, O.; Debouck, D.G.; Gepts, P. Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris). Ann. Bot. 2012, 110, 1573–1580. [Google Scholar] [CrossRef] [PubMed]
  38. Kornegay, J.; White, J.W.; de la Cruz, O.O. Growth habit and gene pool effects on inheritance of yield in common bean. Euphytica 1992, 62, 171–180. [Google Scholar] [CrossRef]
  39. Kwak, M.; Velasco, D.; Gepts, P. Mapping homologous sequences for determinacy and photoperiod sensitivity in common bean (Phaseolus vulgaris). J. Hered. 2008, 99, 283–291. [Google Scholar] [CrossRef] [PubMed]
  40. Torricelli, R.; Lorenzetti, S.; Falcinelli, M. La valorizzazione di una varietà locale di fagiolo attraverso il miglioramento genetico. Sementi Elette 2003, 6, 35–39. [Google Scholar]
  41. Garruti, R.D.S.; Bourne, M.C. Effect of Storage Conditions of Dry Bean Seeds ( Phaseolus on Texture Profile Parameters After Cooking vulgaris L.). J. Food Sci. 1985, 50, 1067–1071. [Google Scholar] [CrossRef]
  42. Yan, W.; Tinker, N.A. Biplot analysis of multi-environment trial data: Principles and applications. Can. J. Plant Sci. 2006, 86, 623–645. [Google Scholar] [CrossRef]
  43. Voysest, O.; Dessert, M. Bean cultivars: Classes and commercial seed types. In Common Beans: Research for Crop Improvement; van Schoonhoven, A., Yovest, O., Eds.; CAB International, in association with CIAT: Wallingford, UK, 1991; pp. 119–162. [Google Scholar]
  44. De Ron, A.M.; Rodiño, A.P.; Santalla, M.; González, A.M.; Lema, M.J.; Martín, I.; Kigel, J. Seedling Emergence and Phenotypic Response of Common Bean Germplasm to Different Temperatures under Controlled Conditions and in Open Field. Front. Plant Sci. 2016, 7, 1087. [Google Scholar] [CrossRef] [PubMed]
  45. Thung, M. Bean Agronomy in Monoculture. In Common Beans: Research for Crop Improvement; van Schoonhoven, A., Yovest, O., Eds.; CAB International, in association with CIAT: Wallingford, UK, 1991; pp. 737–834. [Google Scholar]
  46. Yu, K.; Park, S.J.; Poysa, V.; Gepts, P. Integration of simple sequence repeat (SSR) markers into a molecular linkage map of common bean (Phaseolus vulgaris L.). J. Hered. 2000, 91, 429–434. [Google Scholar] [CrossRef] [PubMed]
  47. Gaitán-Solís, E.; Duque, M.C.C.; Edwards, K.J.J.; Tohme, J. Microsatellite Repeats in Common Bean (Phaseolus vulgaris): Isolation, Characterization, and Cross-Species Amplification in Phaseolus spp. Crop Sci. 2002, 42, 2128–2136. [Google Scholar] [CrossRef]
  48. Guerra-Sanz, J.M. Short Communication—New SSR markers of Phaseolus vulgaris from sequence databases. Plant Breed. 2004, 89, 87–90. [Google Scholar] [CrossRef]
  49. Blair, M.W.; Pedraza, F.; Buendia, H.F.; Gaitán-Solís, E.; Beebe, S.E.; Gepts, P.; Tohme, J. Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.). Theor. Appl. Genet. 2003, 107, 1362–1374. [Google Scholar] [CrossRef] [PubMed]
  50. Blair, M.W.; Muñoz-Torres, M.C.; Pedraza, F.; Giraldo, M.C.; Buendía, H.F.; Hurtado, N. Development of microsatellite markers for common bean (Phaseolus vulgaris L.) based on screening of non-enriched, small-insert genomic libraries. Genome 2009, 52, 772–782. [Google Scholar] [CrossRef] [PubMed]
  51. Córdoba, J.M.; Chavarro, C.; Rojas, F.; Muñoz, C.; Blair, M.W. Identification and Mapping of Simple Sequence Repeat Markers from Common Bean (Phaseolus vulgaris L.) Bacterial Artificial Chromosome End Sequences for Genome Characterization and Genetic–Physical Map Integration. Plant Genome J. 2010, 3, 154. [Google Scholar] [CrossRef]
  52. Raggi, L.; Tissi, C.; Mazzucato, A.; Negri, V. Molecular polymorphism related to flowering trait variation in a Phaseolus vulgaris L. collection. Plant Sci. 2014, 215–216, 180–189. [Google Scholar] [CrossRef] [PubMed]
  53. Martins, W.S.; Lucas, D.C.S.; de Souza Neves, K.F.; Bertioli, D.J. WebSat—A web software for microsatellite marker development. Bioinformation 2009, 3, 282–283. [Google Scholar] [CrossRef] [PubMed]
  54. Rozen, S.; Skaletsky, H. Primer3 on the WWW for general users and for biologist programmers. In Methods in Molecular Biology: Bioinformatics Methods and Protocols; Krawetz, S., Misener, Eds.; Humana Press Inc.: Totowa, NJ, USA, 1999; pp. 365–386. [Google Scholar]
  55. Shen, Z.; Qu, W.; Wang, W.; Lu, Y.; Wu, Y.; Li, Z.; Hang, X.; Wang, X. MPprimer: A program for reliable multiplex PCR primer design. BMC Bioinform. 2010, 11, 143. [Google Scholar] [CrossRef] [PubMed]
  56. Kimura, M.; Crow, J.F. The Number of Alleles That Can Be Maintained in a Finite Population. Genetics 1964, 49, 725–738. [Google Scholar] [PubMed]
  57. Peakall, R.; Smouse, P.E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  58. Liu, K.; Muse, S.V. PowerMarker: An integrated analysis environment for genetic marker analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef] [PubMed]
  59. Excoffier, L.; Lischer, H.E.L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef] [PubMed]
  60. Peakall, R.; Smouse, P.E. GenALEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [PubMed]
  61. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  62. Mcclean, P.E.; Burridge, J.; Beebe, S.E.; Rao, I.M.; Porch, T.G. Crop improvement in the era of climate change: An integrated, multi-disciplinary approach for common bean (Phaseolus vulgaris). Funct. Plant Biol. 2011, 38, 927–933. [Google Scholar] [CrossRef]
  63. Rainey, K.M.; Griffi, P.D. Inheritance of Heat Tolerance during Reproductive Development in Snap Bean (Phaseolus vulgaris L.). Differences 2005, 130, 700–706. [Google Scholar]
  64. Omae, H.; Kumar, A.; Shono, M. Adaptation to high temperature and water deficit in the common bean (Phaseolus vulgaris L.) during the reproductive period. J. Bot. 2012, 2012, 803413. [Google Scholar] [CrossRef]
  65. Porch, T.G.; Smith, J.R.; Beaver, J.S.; Griffiths, P.D.; Station, E.; Street, W.N.; Canaday, C.H. TARS-HT1 and TARS-HT2 Heat-tolerant Dry Bean Germplasm. Hortic. Sci. 2010, 45, 1278–1280. [Google Scholar]
  66. Tolessa, T.T.; Keneni, G.; Sefera, T.; Jarso, M.; Bekele, Y. Genotype × Environment Interaction and Performance Stability for Grain Yield in Field Pea (Pisum sativum L.) Genotypes. Int. J. Plant Breed. 2013, 7, 116–123. [Google Scholar]
  67. Edugbo, R.E.; Nwofia, G.E.; Fayeun, L.S. An Assessment of Soybean (Glycine max, L. Merrill ) Grain Yield in Different Environments Using AMMI and GGE Biplot Models in Humidorest Fringes of Southeast Nigeria. Agric. Trop. Subtrop. 2015, 48, 82–90. [Google Scholar] [CrossRef]
  68. Mehari, M.; Sentayehu, A.; Berhane, L. Genotype X Environment Interaction and Yield Stability of Malt Barley Genotypes Evaluated in Tigray, Ethiopia Using the Ammi Analysis. Asian J. Plant Sci. 2014, 13, 73–79. [Google Scholar]
  69. Yahiaoui, S.; Cuesta-Marcos, A.; Gracia, M.P.; Medina, B.; Lasa, J.M.; Casas, A.M.; Ciudad, F.J.; Montoya, J.L.; Moralejo, M.; Molina-Cano, J.L.; et al. Spanish barley landraces outperform modern cultivars at low-productivity sites. Plant Breed. 2014, 133, 218–226. [Google Scholar] [CrossRef]
  70. Ceccarelli, B.S.; Erskine, W.; Hamblin, J.; Grando, S. Genotype by environment interaction and international breeding programmes. Exp. Agric. 1994, 30, 177–187. [Google Scholar] [CrossRef]
  71. Rao, I.M.; Beebe, S.E.; Polania, J.; Grajales, M.; Cajiao, C.; Ricaurte, J.; García, R.; Rivera, M. Evidence for genotypic differences among elite lines of common bean in the ability to remobilize photosynthate to increase yield under drought. J. Agric. Sci. 2017, 155, 857–875. [Google Scholar] [CrossRef]
Figure 1. A schematic representation of the common bean breeding programme carried out at Dipartimento di Scienze Agrarie, Alimentari e Ambientali (DSA3). “R” and “r” indicate the dominant and recessive alleles controlling the climbing and dwarf trait, respectively.
Figure 1. A schematic representation of the common bean breeding programme carried out at Dipartimento di Scienze Agrarie, Alimentari e Ambientali (DSA3). “R” and “r” indicate the dominant and recessive alleles controlling the climbing and dwarf trait, respectively.
Sustainability 10 00777 g001
Figure 2. Different views of the GGE biplot for yield data of 17 common bean lines and two cultivars in eight environments under organic management. In the biplot, each environment is a location by year combination. Location labels are: L1, San Martino in Campo, Perugia, IT; L2, Rispescia, Grosseto, IT; L3, Fressingfield, Suffolk, UK; L4, Bellegarde du Razes, Limoux, FR. In the figure, numbers from 1 to 17 refer to the UNIPG corresponding breeding lines; 18 to Coco nano (Cn); and 19, to Coco blanc Gautier (Cb). (a) relationships among the tested lines (×) and test environments (+ vectors); (b) ranking of the lines by mean yield potential and stability. In this biplot view, the average environment axes (AEA), and average environment coordinates (AEC) are also shown; (c) The “which-won-where” view of the GGE biplot. ME 1, ME 2 and ME-3 indicate the three meaningful mega-environments. Biplots (a,c) are environment-centred, without scaling, and environment-metric preserving; biplot (b) is environment-centred, without scaling, and genotype-metric preserving.
Figure 2. Different views of the GGE biplot for yield data of 17 common bean lines and two cultivars in eight environments under organic management. In the biplot, each environment is a location by year combination. Location labels are: L1, San Martino in Campo, Perugia, IT; L2, Rispescia, Grosseto, IT; L3, Fressingfield, Suffolk, UK; L4, Bellegarde du Razes, Limoux, FR. In the figure, numbers from 1 to 17 refer to the UNIPG corresponding breeding lines; 18 to Coco nano (Cn); and 19, to Coco blanc Gautier (Cb). (a) relationships among the tested lines (×) and test environments (+ vectors); (b) ranking of the lines by mean yield potential and stability. In this biplot view, the average environment axes (AEA), and average environment coordinates (AEC) are also shown; (c) The “which-won-where” view of the GGE biplot. ME 1, ME 2 and ME-3 indicate the three meaningful mega-environments. Biplots (a,c) are environment-centred, without scaling, and environment-metric preserving; biplot (b) is environment-centred, without scaling, and genotype-metric preserving.
Sustainability 10 00777 g002
Figure 3. Plot of the first two axes of the PCoA. Open circles and open squares refer to breeding lines obtained using the breeding approach A (backcross) and B (F1 self-pollination), respectively; filled triangles refer to the two parents (Gn and Cn). In this genetic assessment Cb was identical to Cn such that they are not distinguishable in the Figure.
Figure 3. Plot of the first two axes of the PCoA. Open circles and open squares refer to breeding lines obtained using the breeding approach A (backcross) and B (F1 self-pollination), respectively; filled triangles refer to the two parents (Gn and Cn). In this genetic assessment Cb was identical to Cn such that they are not distinguishable in the Figure.
Sustainability 10 00777 g003
Figure 4. Neighbor joining tree based on the genetic distance matrix. Open circles and open squares refer to breeding lines obtained using the breeding approach A (backcross) and B (F1 self-pollination) respectively; filled triangles refer to the two parents (Gn and Cn).
Figure 4. Neighbor joining tree based on the genetic distance matrix. Open circles and open squares refer to breeding lines obtained using the breeding approach A (backcross) and B (F1 self-pollination) respectively; filled triangles refer to the two parents (Gn and Cn).
Sustainability 10 00777 g004
Table 1. List of studied materials with their biological status and pedigree notes.
Table 1. List of studied materials with their biological status and pedigree notes.
IDNameBiological StatusPedigree Note
Line-1Solibam-01Breeding lineBC1F1 offspring
Line-2Solibam-02Breeding lineBC1F1 offspring
Line-3Solibam-03Breeding lineBC1F1 offspring
Line-4Solibam-04Breeding lineBC1F1 offspring
Line-5Solibam-05Breeding lineF2 offspring
Line-6Solibam-06Breeding lineF2 offspring
Line-7Solibam-07Breeding lineF2 offspring
Line-8Solibam-08Breeding lineBC1F1 offspring
Line-9Solibam-09Breeding lineF2 offspring
Line-10Solibam-10Breeding lineBC1F1 offspring
Line-11Solibam-11Breeding lineBC1F1 offspring
Line-12Solibam-12Breeding lineBC1F1 offspring
Line-13Solibam-13Breeding lineBC1F1 offspring
Line-14Solibam-14Breeding lineBC1F1 offspring
Line-15Solibam-15Breeding lineBC1F1 offspring
Line-16Solibam-16Breeding lineBC1F1 offspring
Line-17Solibam-17Breeding lineBC1F1 offspring
Gn *GnocchettoLandraceParental line
CnCoco nanoCultivarParental line
CbCoco blanc GautierCultivarControl
* The Gn landrace was not included in the MET trials.
Table 2. Weather data recorded in the different tested environments (i.e., location by year). The environments are designated by a code combining the location code and the year (e.g., L1-10 = San Martino in Campo, 2010).
Table 2. Weather data recorded in the different tested environments (i.e., location by year). The environments are designated by a code combining the location code and the year (e.g., L1-10 = San Martino in Campo, 2010).
San Martino in Campo, Perugia, ITRispescia, Grosseto, ITFressingfield, Suffolk, UKBellegarde du Razes Limoux, FR
Coordinates43°02′35.34′′ N;
12°25′00.94′′ E
42°42′11.81′′N;
11°08′31.33′′E
52°21′14.15′′N;
1°21′16.73′′E
43°06′31.07′′N;
2°02′44.43′′E
Altitude (m a.s.l.)1742650310
EnvironmentL1-10L1-11L1-12L2-11L2-12L3-10L3-11L4-10
Year20102011201220112012201020112010
Tmin a SF b13.015.113.115.214.711.59.014.1
Tmax c SF26.828.428.027.427.120.918.425.1
Rh d SF62.359.456.967.044.868.970.368.1
Rnfl e SF95.886.863.441.824.649.038.081.3
Tmin FM f16.816.717.617.318.611.111.617.2
Tmax FM31.931.534.129.631.518.720.029.2
Rh FM51.351.341.166.842.377.375.859.0
Rnfl FM87.263.653.6116.232.6111.882.065.4
Sowing date20 May24 May16 May20 May21 May2 June12 May21 May
Average dtf g4544474740435341
Average dtm h11210310798107136121101
a Minimum temperature (°C); b From sowing to flowering; c Maximum temperature (°C); d Relative humidity (%); e Rain precipitations (mm); f From flowering to maturity; g Days to flowering; h Days to maturity.
Table 3. Mean yield (g plant−1) of the 17 lines and the two cultivars (controls) tested at eight environments under organic management.
Table 3. Mean yield (g plant−1) of the 17 lines and the two cultivars (controls) tested at eight environments under organic management.
IDEnvironment
L1-10L1-11L1-12L2-11L2-12L3-10L3-11L4-10Mean
Line-165.3AB14.6AB27.7A39.5CD18.8D72.2BCDE119.3A13.0AB46.31ABC
Line-264.4ABC14.7AB15.7AB56.0ABCD28.1ABCD69.6BCDE110.6AB20.5AB47.44ABC
Line-368.4A13.9AB6.4B52.5ABCD18.2D85.5ABCD94.8AB24.6A45.53BC
Line-457.7ABCD19.6AB24.3AB65.3ABC41.7ABC106.5A84.9AB17.6AB52.21AB
Line-523.8DE9.9AB11.1AB43.3BCD29.7ABCD61.7CDE82.4B19.2AB35.12DEF
Line-611.7E3.2B15.0AB35.1D18.9D44.3E101.3AB20.7AB31.26F
Line-715.3E4.6AB9.2AB55.1ABCD29.9ABCD57.8DE48.7C15.3AB29.50F
Line-852.7ABCD11.7AB13.7AB57.6ABCD30.6ABCD79.6ABCD107.9AB20.9AB46.85ABC
Line-916.3DE4.9AB7.8AB35.4D24.8CD48.3E80.3BC11.9B28.71F
Line-1048.4ABCD25.3A4.3B66.2ABC33.1ABCD89.6ABCD85.4AB19.5AB46.48ABC
Line-1159.5ABC21.9AB29.9A74.5A45.8A83.0ABCD109.4AB21.2AB55.66A
Line-1254.2ABCD20.3AB39.0A67.3ABC44.9AB75.2ABCDE103.3AB16.4AB52.58AB
Line-1341.1BCD9.9AB15.5AB46.2BCD33.8ABCD71.4BCDE78.8BC17.1AB39.24CDE
Line-1436.7BCDE11.8AB22.2AB46.6ABCD19.7D66.5BCDE100.9AB13.6AB39.74CDE
Line-1554.4ABCD20.5AB17.3AB53.5ABCD30.9ABCD94.1ABC82.0B17.1AB46.24ABC
Line-1635.5CDE15.0AB19.8AB53.6ABCD36.6ABCD69.7BCDE90.9AB19.1AB42.53BCD
Line-1735.8BCDE20.1AB42.7A59.0ABCD26.7BCD91.2ABC95.5AB19.0AB48.75ABC
Cn46.5ABCD17.2AB21.4AB65.9ABC27.9ABCD96.5AB98.3AB19.1AB49.10AB
Cb28.9DE11.9AB8.9AB67.9AB32.2ABCD89.8ABC103.6AB19.9AB45.38BC
Mean42.9814.2718.5354.7630.1276.4493.6018.1943.61
Table 4. ANOVA results.
Table 4. ANOVA results.
Source of Variationd.f. aSS bSS (%) bF cp-Value
E d7665,932.455.5207.2***
G e1856,264.24.76.8***
G x E f12685,050.67.11.5***
Error1368393,525.032.8
Total15191,200,772.1100.0
a Degree of freedom; b Sum of squares; c Fisher”s exact test; d Environment; e Genotype; f Genotype by environment interaction; *** Significant at p < 0.001.
Table 5. AMOVA results.
Table 5. AMOVA results.
Source of Variationd.f. aSS bv.c. cv d
Among lines16449.92.594.0%
Among individuals within lines7619.20.13.4%
Within Individuals936.50.12.6%
Total185475.72.7100.0%
a Degree of freedom; b Sum of squares; c Variance components; d Percentage of variation.

Share and Cite

MDPI and ACS Style

Caproni, L.; Raggi, L.; Tissi, C.; Howlett, S.; Torricelli, R.; Negri, V. Multi-Environment Evaluation and Genetic Characterisation of Common Bean Breeding Lines for Organic Farming Systems. Sustainability 2018, 10, 777. https://0-doi-org.brum.beds.ac.uk/10.3390/su10030777

AMA Style

Caproni L, Raggi L, Tissi C, Howlett S, Torricelli R, Negri V. Multi-Environment Evaluation and Genetic Characterisation of Common Bean Breeding Lines for Organic Farming Systems. Sustainability. 2018; 10(3):777. https://0-doi-org.brum.beds.ac.uk/10.3390/su10030777

Chicago/Turabian Style

Caproni, Leonardo, Lorenzo Raggi, Carlo Tissi, Sally Howlett, Renzo Torricelli, and Valeria Negri. 2018. "Multi-Environment Evaluation and Genetic Characterisation of Common Bean Breeding Lines for Organic Farming Systems" Sustainability 10, no. 3: 777. https://0-doi-org.brum.beds.ac.uk/10.3390/su10030777

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