A method for obtaining flexible broccoli varieties for sustainable agriculture
BMC Genetics volume 21, Article number: 51 (2020)
The use of high inputs in agriculture resulted in few varieties (hybrids and pure lines) used in all agricultural systems. Also varieties of vegetables, including broccoli, for organic and low-input agriculture, are almost exclusively hybrids, since there are very few specific breeding programs and varieties for sustainable agriculture systems.
A strategy to overcome this issue is the adoption of specific breeding programs for developing heterogeneous varieties (i.e. synthetics, open pollinated varieties, composite cross populations and mixtures). In fact, heterogeneous varieties are able to evolve and adapt to specific agro-climatic conditions.
The aim of this study was to develop a method (an Evolutionary Breeding Program, EBP) for obtaining heterogeneous varieties of broccoli and test its efficiency in developing highly diverse varieties, as needed in sustainable agriculture.
A synthetic variety originated from a landrace was multiplied in different environments for 3 cycles and morpho-phenological and genetic diversity of the derived populations were assessed.
The presented results are the first and unique indication about the efficiency of a short-time EBP for an allogamous species like broccoli.
Few morphological changes were observed among varieties multiplied in different environments with different agro-climatic conditions. This could be probably due to the initial genetic diversity of the landrace from which the populations were selected and also to the great plasticity of the crop. However, SSR data highlighted a genetic differentiation among populations multiplied for two/three years across Europe and in Central Italy, that was not so evident when considering morphological data only.
Few years of multiplication in different environments resulted in genetically differentiated broccoli populations that nonetheless preserved the original genetic diversity and productivity level and appear to evolve in relationship to different environments: the applied EBP is useful for developing heterogeneous materials for sustainable agriculture.
Brassica oleracea L. (2n = 2x = 18) is an economically important vegetable and its world annual commercial production (cauliflower and broccoli) is over 25 Mt, currently cultivated on about 1.4 Mha worldwide . China is the first country in producing broccoli and cauliflower (10 million tons), meanwhile, in Italy, the estimated production of both vegetables is about 390,000 tons  of which 10% is obtained in organic agriculture . Among the different B. oleracea subspecies, broccoli (B. oleracea L. var. italica Plenck, also identified as B. oleracea L. spp. capitata L. (DC.) convar. Botrytis (L.) Alef. var. italica Plenck ) is an important source of vitamins, minerals and antioxidants.
Currently, organic farmers largely depend on broccoli varieties bred for high external input in conventional farming systems [4, 5], and in particular, broccoli varieties are almost exclusively hybrids , due to lack of specific breeding programs and varieties for sustainable (organic and low-input) agriculture.
Suitable varieties for sustainable agriculture should have specific characteristics (as yield stability under different agro-climatic conditions, resistances to biotic and abiotic stresses, competition ability with weeds) in order to avoid the utilisation of external input like pesticides, herbicides and chemical fertilizers [5, 6]. In the context of sustainable agriculture, to allow varieties to evolve and adapt to specific agro-climatic conditions is also important . This is possible only when heterogeneous varieties (i.e. landraces, synthetics, open pollinated varieties, composite cross populations, mixtures) are used and reproduced on-farm for generations.
As a consequence, breeding programs specifically adopted for sustainable agriculture should consider all the above. Among many breeding approaches, the Evolutionary plant Breeding Program (EBP), combining natural and artificial selection in target environments, is suggested to be effective in obtaining heterogeneous varieties for organic and low-input agriculture [8,9,10,11,12,13].
The success of an EBP, like that of any other breeding program, is mainly due the genetic diversity of the initial materials. For example, EBP in cereals has been based on composite cross populations (CCP) or mixtures [8, 10, 13, 14] that can be originated from old varieties or landraces. Also for allogamous species, genetically variable varieties, like open pollinated, synthetic varieties or landraces, which have advantages in terms of yield stability and resistance , have been recommended. In particular, landraces could be the best starting materials for developing varieties suitable for organic and low-input agriculture [5, 16] due their great genetic diversity that provides a buffer against environment fluctuations and adaptation to specific environments across time.
An assessment of genetic diversity is therefore essential for the organization, development and control of a breeding program. The genetic diversity and structure of populations under selection can be efficiently evaluated by molecular markers. Depending on their type (i.e. neutral or related-to-gene markers) they can also provide information about the intensity of both natural and human selection and about the genetic divergence of populations developed across years and environments. To this purpose, different types of molecular markers have been used. For example, genetic changes in maize composite populations under Reciprocal Recurrent Selection (RRS) were analyzed by using Restriction Fragment Length Polymorphism (RFLP) [17, 18], Simple Sequence Repeats (SSR) [19,20,21] and Single Nucleotide Polymorphism (SNP) . SSR were also used to evaluate the genetic changes in different maize populations obtained after about 20 years of stratified mass selection on two historical Portuguese populations . Considering autogamous species, the genetic evolution over time and space of wheat CCP was assessed by RFLP [24, 25] and by SSR  and that of a barley CCP by using SSR .
Agronomic performances of different broccoli heterogeneous varieties or other broccoli breeding materials for sustainable agriculture have been already reported [5, 27,28,29], but few data are available on genetic diversity of broccoli breeding materials [30,31,32] and on their evolution and adaptations under different agro-climatic conditions.
Considering this lack of knowledge, the aims of this paper were to i) describe a EBP method for developing flexible broccoli varieties (i.e. populations, for sustainable agriculture), ii) assess the adaptation potential and the morpho-phenological and genetic changes of the populations developed with this method in different environments and iii) in the adaptation area across years.
Comparisons of populations multiplied for up to two cycles in different environments
Morpho-phenological characterisation and differentiation
The ANOVA analysis showed that the entries differ for all of the morpho-phenological traits recorded with great differences generally observed between the populations and three of the hybrids (HH, HHH, HHHH) (Table 1). For the traits weight of the main head (YLD) and number of secondary heads (HN) the populations derived from the initial synthetic by multiplication for two years were not different from each other. In particular, YLD values ranged from 11.05 g (Syn2-TER) to 324.58 g (HHH), as populations formed a little main head and many (up to 36) secondary heads while most of the hybrids did not form any secondary heads HHH, among the hybrids, and Syn2-FR, among the populations, showed the highest total yield (324.8 and 266.2 g, respectively).
Considering the inflorescence emission and maturation date (DH range = 101–204 dd; DM range = 134–219 dd), generally, the latest were the populations, the earliest the hybrids, with the exception of H. LR, Syn1-PG, Syn2-FR, Syn2-UK and Syn2-GR were the latest for inflorescence emission date, with Syn2-GR also late in maturation. Regarding traits related to head characteristics, generally, hybrids showed higher values than populations: HeH range = 6.41–9.59 cm, (Syn2-TER and H, respectively), HeW range = 3.29–14.53 cm (Syn2-UK and HH, respectively) and SW range = 0.94–3.69 cm (Syn2-TER and HHHH, respectively). The hybrids showed the lowest plant height, while LR and Syn1-PG the highest (PH range = 22.05–50.78 cm, HHHH and LR, respectively). Regarding plant diameter (PD range = 36.31–80.39 cm, in H and Syn2-PG, respectively), all the entries showed similar values, with the exception of H, HH and Syn2-TER showing the lowest values. No differences were found among populations regarding vigour with the only exception of Syn2-TER (4.0) which was the less vigorous one.
The first two PCA axes accounted for 83.57% of total data variability, with the addition of the third axis the 97.06% of the total variability was explained (Fig. 1). The entries were clustered in two groups: the first one included HH, HHH and HHHH, while the second included the populations with Syn2-TER and H as subgroups of this cluster. Most of the morpho-phenological traits had large effects on the ordination of entries (as indicated by the lengths of their vectors along the first two Axis).
The variables were clustered into four groups: the first one included YLD, SW and HeW, the second one included T-YLD and HeH, the third one PD and V, while the last group included all the other variables. Within each group, variables were highly positively correlated, but the variables of the first group were negatively correlated to those belonging to the third group. The variables of the second group were not correlated with those included in the third group, while they were positively correlated with those belonging to the first group (Fig. 1).
Genetic characterisation and differentiation
Averages and standard errors relative to the main genetic parameters calculated by entry over the 22 SSR loci are reported in Table 2.
All the entries were characterised by similar number of alleles (N). In general, Na was higher than Ne in all the entries. In populations, Na ranged from 2.45 (±0.205, Syn2-TER) to 3.50 (±0.327, LR) and Ne ranged from 1.78 (±0.126, Syn2-FR) to 2.38 (±0.188, Syn1-PG), hybrids generally had lower values of Na and Ne than the populations.
Populations did not differ for Ho and He, while obviously differed from hybrids.
Fixation index (F) values were close to 0 for all the populations. The existence of negative values in some populations and in all the hybrids hinted to a heterozygosity excess.
Finally, all Fst values obtained from pairwise comparisons were highly significant (P < 0.01), with the only exception of the comparison between LR and Syn1-PG (Table 3). Considering two years of multiplication, significant pairwise Fst comparisons between the LR and derived populations increased from 0.026 to 0.066. In addition, a crescendo of differentiation was found between Syn1-PG and populations derived by multiplication in other environments (from 0.039 between Syn1-PG and Syn2-GR, to 0.070 between Syn1-PG and Syn2-TER).
Similarly, a crescendo of differentiation was found between Syn2-PG and the same generation of multiplication in different environments (from 0.046 between Syn2-PG and Syn2-GR to 0.103 between Syn2-PG and Syn2-TER). The highest differentiation value was found between Syn2-FR and Syn2-TER, populations which were multiplied in very contrasting environments. As expected, the genetic differentiation between hybrids and between populations and hybrids was very high.
AMOVA (not shown) worked out on the total data showed that among- and within- population diversity accounted for 33 and 67% of total genetic variation, respectively. However, when it was worked out only considering the LR and its derived populations, the proportions were 11 and 82%, respectively.
NJ tree results showed the presence of two main clusters (Additional file 4: Figure S3). A cluster included most of the individuals belonging to populations while the other one to the hybrids. The first cluster showed others, although not completely differentiated, sub-clusters: Syn2-TER, Syn2-FR and Syn2-UK genotypes mostly clustered together. The within-population genetic variation was evident in populations, while hybrids generally showed a high genetic uniformity in accordance with their genetic structure.
The first three components of the PCoA accounted for 85.90% of genetic variation (Fig. 2). PCo1 (52.94%) clearly separated populations from hybrids, while PCo2 (26.47%) separated the hybrid H from the other entries. PCo3 (6.49%) mainly discriminated Syn2-FR and Syn2-TER from Syn2-UK and Syn2-GR.
Comparison of populations multiplied for up to three cycles in Central Italy
Morpho-phenological and genetic characterisation
The ANOVA analysis showed that the assessed entries differ for many of the morpho-phenological traits recorded (Additional file 5: Table S2). All populations were clearly discriminated from HH, HHH and HHHH. However, populations multiplied at the same site, with very few exceptions, did not show much differences and, generally, all populations obtained after three cycles of multiplication showed similar values to those found in the populations obtained after two cycles of multiplication. Worth of note is only that Syn2-GR was different from Syn3-GR and Syn2-TER from Syn3-TER for HeH, Syn2-PG was different from the LR and the other populations multiplied in Perugia for DH, DM and PD and, finally, Syn3-PG was different from LR, Syn1-PG and Syn2-PG for V.
The PCA (with the first three axes accounting for 96.54% of total data variation) and the biplot showed that the three hybrids and hybrid H with populations clustered separately and that the variable clustering and relationship were similar to what previous observed (i.e. considering the morpho-phenological comparisons of populations multiplied for up two cycles in different environments, Additional file 6: Figure S4).
Averages and standard errors relative to the main genetic parameters calculated by entry over the 22 SSR loci are reported in Additional file 7: Table S3. A part differences between entries, populations obtained with three years of multiplication did not significantly differ from those obtained with two years of multiplication for all genetic parameters, with the exception of Na observed in Syn2-GR and Syn3-GR. Finally, pairwise Fst comparisons (Table 3) ranged from 0.005 to 0.159 between populations. Population differentiation increased across locations and generations of multiplication. The highest differentiation was found between Syn2-TER and Syn3-TER.
NJ tree results, beside a clear differentiation between populations and hybrids, showed that Syn2-TER, Syn3-TER, Syn3-GR, and partially Syn3-PG appeared to belong to four different and distinguished clusters (Additional file 8: Figure S5). The first three components of the PCoA accounted for 82.63% of genetic variation (Additional file 9: Figure S6). Populations and H were clearly discriminated from HH, HHH and HHHH. Finally, Syn3-GR, Syn2-TER and Syn3-TER were clearly differentiated from all the other populations.
To the best of our knowledge this is the first insight into the morphological and genetic short-term evolution and local adaptation of broccoli populations developed by a EBP for sustainable agriculture.
The populations here studied could be compared to natural populations, that are mainly represented by genotypes derived from the most fit plants and that mainly evolve in relationship to environmental pressures, because the EBP applied (i.e. the choice of the most vigorous plants for reproduction) favoured plants with the largest heads and then the greatest production of seed in each environment.
Although several experiences of evolutionary breeding for autogamous species have been reported by different authors ([8, 12] and references therein, [33, 34]), a few EBP have been registered for allogamous species like broccoli [23, 35, 36].
As for the comparison of populations developed in different countries, the selection pressures (those applied by different researchers by choosing three plants during multiplication of populations and those linked to different pedo-climatic conditions) did not seem to have heavily changed the morpho-phenological traits and the performances of populations in comparison with LR from which they derived. The populations showed few morphological differences probably due to their intrinsic plasticity, the initial genetic diversity of LR and the mating system of the species, factors that, combined, maintain a large variability . In addition, it is possible that the most vigorous chosen plants per generation had the highest heterozygosity and different alleles were consequently maintained across generations. Similar results were also found by [35, 36] who evaluated the morphological changes in spinach and maize in few years of multiplication in different environments, respectively.
Only populations multiplied in TER were morphologically well distinguished from the other populations. TER is an environment that showed minimum temperatures lower than the others among December and March: these environmental factors could have negatively influenced the vegetative stage of plants and consequently the reproductive stage, since from December to March the heads appear and then proceed to maturation.
Also, the EBP here described (selection and mating of the best tree plants for two years) did not seem to have resulted in any inbreeding depression, suggesting that the genetic diversity remained high.
In spite of that, genetic data showed a differentiation among populations that morphological data did not assess. Only two years of multiplication caused a low to moderate differentiation. Worth of note is the progressive increase in differentiation values observed passing from different population generations and from a multiplication environment to another.
Of course, a longer breeding period or a different selection method than the one implemented in this study could have allowed a greater differentiation among populations, as it happens in the process of broccoli LR differentiation . For example, as suggested by results of , in which the final versions of two maize populations showed low differentiation from the original populations, mass selection resulted less selective than a RRS implemented for a similar period . In addition, in both maize mass and RSS, a long breeding program increased the genetic differentiation among original and final breeding populations [18, 37], probably as a consequence of allele fixation [18, 38]. In fact, if 2–4 cycles of RRS slightly reduce diversity, with an allele loss up to 11% [20, 22, 39], 11–12 cycles of RRS reduce diversity much more, with an allele loss up to 37% [17, 40].
Our genetic results were similar to those found by  (comparing two open-pollinated maize populations undergone to stratified mass selection for 19 and 25 years) and to those found by  who (comparing two maize synthetic varieties after two cycles RRS) found that the genetic differentiation among versions of the same population was very low.
In this study, the short selection period and/or the low pressure breeding method resulted in the maintenance of a high intra-population genetic diversity.
Also in maize a study of molecular changes during intra- and inter- recurrent selection on one adapted and one non adapted populations suggested that natural selection acted stronger than breeding in differentiating the populations obtained in the non adapted population only . Similarly, passing from a multiplication environment to another, we observed an increase of differentiation between the same population generations. Although in this study on broccoli environmental pressures did not strongly affected genetic diversity as in maize , indeed they played a role in differentiating populations.
The comparisons of populations developed in Central Italy for up to three cycles confirmed what observed comparing populations developed in different countries. The SSR data highlighted a genetic differentiation among populations that was not evident with only morphological data. With respect to populations obtained after two cycles of multiplication, one more year of multiplication determined a greater genetic differentiation among populations: genetic distance between Syn1-PG and Syn3-PG was greater than between Syn1-PG and Syn2-PG. Also, the highest Fst values were those between the final versions of populations multiplied in different environments.
In addition, genetic data showed that the environment of adaptation (Perugia, the environment in which LR is grown and the original selection was carried out in order to develop the Syn1-PG) was less selective than the others. In fact, even if the genetic distance among the Syn1-PG and the derived populations was increasing along years, the increase was larger in the TER and GR environments, which are not the adaptation environments.
During an EBP the selective pressure of environments is crucial, but time is also an important factor in driving the oscillations of the allele frequencies and the subsequent adaptation of the genotypes.
One of the aims of this study was to assess the adaptation potential of a heterogeneous variety when multiplied in different environments. Although a longer multiplication period would be necessary to have a final confirmation, the obtained data show that a high genetic diversity is maintained all over the tested multiplication environments and in spite of their pedological and climatic differences in at least three years. This diversity can then continue to evolve in response to different selection pressures in each environment in the future, as it has always happened with any introduction of variable materials in different environments in the past.
Few years of multiplication in different environments resulted in differentiated populations which, nonetheless, preserved the original genetic diversity and productivity level and appear able to evolve in relationship to different environments.
As such the EBP here described also appears to be a good method to develop broccoli materials (populations) for sustainable agricultural systems, those systems where the use of external input (pesticides, irrigation, fertilisers that help a crop facing adversities) is limited or absent.
In sustainable agricultural systems, protection from biotic and abiotic stress (which are variable from year to year and from a situation to another) should rely on the diverse response of the single different genotypes at a certain moment and, consequently, on the ability of the crop, as a whole, to change its composition across time, making the crop community diverse and resistant/tolerant across time and situation. This changing ability is also important because sustainable systems often have very specific characteristics in term of soil fertility, type of occurring biotic and abiotic stress, agronomic and protection techniques applied by farmers (also in relation to their own skills and preferences). In most cases none is similar to another. Finally, it is also important because we are living in a period of rapid and unpredictable climatic changes, with the connected changes of biological forms negatively affecting crops.
For all these reasons, in sustainable agriculture populations like those developed in this study, which are able to evolve across time and locations, appear to be much more suitable materials than hybrids and greater efforts than presently carried out should be applied to develop and diffuse them. Of course, where external inputs can counteract biotic and abiotic constrains and the cultivation environment is uniform, hybrid varieties can be the most useful materials.
Materials and methods
Plant breeding procedure
Aiming to develop materials suitable for organic agriculture (i.e. possessing the adaptation and the intrinsic diversity needed to suit a particular environment) we specifically focused our work on a landrace (LR) instead of on other materials.
We chose a sprouting broccoli LR from Umbria (Central Italy) that is named “Zolfino”, highly appreciated on the local market and generally grown under low-input conditions. The LR is characterised by late flowering and maturation: requiring low winter temperatures for flower induction, it is harvested in a period spanning from beginning of March to mid-April. Initially the main head, then the lateral buds are harvested in the period.
The initial material of LR was obtained by a farmer family which cultivated it since generations. After the farmer abandoned his own population, this Department (DSA3) continued its cultivation by yearly reproducing the seed and applying the selection procedures of the mother plants (MP) suggested by the farmer: the most vigorous 2 or 3 MP, among the 25 usually grown, are chosen and intercrossed, while the field is cleared from other B. oleracea plants, since all the species varieties are interfertile. It should be noted that this procedure is common across Europe for farmers reproducing their B. oleracea LR (V. Negri personal communication).
Starting a breeding program to develop an initial synthetic in 2009, we asked the farmer to choose 17 MP among the 25 grown instead of the usual number (Fig. 3). The following year, 10 plants for each of the 17 MP half sib progenies (MPHS) were grown in a trial (5 plants for each MPHS in each of 2 replicates, for a total of 170 plants). Such a high number of plants was adopted in order to maintain the initial LR diversity as much as possible. Among them, the farmer was asked to select, according to his personal opinion, the best 8 MPHS and then, within each of them, the best 5 plants. Accordingly, a synthetic variety (Syn0) made of 40 plants from the initial 8 MPHS, that mirrors the original LR, was obtained.
The multiplication of Syn0 produced a Syn1-PG (2009–2010), the initial (foundation) material of this study.
To start the EBP the Syn1-PG seed was then multiplied in isolation in different areas of Central Italy, in France and in United Kingdom as follows (Fig. 3):
for two consecutive years at DSA3 (Perugia, Umbria Region) in the same area of its origin giving rise to two populations named Syn2-PG (2010–2011) and Syn3-PG (2011–2012), respectively;
for two consecutive years in two other areas of Central Italy with different pedo-climatic conditions from Perugia: Grosseto (Tuscany) and Terminillo (Lazio), giving the populations named Syn2-GR, Syn3-GR and Syn2-TER and Syn3-TER, respectively;
for one year only (2010–2011) in United Kingdom giving a population named Syn2-UK;
for one year only (2010–2011) in France giving a population named Syn2-FR.
Each multiplication environment (GR, TER, UK and FR) is different from PG for several pedo-climatic traits (Additional file 1: Figure S1, ). Such different environmental conditions were purposely chosen to assess the degree of restriction in genetic diversity due to selective pressure.
As it is usual for this type of LR, in all sites the seeds obtained from each multiplication cycle was sowed in pots at the beginning of July and the plantlets were transplanted in the experimental field at the beginning of September. The seed was harvested in June in Italy, and about a month later in France and United Kingdom.
In each multiplication cycle and in each environment only the three most vigorous plants (in agreement with farmer indications) were intercrossed in order to obtain seed. The inflorescences of the other plants were cut before flowering in order to evaluate agronomic traits of the commercial heads.
Each year, the obtained seed was sent to DSA3 in order to store it in its Genebank (FAO code ITA363) for the successive analyses (see below).
For each environment of multiplication and for the environment of final agronomic trial, climatic data, i.e. average daily temperature, average of maximum daily temperatures, average of minimum daily temperatures (°C) and total rainfall precipitations (mm) (obtained from near-by weather stations), were recorded from sowing to harvest (Additional file 1: Figure S1).
Evaluated plant material
The materials obtained from multiplication in different environments was used to assess the occurrence morpho-phenological and genetic changes.
Firstly, particularly to understand the effects of multiplication in different countries, we assessed the diversity of genotypes belonging to 11 entries of broccoli: seven populations and four hybrids used as controls. Specifically: the Umbrian LR, the Syn1-PG, its five two year multiplications, Syn2-PG, Syn2-GR, Syn2-TER, Syn2-UK and Syn2-FR, and the hybrids Santee (Elsomseed), Ironman (Monsanto), Belstar (Bejo Seeds) and Natura Vallata (Cooperativa Agricola Cesenate), named H, HH, HHH and HHHH, respectively.
Secondly, we assessed other three populations obtained in Central Italy with a further multiplication cycle (Syn3-PG, Syn3-TER and Syn3-GR) in order to obtain information about further evolutionary processes eventually occurred.
Assessment of morpho-phenological changes occurred with multiplications
In order to assess eventual morpho-phenological changes, all the obtained multiplications were evaluated in the same environment.
Seeds from the 14 broccoli entries were sown in pots at the beginning of July and the obtained plantlets were transplanted in a DSA3 experimental field in Perugia (Umbria, Italy) at the beginning of September. A total of 448 individual plants were evaluated. Plants were arranged by using a randomised block design with four replications. In each plot, eight plants of each entry were grown, for a total of 32 plants evaluated per entry. Spacing among plants was set at 1 × 1 m because the plants are very vigorous and this spacing is generally used by farmers growing the “Zolfino” LR (Additional file 2: Figure S2). The management practices applied to the agronomical trial were similar to those used in organic agriculture conditions in the area. Fertilization was ensured by pre-implant livestock manuring, by using 1 t ha− 1 of the commercial product NUTEX LETAME (Agroqualità, 3% N and 3% P2O5 content) that was incorporated by using a rotary cultivator. As above, this is the fertilization level usually applied by farmers. Weed control was always performed by hand-weeding and no chemical control of pathogens and insects was applied. Broccoli were irrigated immediately after transplanting and during the following month, aiming at fully restoring crop evapotranspiration.
Twelve morpho-phenological traits were recorded for each plant on all entries (Table 4).
To determine the significance of the sources of variation, the recorded data were processed by analysis of variance (ANOVA) using a linear model, where an individual trait value Tij of the levels i of the effect “Genotype” (G), and j of the effect “Block” (B), is: Tij = m + Gi + Bj + eij, where m is the grand mean and eij is the experimental error, and means were separated using the Tukey Honest Significant Difference (HSD), as implemented in the HSD.test() function of the “agricolae” package in R [41, 42].
In order to study the correlation pattern of the considered morpho-phenological traits and the differences among entries, a principal components analysis (PCA) was performed by using the function PCA() of the “FactoMineR” package . PCA is a factor analytic technique that is used for the ordination of observations in a reduced rank space. It is used to visualise the weight of variables in determining entry ordination (as shown by vector lengths, by angles between vectors and by the reciprocal positioning of variable vectors and entry symbols). As a preliminary step before performing a PCA on the complete dataset, the missing values were imputed with the Principal Components Analysis model by the function imputePCA() of the “missMDA” package .
Assessment of genetic changes occurred with multiplications
Genomic DNA extraction
Fifty mg of fresh leaf tissue were collected from each of the 448 plants (32 plants for each of 14 entries) and high quality DNA was isolated using the DNeasy 96 Plant Kit (Quiagen). DNA quality and concentration were checked by spectrophotometry using NanoDrop 2000 (Thermo Scientific) and 1.0% (w/v) agarose gel separation.
Broccoli individuals were genotyped using 17 putative neutral microsatellite markers (SSR) [45,46,47] and five microsatellites related to genes involved in flowering control and cold stress response (EST-SSR)  (Additional file 3: Table S1). According to their position on B. oleracea chromosomes, the selected markers are not in linkage (i.e. for each chromosome, the two selected markers are located on the two different arms). All the DNA amplifications were carried out as described in . Obtained amplicons were separated in a 3130xl capillary sequencer (Applied Biosystems), sized according to the internal size standard GeneScanTM23 500LIZ® (Applied Biosystems), visualised and scored using the GENEMAPPER software (Applied Biosystems).
For each SSR marker the number of alleles and allele size range were recorded. The number of successfully analysed (N) and observed (Na) genotypes, the effective (Ne) number of observed alleles, observed (Ho), as well as the expected (He) heterozygosity and the Fixation Index (F), were worked out by entry (i.e. LR, Syn1-PG, its derived populations and hybrid controls) with the use of GENALEX software . For each parameter, the average values were compared by using the Tukey Honest Significant Difference (HSD) as implemented in the HSD.test() function of the “agricolae” package in R [41, 42].
Pairwise comparison Fst values and two Genetic Distance (GD) matrices between individuals following the formula of  were calculated using GENALEX software . Two GD-based Neighbour Joining (NJ) trees were drawn using the on-line software iTOL . To better understand the genetic relationship among entries, two Principal Coordinate Analysis (PCoA) were also worked out by using the same distance matrices and drawn by “scatterplot3d” package in R [42, 54].
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Dipartimento di Scienze Agrarie, Alimentari e Ambientali
Evolutionary Breeding Program
Composite Cross Populations
Reciprocal Recurrent Selection
Simple Sequence Repeats
Restriction Fragment Length Polymorphism
Single Nucleotide Polymorphism
Mother Plant Half Sib
United Kingdom Environment
Days to heading
Days to maturity
Principal components analysis
Principal Coordinate Analysis
Number of successfully analysed genotypes
Number of successfully genotypes
Effective number of observed alleles
FAOSTAT. FAOSTAT crops production data. 2018. http://www.fao.org/faostat/en/#data/QC.
SINAB. Superfici biologiche per anno, prodotto e regione (ettari) - Cavolfiore e broccoli. 2018. http://www.sinab.it/content/superfici.
Hammer K, Gladis T, Laghetti G, Pignone D. The wild and the grown-remarks on Brassica. Int J AgriSci. 2013;3:453–80 http://www.cabdirect.org/abstracts/20133227404.html. Accessed 13 May 2014.
Lammerts van Bueren ET, Struik PC, Jacobsen E. Ecological concepts in organic farming and their consequences for an organic crop ideotype. Neth J Agric Sci. 2002;50:1–26. https://doi.org/10.1016/S1573-5214(02)80001-X.
Torricelli R, Ciancaleoni S, Negri V. Performance and stability of homogeneous and heterogeneous broccoli (Brassica oleracea L. var. italica Plenck) varieties in organic and low-input conditions. Euphytica. 2014;199.
Lammerts van Bueren ET, Jones SS, Tamm L, Murphy KM, Myers JR, Leifert C, et al. The need to breed crop varieties suitable for organic farming, using wheat, tomato and broccoli as examples: A review. NJAS - Wageningen J Life Sci. 2011;58:193–205. https://doi.org/10.1016/j.njas.2010.04.001.
Ceccarelli S. Efficiency of plant breeding. Crop Sci. 2015;55:87. https://doi.org/10.2135/cropsci2014.02.0158.
Dawson J, Goldringer I. Breeding for genetically diverse populations: variety mixtures and evolutionary populations. In: Lammerts Van Bueren ET, Myers J, editors. Chichester: Wiley-Blackwell; 2012. p. 77–98. https://doi.org/10.1002/9781119945932.ch5.
Döring TF, Annicchiarico P, Clarke S, Haigh Z, Jones HE, Pearce H, et al. Comparative analysis of performance and stability among composite cross populations, variety mixtures and pure lines of winter wheat in organic and conventional cropping systems. Field Crop Res. 2015;183:235–45. https://doi.org/10.1016/j.fcr.2015.08.009.
Finckh M, Brumlop S, Goldringer I, Steffan P, Wolfe MS. Maintenance of diversity in naturally evolving composite cross wheat populations in Europe. Collected papers of the 1st IFOAM conference on organic animal and plant breeding. Tholey-Theley. 2009. p 145–52.
Murphy K, Lammer D, Lyon S, Carter B, Jones SS. Breeding for organic and low-input farming systems: an evolutionary–participatory breeding method for inbred cereal grains. Renewable Agric Food Syst. 2005;20:48–55. https://doi.org/10.1079/RAF200486.
Phillips SL, Wolfe MS. Evolutionary plant breeding for low input systems. J Agric Sci. 2005;143:245–54.
Raggi L, Ceccarelli S, Negri V. Evolution of a barley composite-derived population: an insight gained buy molecular markers. J Agric Sci. 2015:1–7.
Goldringer I, Prouin C, Rousset M, Galic N, Bonnin I. Rapid differentiation of experimental populations of wheat for heading time in response to local climatic conditions. Ann Bot. 2006;98:805–17.
Kutka F, Smith M. How many parents give the highest yield in predicted synthetic and composite populations of maize? Crop Sci. 2007;47:1905–13.
Koutis K. Multienvironmental evaluation of wheat landraces by GGE Biplot analysis for organic breeding. Agric Sci. 2012;3:66–74. https://doi.org/10.4236/as.2012.31009.
Labate JA, Lamkey KR, Lee M, Woodman WL. Temporal changes in allele frequencies in two reciprocally selected maize populations. Theor Appl Genet. 1999;99:1166–78.
Labate JA, Lamkey R, Lee M, Woodman WL. Molecular genetic diversity after reciprocal recurrent selection in BSSS and BSCBI maize populations. Crop Sci. 1997;37:416–23.
Rossini Pinto L, Carneiro Vieira ML, Lopes De Souza C, Pereira De Souza A. Genetic-diversity assessed by microsatellites in tropical maize populations submitted to a high-intensity reciprocal recurrent selection. Euphytica. 2003;134:277–86.
Ordas B, Malvar RA, Díaz R, Butron A. Molecular changes in two maize (Zea mays L.) synthetics after reciprocal selection with two alternative methods. Mol Breed. 2015;35(4):111.
Peña-Asin J, Álvarez A, Ordas B. Molecular changes during intra and inter recurrent selection of two populations of maize: one adapted and one non adapted to the selection environment. Euphytica. 2013;193:359–67.
Kolawole AO, Menkir A, Gedil M, Blay E, Ofori K, Kling JG. Genetic divergence in two tropical maize composites after four cycles of reciprocal recurrent selection. Plant Breed. 2017;136:41–9.
Alves ML, Belo M, Carbas B, Brites C, Paulo M, Mendes-Moreira P, et al. Long-term on-farm participatory maize breeding by stratified mass selection retains molecular diversity while improving agronomic performance. Evol Appl. 2018;11:254–70.
Enjalbert J, Goldringer I, Paillard S, Brabant P. Molecular markers to study genetic drift and selection in wheat populations. Comptes Rendus Biologies. 1999;50:283–90.
Enjalbert J, Boeuf C, Belcram H, Leroy P. Use of multiparental inbred populations to determine allelic relationships of molecular markers. Plant Breed. 1999;118:88–90.
Rousselle Y, Thomas M, Galic N, Bonnin I, Goldringer I. Inbreeding depression and low between-population heterosis in recently diverged experimental populations of a selfing species. Heredity. 2011;106:289–99. https://doi.org/10.1038/hdy.2010.72.
Ciancaleoni S, Onofri A, Torricelli R, Negri V. Broccoli yield response to environmental factors in sustainable agriculture. Eur J Agron. 2016;72:1-9.
Sahamishirazi S, Moehring J, Zikeli S, Fleck M, Claupein W, Graeff-hoenninger S. Agronomic performance of new open pollinated experimental lines of broccoli (Brassica oleracea L. var. italica) evaluated under organic farming. PLoS One 2018;13:1–13. https://doi.org/https://doi.org/10.1371/journal.pone.0196775.
Renaud ENC, Lammerts Van Bueren ET, Myers JR, Paulo MJ, Van Eeuwijk FA, Zhu N, et al. Variation in broccoli cultivar phytochemical content under organic and conventional management systems: implications in breeding for nutrition. PLoS One. 2014;9(7).
Ciancaleoni S, Raggi L, Negri V. Genetic outcomes from a farmer-assisted landrace selection programme to develop a synthetic variety of broccoli. Plant Genetic Resources: Characterisation and Utilisation, vol. 12; 2014.
Branca F, Chiarenza GL, Cavallaro C, Honghui G, Zhenqing Z, Tribulato A. Diversity Sicilian broccoli (Brassica oleracea var. italica) and cauliflower (Brassica oleracea var. botrytis) landraces and their distinctive bio-morphological, antioxidant, and genetic traits. Genet Resour Crop Evol 2018; 65:485–502. https://doi.org/https://doi.org/10.1007/s10722-017-0547-8.
Li Z, Mei Y, Liu Y, Fang Z, Yang L, Zhuang M, Zhang Y, Lv H. The evolution of genetic diversity of broccoli cultivars in China since 1980. Sci Hortic 2019;250:69–80. https://doi.org/https://doi.org/10.1016/j.scienta.2019.02.034.
Raggi L, Ciancaleoni S, Torricelli R, Terzi V, Ceccarelli S, Negri V. Evolutionary breeding for sustainable agriculture: selection and multi-environmental evaluation of barley populations and lines. Field Crop Res. 2017;204:76–88.
Suneson C. An evolutionary plant breeding method. Agron J. 1956;48:188–91.
Serpolay E, Schermann N, Dawson J, Lammerts van Bueren ET, Goldringer I, Chable V. Phenotypic changes in different spinach varieties grown and selected under organic conditions. Sustainability. 2011;3:1616–36.
Serpolay-Besson E, Giuliano S, Schermann N, Chable V. Evaluation of evolution and diversity of maize open-pollinated varieties cultivated under contrasted environmental and farmers’ selection pressures: a phenotypical approach. Open J Genet. 2014;4:125–45.
Zou CY, Li LJ, Yang KC, Pan GT, Rong TZ. Effects of mass selection on maize synthetic populations. Acta Agron Sin. 2010;36:76–84. https://doi.org/10.1016/S1875-2780(09)60028-6.
Hagdorn S, Lamkey KR, Frisch M, Guimarães PEO, Melchinger AE. Molecular genetic diversity among progenitors and derived elite lines of BSSS and BSCB1 maize populations. Crop Sci. 2003;43:474–82.
Butrón A, Tarrío R, Revilla P, Ordás A, Malvar RA. Molecular changes in the maize composite EPS12 during selection for resistance to pink stem borer. Theor Appl Genet. 2005;110:1044–51.
Solomon KF, Martin I, Zeppa A. Temporal genetic structure patterns in tropical maize populations under reciprocal recurrent selection. Euphytica. 2010;176:239–49.
De Mendiburu F. Una herramienta de analisis estadistico para la investigacion agricola. PhD Thesis. Universidad Nacional Agraria La Molina (UNI-PERU); 2009.
R Development Core Team. R: a language and environment for statistical computing. Vienna. Retrieved from http://www.R-project.org: R Foundation for Statistical Computing; 2013.
Husson F, Josse J, Le S, Mazet J. FactoMineR, multivariate exploratory data analysis and data mining with R. R package version 12 Retrieved from http,//CRANR project.org/package = FactoMineR. 2013.
Josse J, Husson F. missMDA: a package for handling missing values in multivariate data analysis. J Stat Softw. 2016;70:1–31.
Lowe A, Moule C, Trick M, Edwards KJ. Efficient large-scale development of microsatellites for marker and mapping applications in Brassica crop species. Theor Appl Genet. 2004;108:1103–12. https://doi.org/10.1007/s00122-003-1522-7.
Cheng X, Xu J, Xia S, Gu J, Yang Y, Fu J, et al. Development and genetic mapping of microsatellite markers from genome survey sequences in Brassica napus. Theor Appl Genet. 2009;118:1121–31. https://doi.org/10.1007/s00122-009-0967-8.
Li H, Chen X, Yang Y, Xu J, Gu J, Fu J, et al. Development and genetic mapping of microsatellite markers from whole genome shotgun sequences in Brassica oleracea. Mol Breed. 2011;28:585–96. https://doi.org/10.1007/s11032-010-9509-y.
Aksoy S, Almeida-val VMF, Azevedo VCR, Baucom R, Bazaga P, Beheregaray LB, et al. Permanent Genetic Resources added to Molecular Ecology Resources Database 1 October 2012–30 November 2012. Mol Ecol Resour Note. 2013;13:341–3.
Ciancaleoni S, Chiarenza GL, Raggi L, Branca F, Negri V. Diversity characterisation of broccoli (Brassica oleracea L. var. italica Plenck) landraces for their on-farm (in situ) safeguard and use in breeding programs. Genet Resour Crop Evol. 2013;61:451–64. https://doi.org/10.1007/s10722-013-0049-2.
Peakall R, Smouse PE. Genalex 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95. https://doi.org/10.1111/j.1471-8286.2005.01155.x.
Michalakis Y, Excoffier L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics. 1996;142:1061–4.
Nei M. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics. 1978;89:583–90.
Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44:W242–5.
Ligges U, Machle M. Scatterplot3d - an R package for visualizing multivariate data. J Stat Softw. 2003;8:1–20.
Thanks are due to the personnel of the Centro Appenninico del Terminillo of Università degli Studi di Perugia and in particular Dr. Valeriano Valeriani, of The Organic Research Centre (ORC) and of Istitute Technique de l’Agriculture Biologique (ITAB) for multiplying of populations in Central Italy, United Kingdom and France, respectively.
The research leading to these results was funded by the European Community’s Seventh Framework Programme (FP7/2007–2013) under the grant agreement n° 245058 SOLIBAM (Strategies for Organic and Low-Input Integrated Breeding and Management).
The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Monthly precipitation (mm) and mean of minimum, maximum and average temperature (°C). Data registered during the multiplication seasons in Perugia in 2009/2010 (PG 2009/2010), in Central Italy (Perugia, Grosseto e Terminillo) in 2010/2011 and 2011/2012 (PG 2010/2011, PG 2011/2012, GR 2010/2011, GR 2011/2012, TER 2010/2011 and TER 2011/2012), in United Kingdom in 2010/2011 (UK 2010/2011) and in France in 2010/2011 (FR 2010/2011), and during the final agronomic trial carried out in Perugia (PG 2013/2014). The bar chart refers to monthly precipitation, lines with squares, diamonds and circles refer to maximum, average and minimum temperature, respectively.
Final agronomic trial carried out in the DSA3 experimental field in Perugia (Umbria, Italy). Broccoli multiplications carried out across Europe and in Central Italy were grown in the DSA3 experimental field in Perugia. The populations tested were: LR, Syn1-PG, Syn2-PG, Syn3-PG, Syn2-GR, Syn3-GR, Syn2-TER, Syn3-TER, Syn2-UK, Syn2-FR. H, HH, HHH and HHHH are the hybrid varieties used as controls.
SSR markers used for analyze genetic diversity of broccoli original LR, Syn1-PG, its derived populations and hybrid controls. Bibliographic reference codes, genebank entry, primer sequences, repeated motif, linkage group (LG), and band range (in base pairs) relative to the 22 microsatellites used (1Love et al. 2004; 2Li et al. 2011; 3Cheng at al. 2009; 4Aksoy et al. 2013).
Neighbour-joining tree of genetic distances of the entries (two years of multiplication). B. oleracea var. italica original LR, Syn1-PG, its derived population by two years of multiplication across Europe and hybrid controls are coded with abbreviations and coloured with different colours as follows: LR-black, Syn1-PG-red, Syn2-PG-yellow, Syn2-GR-purple, Syn2-TER-dark green, Syn2-UK-blue, Syn2-FR-orange, H-violet, HH-brown, HHH-pink, HHHH-saddlebrown. Each individual of each entry are code whit a different number.
Average values and standard deviations of morpho-phenological traits for the entries (three years of multiplication). B. oleracea var. italica original LR, Syn1-PG, its derived populations by three years of multiplication in Central Italy and hybrid controls. Averages followed by the same letters are not significantly different at P < 0.05 (Tukey HSD).
PCA based on morpho-phenological traits recorded for each entries (three years of multiplication). B. oleracea var. italica original LR, Syn1-PG, its derived populations by three years of multiplication in Central Italy and hybrid controls. Coloured symbols represent the different entries analysed as from the legend, grey dotted arrow represent the morpho-phenological traits coded as from Table 4.
Averages and standard errors relative to the main genetic parameters calculated by entry (three years of multiplication). Number of successfully analysed genotypes (N), observed (Na), and effective alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), and Fixation Index (F) worked out for original LR, Syn1-PG, its derived populations by three years of multiplication in Central Italy and hybrid controls on the basis of data recorded overall for the 22 markers used.
Neighbour-joining tree of genetic distances of the entries (three years of multiplication). B. oleracea var. italica original LR, Syn1-PG, its derived populations by three years of multiplication in Central Italy and hybrid controls are coded with abbreviations and coloured with different colours as follows: LR-black, Syn1-PG-red, Syn2-PG-yellow, Syn3-PG-light green, Syn2-GR-purple, Syn3-GR-grey, Syn2-TER-dark green, Syn3-TER-light blue, H-violet, HH-brown, HHH-pink, HHHH-saddlebrown. Each individual of each entry are code whit a different number.
PCoA based on average genetic distance of the entries (three years of multiplication). B. oleracea var. italica original LR, Syn1-PG, its derived populations by three years of multiplication in Central Italy and hybrid controls. Coloured symbols represent the different entries analysed as from the legend.
About this article
Cite this article
Ciancaleoni, S., Negri, V. A method for obtaining flexible broccoli varieties for sustainable agriculture. BMC Genet 21, 51 (2020). https://doi.org/10.1186/s12863-020-00846-2