Next Article in Journal
CO2-Induced Changes in Wheat Grain Composition: Meta-Analysis and Response Functions
Next Article in Special Issue
Drought Stress-Related Physiological Changes and Histone Modifications in Barley Primary Leaves at HSP17 Gene
Previous Article in Journal
Efficient Partitioning of Assimilates in Stress-Tolerant Groundnut Genotypes under High-Temperature Stress
Previous Article in Special Issue
Metabolite Profiling for Leaf Senescence in Barley Reveals Decreases in Amino Acids and Glycolysis Intermediates
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Statistical Crop Model to Explain the Relationship between Seed Yield and Phenotypic Diversity within the Brassica napus Genepool

1
Department of Food and Nutritional Sciences and Centre for Food Security, University of Reading, Whiteknights, PO Box 226, Reading, Berkshire RG6 6AP, UK
2
School of Mathematical Sciences, University of Nottingham, University Park, Nottingham NG7 2RD, UK
3
Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, UK
4
Warwick Crop Centre, University of Warwick, Wellesbourne CV35 9EF, UK
5
School of Agriculture, Policy and Development and Centre for Food Security, University of Reading, Earley Gate, Whiteknights Road, PO Box 237, Reading RG6 6AR, UK
6
Rothamsted Research, Harpenden, Hertfordshire AL5 2JQ, UK
7
Southern Cross Plant Science, Southern Cross University, PO Box 157, Lismore, NSW 2480, Australia
8
Office of the Vice-Chancellor, 18 Portland Villas, University of Plymouth, Plymouth, Devon PL4 8AA, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to the work.
Submission received: 30 March 2017 / Revised: 14 April 2017 / Accepted: 19 April 2017 / Published: 22 April 2017
(This article belongs to the Special Issue Senescence of Crop Plants)

Abstract

:
Plants are extremely versatile organisms that respond to the environment in which they find themselves, but a large part of their development is under genetic regulation. The links between developmental parameters and yield are poorly understood in oilseed rape; understanding this relationship will help growers to predict their yields more accurately and breeders to focus on traits that may lead to yield improvements. To determine the relationship between seed yield and other agronomic traits, we investigated the natural variation that already exists with regards to resource allocation in 37 lines of the crop species Brassica napus. Over 130 different traits were assessed; they included seed yield parameters, seed composition, leaf mineral analysis, rates of pod and leaf senescence and plant architecture traits. A stepwise regression analysis was used to model statistically the measured traits with seed yield per plant. Above-ground biomass and protein content together accounted for 94.36% of the recorded variation. The primary raceme area, which was highly correlated with yield parameters (0.65), provides an early indicator of potential yield. The pod and leaf photosynthetic and senescence parameters measured had only a limited influence on seed yield and were not correlated with each other, indicating that reproductive development is not necessarily driving the senescence process within field-grown B. napus. Assessing the diversity that exists within the B. napus gene pool has highlighted architectural, seed and mineral composition traits that should be targeted in breeding programmes through the development of linked markers to improve crop yields.

1. Introduction

Crop yield is a complex trait determined by a number of contributing environmental and genetic factors. Mathematical modelling approaches are one way of synthesising information and simplifying the complex interactions that exist throughout a plant’s life cycle to gain information about the most appropriate target traits that could increase yields. The rationale is that, if used as early predictors of final yield, models have the potential to inform growers of alterations in farming practices that, if implemented early in the season, could increase yields. Models exist for the three major cereal crops grown worldwide (rice, wheat and maize) [1,2,3], and all aim to predict yield in the face of environmental or genetic variation. Crops such as oilseed rape (canola, rapeseed, colza) are of increasing economic importance, yet they have an indeterminate growth habit compared to cereals, having been domesticated for only 4000 years compared to the 10,000 years for wheat, and therefore, require a dedicated approach in order to determine yield components that are of use to growers and breeders. Several models of yield prediction in oilseed rape have been generated previously [4,5,6,7,8]. As with all models, they have limitations; for instance, only applying data from optimum growth conditions, using just one variety of oilseed rape and/or being solely based on parameters derived from the literature. Hence, whilst modelling oilseed rape growth provides a useful tool to evaluate the parameters affecting yield, to date, this technique is still in its infancy and has not had the power of the diversity trial used in the present study to inform the modelling approach.
Oilseed rape (Brassica napus) is a crop species primarily harvested for its oil-containing seeds within temperate regions and is globally ranked as the third leading source of plant oil after soybean (Glycine max) and oil palm (Elaeis guineensis) [9]. The oil is used in the food industry, but the seed and remaining biomass material also have a number of other roles, including use as a protein meal for the feed industry and as a feedstock for biofuels. Its success as an oilseed can be attributed to the fact that its seeds are composed of ~25% protein and ~50% oil (w/w); the oil having a desirable composition of almost entirely unsaturated fatty acids [10], of which the primary components are the unsaturated fatty acids oleic acid (C18:1), linoleic acid (C18:2) and linolenic acid (C18:3; [11]. Brassica napus (AACC) is an amphidiploid species arising from a spontaneous hybridisation between B. rapa (AA) and B. oleracea (CC) [12]. Compared to other crops, such as wheat, oilseed rape has only recently undergone a domestication event [13]; this is highlighted by the increase in its cultivation as a crop between 1961 and 2013 when there was a 481% global increase in the area of oilseed rape harvested, contrasting with a 7% increase in wheat over the same period [14]. This observation could be partly attributed to the fact that oilseed rape is often used as a break crop within a wheat rotation, yet within these few decades, seed quality has already been improved by reducing the concentration of erucic acid and glucosinolates (GSLs), which were believed to be anti-nutritional seed components [15]. This may have been achieved at a cost, as one study linked Quantitative Trait Loci (QTL) for seed yield in winter oilseed rape to high Glucosinolate (GSL) concentration, revealing that reducing the latter has the potential to negatively impact seed yield [16].
As crops become domesticated, there is strong selection against the natural variation in plant development and architecture within the population, leading to increasingly uniform monocultures that are suited to mechanised harvest at a single time point. In general, the growth habit, yield, pest/pathogen defence system and response to environmental conditions throughout the growing season largely depend on the genotype of the accession sown [17].
In order to maximise crop yield, a plant within a field-based canopy has to balance the extent of vegetative growth, the length of the photosynthetic period, the number of reproductive structures in which it invests and the amount of resource (protein, lipid, carbohydrate) that it imports into each seed [18]. It is also important that the plants within a canopy retain the adaptive capacity (plasticity) to respond to environmental changes and exploit any extra resources that may accrue during the growing season [19]. The optimal idiotype of oilseed rape is yet to be defined; the crop currently produces extensive above-ground vegetative biomass and hence has a low harvest index (the ratio between harvestable yield and total plant biomass). Our hypothesis is that plants with high yield in terms of seed mass per plant achieve this as a consequence of more pods per plant rather than a positive increase in resource accumulation within individual seeds or pods. Whilst the vegetative structure of the plant contributes towards increased pod production by providing more sites at which pods can be formed, it can also reduce the amount of incident light reaching photosynthetic components lower down in the canopy. Even the relatively small petals are believed to block out ~60% of the photosynthetically-active radiation (PAR) [20].
The timing of leaf senescence is not predetermined, but is influenced by both environmental [21] and genotypic factors [22,23], which will ultimately affect both the number of fruits that develop and the extent to which the seed will fill. However, the relationship between senescence and yield is a complex one and not fully understood [24], especially in oilseed rape. The role of the pod wall extends far beyond that of a protective organ [18,25], as the pods are themselves photosynthetic units that senesce, contributing 50%–60% of the final dry mass of the mature plant [26]. It has been noted that Arabidopsis accessions that naturally senesce early have a greater number of pods per plant compared to late-senescing accessions (although the nutritional composition and weight of these seeds were not stated), indicating a greater investment in reproductive as opposed to vegetative development [22,23]. Some Arabidopsis mutants that exhibit a delayed senescence phenotype, such as the abscisic acid insensitive (abi3) mutant, develop more pods per plant compared to the wild-type [27]. This delay was predicted to enhance yields by providing a longer photosynthetic period in which to synthesise photoassimilates for subsequent re-allocation into reproductive structures. Whilst this might increase carbon assimilation in the seed, nitrogen (N) remobilisation becomes delayed in late senescing Arabidopsis phenotypes [28], so while delaying senescence might result in higher yields, the trade-off is that the seeds contain less protein [17]. Hence, when choosing ideal traits for crops, there is confusion over whether maximum yield will be achieved from fast or slow senescing lines [29].
The diversity trials performed as part of the Oilseed RapE Genetic Improvement Network (OREGIN) project [30,31] sampled in the current study provided a robust basis to begin modelling, as they represent the genetic diversity that exists across the domesticated B. napus gene pool, therefore providing a useful insight into direct future breeding programmes. The present study developed statistical models of yield across two years that were able to elucidate the most important architectural and quality traits related to yield and that provide a predictive model of yield based on parameters that could be measured prior to pod formation and then acted on to maximize yield.

2. Results

2.1. Development of an Efficient Statistical Model to Explain Yield

Using data for 45 architectural and physiological traits measured in Year 1 of the field trial, forward selection was used to develop a deterministic statistical model to explain log[yield]. Data from only the varieties present in both years of the trial were included. Yield was defined as the weight of seeds harvested per plant. The varieties were grown under two nitrogen treatments in Year 1, but only in residual soil nitrogen in Year 2; residual N was found to be 18–22 kg ha−1 in both years, and high N in Year 1 was determined as 148 kg ha−1. A two-factor ANOVA from Year 1 data to compare the significance of variety and nitrogen on yield showed that both terms are significant factors in predicting yield, with variety explaining 44.5% of the variance and nitrogen explaining 23.1% of the variance, respectively (Table S1). A comparison by ANOVA using variety and year terms showed that variety was significant (explained 37% of the variance), but year was not (explained 0.7% of the variance), thus validating our approach of comparing data across years. Variety was excluded as a significant term from the analysis in order to examine which traits consistently contribute to high or low yield across different accessions (Table 1). The model explained over 50% of the yield variance with the inclusion of five trait terms: chlorophyll a from canopy leaves, seed linolenic acid content, number of seeds per pod, the rate of N applied to the plot, glucosinolate content of the seeds. The effect of N was positive overall and confirmed that addition of N increases yield. However, N treatment (high/low) did not segregate the dataset into two distinct groups according to yield; some varieties grown on low N had higher yield than others grown on high nitrogen and vice versa (Figure 1a). A two-way ANOVA of log[yield] against variety and N showed that the size of the effect differs between varieties, and for most varieties in the present study, the effect of N on seed yield per plant is not significant (Figure 1b), although if more years/blocks were available, other varieties in the middle of the graph may also show a significant response.

2.2. Genetic Factors Are Stronger Predictors of Yield than Environmental Factors within Year

In order to explore models that could predict yield independently of N application, and in preparation for characterising Year 2 data where only one N treatment was included in the trial, statistical models were developed independently for low and high N treatments (Figure 2 and Table 2). Both models performed slightly less well than the combined model (Figure 1a). Each shared some covariates with the combined model, with seeds per pod and seed linolenic acid common between the combined and low N models and log[seed glucosinolate] common between the combined and high N models. There were no common covariates between the low and high N models. Having seed compositional traits as core terms in both the combined and separated models, as opposed to the terms being entirely of growth and development traits, suggests that the combined model is applicable to a range of growing environments. The six terms in the low N model accounted for 52.67% of the variance in yield. The low N treatment demonstrated the importance of senescence when resources are limited; a higher percentage of green pods at the stages examined showed a positive and significant (p < 0.05) relationship with yield (Table 2). Variety was excluded as a candidate for this model; however, in a separate model, where only variety was considered, it explained 44.50% of the variation in yield (across all plots), rising to 60.92% in high N plots and 78.43% in low N plots. Therefore, although N does have an impact on yield, the genetic basis of yield (i.e., variety) is a much more powerful predictor.

2.3. Developing an Improved Model to Explain Yield

Application of the Year 1 low N model to Year 2 data (the crop was only grown under low N in Year 2) showed that the model was poor at explaining yield, being able to only account for 25.48% of yield variation (Table 3, Figure 3a). This is half of the yield variation explained by the same model in Year 1 (52.67%) and demonstrates that within-year variation in yield is dominated by genetic factors, whilst between-year variation is strongly influenced by the environmental conditions to which the plants were exposed. However, although the estimates of variation for the terms used in the model are different between years, we see the same directional effects, i.e., there is a negative estimate of variation for seed linolenic acid and a positive estimate of variation for seeds per pod.
To develop the model further, a second forward selection process was performed to determine which traits best explained log[yield] in Year 2. Plants in Year 2 were grown on residual N only, and additional traits were recorded. As predicted, variety was highly significant (p = 0.009) in predicting log[yield], but as with Year 1, this was not included as a term in the following analyses in order to determine traits that contribute to the yield difference observed between varieties. Seed compositional traits, such as linolenic acid, total seed oil and seed protein content, were important factors; however, total vegetative mass (excluding seeds) was the most significant term, explaining 73.70% of log[yield] (Table 3, Figure 3b). Seed oil was shown to be a positive term, whereas protein content had a negative impact on yield, demonstrating that the oil:protein ratio is important in determining total yield. In total, the forward selection model accounted for 95.12% of the variation contributing to log[yield].
Given that total above-ground biomass at harvest was such an important predictor of log[yield], we investigated this trait further and used forward selection to identify the traits and variables that explained above-ground vegetative biomass at harvest (Table 3, Figure 3c). A consequence of including a trait in the model is that other traits correlated with it are less likely to be significant themselves and are not included in the model; hence, none of the terms explaining biomass appear in the yield model, as the biomass terms acts as a proxy for them all. The first three predictors of vegetative biomass were all architectural traits with positive effects, but these three traits only explain 77.74% of the vegetative biomass, suggesting that other traits listed in both models make a small, but significant contribution to yield. The first four terms in the biomass model (stem area, number of branches on the primary stem, plant height and pod carotenoid content at Week 41) were all measured pre-harvest, thus providing a means of predicting yield.

2.4. Ten-Fold Cross-Validation of the Year 2 Models

In order to investigate how robust these models are to new data, we divided the plots from Year 2 at random into ten groups (folds). We then fitted the models (estimated the parameters) using nine of the folds and then applied the model to the remaining fold and measured the percentage of variance explained. We repeated this procedure for each of the ten folds in turn. Over the 10 folds, the percentage of variance explained by the yield model ranged from 83.28%–97.90%, with an average of 93.61%. The same approach was applied to the “Vegmass” model, where the average percentage of variance explained was 67.57%. As expected, there is a slight drop in the percentage of variation explained compared to the original model that used all of the available data, but the decrease is not substantial, suggesting that the models are robust to new data obtained from plants grown under similar environmental conditions.

2.5. Seed Oil Content Rather Than Protein Drives Seed Yield

A correlation matrix was used to visualise positive and negative relationships between different traits (Figure 4). The matrix showed (factors mentioned in the text below are boxed in blue in Figure 4) that the oil:protein ratio was correlated positively with seed number per pod (0.756) and negatively with seed packing density per pod (−0.702). This appeared to be driven by the lipid component, since there was a negative relationship between protein content and seed weight per pod (−0.627). This suggests that seed yield gain for an individual plant was driven by increased oil content. There was also a strong negative relationship between pigment (chlorophyll and carotenoid) content of the pod and protein/sugars (range from −0.239—−0.579). This implies that senescence may drive the remobilisation of sugars and proteins into the developing seed and that retention of pigment in the pod during development does not necessarily imply that resources will be translocated into seeds at a later stage. Flowering window has a non-significant relationship with yield and biomass. The relationship between the majority of the leaf minerals analysed and yield was also non-significant, the exceptions being canopy leaf potassium, which was strongly negatively correlated with yield (−0.432), total biomass (−0.403) and vegetative biomass (−0.375) and early leaf phosphorus, which was positively correlated with yield (0.281), total biomass (0.384) and vegetative biomass (0.408).
Multi-Dimensional Scaling (MDS) was used to project the distances between traits in a schematic diagram (Figure 5). The correlation between each pair of traits was converted into a distance, where a correlation of ±1 was given a distance of zero and a correlation of zero was given a distance of one. We then used MDS to show the best two-dimensional projection of these distances; traits that are positioned close together are highly correlated. Traits that the models (Table 3) deem useful for predicting yield and/or biomass are circled, and colour is used to group traits according to their type: yield, architecture, leaf chlorophyll, pod chlorophyll, progression of plant development, oil/protein content and minerals. Biomass, the most important predictor of yield, was located relatively close to the yield term, and biomass itself is positioned closely to the architectural traits nearby. The next two predictors from the linear yield model, seed protein and canopy leaf chlorophyll b content, are located further away from the yield term. Leaf mineral analysis makes no contribution to the yield models. Perhaps surprisingly, traits relating to floral development are not important for the yield model, so the date of first flower opening, photoperiod, flowering duration and time to pod shatter had no significant impact on yield.

2.6. Analysis of Varieties to Determine Those Produce High Yields and Highlight Gaps in Breeding Potential

Having established which traits are important for predicting yield (Table 3), next we investigated which varieties displayed the traits contributing to high yield. For this analysis, we used the mean of both blocks for each variety for each trait contributing to the models for yield and above-ground vegetative biomass at harvest (Table 3). Each of the traits was scaled to have a mean = 0 and variance = 1. For the cluster analysis, we took the Euclidean distance between each pair of varieties and combined varieties using Ward’s method, creating three groups from the whole dataset (Figure 6). A combined approach was then taken to include the other traits shown to be important in our earlier models and to cluster them all against the varieties studied (Figure 6b). Varieties are ordered according to the cluster analysis shown in Figure 6a. The top cluster contains all of the highest yielding varieties and the bottom cluster all the lowest yielding varieties. High seed protein, seed linolenic acid content, chlorophyll b content of the canopy leaf and early leaf calcium content were all negatively correlated with yield; however, it should be noted that canopy leaf chlorophyll b was a positive term within the statistical model for log[yield]. This is because the correlation is relatively weak with yield (0.11), biomass (0.03) and protein (0.09), so canopy leaf chlorophyll b only becomes an important term once biomass and protein are taken into account. Conversely, the architectural plant traits of raceme area and plant height were positively correlated with yield, together with lipid content and pod weight for all varieties, except Darmor, Rameses and Canard.
In the model for explaining yield, the top two terms were vegetative biomass (positive) and total protein (negative). Total seed oil also (positively) contributed to yield. Therefore, the oil:protein ratio and biomass traits were used to illustrate the clustering of different varieties in two-dimensional space (Figure 7). Canard produced the highest yield, but had a comparatively low oil:protein ratio. It also had a very poor score for raceme stiffness and has been observed to lodge readily in the field; therefore, although Canard had a high vegetative biomass, the architecture was not sufficiently robust. The varieties from the cluster coloured grey also had high yield, but show a higher oil:protein ratio than other clusters, meaning that the yield increase is achieved by assimilating more oil in the seeds. The other groups had lower vegetative biomass, lower oil:protein ratio and lower yield, with Ningyou having the poorest phenotype of all of the varieties. Yield gaps are apparent (green squares on Figure 7) where breeding and selection, or agronomy, could be used to target plants with high oil content from a smaller vegetative biomass (by changing the harvest index to put more resources into seeds), or large plants with high yields of oil-rich seeds, or large plants with low oil content that could potentially be used as biomass crops.

2.7. Traits Measured Early in Crop Development Are Also Good Predictors of Yield

Our earlier model for explaining yield (Table 3) relied on many traits that can only be observed at, or close to, the time of harvest. While that model is useful for selecting the variety for growing in subsequent seasons, a model that forecasts yield during the growing period would also be more useful. Traditional early observations, such as date of flowering and leaf chlorophyll retention, had very poor correlation with seed yield (Figure 8).
Observations in the field of the highest and lowest yielding varieties showed that there was very little visual difference in colour of pods on the main raceme from May–July (Figure 9a), although the absolute leaf chlorophyll levels in April demonstrated clear differences between the varieties (Figure 8). However, when the canopy was viewed as an entirety, there were clear differences in senescence between the high and low yielding varieties (Figure 9b) with high yielding plants having a slower senescing canopy. When the extreme lines (defined as the top ten lines for seed yield, retention of pod chlorophyll and biomass) were analysed, only two high yielding cultivars, Canard and Victor, also had high biomass and slow pod senescence. In other cases, high biomass lines overlapped with high yielding lines, but retention of pod chlorophyll did not correlate well (Figure 10a).
It was essential to develop a model that can predict final yield prior to visual indicators of yield differences. The forward selection procedure was repeated, but only allowing covariates that were unrelated to flowers or pods in the model. The fitted model (Table 4, Figure 10b) explains 63.69% of the variance in yield, compared to 95.12% in the earlier model (Table 3). This reduction was due to the omission of information regarding the later stages of plant development. However, the model still offers a good method for predicting yield approximately 15 weeks prior to harvest, i.e., approximately April of the same year. The model comprises two of the best architectural variables (stem area and numbers of branches on the primary stem) that were used to explain vegetative biomass in the earlier model. Canopy leaf chlorophyll b content, which was in the earlier model for predicting yield, was also included in the predictive model (Table 4). Interestingly, two leaf mineral concentration traits were also selected, canopy leaf potassium and canopy leaf magnesium.

3. Discussion

Seed yield is a complex trait influenced by numerous interacting variables that are governed by both genetic and environmental factors. Several studies have sought to understand the contribution of individual traits to oilseed rape yield [4,6,32,33,34,35] or to model yield such that it can be predicted over future growing seasons [5,7,8,36]. However, there is some disagreement between the studies as to what are the most important traits that influence yield, be it the number of pods per plant [34], the time to maturity [6] or the duration of flowering [32]. One of the reasons for the discrepancies observed between the different studies is the limitations in the number or type of traits that were measured or the use of a limited number of cultivars, which constrained the genotypic variation that could be observed.

3.1. Greater Branching Density Drives Yield Improvement

Oilseed rape architectural traits accounted for the greatest percentage of variation in above-ground biomass at harvest. Significant positive correlations were observed between both total above-ground biomass and vegetative biomass and seed yield. Thus, although partitioning more resources into vegetative development could be seen to divert resources away from reproductive development, a larger above-ground biomass with increased branching, supported by a larger stem and taller plants, gives rise to more available positions upon which pods can develop. In this study all plots were planted at the same seed rate within year, but it is well established that lower seed rates can be used to produce plants with more branches than those in tightly-packed canopies [37].
Furthermore, an increase in seed yield per plant does not seem to be the result of heavier seeds, as judged by the lack of correlation between total seed yield and Thousand Grain Weight (TGW), a result that was also observed in Arabidopsis [17] and kale [38]. This also supports our hypothesis that increased yield in terms of seed mass per plant is mainly a consequence of more pods rather than a positive increase in resource accumulation within individual seeds or pods. Whilst seed yield per plant is positively correlated with plant height, this is not a trait that breeders would wish to emphasise, as tall plants are generally more susceptible to lodging [39], which can decrease seed yields by up to 16% [40]. For instance, transgenic dwarf oilseed rape lines have been found to have increased yields compared to their taller counterparts [41], and therefore, the breeding target would be to shorten the internode length whilst maintaining or increasing branch number on each plant. This is most likely to be achieved by targeting the plant hormone signalling network, which mediates the genetic regulation of responses to the environment, facilitating individual plants to change their architecture over time in response to the external conditions [42]. Cytokinins promote the emergence of axillary buds into fully-formed branches, whereas strigolactones inhibit branching. Plants are sufficiently sophisticated so that branching can be temporally or spatially regulated within an individual plant; indeed several authors regard plants as a system of competing populations of redundant organs since no individual pod, leaf or flower can be regarded as indispensable [43], and it is known that oilseed rape produces more flowers than are required to obtain maximum yield [35].

3.2. Increasing Seed Oil:Protein Ratio

The model also showed that seed oil and protein were important drivers of yield; indeed, the oil:protein ratio (derived from the seed protein and oil terms used in Model 3b) explains a further 8.1% of the variation within yield, once biomass is taken into account. Grami et al. [44] also established that seed oil content is negatively correlated with protein content; therefore, increasing the oil fraction of the seed would naturally enhance the oil:protein ratio and, according to our model, potentially help increase yields independently of changing architectural traits, although the mechanism by which this conversion occurs is presently not clear. A recent study using Arabidopsis populations established that seed carbon and N content were antagonistic and that % N was negatively correlated with yield [45], as it is in the present study. In Arabidopsis, the number of seeds and pods per plant were closely related to the concentration of seed oil and protein within the seeds [19]. It was shown for soybean that a 1-kg−1 increase in oil content will usually lead to a 2-kg−1 decrease in protein due to a negative genetic correlation between the two yield components [46]; therefore, changing the oil:protein ratio towards oil production would naturally influence seed yield when defined as the total weight of seeds per plant. None of the lines studied have both a high above-ground biomass and a high oil:protein ratio, a combination that would be predicted to produce even greater yields than those recorded. The lines that show the greatest above-ground biomass (and highest seed yields) could also have their vegetative material utilised for the production of bioethanol and biogas, easing competition with food crops for land and resources. Whilst this technology is still being developed, especially in the case of appropriate pre-treatments to extract sugars from the lignocellulosic plant cell walls, yields of 20.4 g methane 100 g−1 DM and 10.9 g ethanol 100 g−1 DM have been reported for oilseed rape [47]. Although these values are much lower than those for sugarcane (37 g ethanol 100 g−1 DM) [48], it is still a useful by-product of oilseed rape production and a technology under development.
In the current study, there was only limited variation between the lines in terms of seed total oil content (45%–58%), but the range of individual fatty acids was highly variable between varieties (Table S1). Many of the oilseed rape varieties grown commercially are High Oleic acid, Low Linolenic acid (HOLL) and the current population contained HOLL and non-HOLL varieties. Seed linolenic acid was consistently and negatively associated with yield; a reflection of the breeding programme that has already taken place to increase yield simultaneously with developing varieties that produce healthy, low trans-fat and stable frying oils and have good disease resistance [49]. It also demonstrates that our model is robust and a true reflection of the genetic basis of yield. Transgenic or mutation (TILLING) technologies would appear to be the most effective method of manipulating the composition of fatty acids within the seed such that it fits the purpose of use, be it for the production of lubricants in industry, or of edible vegetable oils [50], or as a way of producing omega 3-rich oils to enter the human food chain [51]. For example, expressing the yeast gene glycerol-3-phosphate dehydrogenase (gpd1) under the control of a seed-specific promoter increased B. napus seed oil content by 40% [52].

3.3. Early Leaf P Status Correlates with Final Yield

The predictive yield model showed that the traits that explain most of the variation associated with yield are architectural traits; beyond this, the concentration of minerals such as canopy leaf potassium (K), early leaf phosphorus (P) and early leaf N significantly influenced yield. Previous work surrounding the effect of P fertilisation on seed yields was inconclusive, with reports of increasing P having no effect on total seed oil and protein concentration [53,54], decreasing seed linoleic acid concentrations [55] or increasing total seed oil concentrations [56]. Therefore, more research is necessary to clarify the effect of P on yield. One hypothesis is that early leaf P reflects good access to this poorly mobile nutrient in the soil and reflects a strong root system capable of supplying resources during crop development. Given that early leaf P status is a trait that can be measured before the pods have begun to develop, it could provide a useful early indicator of seed yield. Equally, raceme width at the time of flowering was a good indicator of biomass and capable of predicting 36.64% of yield. Similarly, a study in maize found that measurements of raceme width early in the growing season also provided a good predictor of final yield [57]; these early indicators of yield could be used by growers to modify inputs and optimise yield in order to maximise profit margins.

3.4. Number of Pods Determines Yield

The model presented here used the weight of seeds per plant as a measure of yield; however, we also explored whether other parameters of yield, which are used commercially, such as TGW or the number of seeds per pod, could be used as a measure of final seed yield. From the correlation matrix, it became evident that traits previously used as a proxy for seed yield per plant, such as the number of seeds per pod [58], seed weight per pod [59], pod length [60] and even TGW, were not significantly correlated with seed yield per plant in this study (Figure 4). This finding indicates that current commercially-used measures, such as TGW, may not be the most appropriate indicator of yield for breeding and selection purposes. Only the numbers of seeds per plant and pods per plant were significantly correlated with seed yield per plant. This indicates that the pod in B. napus is a fairly conserved organ and that a plant will preferentially invest in the production of more pods per plant as opposed to increasing the number or mass of seeds within an individual pod. Once a seed contains the minimum amount of resources necessary to ensure germination viability, then there is no advantage for a plant in investing more resources beyond this in an individual seed, and a better survival strategy is to make more seeds in different pods, a view also reported in a study on Arabidopsis resource allocation [19]. Similarly, the addition of N was found to increase oilseed rape yields through the production of more pods as opposed to affecting individual seed or pod weight [33], something also observed in Year 1 of the current trial where plants were grown on high and low N.

3.5. Flowering Start Date and Duration Do Not Affect Seed Yield

Oilseed rape has been bred to be self-fertile, but it can also be pollinated by wind and insects; the latter was found to enhance seed oil content quality, potentially due to the ability of insects to optimise the timing of fertilisation for the plant [61]. In Arabidopsis grown under controlled conditions, initially under short days and then switched to long days to induce flowering, the plants that flowered later had a lower yield [17]. In contrast, in the current field-based study, it was found that the flowering start date and flowering window were related neither to the number of pods per plant, nor to seed yield. Thus, increasing the flowering duration does not appear to increase pollination efficiency, perhaps because the general decline in pollinators within agricultural systems [62] means insufficient exploitation of the additional flowers. The increased early flowering nature of lines with long flowering duration may also have increased the unwanted presence of pollen beetle pests (Meligethes spp.), which are attracted by the visual and olfactory properties of oilseed rape flowers [63]. Thus, within the current field trial, the unopened buds on plants of those lines that flowered earlier may have been at increased risk from feeding damage than those on plants not yet flowering, especially if pollen resources were scarce due to high intraspecific competition, and this may partly explain why it did not afford a yield advantage to be flowering earlier or over a longer period.

3.6. The Relationship between Plant Senescence and Seed Yield

In this study, neither chlorophyll concentration nor senescence of the pods were related to seed yield per plant or other yield parameters. This is in contrast to previous work on Arabidopsis [64] and wheat [65] where delaying leaf senescence led to an increase in seed yield. The majority of photosynthetic tissue in an oilseed rape plant is raceme material; therefore, a greater biomass would also provide a larger photosynthetic surface area through an increased number of stems. In the present study, canopy leaf chlorophyll b content (a proxy measure for the stage of senescence) had a positive relationship with final yield, suggesting that the leaf resources are mobilised into the developing plant prior to abscission via a route that contributes to seed yield; hence, the parameter was included in the predictive model. Chlorophyll b content is a measure of light-harvesting complex status rather than of reaction centre complexes, and light harvesting complexes (LHC) are second only to Rubisco as sources of recyclable protein [66]. LHC also tend to be in excess, so the chloroplast might be able to lose significant quantities of chlorophyll b and associated proteins without it having a negative effect on photosynthesis. In cereals, photosynthesis of the canopy is largely responsible for the carbon yield of seeds, whereas seed N is mobilized from senescing vegetative tissues [67]. Stems have also been shown to play a pivotal role in B. napus seed filling by contributing 31% of the carbon contained within seeds [68]. Experiments in which soy bean pods were removed from mature plants showed that N fixation continued to take place in the roots [69]; elevated levels of N represent a ‘metabolic sink’ for photosynthate, thus altering carbohydrate metabolism and sugar signalling in source leaves. Leaves therefore accumulated starch products that would normally be mobilised to the pods. De-podding Arabidopsis showed that restricting the number of sinks increased the concentration of lipids in the remaining pods [19]; since B. napus loses leaves relatively early in its reproductive cycle, it seems reasonable to conclude that canopy photosynthesis is unlikely to contribute to oil/protein content in the seed, but it may make a contribution to carbon accumulation. However, as discussed in Bennett et al. [18], Arabidopsis is a weed that treats each pod as an individual unit, and local remobilisation of photosynthates into the seeds makes sense as it is important that each pod shatters as soon as the seed attains viability in order to stand the best chance of reproductive success. In contrast, the artificial selection that occurs within crop species has sought to synchronise seed maturation across the plant; therefore, the photosynthetic window of a pod is less closely related to yield and shatter since the pod sits in ‘suspended animation’ waiting for the whole plant to be ready to shatter. The extent of developmental coordination is poorer in oilseed rape compared to other crop species that have been domesticated for longer, such as cereal crops [70,71], but artificial selection has started to influence the loss of weedy reproductive traits. Other authors have found that seed maturity, as opposed to pod wall maturity, drives pod shatter in oilseed rape [72], highlighting the uncoupling of pod wall senescence in modern varieties from shattering and seed yield.
An increase in the number of pods per plant does not affect the rate of pod senescence or pod development in oilseed rape, giving credence to the idea that pod reproductive development is not necessarily driving the senescence process, but it is temporally fixed under growing conditions that impart only minimal stress on the plant. Just as leaf senescence has been shown to naturally occur in an age-dependent manner [73,74], so whole plant senescence also appears to be controlled by the age of the plant, an assumption borne out by the fact that there was little variation in the lifespan of different oilseed rape lines despite their genetic diversity. Plant longevity appears to be connected to plant mass, with small plants having a shorter lifespan than larger ones; however, within-species lifespan is very similar to keep the population stable [75]. An as yet unknown biological system in B. napus controls the timing of whole plant senescence and accordingly lifespan. Thus, whilst senescence is the visible output of an internal biological clock, the number of pods, plant size and seed yield are the consequences of resource availability.

4. Materials and Methods

In an effort to understand better the physiological basis governing seed yield in oilseed rape, 35 cultivar lines of B. napus, showing substantive divergence in seed yield per plant were assessed over two seasons for 133 different traits, covering developmental, architectural and nutritional parameters, generating a powerful dataset for exploring the factors underpinning yield.

4.1. Site Description and Meteorological Conditions

The experiment was carried out over two growing seasons, Year 1: August 2009–July 2010; Year 2: September 2010–August 2011, at Rothamsted Research, Hertfordshire, UK (51°48′32.6376″ N, 0°21′22.5432″ W). The trials in each year were conducted in different fields, but in both cases, the previous crop was spring barley. The soil type for the Year 1 field was flinty silt clay loam (Batcombe–Carstens), and in Year 2, it was clay with flints (Batcombe). Soil samples (30-cm cores) were taken prior to the application of fertilisers. The soil in Year 1 contained a total residual N concentration of 18 kg N ha−1, which was comprised of 9 kg ha−1 N-NO3 and 9 kg ha−1 N-NH4+, whereas that in Year 2 was 22 Kg N ha−1, comprising 9 kg ha−1 N-NO3 and 13 kg ha−1 N-NH4+.
Weather parameters (rainfall, maximum and minimum temperature) for the experimental period (2009–2011) and the 30-year average were taken from the Rothamsted Research meteorological station (Figure S1). There was a large variation in the rainfall pattern between the two years; e.g., June rainfall in Year 2 was approximately four-fold higher than the same period in Year 1.

4.2. Crop Management and Field Trial Design

The field trials were arranged in a randomised complete block design and used two replicate plots per treatment. In Year 1 (2009–2010), two soil N supplies were investigated, low (= residual) N (~18 kg ha−1) and high N (148 kg ha−1). The high N treatment was applied in a split dose application; 30 kg ha−1 N were applied in the autumn (29 September 2009), and 50 kg ha−1 N were subsequently applied in the spring (12 March 2010 and 07 April 2010) to a maximum concentration of 148 kg ha−1 total available N. The Year 2 trial used only residual nitrogen (22 kg ha−1) across the whole plot. In Year 1, seeds were drilled on 28 August 2009, and each replicate plot (1 × 1.5 m) was drilled in four rows with 120 seeds per plot. In Year 2, 125 plots were drilled on 6 September 2010, but due to bad weather conditions, the drilling of the remaining 307 plots had to be postponed by one week; ANOVA showed that no significant differences were found in yield parameters as a result of this split drilling time Year 2 plots were 1.8 × 1.5 m, drilled with 160 seeds per plot in 14 rows. Pre-germination herbicides and fungicides were applied to both trials at the manufacturers’ recommended rates with pests, diseases and volunteer cereals/weeds controlled according to standard agronomic practice. In Year 1, all of the 61 winter cultivars planted were sampled, and this included 48 winter oilseed rape (WOSR) lines, 8 fodder/forage/salad kales, 4 swedes and 1 synthetic line. In Year 2, of the 84 WOSR varieties planted, a subset of 35 lines was sampled; being selected if they were also present in Year 1 of the trial and had transcriptome sequencing information available.

4.3. Chlorophyll Measurements

Chlorophyll and carotenoid concentration were analysed in early and canopy leaf samples as a proxy measure of senescence, and in Year 2, additional sampling was made from leaf disks from upper and lower leaves and whole pods sampled at set intervals throughout the growing season, March–June and May–July, respectively. The early leaves represented the earliest not fully-expanded leaves, and the canopy leaves were the youngest true leaves on the plant when sampled on 2 February 2011. The lower leaf (sampled 11 March, 6 April and 3 May 2011) was the 5th leaf up from the base of the main raceme, and the upper leaf (sampled 3 May 2011 and 6 June 2011) was that immediately subtending the lowermost pod. All pods were sampled from the main raceme. Chlorophyll was extracted and analysed in the leaf and pod material using the methods described in Bennett et al. [29]. The progression of chlorophyll and carotenoid loss from the leaves and pod walls was calculated by determining the % difference from the previous sampling date. The pod photosynthetic period was calculated as the time from the start of flowering until pod wall chlorophyll fell below 10% of its maximum.

4.4. Leaf Mineral Analysis

Early and canopy leaves from ten plants per plot were dried overnight at 80 °C and stored in a controlled environment room at 15 °C, 15% humidity. Before milling, the leaf samples were re-dried overnight at 60 °C. Plots were analysed in duplicate for their mineral content, and a Kjeldahl digestion was used to analyse total calcium (Ca), potassium (K), magnesium (Mg), manganese (Mn), nitrogen (N) and phosphorus (P) using the methods described in Broadley et al. [76]. Aqueous extracts were subsequently used to measure the nitrate concentration [76].

4.5. Plant Development

The plots were walked periodically between March–July. Growth stage [77] was visually scored over the whole plot and graded accordingly based on 25% of the plants in a plot reaching a particular stage. Flowering duration was defined by the difference between the start date of flowering (when 25% of the plot had reached GS60) and the date when >95% of the plants within the plot had no more flowers.
Pod development was also visually scored across the whole plot and categorized per plot on the following scale: 0: plot still in flower; 0.5: immature pods; 1: mature green pods fully expanded; 1.5: 25% yellow pods; 2: 50% yellow pods (yellow/green); 2.5: 75% yellow pods; 3: 100% yellow pods; 3.5: 50% brown pods; 4: 100% brown pods just before shatter. A plot was recorded as shattering when 25% of the pods in the whole plot were seen to have dehisced.

4.6. Plant Architecture

Several plant architectural traits were measured at the time of harvest. Five plants per plot were selected at random from each plot and were returned to the laboratory for assessment of the following traits: plant height, length of the raceme containing pods, the total number of branches and number of branches arising from the main raceme only. Pod spacing on the main raceme was calculated by dividing the length of the raceme containing pods by the number of pod sites on the raceme, i.e., including sites where pods had subsequently aborted. Raceme area was calculated using a calliper to measure two raceme diameters 10 cm above the ground at 90° to each other and multiplying these values together. Raceme stiffness was visually scored on a scale of 1–9 where: 1, all plants prostrate; 2, main stems leaning at 20°–30° to the horizontal; 3, main stems leaning at 45° to the horizontal; 4, main stems leaning at 70°–80° to the horizontal; 5, main stems leaning slightly, but branched canopy bent over at 10° to the horizontal; 6, main stems leaning slightly, but branched canopy leaning at 20°–30° to the horizontal; 7, main stems leaning slightly, but branched canopy leaning at 45° to the horizontal; 8, main raceme erect with branched canopy leaning at 70°–80° to the horizontal; and 9, all plants completely erect. The above-ground biomass was calculated from the dry weight of 5 plants per plot at harvest; the vegetative mass was calculated from the above-ground biomass minus the seed yield per plant.

4.7. Pod Physiology

Several measures of pod physiology were taken at harvest from pods on the main raceme at Stage 4 (brown and ready to dehisce) just before the plots were harvested. In each case, ten replicate pod samples were taken per plot, and there were four plots per line. Pod length was determined by measuring from the start of the bivalve to the end of the beak, and the seed cavity length was also determined (the length of the pod minus the beak). The seeds were removed from the pod and aborted seeds discarded before the remainder was counted. The seed packing density was determined by dividing the pod cavity length by the number of seeds per pod. The weight of seeds per pod was determined by placing the seeds in a controlled environment room at 10 °C, 15% humidity until they had reached a constant weight. Seed area was calculated using a MARVIN seed scanner (Gta Sensors Gmbh-Germany).

4.8. Yield

At harvest, ten plants per plot were bagged and dried at 20 °C. The plants were then threshed and winnowed before the seeds were placed in a controlled environment room at 15 °C, 15% humidity until they had reached a constant weight. Seed yield was subsequently determined as total seed weight per plant. To calculate the MARVIN Thousand Grain Weight (TGW), 100 harvested seeds per plot (post threshing, winnowing and storage) were placed into the MARVIN seed scanner with a ±10 seed tolerance; any seed that was not representative of the batch (i.e., damaged or discoloured) was removed by hand prior to the analysis, and the seeds were subsequently weighed. The following calculations were used to determine the MARVIN TGW, TGW, number of seeds per plant, number of pods per plant and harvest index.
MARVIN   TGW :   [ Weight   of   seeds   ( g ) Number   of   seeds ]   × 1000   ( g )
TGW :   [ Seed   weight   per   pod   ( g ) Number   of   seeds   per   pod ]   × 1000   ( g )
Number   of   seeds   per   plant :   Seed   yield   per   plant   ( g )   ÷   [ TGW   ( g ) 1000 ]
Number   of   pods   per   plant :   [ Seed   yield   per   plant   ( g ) Seed   weight   per   pod   ( g ) ]
Harvest   index :   [ Seed   yield   per   plant   ( g ) Total   above   ground   biomass   ( g ) ]

4.9. Seed Composition

Seed oil, fatty acid, protein, moisture and glucosinolate content were determined by Near Infrared Spectroscopy (NIRS) using the methods described in Kelly et al. [78]. In the current study, seed samples from five plants per plot were analysed.

4.10. Statistics

Statistical analyses were performed using the statistical software package R (R Core Team, Florham Park, New Jersey, USA, 2013), except for the yield mean per plant and standard error of the mean, which were calculated using IBM SPSS version 19 (IBM Corp., Armonk, NY, USA). For the Pearson’s correlation analysis of all traits, with 9180 pairwise tests, a Bonferroni multiple comparison correction was applied making the cut-off value for significant correlations ±0.496. To determine the main traits influencing seed yield per plant and above-ground biomass, linear regression models were developed using forward selection based on the 133 traits listed in Supplementary Table S1. Starting with the null model (no traits included), forward selection identifies the most significant trait to include in the model based on an F-test. The process was repeated multiple times, at each stage selecting the most significant model with the n + 1 variable compared to the current model with n variables based on the F-test. The process was stopped when none of the remaining traits significantly improved the estimation of the response variable at the 5% significance level. Some traits, including yield, underwent a log transformation prior to inclusion in the model to improve the linear relationship between the response variable and the other traits. After selecting the traits identified in Table 3 and standardising, the traits and lines underwent hierarchical clustering using Ward linkage, to compile the dendrogram, heat map and scatter diagram (Figure 4 and Figure 6).

5. Conclusions

We have developed a forward stepwise multiple linear regression model to determine the traits most important in influencing seed yield per plant in B. napus. From the 133 traits recorded, that could have potentially influenced seed yield, above-ground biomass and protein were found to be the most important factors, together accounting for 94.36% of the variation in seed yield. Increasing our understanding of resource allocation strategies could lead to improved plant breeding via marker-assisted selection, and this would have important consequences for crop breeding targets and meeting future food security demands. It would appear that year-on-year environmental weather conditions have a bigger impact than soil conditions, since the Year 2 model fitted to Year 1 data was less good than the same model fitted to different soil nitrogen regimes in Year 1. However, the refined model developed in Year 2 included more informative parameters that provided more robust indicators of yield. We therefore propose that this model would require adjustment for climate and farm-specific soil mineral content, but that it could provide the basis of a model for oilseed rape yield across the country.
Overall, our findings demonstrate that there is a diverse array of resource allocation strategies within the B. napus gene pool, and different approaches will need to be applied depending on which yield components are to be targeted. Whilst biomass and protein are the main contributors towards the variation observed in seed yield, the traits that underpin these variables, such as raceme width for biomass and early leaf mineral concentration, could provide useful early indicators of seed yield, as well as potential indirect targets for improving yields.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4395/7/2/31/s1: Figure S1: Monthly figures for (a) total rainfall and (b) mean maximum and minimum temperature at Rothamsted Research; Table S1: Analysis of variance table for (a) variety and nitrogen terms on Year 1 data and (b) variety and year terms on data from both years; Table S2: Mean, range (min and max), standard deviation (SD) and coefficient of variation (CV) of the physiological and biochemical traits measured in 35 B. napus lines.

Acknowledgments

This work was supported by funding from AL Tozer Ltd., Pyports, Downside Bridge Road, Cobham, Surrey, KT11 3EH, U.K., and The University of Reading Endowment Trust Fund to enable Emma Bennett to undertake her PhD. The OREGIN Diversity Trials at Rothamsted Research were supported by the U.K. Defra Project IF0144. Rothamsted Research is a national institute of bioscience strategically funded by the U.K. Biotechnology and Biological Sciences Research Council (BBSRC). The authors are grateful to Helen Ougham (IBERS, Aberystwyth) for critical reading of the manuscript. The collection of data described in this paper would not have been possible without the assistance of several dedicated undergraduate students from the University of Reading: Jason Fitzgerald, Manon Morin and Leo Shum, supported by Ian Bennett; and from Rothamsted Research, extensive assistance was provided by Scott Pilkington-Bennett, Stephanie Cooper, Carl Halford, Matthew Skellern, Nigel Watts and the farm staff who grew and maintained the trials. The OREGIN project was managed by Jacqueline Barker, and Sue Welham provided statistical rigour in the experimental design of the field trials.

Author Contributions

E.J.B. performed the majority of trait screening and preliminary analyses and drafted the text. C.J.B. performed the statistical modelling. P.W.C.C. curated all OREGIN data. S.M.C., P.E., J.H., C.L. contributed trait data and to the writing process. G.T. and G.K. obtained funding for the OREGIN trials and managed the project. J.A.R. and C.W. supervised E.J.B., conceived of the experiments, obtained funding for the work and wrote the finished article.

Conflicts of Interest

The authors declare no conflict of interest. The founding sponsors had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript, nor in the decision to publish the results.

References

  1. Wang, E.; Engel, T. Simulation of phenological development of wheat crops. Agric. Syst. 1998, 58, 1–24. [Google Scholar] [CrossRef]
  2. Yang, H.S.; Dobermann, A.; Lindquist, J.L.; Walters, D.T.; Arkebauer, T.J.; Cassman, K.G. Hybrid-maize—A maize simulation model that combines two crop modeling approaches. Field Crops Res. 2004, 87, 131–154. [Google Scholar] [CrossRef]
  3. Son, N.T.; Chen, C.F.; Chen, C.R.; Minh, V.Q.; Trung, N.H. A comparative analysis of multitemporal MODIS EVI and NDVI data for large-scale rice yield estimation. Agric. For. Meteorol. 2014, 197, 52–64. [Google Scholar] [CrossRef]
  4. Habekotte, B. Evaluation of seed yield determining factors of winter oilseed rape (Brassica napus L.) by means of crop growth modelling. Field Crops Res. 1997, 54, 137–151. [Google Scholar] [CrossRef]
  5. Habekotte, B. A model of the phenological development of winter oilseed rape (Brassica napus L.). Field Crops Res. 1997, 54, 127–136. [Google Scholar] [CrossRef]
  6. Habekotte, B. Options for increasing seed yield of winter oilseed rape (Brassica napus L.): A simulation study. Field Crops Res. 1997, 54, 109–126. [Google Scholar] [CrossRef]
  7. Deligios, P.A.; Farci, R.; Sulas, L.; Hoogenboom, G.; Ledda, L. Predicting growth and yield of winter rapeseed in a Mediterranean environment: Model adaptation at a field scale. Field Crops Res. 2013, 144, 100–112. [Google Scholar] [CrossRef]
  8. Husson, F.; Wallach, D.; Vandeputte, B. Evaluation of CECOL, a model of winter rape (Brassica napus L.). Eur. J. Agron. 1998, 8, 205–214. [Google Scholar] [CrossRef]
  9. Peltonen-Sainio, P.; Jauhiainen, L.; Hyovela, M.; Nissila, E. Trade-off between oil and protein in rapeseed at high latitudes: Means to consolidate protein crop status? Field Crops Res. 2011, 121, 248–255. [Google Scholar] [CrossRef]
  10. Nesi, N.; Delourme, R.; Bregeon, M.; Falentin, C.; Renard, M. Genetic and molecular approaches to improve nutritional value of Brassica napus L. seed. Comptes Rendus Biol. 2008, 331, 763–771. [Google Scholar] [CrossRef] [PubMed]
  11. Smooker, A.M.; Wells, R.; Morgan, C.; Beaudoin, F.; Cho, K.; Fraser, F.; et al. The identification and mapping of candidate genes and QTL involved in the fatty acid desaturation pathway in Brassica napus. Theor. Appl. Genet. 2011, 122, 1075–1090. [Google Scholar] [CrossRef] [PubMed]
  12. Chalhoub, B.; Denoeud, F.; Liu, S.; Parkin, I.A.P.; Tang, H.; Wang, X.; Chiquet, J.; Belcram, H.; Tong, C.; Samans, B.; et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 2014, 345, 950–953. [Google Scholar] [CrossRef] [PubMed]
  13. Gressel, J. Genetic Glass Ceilings: Transgenics for Crop Biodiversity, 1st ed.; Johns Hopkins University Press: Baltimore, MD, USA, 2008; p. 488. [Google Scholar]
  14. FAOSTAT. Available online: http://faostat3.fao.org/faostat-gateway/go/to/download/Q/QC/E (accessed on 31 October 2015).
  15. Mawson, R.; Heaney, R.K.; Zdunczyk, Z.; Kozlowska, H. Rapeseed meal-glucosinolates and their antinutritional effects. Part II. Flavor and palatability. Mol. Nutr. Food Res. 1993, 37, 336–344. [Google Scholar]
  16. Quijada, P.A.; Udall, J.A.; Lambert, B.; Osborn, T.C. Quantitative trait analysis of seed yield and other complex traits in hybrid spring rapeseed (Brassica napus L.): 1. Identification of genomic regions from winter germplasm. Theor. Appl. Genet. 2006, 113, 549–561. [Google Scholar] [CrossRef] [PubMed]
  17. Masclaux-Daubresse, C.; Chardon, F. Exploring nitrogen remobilization for seed filling using natural variation in Arabidopsis thaliana. J. Exp. Bot. 2011, 62, 2131–2142. [Google Scholar] [CrossRef] [PubMed]
  18. Bennett, E.J.; Roberts, J.A.; Wagstaff, C. The role of the pod in seed development: Strategies for manipulating yield. New Phytol. 2011, 190, 838–853. [Google Scholar] [CrossRef] [PubMed]
  19. Bennett, E.J.; Roberts, J.A.; Wagstaff, C. Manipulating resource allocation in plants. J. Exp. Bot. 2012, 63, 3391–3400. [Google Scholar] [CrossRef] [PubMed]
  20. Byzova, M.; Verduyn, C.; Brouwer, D.; Block, M. Transforming petals into sepaloid organs in Arabidopsis and oilseed rape: Implementation of the hairpin RNA-mediated gene silencing technology in an organ-specific manner. Planta 2004, 218, 379–387. [Google Scholar] [CrossRef] [PubMed]
  21. Buchanan-Wollaston, V.; Earl, S.; Harrison, E.; Mathas, E.; Navabpour, S.; Page, T.; Pink, D. The molecular analysis of leaf senescence—A genomics approach. Plant Biotechnol. J. 2003, 1, 3–22. [Google Scholar] [CrossRef] [PubMed]
  22. Balazadeh, S.; Parlitz, S.; Mueller-Roeber, B.; Meyer, R.C. Natural developmental variations in leaf and plant senescence in Arabidopsis thaliana. Plant Biol. 2008, 10, 136–147. [Google Scholar] [CrossRef] [PubMed]
  23. Levey, S.; Wingler, A. Natural variation in the regulation of leaf senescence and relation to other traits in Arabidopsis. Plant Cell Environ. 2005, 28, 223–231. [Google Scholar] [CrossRef]
  24. Flood, P.J.; Harbinson, J.; Aarts, M.G.M. Natural genetic variation in plant photosynthesis. Trends Plant Sci. 2011, 16, 327–335. [Google Scholar] [CrossRef] [PubMed]
  25. Wagstaff, C.; Yang, T.J.W.; Stead, A.D.; Buchanan-Wollaston, V.; Roberts, J.A. A molecular and structural characterization of senescing Arabidopsis siliques and comparison of transcriptional profiles with senescing petals and leaves. Plant J. 2009, 57, 690–705. [Google Scholar] [CrossRef] [PubMed]
  26. Lewis, G.J.; Thurling, N. Growth, development, and yield of three oilseed Brassica species in a water-limited environment. Aust. J. Exp. Agric. 1994, 34, 93–103. [Google Scholar] [CrossRef]
  27. Robinson, C.K.; Hill, S.A. Altered resource allocation during seed development in Arabidopsis caused by the abi3 mutation. Plant Cell Environ. 1999, 22, 117–123. [Google Scholar] [CrossRef]
  28. Diaz, C.; Lemaitre, T.; Christ, A.; Azzopardi, M.; Kato, Y.; Sato, F.; Morot-Gaudry, J.-F.; Le Dily, F.; Masclaux-Daubresse, C. Nitrogen recycling and remobilization are differentially controlled by leaf senescence and development stage in Arabidopsis under low nitrogen nutrition. Plant Physiol. 2008, 147, 1437–1449. [Google Scholar] [CrossRef] [PubMed]
  29. Bennett, E.J.; Roberts, J.A.; Wagstaff, C. Use of mutants to dissect the role of ethylene signalling in organ senescence and the regulation of yield in Arabidopsis thaliana. J. Plant Growth Regul. 2014, 33, 56–65. [Google Scholar] [CrossRef]
  30. Hopkins, C.J.; Welham, S.J.; Teakle, G.R.; Peplow, K.-S.; Pink, D.; Carion, P.W.C.; King, G.J.; Barker, J. Oilseed Rape Genetic Improvement Network Diversity Demonstration Trial—Year 1; Rothamsted Research: Harpenden, UK, 2009–2010. [Google Scholar]
  31. Hopkins, C.J.; Welham, S.J.; Teakle, G.R.; Peplow, K.-S.; Pink, D.; Carion, P.W.C.; King, G.J.; Barker, J. Oilseed Rape Genetic Improvement Network Diversity Demonstration Trial—Year 2; Rothamsted Research: Harpenden, UK, 2010–2011. [Google Scholar]
  32. Ali, N.; Javidfar, F.; Elmira, J.Y.; Mirza, A.M.Y. Relationship among yield components and selection criteria for yield improvement in winter rapeseed (Brassica napus L.). Pak. J. Bot. 2003, 35, 167–174. [Google Scholar]
  33. Allen, E.J.; Morgan, D.G. A quantitative analysis of the effects of nitrogen on the growth, development and yield of oilseed rape. J. Agric. Sci. 1972, 78, 315–324. [Google Scholar] [CrossRef]
  34. Lu, G.-Y.; Zhang, F.; Zheng, P.-Y.; Cheng, Y.; Liu, F.-I.; Fu, G.-P.; Zhang, X.-K. Relationship among yield components and selection criteria for yield improvement in early rapeseed (Brassica napus L.). Agric. Sci. China 2011, 10, 997–1003. [Google Scholar] [CrossRef]
  35. Berry, P.M.; Spink, J.H. A physiological analysis of oilseed rape yields: Past and future. J. Agric. Sci. 2006, 144, 381–392. [Google Scholar] [CrossRef]
  36. Gabrielle, B.; Denoroy, P.; Gosse, G.; Justes, E.; Andersen, M.N. A model of leaf area development and senescence for winter oilseed rape. Field Crops Res. 1998, 57, 209–222. [Google Scholar] [CrossRef]
  37. Diepenbrock, W. Yield analysis of winter oilseed rape (Brassica napus L.): A review. Field Crops Res. 2000, 67, 35–49. [Google Scholar] [CrossRef]
  38. Bennett, E.J.; Gawthrop, F.; Yao, C.; Boniface, C.; Ishihara, H.; Roberts, J.A.; Wagstaff, C. Effects of planting density and nitrogen application on seed yield and other morphological traits of the leafy vegetable kale (Brassica oleracea). Ann. Appl. Biol. 2013, 119, 201–216. [Google Scholar]
  39. Navabi, A.; Iqbal, M.; Strenzke, K.; Spaner, D. The relationship between lodging and plant height in a diverse wheat population. Can. J. Plant Sci. 2006, 86, 723–726. [Google Scholar] [CrossRef]
  40. Islam, N.; Evans, E.J. Influence of lodging and nitrogen rate on the yield and yield attributes of oilseed rape (Brassica napus L.). Theor. Appl. Genet. 1994, 88, 530–534. [Google Scholar] [CrossRef] [PubMed]
  41. Al-Ahmad, H.; Dwyer, J.; Moloney, M.; Gressel, J. Mitigation of establishment of Brassica napus transgenes in volunteers using a tandem construct containing a selectively unfit gene. Plant Biotechnol. J. 2006, 4, 7–21. [Google Scholar] [CrossRef] [PubMed]
  42. Leyser, O. Auxin, self-organisation, and the colonial nature of plants. Curr. Biol. 2011, 21, R331–R337. [Google Scholar] [CrossRef] [PubMed]
  43. Sachs, T.; Novoplansky, A.; Cohen, D. Plants as competing populations of redundant organs. Plant Cell Environ. 1993, 16, 765–770. [Google Scholar] [CrossRef]
  44. Grami, B.; Stefansson, B.R.; Baker, R.J. Genetics of protein and oil content in summer rape: Heritability, number of effective factors, and correlations. Can. J. Plant Sci. 1977, 57, 937–943. [Google Scholar] [CrossRef]
  45. Chardon, F.; Jasinski, S.; Durandet, M.; Lecureuil, A.; Soulay, F.; Bedu, M.; Guerche, P.; Masclaux-Daubresse, C. QTL meta-analysis in Arabidopsis reveals an interaction between leaf senescence and resource allocation to seeds. J. Exp. Bot. 2014, 65, 3949–3962. [Google Scholar] [CrossRef] [PubMed]
  46. Specht, J.E.; Hume, D.J.; Kumudini, S.V. Soybean yield potential—A genetic and physiological perspective. Crop Sci. 1999, 39, 1560–1570. [Google Scholar] [CrossRef]
  47. Petersson, A.; Thomsen, M.H.; Hauggaard-Nielsen, H.; Thomsen, A.-B. Potential bioethanol and biogas production using lignocellulosic biomass from winter rye, oilseed rape and faba bean. Biomass Bioenergy 2007, 31, 812–819. [Google Scholar] [CrossRef]
  48. Pandey, A.; Soccol, C.R.; Nigam, P.; Soccol, V.T. Biotechnological potential of agro-industrial residues. I: Sugarcane bagasse. Bioresour. Technol. 2000, 74, 69–80. [Google Scholar] [CrossRef]
  49. Orson, J.; Booth, E.; Merritt, C.; Lea, C. Growing High Oleic Low Linolenic (HOLL) Oilseed Rape for Specialised Markets; HGCA Project Report 442; AHDB Cereals & Oilseeds: Kenilworth, UK, 2008. [Google Scholar]
  50. Venegas-Caleron, M.; Sayanova, O.; Napier, J.A. An alternative to fish oils: Metabolic engineering of oil-seed crops to produce omega-3 long chain polyunsaturated fatty acids. Prog. Lipid Res. 2010, 49, 108–119. [Google Scholar] [CrossRef] [PubMed]
  51. Betancor, M.B.; Sprague, M.; Usher, S.; Sayanova, O.; Campbell, P.J.; Napier, J.A.; Tocher, D.R. A nutritionally-enhanced oil from transgenic Camelina sativa effectively replaces fish oil as a source of eicosapentaenoic acid for fish. Sci. Rep. 2015, 5, 8104. [Google Scholar] [CrossRef] [PubMed]
  52. Vigeolas, H.; Waldeck, P.; Zank, T.; Geigenberger, P. Increasing seed oil content in oil-seed rape (Brassica napus L.) by over-expression of a yeast glycerol-3-phosphate dehydrogenase under the control of a seed-specific promoter. Plant Biotechnol. J. 2007, 5, 431–441. [Google Scholar] [CrossRef] [PubMed]
  53. Brennan, R.F.; Bolland, M.D.A. Effect of fertiliser phosphorus and nitrogen on the concentrations of oil and protein in grain and the grain yield of canola (Brassica napus L.) grown in south-western Australia. Aust. J. Exp. Agric. 2007, 47, 984–991. [Google Scholar] [CrossRef]
  54. Brennan, R.F.; Bolland, M.D.A. Comparing the nitrogen and phosphorus requirements of canola and wheat for grain yield and quality. Crop Pasture Sci. 2009, 60, 566–577. [Google Scholar] [CrossRef]
  55. Krueger, K.; Goggi, A.S.; Mallarino, A.P.; Mullen, R.E. Phosphorus and potassium fertilization effects on soybean seed quality and composition. Crop Sci. 2013, 53, 602–610. [Google Scholar] [CrossRef]
  56. Lickfett, T.; Matthäus, B.; Velasco, L.; Möllers, C. Seed yield, oil and phytate concentration in the seeds of two oilseed rape cultivars as affected by different phosphorus supply. Eur. J. Agron. 1999, 11, 293–299. [Google Scholar] [CrossRef]
  57. Mourtzinis, S.; Arriaga, F.J.; Balkcom, K.S.; Ortiz, B.V. Corn grain and stover yield prediction at R1 growth stage. Agron. J. 2013, 105, 1045–1050. [Google Scholar] [CrossRef]
  58. Zhang, L.; Liu, P.; Hong, D.; Huang, A.; Li, S.; He, Q.; Yang, G. Inheritance of seeds per silique in Brassica napus L. using joint segregation analysis. Field Crops Res. 2010, 116, 58–67. [Google Scholar] [CrossRef]
  59. Yang, P.; Shu, C.; Chen, L.; Xu, J.; Wu, J.; Liu, K. Identification of a major QTL for silique length and seed weight in oilseed rape (Brassica napus L.). Theor. Appl. Genet. 2012, 125, 285–296. [Google Scholar] [CrossRef] [PubMed]
  60. Chay, P.; Thurling, N. Identification of genes controlling pod length in spring rapeseed, Brassica napus L.; and their utilization for yield improvement. Plant Breed. 1989, 103, 54–62. [Google Scholar] [CrossRef]
  61. Bommarco, R.; Marini, L.; Vaissière, B. Insect pollination enhances seed yield, quality, and market value in oilseed rape. Oecologia 2012, 169, 1025–1032. [Google Scholar] [CrossRef] [PubMed]
  62. Biesmeijer, J.C.; Roberts, S.P.M.; Reemer, M.; Ohlemüller, R.; Edwards, M.; Peeters, T.; Schaffers, A.P.; Potts, S.G.; Kleukers, R.; Thomas, C.D.; et al. Parallel declines in pollinators and insect-pollinated plants in Britain and the Netherlands. Science 2006, 313, 351–354. [Google Scholar] [CrossRef] [PubMed]
  63. Cook, S.M.; Rasmussen, H.B.; Birkett, M.A.; Murray, D.A.; Pye, B.J.; Watts, N.P. Behavioural and chemical ecology underlying the success of turnip rape (Brassica rapa) trap crops in protecting oilseed rape (Brassica napus) from the pollen beetle (Meligethes aeneus). Arthropod-Plant Interact. 2007, 1, 57–67. [Google Scholar] [CrossRef]
  64. Gan, S.S.; Amasino, R.M. Inhibition of leaf senescence by autoregulated production of cytokinin. Science 1995, 270, 1986–1988. [Google Scholar] [CrossRef] [PubMed]
  65. Verma, V.; Foulkes, M.J.; Worland, A.J.; Sylvester-Bradley, R.; Caligari, P.D.S.; Snape, J.W. Mapping quantitative trait loci for flag leaf senescence as a yield determinant in winter wheat under optimal and drought-stressed environments. Euphytica 2004, 135, 255–263. [Google Scholar] [CrossRef]
  66. Makino, A.; Osmond, B. Effects of nitrogen nutrition on nitrogen partitioning between chloroplasts and mitochondria in pea and wheat. Plant Physiol. 1991, 96, 355–362. [Google Scholar] [CrossRef] [PubMed]
  67. Gifford, R.M.; Evans, L.T. Photosynthesis, carbon partitioning, and yield. Ann. Rev. Plant Physiol. 1981, 32, 485–509. [Google Scholar] [CrossRef]
  68. Brar, G.; Thies, W. Contribution of leaves, stem, siliques and seeds to dry matter accumulation in ripening seeds of rapeseed, Brassica napus L. Z. Pflanzenphysiol. 1977, 82, 1–13. [Google Scholar] [CrossRef]
  69. Nakano, H.; Muramatsu, S.; Makino, A.; Mae, T. Relationship between the suppression of photosynthesis and starch accumulation in the pod-removed bean. Aust. J. Plant Physiol. 2000, 27, 167–173. [Google Scholar]
  70. Hay, R.; Kirby, E. Convergence and synchrony-a review of the coordination of development in wheat. Aust. J. Agric. Res. 1991, 42, 661–700. [Google Scholar] [CrossRef]
  71. Price, J.S.; Hobson, R.N.; Neale, M.A.; Bruce, D.M. Seed losses in commercial harvesting of oilseed rape. J. Agric. Eng. Res. 1996, 65, 183–191. [Google Scholar] [CrossRef]
  72. Meakin, P.J.; Roberts, J.A. Anatomical and biochemical changes associated with the induction of oilseed rape (Brassica napus) pod dehiscence by Dasineura brassicae (Winn.). Ann. Bot. 1991, 67, 193–197. [Google Scholar] [CrossRef]
  73. Hensel, L.L.; Grbic, V.; Baumgarten, D.A.; Bleecker, A.B. Developmental and age-related processes that influence the longevity and senescence of photosynthetic tissues in Arabidopsis. Plant Cell 1993, 5, 553–564. [Google Scholar] [CrossRef] [PubMed]
  74. Nooden, L.D.; Penney, J.P. Correlative controls of senescence and plant death in Arabidopsis thaliana (Brassicaceae). J. Exp. Bot. 2001, 52, 2151–2159. [Google Scholar] [CrossRef] [PubMed]
  75. Marbà, N.; Duarte, C.M.; Agustí, S. Allometric scaling of plant life history. Proc. Natl. Acad. Sci. USA 2007, 104, 15777–15780. [Google Scholar] [CrossRef] [PubMed]
  76. Broadley, M.R.; Hammond, J.P.; King, G.J.; Astley, D.; Bowen, H.C.; Meacham, M.C.; Mead, A.; Pink, D.A.; Teakle, G.R.; Hayden, R.M.; et al. Shoot calcium and magnesium concentrations differ between subtaxa, are highly heritable, and associate with potentially pleiotropic loci in Brassica oleracea. Plant Physiol. 2008, 146, 1707–1720. [Google Scholar] [CrossRef] [PubMed]
  77. Lancashire, P.; Bleiholder, H.; Boom, T.; van den Bout-van den Beukel, C.J.P.; Langelüddeke, P.; Stauss, R.; Weber, E.; Velde, T.; Lancashire, P.D.; Boom, T.V.D.; et al. A uniform decimal code for growth stages of crops and weeds. Ann. Appl. Biol. 1991, 119, 561–601. [Google Scholar] [CrossRef]
  78. Kelly, A.A.; Shaw, E.; Powers, S.J.; Kurup, S.; Eastmond, P.J. Suppression of the SUGAR-DEPENDENT1 triacylglycerol lipase family during seed development enhances oil yield in oilseed rape (Brassica napus L.). Plant Biotechnol. J. 2013, 11, 355–361. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Responses of 35 varieties of Brassica napus to nitrogen fertilizers in Year 1. (a) Comparison of observed yield to that predicted by the statistical model for Year 1 data from low (18 kg ha−1; open circle) and high (148 kg ha−1; closed circle) nitrogen plots (R2 = 0.5404). (b) Size of the nitrogen effect for each variety. Bars = SEM. Varieties for which the addition of nitrogen had a significant impact on yield are marked with *.
Figure 1. Responses of 35 varieties of Brassica napus to nitrogen fertilizers in Year 1. (a) Comparison of observed yield to that predicted by the statistical model for Year 1 data from low (18 kg ha−1; open circle) and high (148 kg ha−1; closed circle) nitrogen plots (R2 = 0.5404). (b) Size of the nitrogen effect for each variety. Bars = SEM. Varieties for which the addition of nitrogen had a significant impact on yield are marked with *.
Agronomy 07 00031 g001
Figure 2. Comparison of observed yield of 35 varieties of Brassica napus to that predicted by the statistical model for Year 1 data using data from (a) high nitrogen (148 kg ha−1) plots (R2 = 0.4922) and (b) low nitrogen (18 kg ha−1) plots (R2 = 0.5221).
Figure 2. Comparison of observed yield of 35 varieties of Brassica napus to that predicted by the statistical model for Year 1 data using data from (a) high nitrogen (148 kg ha−1) plots (R2 = 0.4922) and (b) low nitrogen (18 kg ha−1) plots (R2 = 0.5221).
Agronomy 07 00031 g002
Figure 3. Models: (a) Comparing observed yield of Brassica napus in Year 2 to that predicted by the statistical model developed on Year 1 data using plants grown on low nitrogen (18 kg ha−1) plots. R2 = 0.2548. (b) Explaining yield in Year 2 using forward selection from 133 possible terms. The final model used nine terms to explain 95.12% of the yield variance using plants grown on low nitrogen (18 kg ha−1) plots in a randomised block design. R2 = 0.951. (c) Explaining vegetative biomass in Year 2 using forward selection from 133 possible terms. The final model used four terms to explain 80.04% of the biomass variance using plants grown on low nitrogen (18 kg ha−1) plots in a randomised block design. R2 = 0.8004.
Figure 3. Models: (a) Comparing observed yield of Brassica napus in Year 2 to that predicted by the statistical model developed on Year 1 data using plants grown on low nitrogen (18 kg ha−1) plots. R2 = 0.2548. (b) Explaining yield in Year 2 using forward selection from 133 possible terms. The final model used nine terms to explain 95.12% of the yield variance using plants grown on low nitrogen (18 kg ha−1) plots in a randomised block design. R2 = 0.951. (c) Explaining vegetative biomass in Year 2 using forward selection from 133 possible terms. The final model used four terms to explain 80.04% of the biomass variance using plants grown on low nitrogen (18 kg ha−1) plots in a randomised block design. R2 = 0.8004.
Agronomy 07 00031 g003
Figure 4. Correlation matrix showing the relationship between all traits measured from Brassica napus plants in Year 2. White squares = positive correlation of +1; red squares indicate a strong negative correlation. Note that where a measurement was taken on multiple occasions, some time points are not shown for clarity. Blue squares indicate regions of the matrix that are discussed in the main text.
Figure 4. Correlation matrix showing the relationship between all traits measured from Brassica napus plants in Year 2. White squares = positive correlation of +1; red squares indicate a strong negative correlation. Note that where a measurement was taken on multiple occasions, some time points are not shown for clarity. Blue squares indicate regions of the matrix that are discussed in the main text.
Agronomy 07 00031 g004
Figure 5. Multi-dimensional scaling of all traits measured from Brassica napus. The correlation between each pair of traits was converted into a distance, where a correlation of ±1 was given a distance of zero and a correlation of zero was given a distance of one; a two-dimensional projection of these distances is shown in the figure. Traits that are close together in this space are highly correlated. Traits are coloured according to the type of trait: measures of yield (black), architecture (red), lower leaf and upper leaf chlorophyll (green), pod chlorophyll (blue), timing of plant development (cyan), NIRS data (pink), early leaf and canopy leaf chlorophyll (yellow) and early leaf and canopy leaf minerals (grey). Traits used in the models for predicting yield and/or vegetative biomass are circled.
Figure 5. Multi-dimensional scaling of all traits measured from Brassica napus. The correlation between each pair of traits was converted into a distance, where a correlation of ±1 was given a distance of zero and a correlation of zero was given a distance of one; a two-dimensional projection of these distances is shown in the figure. Traits that are close together in this space are highly correlated. Traits are coloured according to the type of trait: measures of yield (black), architecture (red), lower leaf and upper leaf chlorophyll (green), pod chlorophyll (blue), timing of plant development (cyan), NIRS data (pink), early leaf and canopy leaf chlorophyll (yellow) and early leaf and canopy leaf minerals (grey). Traits used in the models for predicting yield and/or vegetative biomass are circled.
Agronomy 07 00031 g005
Figure 6. Cluster analysis of the Brassica napus traits and varieties. (a) Cluster dendrogram of the variety set; each cluster is blocked in green, black (grey), red or blue to correspond with Figure 6b and the colours used to circle varieties in Figure 7. The study included 35 different plant varieties. The mean of both blocks was used for each variety for each trait; each of the traits was scaled to have a mean = 0 and variance = 1, and then cluster analysis was performed by taking the Euclidean distance between each pair of varieties and combining varieties using Ward’s method. The distance between varieties is defined as equal to the sum of the squares of differences across all traits analysed in Figure 5. (b) Analysis of the traits contributing to the yield and vegetative biomass models clustered by variety.
Figure 6. Cluster analysis of the Brassica napus traits and varieties. (a) Cluster dendrogram of the variety set; each cluster is blocked in green, black (grey), red or blue to correspond with Figure 6b and the colours used to circle varieties in Figure 7. The study included 35 different plant varieties. The mean of both blocks was used for each variety for each trait; each of the traits was scaled to have a mean = 0 and variance = 1, and then cluster analysis was performed by taking the Euclidean distance between each pair of varieties and combining varieties using Ward’s method. The distance between varieties is defined as equal to the sum of the squares of differences across all traits analysed in Figure 5. (b) Analysis of the traits contributing to the yield and vegetative biomass models clustered by variety.
Agronomy 07 00031 g006
Figure 7. Varietal distribution across the Brassica napus diversity trial in two-dimensional space using vegetative biomass and oil:protein ratio. The yield of each variety is shown by the radius of the circle; colours are used to identify the varieties within each cluster as determined in Figure 6. Green boxes indicate yield gaps where no varieties currently exist.
Figure 7. Varietal distribution across the Brassica napus diversity trial in two-dimensional space using vegetative biomass and oil:protein ratio. The yield of each variety is shown by the radius of the circle; colours are used to identify the varieties within each cluster as determined in Figure 6. Green boxes indicate yield gaps where no varieties currently exist.
Agronomy 07 00031 g007
Figure 8. Scores for nine traits associated with resource allocation and yield for the thirty-five varieties used in this study, ordered by yield (outer ring) and coloured from largest (red) to smallest (green) in each ring. Averages of all replicates over both blocks were used. White spaces = missing data.
Figure 8. Scores for nine traits associated with resource allocation and yield for the thirty-five varieties used in this study, ordered by yield (outer ring) and coloured from largest (red) to smallest (green) in each ring. Averages of all replicates over both blocks were used. White spaces = missing data.
Agronomy 07 00031 g008
Figure 9. Contrast in uniformity of Brassica napus. (a) Pod senescence on primary stems and (b) Canopy senescence. Canard, Darmor and Victor TVSL are the three highest yielding lines; Ningyou7, POH285 and Shen-li Jutsai are the three lowest yielding in terms of seed mass per plant. Dates in (a) are the dates the images were taken. All images in (b) were taken on 13 June 2011, with the dates shown in brackets being the final harvest date for each line.
Figure 9. Contrast in uniformity of Brassica napus. (a) Pod senescence on primary stems and (b) Canopy senescence. Canard, Darmor and Victor TVSL are the three highest yielding lines; Ningyou7, POH285 and Shen-li Jutsai are the three lowest yielding in terms of seed mass per plant. Dates in (a) are the dates the images were taken. All images in (b) were taken on 13 June 2011, with the dates shown in brackets being the final harvest date for each line.
Agronomy 07 00031 g009
Figure 10. Plant traits as drivers for final yield in Brassica napus. (a) Top 10 lines in each category showing overlap between them. Temple is one of the most commonly-grown varieties, yet it is atypical in linking chlorophyll retention to high yield. In general, biomass and yield overlap (big plants = more seed) and varieties with a delayed senescence phenotype do not yield highly. (b) Model for predicting yield using traits that were observed at or before the time of flowering. The model was developed using nine terms and applied to Year 2 data collected from plants grown on low nitrogen (22 kg ha−1) plots in a randomised block design. R2 = 0.6369.
Figure 10. Plant traits as drivers for final yield in Brassica napus. (a) Top 10 lines in each category showing overlap between them. Temple is one of the most commonly-grown varieties, yet it is atypical in linking chlorophyll retention to high yield. In general, biomass and yield overlap (big plants = more seed) and varieties with a delayed senescence phenotype do not yield highly. (b) Model for predicting yield using traits that were observed at or before the time of flowering. The model was developed using nine terms and applied to Year 2 data collected from plants grown on low nitrogen (22 kg ha−1) plots in a randomised block design. R2 = 0.6369.
Agronomy 07 00031 g010
Table 1. Statistical model to explain log[yield].
Table 1. Statistical model to explain log[yield].
CoefficientsEstimateStandard Errort-ValuePr (>|t|)SignificanceCumulative % Variance Explained
(Intercept)2.4720.2858.6721.42 × 10−14***-
Log(canopy leaf Chl a)0.4910.1792.7400.0070**28.81
Seed linolenic acid−0.0580.016−3.6350.0004***44.05
Seeds per pod0.0190.0082.2870.0238*48.70
Nitrogen applied to plot0.3290.0883.7220.0003***52.91
Log(seed glucosinolate)−0.0670.031−2.2000.0296*54.59
Forty-five traits were considered as candidates for terms in the model. Starting with no terms in the model, terms were added using forward selection (see the Materials and Methods). Traits that do not appear in the model did not significantly improve the estimation of log[yield]. The data modelled represented 35 varieties grown at residual N (~18 kg ha−1) and high N (148 kg ha−1) concentrations in duplicate plots in a randomised block design. Residual standard error: 0.3629 on 131 degrees of freedom; multiple R-squared: 0.5459; adjusted R-squared: 0.5285; significance codes: 0.001 ‘***’; 0.01 ‘**’; 0.05 ‘*’.
Table 2. Statistical model to explain log[yield] under high and low nitrogen.
Table 2. Statistical model to explain log[yield] under high and low nitrogen.
High Nitrogen CoefficientsEstimateStandard Errort-ValuePr (>|t|)SignificanceCumulative % Variance Explained
(Intercept)1.6490.8331.9800.0518.-
Seed protein content−0.0510.018−2.7670.0073**33.77
Log[early leaf Mn]0.6320.2232.8360.0060**40.51
Pod dehiscence Low0.6720.2822.3790.0202*-
Pod dehiscence Medium0.7310.2780.2680.0107*46.17
Log[seed glucosinolate]−0.1080.054−2.0190.0475*49.26
Low Nitrogen Coefficients
(Intercept)2.1130.3426.1807.28 × 10−8***-
Seed linolenic acid−0.0330.021−1.5970.1159-35.59
Seeds per pod0.0280.0092.9950.0041**41.66
Pod development: Yellow stage−0.2760.243−1.1370.2602--
Pod development: Yellow-green stage−0.0480.216−0.2220.8253-52.67
Forty-five traits were considered as candidates for terms in the model. Starting with no terms in the model, terms were added using forward selection (see Methods). Traits that do not appear in the model did not significantly improve the estimation of log[yield]. The data modelled represented 35 varieties grown at high N (148 kg ha−1) and low (18 kg ha−1) concentrations in duplicate plots in a randomised block design. Residual standard error: 0.3795 on 67° of freedom; multiple R-squared: 0.4926; adjusted R-squared: 0.4547; residual standard error: 0.2931 on 57° of freedom; multiple R-squared: 0.5267; adjusted R-squared: 0.4769; significance codes: 0.001 ‘***’; 0.01 ‘**’; 0.05 ‘*’; 0.1 ‘.’.
Table 3. Statistical model to explain log[yield] using the Year 1 model on Year 2 data, the forward selection statistical model on Year 2 data and the forward selection statistical model to explain the log[biomass] term from Year 2 data.
Table 3. Statistical model to explain log[yield] using the Year 1 model on Year 2 data, the forward selection statistical model on Year 2 data and the forward selection statistical model to explain the log[biomass] term from Year 2 data.
Year 1 Model On Year 2 Data
CoefficientsEstimateStandard Errort-ValuePr (>|t|)SignificanceCumulative % Variance Explained
(Intercept)2.1860.3256.7303.65 × 10−9***-
Linolenic acid−0.0080.023−0.3570.7224--
Seeds per pod0.0120.0071.7770.0798.-
Pod development: Yellow stage0.1920.1031.8630.0666.-
Pod development: Yellow-green stage0.3320.1033.2190.0019**-
Year 2 Model on Year 2 Data
(Intercept)−2.8440.494−5.7572.42 × 10−7***-
Log[Vegmass]1.2620.04230.377<2 × 10−16***73.70%
Total Protein−0.0190.010−1.9520.0551.89.16%
Canopy leaf Chlorophyll b0.3370.0724.6661.55 × 10−5***91.41%
Total Lipid0.0150.0062.7360.0204*92.04%
Seed weight per pod0.0010.0004.0360.0001***92.80%
No. pods on primary raceme0.0070.0014.4273.67 × 10−5***93.50%
Canopy leaf calcium−0.1250.029−4.3365.07 × 10−5***94.40%
Linolenic acid−0.0210.008−2.5170.0143*94.81%
log(Pod Chlorophyll a content Wk43)−0.0240.012−2.0460.0448*95.12%
Log[Vegmass] Model on Year 2 Data
(Intercept)2.4230.10822.435<2 × 10−16***-
Stem area0.0010.0003.4070.0011**47.16%
Log[No. branches on primary stem]0.0290.0039.5742.06 × 10−14***59.93%
Height0.0050.0017.3732.42 × 10−10***77.74%
Total branch No.0.0340.0049.7932.04 × 10−14***75.63%
Pod carotenoid content (Wk41)0.0210.0072.8570.005**80.04%
The data modelled represented 35 varieties grown at residual N (~18 kg ha−1) in duplicate plots in a randomised block design. Residual standard error: 0.2893 on 71° of freedom; multiple R-squared: 0.2548; adjusted R-squared: 0.2128; residual standard error: 0.07675 on 66° of freedom; multiple R-squared: 0.9512; adjusted R-squared: 0.9446; residual standard error: 0.102 on 71° of freedom; multiple R-squared: 0.8004; adjusted R-squared: 0.7891; significance codes: 0.001 ‘***’; 0.01 ‘**’; 0.05 ‘*’; 0.1 ‘.’.
Table 4. Forward selection statistical model to predict the yield from traits observable at or prior to flowering using Year 2 data. The data modelled represented 35 varieties grown at residual N (~18 kg ha−1) in duplicate plots in a randomised block design.
Table 4. Forward selection statistical model to predict the yield from traits observable at or prior to flowering using Year 2 data. The data modelled represented 35 varieties grown at residual N (~18 kg ha−1) in duplicate plots in a randomised block design.
CoefficientsEstimateStandard Errort-ValuePr (>|t|)SignificanceCumulative % Variance Explained
(Intercept)−0.3770.3591.0500.2975--
Stem area0.0020.0013.1500.00241**36.64%
Log[No. branches on primary stem]0.7240.1514.8118.51 × 10−6***44.22%
Length of flower set0.0090.0033.2940.0016**51.73%
Canopy leaf potassium−0.1880.057−3.2930.0016**56.92%
Canopy leaf chlorophyll b0.5500.1852.9720.0041**60.19%
Canopy leaf Mg1.5940.6172.5820.0120*63.69%
Residual standard error: 0.2048 on 69° of freedom; multiple R-squared: 0.6369; adjusted R-squared: 0.6054; significance codes: 0.001 ‘***’; 0.01 ‘**’; 0.05 ‘*’.

Share and Cite

MDPI and ACS Style

Bennett, E.J.; Brignell, C.J.; Carion, P.W.C.; Cook, S.M.; Eastmond, P.J.; Teakle, G.R.; Hammond, J.P.; Love, C.; King, G.J.; Roberts, J.A.; et al. Development of a Statistical Crop Model to Explain the Relationship between Seed Yield and Phenotypic Diversity within the Brassica napus Genepool. Agronomy 2017, 7, 31. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy7020031

AMA Style

Bennett EJ, Brignell CJ, Carion PWC, Cook SM, Eastmond PJ, Teakle GR, Hammond JP, Love C, King GJ, Roberts JA, et al. Development of a Statistical Crop Model to Explain the Relationship between Seed Yield and Phenotypic Diversity within the Brassica napus Genepool. Agronomy. 2017; 7(2):31. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy7020031

Chicago/Turabian Style

Bennett, Emma J., Christopher J. Brignell, Pierre W. C. Carion, Samantha M. Cook, Peter J. Eastmond, Graham R. Teakle, John P. Hammond, Clare Love, Graham J. King, Jeremy A. Roberts, and et al. 2017. "Development of a Statistical Crop Model to Explain the Relationship between Seed Yield and Phenotypic Diversity within the Brassica napus Genepool" Agronomy 7, no. 2: 31. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy7020031

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

Article Metrics

Back to TopTop