Genetic effects and genotype × environment interactions govern seed oil content in Brassica napus L.

As seed oil content (OC) is a key measure of rapeseed quality, better understanding the genetic basis of OC would greatly facilitate the breeding of high-oil cultivars. Here, we investigated the components of genetic effects and genotype × environment interactions (GE) that govern OC using a full diallel set of nine parents, which represented a wide range of the Chinese rapeseed cultivars and pure lines with various OCs. Our results from an embryo-cytoplasm-maternal (GoCGm) model for diploid seeds showed that OC was primarily determined by genetic effects (VG) and GE (VGE), which together accounted for 86.19% of the phenotypic variance (VP). GE (VGE) alone accounted for 51.68% of the total genetic variance, indicating the importance of GE interaction for OC. Furthermore, maternal variance explained 75.03% of the total genetic variance, embryo and cytoplasmic effects accounted for 21.02% and 3.95%, respectively. We also found that the OC of F1 seeds was mainly determined by maternal effect and slightly affected by xenia. Thus, the OC of rapeseed was simultaneously affected by various genetic components, including maternal, embryo, cytoplasm, xenia and GE effects. In addition, general combining ability (GCA), specific combining ability (SCA), and maternal variance had significant influence on OC. The lines H2 and H1 were good general combiners, suggesting that they would be the best parental candidates for OC improvement. Crosses H3 × M2 and H1 × M3 exhibited significant SCA, suggesting their potentials in hybrid development. Our study thoroughly investigated and reliably quantified various genetic factors associated with OC of rapeseed by using a full diallel and backcross and reciprocal backcross. This findings lay a foundation for future genetic studies of OC and provide guidance for breeding of high-oil rapeseed cultivars.


Background
The seed oil content (OC) is a key measure of rapeseed quality and is also a complicated quantitative trait easily affected by the environment and difficult to investigate [1][2][3]. Previous studies have demonstrated that the OC of rapeseed is mainly controlled by genotype and genotype × environment interactions (GE) [4][5][6]; in addition, it is governed by multiple genes mainly through additive effect, and thus can be altered through breeding and selection [7][8][9][10].
A previous study on summer rapeseed has suggested that OC might be primarily controlled by maternal factors or embryo genotype or xenia [4]. The strong influence of maternal effect on the OC of F 1 seeds is usually accompanied by weak xenia [11][12][13]. For maternal effect, several forms have been proposed, including the maternal inheritance of plastid, endosperm, seed coat, and maternal provision of nutrients [14,15]. Seed lipid synthesis is independent of the leaf photosynthesis and the phloem transport of photosynthate [16], but mainly requires the supply of photosynthate from the silique wall [6,12]. Photosynthesis of the silique wall, sugar transport in the seed coat, and the expression of fatty acid synthesis-related genes in the embryo can significantly influence the OC [17,18]. In addition, the storage substance in seed is determined not only by the availability of assimilates (source strength), but also by the intrinsic traits of the seed (sink strength), which are controlled by the embryo genotype [19]. Therefore, variation in OC of rapeseed may be governed by multiple genetic components, including embryo, cytoplasmic, xenia, maternal, and GE effects [11-13, 20, 21].
Xenia, which represents the effect of pollen on the development and characters of seed, can be demonstrated by analyzing the differences between seeds fertilized using different pollen sources from the same plant [22,23]. The maternal effect on OC can be investigated using reciprocal crosses [24,25], and ancillary data from backcross progeny can be used to separate cytoplasmic effect from maternal effect [26,27]. A significant difference between reciprocal backcrosses is strongly indicative of cytoplasmic effect. Cytoplasmic effect and maternal effect may also be distinguished by comparing reciprocal F 2 seeds [28]. In contrast, a difference between reciprocal F 1 hybrids that does not exist between reciprocal F 2 hybrids would indicate that OC is determined by maternal effect without cytoplasmic effect [11].
Diallel mating designs have been frequently used to investigate the genetic effects of parents or to determine which cultivars are the best combiners for favorable alleles in hybrids. Diallel analysis provides information on the genetic behaviors of these attributes in the F 1 generation [29]. The four methods of Griffing have usually been used to obtain genetic information on the basis of data from only one year or one location [30]; however, it has been suggested that considering multiple types of environment data could provide more reliable genetic information on the material tested [31]. The variance of general combining ability (GCA) incorporates additive epistasis, whereas that of specific combining ability (SCA) incorporates dominance epistasis [30,32]. The observation of high additive effect for a specific trait indicates higher heritability and less environmental influence and will facilitate the selection of this trait [33]. Additive effect efficiently responds to selection, whereas non-additive effects, such as dominance and epistatic components, increase hybrid vigor in the cross combinations of cultivars. These facts suggest that the evaluation of GCA of a specific trait can guide the breeders to select the parents that can be used in breeding program for that trait, while SCA indicates the heterosis of a specific trait, and significant SCA of a cross suggests the presence of non-additive gene action [34].
In the present study, we investigated whether OC is determined by embryo, maternal, cytoplasmic, xenia or GE effects, or a combination of these factors. To this end, we estimated the components of different genetic systems and their corresponding GE. In addition, we investigated the roles of GCA, SCA, and reciprocal effect in the inheritance of OC. An improved understanding of the genetic components of OC and the combining abilities of inbred lines will help breeders to develop high-oil rapeseed cultivars and hybrids for particular geographic locations.

Materials and field experiments
The experiments were carried out from September of 2009 to May of 2012. Nine semi-winter rapeseed lines with differences in seed oil content (OC) were selected as parental lines from the Rapeseed Laboratory of Huazhong Agricultural University (Table 1). Three high-oil lines (HO; H1, H2, H3), three medium-oil lines (MO; M1, M2, M3), and three low-oil lines (LO; L1, L2, L3), were crossed in a 9 × 9 diallel mating scheme to produce 72 F 1 hybrids (F 1 seeds from maternal plants), including 36 crosses and 36 of their reciprocals in March of 2010. These 72 crosses were performed again in March of 2012, and these F 1 seeds were only analyzed for maternal effect and xenia on the OC of rapeseed.
A total of 81 genetic entries (72 F 1 hybrids from maternal plants and 9 parents) obtained from the complete diallel cross of the nine parents in May of 2010, were directly sown in the field at Huazhong Agricultural University, Wuhan, China around September 27th of both 2010 and 2011. The parental and F 1 seeds were planted in a randomized block design with three replications over two consecutive growing seasons, and each block contained two rows with 10 plants in each row at a space of 30 cm × 15 cm. The OC of self-pollinated Environmental factors, such as temperature, light, and weather during seed development and difference in flowering time, could influence OC and the results of the study. Therefore, since F 1 plants have the same developmental timing, backcross using an F 1 plant as the maternal parent could be used to solve these problems [27]. HO (H2, H3) and LO (L1, L2) were selected for producing backcross (BC) and reciprocal backcross (RBC) generations. The four parents and their eight F 1 hybrids produced in March of 2010 were directly sown in the field at Huazhong Agricultural University, Wuhan, China around September 27th of both 2010 and 2011. In the spring of both 2011 and 2012, a total of 32 BC and RBC combinations were produced, including eight BC 1 s, eight BC 2 s (with F 1 plants as female parents), eight RBC 1 s, and eight RBC 2 s (with F 1 plants as male parents), and each of the combinations was repeated three times and used to study the maternal, xenia and cytoplasmic effects on OC.
All plants were isolated by paper bags at the beginning of flowering to ensure self-pollination of parents and F 2 seeds, and the bags were removed at the end of full flowering. All plants were harvested and threshed at early May in 2010, 2011 and 2012. The total OC was measured using the Foss NIR-Systems 5000 near-infrared reflectance spectroscope (NIR-Systems, Inc., Silver Spring, MD, USA) [35], with the parameters described by Gan et al. [36].

Statistical analysis
Data were analyzed using a general linear model (GLM) [37]. Based on the GLM, least-square means were used to compute the combining ability based on Griffing's diallel analysis with Method 1 (full diallel set), Model 1 [30]. All variance analyses and combining ability (GCA, SCA, reciprocal) estimation were performed using the SAS codes published by Zhang and Kang [31]. On the basis of Cong [38] and Duan et al. [39], the method proposed by Wang et al. [11] was used to estimate the genetic components of maternal effect and xenia in rapeseed. The OC of F 1 seeds was calculated as F 1 = MP 1 + XP 2 , where M was the value of maternal effect and X was the value of xenia (X = 1-M).
The partitions of embryo, maternal, cytoplasmic effects as well as the corresponding GE were estimated using the embryo-cytoplasm-maternal (GoCGm) model [40,41] for diploid seeds in QGAStation1.0 (http://ibi. zju.edu.cn/software/qga/index.htm). The genetic variance components were estimated using the minimum norm quadratic unbiased estimation (MINQUE) (0/1) method [42], and the genetic effects of each parent were investigated using the adjusted unbiased prediction (AUP) method [43]. Standard errors of the estimated variance and predicted genetic effects were then analyzed using a Jackknife procedure [44], and t-tests were used to test for significant differences in the traits examined.

Phenotypic variation
We found that the self-pollinated seeds from nine parents, including HO (H1, H2, and H3), MO (M1, M2, and M3), and LO (L1, L2, and L3), exhibited significant differences in OC, and that the differences between the parents with the highest and lowest OC in each year was about 12.00% (Tables 1, 2). The OC of the 72 F 1 hybrids was strongly influenced by the maternal parent. In 2012, the mean OC of the eight F 1 lines with H1 as the female parent was 42.01%, which was comparable to that of H1 (42.79%), whereas that of the eight F 1 lines with H1 as Table 2 Oil contents (%) of the complete diallel crosses F 1 seeds harvested in 2010 and 2012 (Mean ± SE) Lines 2010 2012 Mean value of the F 1 derived from the line as male crossed with other 8 lines the male parent was 38.04% (Table 2). Similarly, the mean OC of the eight F 1 lines with L3 as the female parent was 32.44%, which was comparable to that of L3 (29.63%), and that of the eight F 1 lines with L3 as the male parent was 37.10%, which was much higher than that of L3. Similar trends were observed in the other sets of crosses in both 2010 and 2012. Thus, the OC of F 1 hybrids was similar to that of maternal parent, even with different male parents. The results in Table 3 showed that the mean OC of open-pollinated F 2 seeds was about 2% higher than that of self-pollinated F 2 seeds, and the mean values for most of the F 2 seeds from reciprocal crosses were more or less similar, despite a small amount of inbreeding depression. In addition, differences between the F 1 reciprocal crosses disappeared in the F 2 generation for most combinations, thus demonstrating a maternal effect and little or no cytoplasmic effect (Tables 2, 3). The analysis of backcrosses confirmed the predominant influence of maternal parent on OC, even with the elimination of differences in flowering time (Fig. 1). Furthermore, the mean OCs of (H2/L1)H2 and (H2/L1)L1 were similar (Fig. 1a, e), and it was the same case for (L1/H2)H2 and (L1/H2)L1 (Fig. 1a, e). However, when H2 was used as the maternal parent to cross with four F 1 hybrid lines (L1/H2, H2/L1, H2/L2, L2/H2), the various backcrosses exhibited different OCs, indicating the presence of a slight xenia (Fig. 1a, b, e, f ). Similar trends were observed in other three sets of backcrosses, confirming that maternal effect was the main determinant of OC and that the influence of pollen source was minor (Fig. 1).

Variance analysis of means and gene action
A combined analysis of variance also indicated significant differences in OC among genotypes ( Table 4). The significant mean squares of GE for OC indicated that the magnitude of the trait in different genotypes varied over the two years, and the differences between plants in different blocks were also significant, indicating that OC was easily affected by the environment. GCA and SCA also significantly contributed to the OC variation over the years, indicating the importance of both additive and non-additive effects on OC (Table 4). Further partitioning of reciprocal sum squares indicated that maternal effect was significant, whereas non-maternal effect was not, suggesting that OC was not under strict nuclear control and could also be influenced by cytoplasm. The GCA × Environment showed significant effect on OC, suggesting that the additive variances were influenced by the environment. However, the SCA × Environment did not show any effect. Therefore, the effect of environment on non-additive components of genetic variance was not significant.

Combining ability
The GCA was significant for all parents even though some parents had positive values and other had negative values, whereas no parental line exhibited significant maternal effect (Table 5). H2 (3.21 ** ) and H1 (2.87 ** ) had the highest GCA for OC, whereas L3 (−3.06 ** ) and L2 (−2.06 ** ) were negative combiners with reduced OC in F 1 . The SCA in several crosses was significant, and reciprocal effect was only significant in two crosses. Lines H3 and M2 showed significantly positive GCA, whereas M1, M3 and L1 exhibited significantly negative GCA. Crosses H1 × M3, H1 × L1, H1 × L3, H2 × H3, H2 × M2 and H3 × M2 exhibited significantly positive SCA and were the best specific combiners for improvement of OC (Table 5), The negative SCA of H1 × H2, H1 × H3, H1 × M2, H3 × M1 and L2 × L3 observed with lower OC suggested that they were poor parental combinations for breeding. The reciprocal effect was estimated to be not significant in most reciprocal crosses, except for the reciprocals of L3 × M3 and L1 × H3. The cross L3 × M3 showed significant positive reciprocal effect indicating a negative cytoplasmic effect from the maternal parent L3. The negative reciprocal effect of L1 × H3 indicated a positive cytoplasmic effect from maternal parent L1. Therefore, L3 and L1 should be respectively exploited as male and maternal parent in crosses for OC improvement.  (Table 6). On average across the two years, maternal effect accounted for 0.87 of the observed variation, whereas xenia accounted for only 0.13 of it, confirming a very strong maternal effect on the OC of F 1 seeds and a weak effect of xenia. Some crosses (H1 × MO and M3 × LO) exhibited values of maternal effect around 1 consistently across the two years, suggesting 100% maternal influence. Some crosses (H3 × LO and M1 × LO) showed maternal effect and xenia varying markedly between two years and also had large standard errors (SE), suggesting that these crosses had relatively weak maternal effects.

Cytoplasmic effect
Analyses of backcrosses revealed the minor or negligible influence of cytoplasmic effect, although the backcrosses also displayed some variations between 2011 and 2012 ( Fig. 1), which demonstrated that environment might have some influence on cytoplasmic effect. In 2011, for example, the average OC of (H3/L2)H3 was higher than that of (L2/H3)H3, and that of (H3/L2)L2 was higher than that of (L2/H3)L2 as well. However, in 2012, the mean OC of (H3/L2)L2 was higher than that of (L2/H3)L2, which might be owing to positive cytoplasmic effect, whereas that of (H3/L2)H3 was lower than that of (L2/H3)H3, which might be due to negative cytoplasmic effect. The mean OC of (H2/L2)L2 was higher than that of (L2/H2)L2, and that of (H2/L2)H2 was higher than that of (L2/H2)H2, which could be attributed to the positive cytoplasmic effect (Fig. 1b, f). The average OC of (H2/L1)H2 was higher than that of (L1/H2)H2, whereas that of (H2/L1)L1 was similar to that of

Components of total genetic variance
The results from GoCGm model showed that 13.81% of the phenotypic variance (V P ) in OC could be attributed to environmental variations and experimental errors, whereas the rest (86.19%) was attributable to  (Table 7). Genetic and GE variance components respectively accounted for 48.32% and 51.68% of the total genetic variance (V G + V GE ). Embryo additive (V A ), embryo additive interaction (V AE ), maternal additive (V Am ) and maternal dominant interaction (V DmE ) variance components were significant, whereas embryo dominant interaction (V DE ), cytoplasm interaction (V CE ), maternal additive (V AmE ), and maternal dominant (V Dm ) variance components were not. Maternal variances (V Am + V Dm + V AmE + V DmE ) explained 75.03% of the total genetic variance, whereas embryo and cytoplasmic effects accounted for only 21.02% and 3.95% of it, respectively. Overall, these results indicated that OC was predominantly influenced by maternal effect, followed by embryo and cytoplasmic effects.

Estimation of genetic components of parents
H2 showed the highest additive effect, whereas L3 displayed the lowest additive effect, and we also found that H1, H2, H3, M1 and L1 exhibited significant positive maternal additive effect, whereas other parents showed significant negative maternal additive effect ( Table 8). L2 exhibited the lowest maternal additive effect, whereas H1 showed the highest maternal additive effect. Cytoplasmic effects were positive in H1, H2, M1, M3 and L1, whereas the cytoplasmic effects of the remaining four lines were insignificantly or significantly negative. According to our estimated values, H1 exhibited the highest cytoplasmic effect (1.79 ± 0.86) among the nine parents, indicating the high probability of generating high-oil rapeseed with this line being used as the maternal parent. In addition, H2, M1, M2, and L2 displayed negative homozygous dominance, whereas the remaining five lines exhibited positive homozygous dominance (Table 8). L3 showed the highest homozygous dominance effect, indicating that it is more likely produce high heterosis. Moreover, all of the genetic effects of H1 were positive, suggesting that it might be an ideal maternal parent for breeding high-oil rapeseed.

Discussion
Investigation of the influence of different genetic systems can facilitate a better understanding of the nature of gene interactions that could influence OC. Our study demonstrates that the variation of OC is mainly determined by genetic and GE components (Table 7), though the environment also has a significant influence [7,8]. However, the influence of GE should not be neglected, since GE accounted for 51.68% of the total genetic variance. Significant interactions have been observed between  male or female and the environment [13]. GE is a critical factor for developing varieties with wide geographical adaptability and should be taken into account in genetic model and breeding process. Data from the genetic model clearly demonstrated that OC was mainly controlled by the maternal parent, with the maternal effect accounting for 75% of the genetic variance. Embryo effect accounted for 21% of the genetic variance and cytoplasmic effect was detectable only at a very low level ( Table 7). The present study confirms the earlier findings that the OC of F 1 seeds is mainly controlled by maternal effect, and xenia is weak [11,13,45]. Analysis of reciprocal backcrosses demonstrates that maternal, xenia, and cytoplasmic effect can all influence OC (Fig. 1). Therefore, our genetic analysis confirms that OC is determined by maternal, embryo, cytoplasm, xenia and GE effects. Previous studies have shown that the photosynthesis in silique wall makes a crucial contribution to OC [6,12,16,17], providing an explanation for the influence of maternal parent on OC. Therefore, the selection for high-oil rapeseed would be more effective based on the photosynthetic activity of silique wall than based on the performance of maternal parent. It should be noted that xenia has a direct genetic effect on OC in many crops [11,12,46,47], and is applied not only in genetic and physiological research but also in crop breeding and production. Similar to the case of corn and soybean, controlling the pollen source and considering xenia are required in both breeding programs and investigation of OC in rapeseed. Therefore, decision on parents and cross direction are very important for hybrid rapeseed breeding and production.
The OC of self-pollinated seeds was lower than that of open-pollinated seeds harvested from the same plant in the F 2 generation (Table 3). This result is similar to the finding of Hom et al. [48]. Earlier studies reported that the OC of F 1 seeds was significantly lower than that of  self-pollinated seeds when a HO female was crossed with a LO male, but it was significantly higher than that of self-pollinated seeds when a LO female was crossed with a HO male [11,12]. Similar results were obtained for soybean [49]. Although xenia has a direct genetic effect on OC of rapeseed [11,12,45], different pollens may not actually be responsible for the difference in OC between self-pollinated and open-pollinated seeds. During seed development, self-pollinated seeds were isolated using paper bags, whereas open-pollinated seeds were not. Thus, the difference in microclimate may have influenced the OC. Concurrently elevated CO 2 and temperature could reduce seed biomass by half and further reduce losses in fatty acids and OC [50]. Seeds grown in high-light environments tend to have higher OC [51]. Hua et al. [12] suggested that local and tissue-specific photosynthetic activity in the silique wall were the main determinant for OC, and factors of weather and temperature could also influence photosynthesis of the silique wall, thereby affecting OC. Thus, the influence of microclimate on the photosynthesis of silique wall might explain the lower OC of self-pollinated seeds. Our analysis of combining ability indicated the presence of both additive and non-additive gene actions in the parental lines (Table 4). Additive effect is equivalent to GCA, parental lines with high GCA for a specific trait may be better candidates as parental lines [30,32,34,52]. Lines that had positive additive and maternal additive effects, such as H1, H2 and M1, are desirable general combiners that can be used for OC improvement. Cytoplasmic genes can persist through generations and are also expressed as additive effect [52]. Therefore, lines with positive cytoplasmic effect, such as H1, H2 and L1, could be used as maternal parents. It is evident that H1 and H2 are good combiners and that the genotypes might be the best candidates as maternal parents for improvement of OC. Hybrids in both pol and mur cytoplasmic male sterility systems exhibited significantly lower OC (by 1.3% and 1.4%, respectively) than identical hybrids constructed in the cytoplasm (nap) common in rapeseed [53,54]. Line H1, which had positive additive, maternal additive, homozygous maternal dominance, and cytoplasmic effects, would be a more effective parent for high-oil rapeseed breeding. Crosses such as H3 × M2 and H1 × M3, which exhibited significant SCA, could be used in the development of hybrid varieties.

Conclusion
Based on analyses of a 9 × 9 full diallel scheme and selected backcross and reciprocal backcross, we concluded that the OC is simultaneously controlled by genetic components of maternal, embryo, xenia, cytoplasmic, and GE effects in rapeseed. Maternal effect is the most important factor, accounting for 75% of the genetic variance, followed by embryo effect, which accounts for 21% of the genetic variance. Therefore, selection of maternal parents is paramount in the genetic improvement of OC in rapeseed. Together with previous studies, additional information regarding the role of genetic components in determining OC could help breeder to better manipulate the OC of rapeseed.