A method for obtaining flexible broccoli varieties for sustainable agriculture

doi:10.21203/rs.2.12167/v1

Download PDF

Research article

A method for obtaining flexible broccoli varieties for sustainable agriculture

https://doi.org/10.21203/rs.2.12167/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 07 May, 2020

Read the published version in BMC Genetics →

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

Background

The use of high inputs in agriculture has resulted in few varieties (hybrids and pure lines) used in all agricultural systems. Also varieties of vegetables, including broccoli (Brassica oleracea L.), for organic and low-input agriculture, are almost exclusively hybrids, due to lack of specific breeding programs and varieties for sustainable agriculture.

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). The heterogeneous varieties, in fact, are able to evolve and adapt to specific agro-climatic conditions.

Aim of this study was to develop a method (an Evolutionary Breeding Program, EBP) for obtaining heterogeneous varieties and test its efficiency in giving rise to 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 was assessed.

Results

The resulted presented 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 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 with morphological data.

Conclusions

Few years of multiplication in different environments resulted in genetically differentiated broccoli populations that none the less 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.

Population Genetics

Molecular Genetics

Medical Genetics

Brassica oleracea L.

Evolutionary Breeding Program

populations

sustainable agriculture

genetic diversity

adaptation

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 [1]. 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 [1] of which 10% is obtained in organic agriculture [2]. Among the different B. oleracea subspecies, broccoli (spp. capitata L. (DC.) convar. botrytis (L.) Alef. var. italica Plenck [3]) is an important source of vitamins, minerals and antioxidants.

Presently, organic farmers depend largely on broccoli varieties bred for high external input in conventional farming systems [4, 5], and in particular, broccoli varieties are almost exclusively hybrids [5], 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 (like 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 [7]. 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 breeding programs, 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–13].

The success of an EBP, like 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 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 [15], 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 which 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–21] and Single Nucleotide Polymorphism (SNP) [22]. SSR were also used in order to evaluate the genetic changes in different maize populations obtained after about 20 years of stratified mass selection on two historical Portuguese populations [23]. Considering autogamous species, the genetic evolution over time and space of wheat composite cross populations (CCP) was assessed by RFLP [24, 25] and by SSR [26] and that of a barley CCP by using SSR [13].

Agronomic performances of different broccoli heterogeneous varieties or broccoli breeding materials for sustainable agriculture have been already reported [5, 27–29], but few data are available on genetic diversity of broccoli breeding materials [30] 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.

Plant breeding procedure

The populations examined were originated from a Umbrian landrace (LR) that has been cultivated by a farmer named Vento in the Umbria Region. After the farmer abandoned the landrace, 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 MPs, among the 25 usually grown, are chosen and intercrossed, while the field is cleared from other B. oleracea plants. Starting a breeding program to developed an initial synthetic in 2009, we asked the farmer to choose 17 MP among the 25 grown instead of the usual numbers (Figure 1). 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). 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. Syn1-PG seed was then multiplied in isolation in different areas of Central Italy, in France and in United Kingdom as follows (Figure 1):

i) for two consecutive years at dsa3 (Perugia, Umbria Region) in the same area of its origin giving populations named Syn2-PG (2010-2011) and Syn3-PG (2011-2012), respectively;

ii) for two consecutive years in two other areas of Central Italy with different pedo-climatic conditions from Perugia named Grosseto (Tuscany) and Terminillo (Lazio), giving the populations named Syn2-GR, Syn3-GR and Syn2-TER and Syn3-TER, respectively;

iii) for only one year (2010-2011) in United Kingdom giving a population named Syn2-UK;

iv) for only one year (2010-2011) in France giving a population named Syn2-FR.

In each multiplication cycle and in each environment only the three best 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 head.

Each year, the seed obtained was sent to DSA3 that stored it in its Genebank for the successive analyses (see below).

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) for each environment of multiplication and for the environment of final agronomic trial (Additional File1: Figure S1).

Plant materials

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

In addition, we also assessed three other populations obtained in Central Italy with a further multiplication cycle (Syn3-PG, Syn3-TER and Syn3-GR) in order to obtain information about the further evolutionary processes eventually occurred.

Morpho-phenological characterisation

Seeds from the 14 B. oleracea 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 x 1 m because the plants are very vigorous and this spacing is generally used by farmers growing them. The management practices applied to the agronomical trial were similar to those used in organic agriculture. 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% P₂O₅ 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 1).

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 T_ijof the levels i of the effect “Genotype” (G), and j of the effect “Block” (B), is: T_ij=m + G_i + B_j + e_ij, where m is the grand mean and e_ij is the experimental error, and means were separated using the Tukey Honest Significant Difference (HSD), as implemented in the HSD.test() function of the “agricolare” package in R [31, 32].

In order to study the correlation pattern of 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 [33]. As a preliminary step before performing a PCA on an completed dataset, the missing values of the dataset were imputed with the Principal Components Analysis model by the function imputePCA() of the “missMDA” package [34].

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.

Microsatellite analysis

B.oleracea individuals were genotyped using 17 putative neutral microsatellite markers (SSR) [35–37] and five microsatellites related to genes involved in flowering control and cold stress response (EST-SSR) [38] (Additional file 2: 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 [39]. 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 [40]. For each parameter, the averages were compared by using the Tukey Honest Significant Difference (HSD) as implemented in the HSD.test() function of the “agricolae” package in R [31, 32].

In order to assess the components of molecular variation, an Analysis of Molecular Variance (AMOVA) was worked out following the method of [41] as implemented in GENALEX software [40].

Pairwise comparison Fst values and two Genetic Distance (GD) matrices between individuals following the formula of [42] were calculated using GENALEX software [40]. Two GD-based Neighbour Joining trees were drawn using the on-line software iTOL [43]. 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 [32, 44].

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 2). 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 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). (see Table 2)

Considering the inflorescence emission and maturation date (DH range = 101 - 204 dd; DM range = 1341 - 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 (Figure 2). 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 second group were not correlated with those included in the third group, while they were positively correlated with those belonging to the first group (Figure 2).

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 3.

All the entries were characterised by similar number of alleles (N). In general, Na is 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 4). 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).

All the same 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. (see Table 4)

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 3: Figure S2). A cluster included most of the individuals belonging to populations while the other one 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 (Figure 3). 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 4: 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 values similar to those found in 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 and 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 5: Figure S3)

Averages and standard errors relative to the main genetic parameters calculated by entry over the 22 SSR loci are reported in Additional file 6: Table S3. A part differences between entries, population 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 4) 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 7: Figure S4). The first three components of the PCoA accounted for 82.63% of genetic variation (Additional file 8: Figure S5). 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 which mainly evolve in relationship to environmental pressures because the EBP applied favoured the plants with largest heads and then greater 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, 45, 46], a few EBP have been registered for allogamous species like broccoli [23, 47, 48].

As for comparison of populations developed in different countries, the selection pressure applied by different researchers during multiplication of populations, by choosing and intermating three plants (following the indications of the Vento farmer) 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 to the mating system of the species that probably maintains large variability [30]. Similar results were also found by [47, 48] 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 proceed to maturation.

Also, the EBP described here (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 implemented in this study could have allowed a greater differentiation among populations. For example, as suggested by results of [23], 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 [18]. In addition, in both maize mass and RSS, a long breeding program increased the genetic differentiation among original and final breeding populations [18, 49], probably as a consequence of allele fixation [18, 50]. In fact, if 2 - 4 cycles of RRS slightly reduce diversity, with an allele loss up to 11% [20, 22, 51], 11 – 12 cycles of RRS reduce diversity much more, with an allele loss up to 37% [17, 52].

Our genetic results were similar to those found by [23] (comparing two open-pollinated maize populations undergone to stratified mass selection for 19 and 25 years) and to those found by [20] (comparing two maize synthetic varieties after two cycles RRS) who 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 non adapted population only [21]. 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 [21], indeed they played a role in differentiating population.

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 found among 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 Terminillo and Grosseto 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.

Few years of multiplication in different environments resulted in differentiated populations which, non the less, preserved the original genetic diversity and productivity level and appear able to evolve in relationship to different environments. AS such the EBP here described appears to be a good method to develop broccoli materials (populations) for sustainable agricultural systems, those systems where the use of external input (pesticides, water, fertilisers) that help facing crop advertise, is limited or absent.

Sustainable systems, like the organic 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. In addition to that, it should be noted that we are living in a period of rapid and unpredictable climatic changes with connected changes of biological forms negatively affecting crops.

As such, in sustainable agricultural systems, protection from biotic and abiotic stress (which are variable from year to year and from situation to another ) should rely on the diverse response of the single different genotypes at a certain moment and the ability of the crop as a whole to changes its composition across time, making the crop community diverse and resistant/tolerant across time and situation.

For all these reasons, 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. Hybrids can perform better than populations in “conventional” agricultural systems where external input counteract stress, but do not possess the flexibility needed in the, often so particular, sustainable systems. Neither is feasible to breed specific hybrids for each different system. Finally, populations developed through an EBP also allow farmers yearly reproducing them to have materials always tuned with their own cultivation environment.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of data and material

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare that they have no competing interests.

Funding

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).

Authors' contributions

VN conceived the research; SC performed experiments and analyzed the data; SC and VN wrote the manuscript. All authors read and approved the final manuscript.

Acknowledgment

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 multiplications of populations in Central Italy, United Kingdom and France, respectively.

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. International Journal of AgriScience. 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. Netherlands Journal of Agricultural Science. 2002;50:1–26. doi: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 Journal of Life Sciences. 2011;58:193–205. doi:10.1016/j.njas.2010.04.001.
Ceccarelli S. Efficiency of plant breeding. Crop Science. 2015;55:87. doi:10.2135/cropsci2014.02.0158.
Dawson J, Goldringer I. Breeding for genetically diverse populations: variety mixture and evolutionary populations. In: Lammerts Van Bueren ET, Myers J, editors. Organic crop breeding. John Wiley and Sons, Inc.; 2012. p. 77–98.
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 Crops Research. 2015;183:235–45. doi:10.1016/j.fcr.2015.08.009.
Finckh M, Brumlop S, Goldringer I, Steffan P, Wolfe M. Maintenance of diversity in naturally evolving composite cross wheat populations in Europe. In: Collected Papers of the 1st IFOAM International Conferences on Organic Animal and Plant Breeding. Breeding Diversity. 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 Agriculture and Food Systems. 2005;20:48–55. doi:10.1079/RAF200486.
Phillips SL, Wolfe MS. Evolutionary plant breeding for low input systems. Journal of Agricultural Science. 2005;143:245–54.
Raggi L, Ceccarelli S, Negri V. Evolution of a barley composite-derived population: an insight gained buy molecular markers. Journal of Agricultural Science. 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. Annals of Botany. 2006;98:805–17.
Kutka F, Smith M. How many parents give the highest yield in predicted synthetic and composite populations of maize? Crop Science. 2007;47:1905–13.
Koutis K. Multienvironmental evaluation of wheat landraces by GGE Biplot Analysis for organic breeding. Agricultural Sciences. 2012;3:66–74. doi:10.4236/as.2012.31009.
Labate JA, Lamkey KR, Lee M, Woodman WL. Temporal changes in allele frequencies in two reciprocally selected maize populations. Theoretical and Applied Genetics. 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 Science. 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. Molecular Breeding. 2015;35.
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 Breeding. 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. Evolutionary Applications. 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 Breeding. 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. doi:10.1038/hdy.2010.72.
Ciancaleoni S, Onofri A, Torricelli R, Negri V. Broccoli yield response to environmental factors in sustainable agriculture. European Journal of Agronomy. 2016;72.
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/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.
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. 2014;12.
De Mendiburu F. Una herramienta de analisis estadistico para la investigacion agricola. Universidad Nacional de Ingenieria (UNI-PERU); 2009.
R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Retrieved from http://www. R-project.org. 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. Journal of Statistical Software. 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. Theoretical and Applied Genetics. 2004;108:1103–12. doi: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. Theoretical and Applied Genetics. 2009;118:1121–31. doi: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. Molecular Breeding. 2011;28:585–96. doi: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. Molecular Ecology Resources 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. Genetic Resources and Crop Evolution. 2013;61:451–64. doi:10.1007/s10722-013-0049-2.
Peakall R, Smouse PE. Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006;6:288–95. doi: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 research. 2016;44:W242–5.
Ligges U, Machle M. Scatterplot3d - an R Package for visualizing multivariate data. Journal of Statistical Software. 2003;8:1–20.
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 Crops Research. 2017;204:76–88.
Suneson C. An evolutionary plant breeding method. Agronomy Journal. 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 Journal of Genetics. 2014;4:125–45.
Zou CY, Li LJ, Yang KC, Pan GT, Rong TZ. Effects of mass selection on maize synthetic populations. Acta Agronomica Sinica. 2010;36:76–84. doi: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 Science. 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. Theoretical and Applied Genetics. 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.

Additional file 1: Figure S1.pdf. 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).

Additional file 2: Table S1.pdf. 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 (¹Love et al 2004; ²Li et al. 2011; ³Cheng at al. 2009; ⁴Aksoy et al. 2013).

Additional file 3: Figure S2.tiff. Neighbour-joining tree of genetic distances of the entries (two years of multiplication).

B. oleracea 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 follow: 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.

Additional file 4: Table S2.pdf. Average morpho-phenological traits values for the entries (three years of multiplication).

B.oleracea 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).

Additional file 5: Figure S3.pdf. PCA based on morpho-phenological traits recorded for each entries (three years of multiplication).

B.oleracea 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 1.

Additional file 6: Table S3.pdf. 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.

Additional file 7: Figure S4.tiff. Neighbour-joining tree of genetic distances of the entries (three years of multiplication).

B.oleracea 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 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.

Additional file 8: Figure S5.pdf. PCoA based on average genetic distance of the entries (three years of multiplication)

Due to technical limitations, tables are only available as a download in the supplemental files section.

Download PDF

Journal Publication

published 07 May, 2020

Read the published version in BMC Genetics →

Version 1

posted

You are reading this older preprint version

Read the latest preprint version →

A method for obtaining flexible broccoli varieties for sustainable agriculture

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Materials and Methods

Results

Discussion

Conclusions

Declarations

References

Additional Files

Tables

Supplementary Files

Status:

Journal Publication

Version 1