Pollen Grain Number and Viability Contribute to Variation in Effective Ovule Number Among Oilseed Rape Genotypes

Pollen Grain Number and Viability Contribute to Variation in Effective Ovule Number Among Oilseed Rape Genotypes | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Pollen Grain Number and Viability Contribute to Variation in Effective Ovule Number Among Oilseed Rape Genotypes Xiao Guo, Bao-Luo Ma, Neil B. McLaughlin, Yajun gao, Xiaoming Wu, and 1 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-830636/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 21 Dec, 2021 Read the published version in Euphytica → Version 1 posted 4 You are reading this latest preprint version Abstract While no significant differences in initial ovule number were found among oilseed rape genotypes, there was a large variation in effective ovule number (EON), which determines the final seeds per silique (SPS), a critical component of yield. In this study, we selected 18 oilseed rape genotypes with contrasting nitrogen utilization efficiency (NUtE) to unravel the main factors responsible for different EON and determine the critical period of EON formation under both a field and a pot experiments from 2016-2018. The high NUtE genotypes displayed significantly higher NUtE by 14.3%, along with greater yield per plant (29.4%) and SPS (21.1%) than the low NUtE genotypes. The greater productivity of the high NUtE genotypes was associated with 44.1% higher pollen grain number, 23.5% greater pollen viability, and 39.3% lower ovule abortion rate, compared to the low NUtE genotypes. In addition, the heart stage was the critical ovule development period for delineating the variability of EON among contrasting NUtE oilseed rape genotypes, when the high NUtE genotypes displayed higher silique net photosynthetic rate, surface area, biomass, and RNA expression levels. Taken together, this study indicated the pollen grain number, pollen viability and ovule abortion rate contributed to the final variation in EON and the heart stage was the critical period of determining the EON differences among contrasting NUtE genotypes. Increasing pollen grain number and pollen viability, and decreasing ovule abortion rate before heart stage should be the prerequisite for breeders to improve yield and NUtE of oilseed rape genotypes. Agronomy Molecular Biology General Biochemistry effective ovule number pollen grain number pollen viability ovule abortion rate heart stage Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Introduction Oilseed rape ( Brassica napus L.) is the world’s second largest crop source of vegetable oil following soybean, with increasing significance in the international market. Yield, one of the most important traits in oilseed rape breeding, is the product of seeds per silique (SPS), siliques per plant (SPP) and thousand seed weight (TSW)(Ma et al. 2015 ). As with other crops, yield components are critical to yield formation and often compensate each other, depending on genotypes and interaction with environment (Yang et al. 2017a ; Shi et al. 2015 ). Additionally, the three yield components (SPS, SPP and TSW) of oilseed rape display different degree of negative correlation, which indicates that the final seed yield can be improved by increasing the individual yield components (such as SPS) (Yang et al. 2017b ). SPS, which is usually negatively correlated with the other two yield components (SPP and TSW) for competition among carbohydrate sinks, has received much attention (Zhang et al. 2012 ; Yang et al. 2016 ; Wang et al. 2016b ; Yang et al. 2017a ; Gardner and Latta 2007 ). On the other hand, in the study conducted by He et al. ( 2017 ), the SPS was positively correlated with nitrogen utilization efficiency (NUtE; seed dry weight / shoot N accumulation), which indicated that SPS can be used as a simple and applicable indicator for the evaluation of NUtE in oilseed rape. During the pollination process, pollen grains have to reach and germinate on the stigma, and then form pollen tubes to accurately reach the female gametes. After that, fertilization is initialized once the male gametes meet the female gametophyte within the ovule (Palanivelu and Tsukamoto 2012 ). Thus, successful reproduction relies on the pollen grains number, grain viability and female ovule number (Gao et al. 2017 ; Li et al. 2015 ). Following fertilization, the developmental process from embryogenesis to final seed can be divided into five distinct stages: pre-embryo, globular, heart, torpedo and maturation (Andriotis et al. 2010 ; Tan et al. 2011 ; Hehenberger et al. 2012 ). At the same time, the fertilized ovules may fail at any one of the stages due to the internal and external environmental conditions, and the remaining ovules (effective ovule number, EON) eventually develop into viable seeds (SPS). Therefore, to some degree, the EON determines the SPS. While no significant difference in the initial ovule number was observed among the diverse oilseed rape genotypes (Li et al. 2015 ; Yang et al. 2016 ), large variations (from 5 to 35) in the SPS have been reported (Yang et al. 2016 ). Based on this, we hypothesized that there is a critical period from the pre-embryo to maturation that determines the variation in EON among contrasting NUtE genotypes. Biologically, the EON of oilseed rape genotypes is primarily determined by the number of flowers initiated, the number of ovules per placenta, the proportion of ovules to be fertilized (the proportion of fertile ovules × the proportion of fertile ovules to be fertilized), the proportion of fertilized ovules to develop into seeds as well as the frequency of embryonic abortions (the proportion of fertilized ovules to develop into seeds) (Jeong et al. 2012 ; Yang et al. 2016 ; Yang et al. 2017b ). Embryonic abortion, which occurs at different developmental stages within the fruit, has been shown to be non-random and based on the genotype of the embryo and environmental factors as well as their interaction (Calvino 2014 ), such as pollen grain germination (Wang et al. 2014 ), silique photosynthesis (Khan et al. 2018 ; Labra et al. 2017 ), nutrient accumulation (Li et al. 2020 ) and the number of flowers (Wang et al. 2011 ; Brookes et al. 2010 ). Although the developmental processes from the ovule fertilization to seed maturity is well known among oilseed rape genotypes, there is a paucity of knowledge that pinpointed the main factors that determine the EON, especially in pollen indices. Thus, establishment of the link between EON and pollen indices will provide basic information for breeders to take steps to reduce ovule abortion and increase yield of oilseed rape genotypes. It has been reported that the EON is strongly and positively correlated with NUtE, and nitrogen utilization-efficient genotypes had higher EON than the inefficient genotypes (He et al. 2017 ). Therefore, 18 highly pure rapeseed lines with contrasting NUtE and broad genetic background were used to systematically dissect the main causes and critical period of the natural variation of EON in rapeseed under both a field and a pot experiments from 2016–2017 and 2017–2018. In detail, the main objectives of the current study were to: i) unravel the main factors that restrict EON, and ii) determine the critical period of EON formation. It was hoped that our results would be helpful for the researchers that work in breeding for high yielding and NUtE oilseed genotypes for sustainable agriculture in the future. Materials And Methods Plant materials The results from a previous study, in which 50 different oilseed rape genotypes were examined under both pot and field experiments, allowed us to categorize the tested genotypes into four groups, based on their NUtE values: Nt-responder, Nt-nonresponder, Nt-efficient, and Nt-inefficient. Nt-responder referred to genotypes with an NUtE value above the mean at high N (0.3 g N kg − 1 dry soil or 180 kg N ha − 1 ), while genotypes with a NUtE value below the mean were categorized into the Nt-nonresponder group. At the low N supply (0.1 g N kg − 1 dry soil or 0 kg N ha − 1 ), genotypes displayed a NUtE value above the mean were called Nt-efficient, and Nt-inefficient genotypes were those genotypes with a NUtE value below the mean (He et al. 2017 ). According to the rank of responses reported previously, we chose 18 genotypes with distinctive NUtE for this study. These included 5 Nt-responder, 5 Nt-nonresponder, 4 Nt-efficient and 4 Nt-inefficient and are described in Table S1. Experimental design In this study, a field and a pot experiments were conducted at Yangling district, Shaanxi province, P.R. China (34⁰ 24’ N, 108⁰ 08’ E) from 2016–2018. In the field trail, the soil was with a pH of 7.65, containing 13.7 g kg − 1 organic matter, 1.19 g of total N kg − 1 , 24.7 mg of available N kg − 1 , 15.7 mg of Olsen-P kg − 1 and 76.9 mg of available K kg − 1 . In the pot experiment, the soil had the following properties: 9.21 g kg − 1 organic matter, 0.92 g of total N kg − 1 , 22.0 mg of available N kg − 1 , 10.7 mg of Olsen-P kg − 1 and 73.6 mg of available K kg − 1 , and a pH of 7.5. Our previous study indicated that the Nt-responder and Nt-nonresponder genotypes were expressed only under high N supply conditions, while the Nt-efficient and Nt-inefficient were expressed under low N supply conditions. Therefore, the 10 high N genotypes were planted with high N level (150 kg N ha − 1 in the field experiment and 0.30 g of N kg − 1 dry soil in the pot experiment), while the 8 low N genotypes were planted with low N level (0 kg N ha − 1 in field experiment and 0.10 g of N kg − 1 dry soil in pot experiment). In the field experiment, the high or low plot of 10 m in length × 2 m wide or 8 m length × 2 m width, consisted of forty or thirty-two rows of oilseed rape with row spacing of 50 cm. Each experiment was arranged in a randomized complete block design with four replicates. In both experiments, sufficient phosphorus as superphosphate (P 2 O 5 135 kg ha − 1 and 0.20 g kg − 1 dry soil) and potassium (K 2 O 150 kg ha − 1 and 0.30 g kg − 1 dry soil) were supplied. Observations and measurements In order to monitor the dynamics of silique development, flowering and fertilization process of the genotypes were categorized into: pre-embryo (1–4 days after flowering), globular (5–8 days after flowering), heart (9–14 days after flowering), torpedo (15–22 days after flowering), and maturation (23–30 days after flowering) stages (Andriotis et al. 2010 ; Tan et al. 2011 ; Hehenberger et al. 2012 ). For accurate monitoring, fully-opened flowers on the main inflorescence of each plant were marked with colored strings each day. In the pot experiment, five marked siliques of each plant at each developmental stage were chosen for the determination of net photosynthetic rate, silique length, width, surface area, EON, biomass and RNA expression levels. At the flowering stage, we removed the petals from the sampled buds with tweezers, and recorded the stamen and anther number. Then, the pollen grain number and pollen viability were measured according to Lankinen et al. ( 2018 ): Pollen grains were sprinkled on a microscope slide in Hoekstra medium with 16% sucrose. The microscope slides were placed in a dark constant temperature incubator with a temperature of 25˚C for 3 h. Pollen germination was terminated by adding 100% glycerol. Pollen grain number and germination rate (pollen viability) were determined under a light microscope (Axioplan 2; Zeiss) as the percentage of germinated pollen grains from 100 pollen grains in a randomly chosen area with the following equation: Pollen viability = the total number of germinated pollen grains / the total number of pollen grains × 100%. The initiation ovule number and abortion rate were determined according to Wang et al. ( 2011 ): Ovule abortion rate = the aborted ovule number / the initial ovule number. Measurement of agronomic traits At maturation stage, plant height was measured from the base of the stem to the tip of the main inflorescence (Sun et al. 2016 ). The point of measurement of stem diameter was set at 10 cm from the base of the main stem (Ohashi et al. 2006 ). First, valid branch height was measured as the height from the base of the stem to the effective primary branches at the bottom of the main stem. Then, the number of first valid branches was measured according to Xu et al. ( 2014a ). Measurement of yield and N use efficiency Number of siliques per plant was measured as the number of effective pods on the main inflorescence, branch inflorescence and the whole plant (Shi et al. 2015 ). Seed weight of each plant was measured by weighing 500 fully developed seeds with four replications; the weight of 500 seeds was then converted to 1,000-seed weight for easy comparison with other studies (Fan et al. 2010 ). Seed yield per plant was measured as the average dry weight of seeds of the four randomly selected plants from each genotype (Shi et al. 2011 ). Each of the ground seed and straw samples was weighed, and then digested and determined for N concentration by following the Kjeldahl procedure. Seed and shoot (seed + straw) N accumulations were calculated as the product of dry matter by the respective N concentration. According to He et al. ( 2017 ), N utilization efficiency (NUtE) and N uptake efficiency (NUpE) were estimated with the following equations: NUtE = seed yield / shoot N accumulation; NUpE = shoot N accumulation / N supply. Measurement of siliques indices At each stage (pre-embryo, globular, heart, torpedo, and maturation), silique net photosynthetic rate was measured on the marked siliques between 09:00 and 11:00, using a Portable Photosynthesis System (Li-6400, LI-COR, Lincoln, NE, USA). Afterwards, the marked siliques were removed from the main inflorescence. Half of each of the sampled siliques was rapidly frozen in liquid nitrogen, and stored at − 80°C for later RNA transcript analysis, and the other half was used for the determination of surface area (Hua et al. 2012 ) and normal ovule numbers which were counted under a light microscope. The sampled siliques were oven-dried at 105°C for 30 min and then at 70°C until a constant dry weight was reached and then were weighed. Quantitative RT-PCR analysis Silique samples of the 18 diverse NUtE genotypes were frozen in liquid nitrogen immediately after collection and stored at -80℃. Approximately 100 mg of tissue was ground in liquid nitrogen and total RNA was extracted using an EZNATM Plant RNA Kit (Omega Bio-Tek Inc., Norcross, GA, USA). The sample was used for cDNA synthesis using TransScript® First-Strand cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China) according to the manufacturer’s instructions. Quantitative PCR analyses were performed in an QuantStudio® Design and Analysis (QuantStudio 5, Life Technologies, CA, USA), using a TransStart Tip Green qPCR SuperMix (TransGen Biotech, China). All of the primers used in the analysis are listed in Table S2. The low NUtE genotypes (Nt-nonresponder and Nt-inefficient) were used as references (Nt-nonresponder for Nt-responder and Nt-inefficient for Nt-efficient, respectively). Data were expressed as the mean of four biological replicates ± SD. Statistical Analysis All the data from both the field and pot experiments, were subjected to analysis of variance. Estimates of correlation coefficients and tests of significance, were performed using the SPSS version 17.0. Principal component analysis was performed using the Canoco 5.0 software. When the analysis showed significance, mean comparisons were made according to the least-significant difference test at P ≤ 0.05 (LSD 0.05 ). All the figures were created using Origin 9.0 software, Microsoft Excel and PowerPoint 2016. Results Stamen and ovule characteristics Our results showed that the stamens of the contrasting NUtE genotypes were identical and normal in the field and pot experiments (Fig. 1A-D), and no significant differences in anther and filament number were found (Fig. 1E-H, Supplementary Fig. 1). On average, the high NUtE (Nt-responder and Nt-efficient) genotypes showed 44.1% (52.7% and 33.1%) greater pollen grain number and 23.5% (21.2% and 26.2%) higher pollen viability than the low NUtE (Nt-nonresponder and Nt-inefficient) genotypes (Fig. 1I-P). However, there were no significant differences in the initial ovule numbers among the high and low NUtE genotypes under both the field and pot experiments (Fig. 2A-H). Silique physical indices We also observed the silique net photosynthetic rate of Nt-efficient genotypes was greater than the Nt-inefficient genotypes from the heart to maturation stages (Fig. 3B), and silique surface area for the high NUtE genotypes was higher than for the low NUtE genotypes from the heart to torpedo stages. The same trend was found in the silique length and width (Supplementary Fig. 2). Finally, higher effective ovule number and silique biomass were recorded for the high NUtE than for the low NUtE genotypes from the heart to maturation stages (Fig. 3E-H). Silique RNA transcript levels In order to understand the observed differences in silique development, 18 of the contrasting NUtE genotypes were selected for measurement of the expression profile of silique development related-genes ( BnARF18 , BnaC9.SMG7b , BnCLV3 and BnDA1 ). Our results found that there was no significant difference in the expression level of the four genes among the high and low NUtE genotypes at the globular stage (Fig. 4). However, from the heart to torpedo stages, the RNA expression levels were higher by 15-40 fold in BnARF18 , 5-21 fold in BnaC9.SMG7b , 5-15 fold in BnCLV3 and 4-18 fold in BnDA1 for the high NUtE than for the low NUtE genotypes. Ovule abortion rate Differences in ovule abortion rate among the contrasting NUtE genotypes also existed in both field and pot experiments, with 39.3% (45.4% in Nt-responder and 34.3% in Nt-efficient) lower ovule abortion rate for the high NUtE than for the low NUtE genotypes (Fig. 5A-H). Agronomic traits, yield and NUE No significant difference of agronomic traits (plant height, stem diameter, first valid branch height and number of first valid branches) among the contrasting NUtE genotypes was recorded under field or pot experiments (Table S3). However, the number of siliques per plant, seeds per silique and yield per plant were respectively 12.4%, 21.1% and 29.4% greater for the high NUtE genotypes than for the low NUtE genotypes. In contrast, the low NUtE genotypes displayed 14.3% higher 1000-seed weight than the high NUtE genotypes in both experiments (Table 1).Compared to the low NUtE, the high NUtE genotypes showed 14.3% and 27.7% larger NUtE and NUpE, respectively (Fig. 6). Correlations between effective ovule number and the agronomic traits Principal component analysis showed that the first two principal components explained about 66%-83% of the total variation (Fig. 7A-B). The first principal component was mainly spanned by traits of EON, anther number, filament number, pollen grain number, pollen viability, initial ovule number, stem diameter, first valid branch height, 1000-seed weight, yield per plant, and NUtE on the positive side of the scale, while ovule abortion rate spanned on the negative side. The vectors indicated NUtE, the yield per plant, pollen grain number and pollen viability exhibited obviously positive correlations with EON, while there was a strong and negative correlation between ovule abortion rate and EON in both the field and the pot experiments. Using the path coefficient analysis, we illustrated the direct and indirect effects of the three traits (pollen grain number, pollen viability and ovule abortion rate) on EON and their interrelationships (Fig. 7C-D). The pollen grain number and pollen viability displayed a positive direct effect on EON, while the ovule abortion rate displayed a negative direct effect on EON. A positive indirect effect of pollen grain number (or pollen viability) was evident on EON via the pollen viability (pollen grain number), however a negative direct effect occurred between EON and ovule abortion rate. Discussion Our study showed that although no significant differences in anther number and initial ovule number between the contrasting NUtE genotypes were found, the high NUtE genotypes displayed the greater pollen grain number, higher pollen viability and lower ovule abortion rate, leading in larger EON, yield and NUtE, compared to the low NUtE genotypes. After flowering, the fertilized ovules went through proembryo, globular, heart, torpedo and maturation stages. The results of silique net photosynthetic rate, surface area, biomass, and RNA expression levels revealed that there was no significant difference of EON from the proembryo to globular stages, but the high NUtE oilseed rape genotypes displayed higher EON than the low NUtE genotypes at the heart stage, until the maturation stage. We suggested that the crop physiologists and plant breeders should take effective measures to increase the pollen number and vigor, and decrease ovule abortion rate before the heart stage, so as to boost EON, seed yield and NUE. A simplified model (Fig. 8 ) was drawn to illustrate the mechanistic understanding of the main reasons and critical period of EON formation among the contrasting NUtE oilseed rape genotypes. Greater pollen grain number, higher pollen viability and lower ovule abortion rate contribute to the larger EON in the high NUtE genotypes The main biological processes that determine SPS are basically known in oilseed rape (Yang et al. 2017b ), however the pollen indices and ovule abortion rate responsible for its natural variation among contrasting NUtE genotypes as well as their relative contributions are poorly known. The specific aspect of present study was to report that the high NUtE genotypes produced greater seed yields and higher NUtE (Table 1 and Fig. 6 ). These high NUtE genotypes displayed a larger EON, which was associated with greater pollen grain number, higher pollen viability, and lower ovule abortion rate, compared to the low NUtE genotypes. The process from the initial ovule to final seeds is affected by many factors, such as pollen sterility (Yang et al. 2016 ), floral bud quantity (Luo et al. 2018 ), inflorescence position (Wang et al. 2014 ), planting density (Khan et al. 2018 ), shading (Labra et al. 2017 ) and silique photosynthesis (Wang et al. 2016a ). In this study, we explored some factors that may have caused variation in EON among the contrasting NUtE oilseed rape genotypes. After ovule fertilization, the ovule abortion is frequently associated with seed yield (Mena-Ali and Rocha 2005 ). In agreement with this study, we observed a negative correlation between EON and the ovule abortion rate (Fig. 7 C-D), and the ovule abortion rate of the high NUtE genotypes was significantly lower than the low NUtE genotypes (Fig. 5 ), resulting in greater seed yield and NUtE (Table 1 and Fig. 6 ). This suggested that breeding towards high NUtE genotypes can be a relatively quick trait-based strategy to reduce ovule abortion and improve seed yield and NUtE. Isoda et al. ( 2010 ) demonstrated the lower supply of photosynthates can cause more pod abortion, thereby reducing seed yields. Therefore, further increase in potential yield will likely rely on improving photosynthesis (Zhu et al. 2018 ). According to our results (Fig. 3 A-B), the greater silique net photosynthetic rate of the high NUtE genotypes than that of the low NUtE genotypes leads to higher photoassimilates supply to the developing pods, and ultimately greater crop yield and NUtE (Table 1 and Fig. 6 ). Additionally, the majority of reproductive abortion occurs early in development (Coello and Martinez-Barajas 2016 ), and negative correlations between EON and the abortion rate were recorded, suggesting ovule abortion rate may have acted as the “upstream” of EON. Therefore, implementation of strategies targeted on reducing ovule abortion rate, such as sufficient supply of water (Gusmao et al. 2012 ), appropriate temperature (Zinn et al. 2010 ) and light (Wang et al. 2011 ) during the early ovule developmental period, may increase EON, yield and NUtE in oilseed rape genotypes. It is generally believed that there were no significant differences in initial ovule number among different oilseed rape genotypes (Yang et al. 2016 ; Li et al. 2015 ). These previous studies are fully in agreement with our present results (Fig. 2 ), indicating the initial ovule number could not contribute to EON variation among oilseed rape genotypes. While no significant difference in the number of anthers among the high and low NUtE genotypes was found (Fig. 1 A-H), the high NUtE genotypes exhibited 44.1% greater pollen grain number and 23.5% higher pollen viability than the low NUtE genotypes (Fig. 1 I-P). It suggested that the gametophytic competition is intensified and the chance of fertilization yielding vigorous zygotes increases when larger pollen grain number and viability are selected in breeding programs (Zhang et al. 2010 ). Chen et al. ( 2014 ) reported that the decreased pollen viability resulted in lower grain yield in rice lines, and lower seed-set in sorghum was due to decreased pollen grain number and viability. Similar result was observed in Pyrus pyrifolia , where the higher number of pollen grains improved pollen germination rate, seed viability, fruit yield and fruit quality (Zhang et al. 2010 ). The present study demonstrated significant and positive correlations of EON with both the pollen grain number and pollen viability (Fig. 7 C-D). It indicated that the greater pollen grain number and pollen viability in the high NUtE genotypes led to the higher EON, and the latter ultimately resulted in greater seed yield and NUtE (Table 1 and Fig. 6 ). In oilseed rape plants, ovules fail to develop into mature seeds because of fewer viable pollens available to fertilize all the ovules, subsequently resulting in aborted ovules (Fu et al. 2018 ; Dai et al. 2018 ). That means the low viable pollen grain number and low pollen viability could therefore, be the “upstream” cause of ovule abortion, and those two pollen indices ultimately regulate the variation in EON by controlling ovule abortion rate. Pollen grains are produced within the anther and released to the environment after maturation (Selinski and Scheibe 2014 ). During pollen development, nitrogen is essential to ensure efficient growth. Inadequate N supply significantly reduced pollen grain number and pollen viability in maize (Zheng et al. 2017 ; Soto et al. 2010 ). Nikolic et al. ( 2012 ) demonstrated a significant correlation between the plant N accumulation and grain yield in wheat. In our study, we found 19.2% greater shoot N accumulation at pre-flowering for the high NUtE than for the low NUtE genotypes (data not shown). Accordingly, we speculated that the pre-flowering shoot N accumulation in the high NUtE genotypes may have served as the N reserve that allowed the plant to produce larger pollen grain number with higher pollen viability, leading to less ovule abortion, greater seed yield and larger NUtE, compared to the low NUtE oilseed rape genotypes. Larger effective ovule number for the high NUtE than for the low NUtE genotypes at the heart stage Because of the complex conditions, especially germplasm resources, the EON formation of diverse oilseed rape genotypes showed great variation during the reproductive stage, the most important process in the life cycle of plants (Yang et al. 2016 ; Pereira et al. 2014 ). In the present study, we first tested the anther and filament number, as well as initial ovule number, and found there was no significant difference among contrasting NUtE genotypes (Fig. 2 ). Based on this, we surmised that the critical period for EON variation occurred after fertilization during one or more stages of the pre-embryo, globular, heart, torpedo or maturation (Andriotis et al. 2010 ; Tan et al. 2011 ; Hehenberger et al. 2012 ). Our results clearly identified the heart stage as the critical stage at which the difference in dynamic change of EON formation became significant among the contrasting NUtE genotypes and remained until the maturation (Fig. 3 – 4 ). Xu et al. ( 2014b ) and Yang et al. ( 2016 ) respectively claimed the globular and torpedo stage as the critical period for EON formation among diverse oilseed rape genotypes, which was inconsistent with our results. Our fine turning procedure with more frequent sampling in the current study allowed us to pinpoint the critical stage which may explain the inconsistency with the literature. The rapeseed silique is one of the main photosynthetic organs in oilseed rape, and it also serves as an important sink organ that stores carbohydrates for filling the developing seeds (Wang et al. 2016a ; Chen et al. 2018 ). In this study, the net photosynthetic rate and biomass of the silique were both greater for the high NUtE than for the low NUtE genotypes at the heart stage (Fig. 3 A-B and 3 G-H), leading to an increase in silique biomass accumulation and sustained photosynthesis which are the major physiological determinants of EON increases (Makino 2011 ). It was expected that plants with greater silique surface area would have higher silique photosynthesis and thus would increase seed yield (Kebede et al. 2014 ). Compared to the low NUtE oilseed rape genotypes, our results showed significantly larger silique surface area, silique length and width in high NUtE genotypes at heart stage, implying that the heart stage is the critical period of EON development among contrasting NUtE oilseed rape genotypes. This was consistent with the results of silique net photosynthetic ratio, biomass and EON. Fertilized ovules development is a complex and important process in the development of flowering plants, and it is indispensable for the next generation production (Gao et al. 2017 ; Wei et al. 2015 ). The silique is the main organ in which fertilized ovules are formed and filled, and this is closely correlated to the growth conditions of plants and is all complex quantitative trait controlled by polygenes (Wang et al. 2016b ). Recent literatures have reported the identification of genes that regulate fertilized ovules or silique development in oilseed rape plants. For example, BnARF18 gene simultaneously regulated silique development and seed weight (Liu et al. 2015 ). BnaC9.SMG7b gene was found to alter fruit development and seed formation with regard to seed morphology (Li et al. 2015 ). The mutation of BnCLV3 produced more multilocular siliques with a significantly higher number of seeds per silique and a higher seed weight (Yang et al. 2018 ), and BnDA1 controlled the seeds weight and organ size (Wang et al. 2017 ). Here, we selected the above four silique development genes to understand the relationship between the EON development and the expression level of genes obtained from the high and low NUtE genotypes. Our results clearly showed higher silique RNA transcript levels of the high NUtE than of the low NUtE genotypes from the heart to torpedo stages (Fig. 4 ). It supported the conjecture that green siliques of oilseed rape functioned as a photosynthesis center, and may contribute assimilates and nutrients into the silique biomass and developing seeds, generating larger seeds (Wang et al. 2016a ; Hua et al. 2012 ; Bennett et al. 2011 ). This may be one of the reasons that the silique indices of the high NUtE genotypes were higher than the low NUtE genotypes, and further research on selection the silique development related-genes may be one of the prerequisites to enhance EON, and ultimately higher yield and NUtE (Table 1 and Fig. 6 ). Taken together, the significant differences in silique development related-genes transcript levels among the high and low NUtE genotypes were found at the heart stage, which led to the larger silique surface area, greater biomass and higher EON, accompanied with higher silique net photosynthetic rate for the high NUtE genotypes, and difference in the dynamic silique indices maintained until maturation. With these results as the supporting evidence, it was concluded that the heart stage is the critical period of determining the EON among contrasting NUtE genotypes in silique development. Conclusion Compared to the low NUtE oilseed rape genotypes, the high NUtE genotypes produced a larger number of pollen grains with higher pollen viability before flowering, and displayed a lower ovule abortion rate after ovule fertilization. At the heart stage, the high NUtE genotypes formed a larger EON than the low NUtE genotypes, and the dynamic difference in EON between the genotypes was maintained until the maturation. Consequently, the high NUtE genotypes produced more SPS and greater yield with larger NUtE than the low NUtE genotypes. With this evidence, we suggest that breeding programs on oilseed rape improvement to achieve higher NUtE and seed yield should be targeted at increasing pollen grain number and pollen viability, and decreasing the ovule abortion rate before the heart stage. Declarations Acknowledgments This study was financially supported, in part by the grants from the National Key R&D Program of China (2018YFD0200907), the Special Fund for Agro-scientific Research in the Public Interest (201503124), and the Innovative Research Team Plan of the Agriculture Ministry. The senior author was sponsored by the MOE-AAFC PhD Training Program. This is a joint contribution of Northwest A & F University and Agriculture and Agri-Food Canada (AAFC). AAFC-ORDC contribution no. 20-xxx. Conflict of interest The authors declare no competing financial interests. Author contributions XG conducted the two experiments, and analyzed the data and wrote the manuscript; YX and YN participated in pot experiment; XW and BC provided the experimental materials; BM, MN and YG drafted the work and revised the manuscript critically. All authors read and approved the final manuscript. References Andriotis VM, Pike MJ, Kular B, Rawsthorne S, Smith AM (2010) Starch turnover in developing oilseed embryos. New Phytologist 187 (3):791-804. doi:10.1111/j.1469-8137.2010.03311.x Bennett EJ, Roberts JA, Wagstaff C (2011) The role of the pod in seed development: strategies for manipulating yield. New Phytologist 190 (4):838-853. doi:10.1111/j.1469-8137.2011.03714.x Brookes RH, Jesson LK, Burd M (2010) Reproductive investment within inflorescences of Stylidium armeria varies with the strength of early resource commitment. Annals of botany 105 (5):697-705. doi:10.1093/aob/mcq026 Calvino A (2014) Effects of ovule and seed abortion on brood size and fruit costs in the leguminous shrub Caesalpinia gilliesii (Wall. ex Hook.) D. Dietr. Acta Botanica Brasilica 28 (1):59-67 Chen L, Chen L, Zhang X, Liu T, Niu S, Wen J, Yi B, Ma C, Tu J, Fu T, Shen J (2018) Identification of miRNAs that regulate silique development in Brassica napus . Plant Science 269:106-117. doi:10.1016/j.plantsci.2018.01.010 Chen M, Chen JJ, Fang JY, Guo ZF, Lu SY (2014) Down-regulation of S-adenosylmethionine decarboxylase genes results in reduced plant length, pollen viability, and abiotic stress tolerance. Plant Cell, Tissue and Organ Culture 116:311-322. doi:10.1007/s11240-013-0405-0) Coello P, Martinez-Barajas E (2016) Changes in nutrient distribution are part of the mechanism that promotes seed development under severe nutrient restriction. Plant Physiology and Biochemistry 99:21-26. doi:10.1016/j.plaphy.2015.11.022 Dai S, Kai W, Liang B, Wang J, Jiang L, Du Y, Sun Y, Leng P (2018) The functional analysis of SlNCED1 in tomato pollen development. Cellular and Molecular Life Sciences 71 (18):3457-3472. doi:10.1007/s00018-018-2809-9 Fan CC, Cai GQ, Qin J, Li QY, Yang MG, Wu JZ, Fu TD, Liu KD, Zhou YM (2010) Mapping of quantitative trait loci and development of allele-specific markers for seed weight in Brassica napus . Theoretical and Applied Genetics 121 (7):1289-1301. doi:10.1007/s00122-010-1388-4) Fu S, Yin L, Xu M, Li Y, Wang M, Yang J, Fu T, Wang J, Shen J, Ali A, Zou Q, Yi B, Wen J, Tao L, Kang Z, Tang R (2018) Maternal doubled haploid production in interploidy hybridization between Brassica napus and Brassica allooctaploids. Planta 247 (1):113-125. doi:10.1007/s00425-017-2772-y Gao H, Zhang Y, Wang W, Zhao K, Liu C, Bai L, Li R, Guo Y (2017) Two membrane-anchored aspartic proteases contribute to pollen and ovule development. Plant Physiology 173 (1):219-239. doi:10.1104/pp.16.01719 Gardner KM, Latta RG (2007) Shared quantitative trait loci underlying the genetic correlation between continuous traits. Molecular Ecology 16 (20):4195-4209. doi:10.1111/j.1365-294X.2007.03499.x Gusmao M, Siddique KHM, Flower K, Nesbitt H, Veneklaas EJ (2012) Water deficit during the reproductive period of grass pea ( Lathyrus sativus L.) reduced grain yield but maintained seed size. Journal of Agronomy and Crop Science 198 (6):430-441. doi:10.1111/j.1439-037X.2012.00513.x He HY, Yang R, Li YJ, Ma A, Cao LQ, Wu XM, Chen BY, Tian H, Gao YJ (2017) Genotypic variation in nitrogen utilization efficiency of oilseed rape ( Brassica napus ) under contrasting N supply in pot and field experiments. Frontiers in plant science 8:1825. doi:10.3389/fpls.2017.01825 Hehenberger E, Kradolfer D, Kohler C (2012) Endosperm cellularization defines an important developmental transition for embryo development. Development 139 (11):2031-2039. doi:10.1242/dev.077057 Hua W, Li RJ, Zhan GM, Liu J, Li J, Wang XF, Liu GH, Wang HZ (2012) Maternal control of seed oil content in Brassica napus : the role of silique wall photosynthesis. The Plant Journal 69 (3):432-444. doi:10.1111/j.1365-313X.2011.04802.x Isoda A, Mao HX, Li ZY, Wang PW (2010) Growth of high-yielding soybeans and its relation to air temperature in Xinjiang, China. Plant Production Science 13 (2):209-217. doi:10.1626/pps.13.209 Jeong N, Suh SJ, Kim MH, Lee S, Moon JK, Kim HS, Jeong SC (2012) Ln is a key regulator of leaflet shape and number of seeds per pod in soybean. The Plant cell 24 (12):4807-4818. doi:10.1105/tpc.112.104968 Kebede B, Rahman H, Chevre AM (2014) Quantitative trait loci (QTL) mapping of silique length and petal colour in Brassica rapa . Plant Breeding 133 (5):609-614. doi:10.1111/pbr.12193 Khan S, Anwar S, Kuai J, Noman A, Shahid M, Din M, Ali A, Zhou GS (2018) Alteration in yield and oil quality traits of winter rapeseed by lodging at different planting density and nitrogen rates. Scientific Reports 8 (1):634. doi:10.1038/s41598-017-18734-8 Labra MH, Struik PC, Evers JB, Calderini DF (2017) Plasticity of seed weight compensates reductions in seed number of oilseed rape in response to shading at flowering. European Journal of Agronomy 84:113-124. doi:10.1016/j.eja.2016.12.011 Lankinen A, Lindstrom SAM, D'Hertefeldt T (2018) Variable pollen viability and effects of pollen load size on components of seed set in cultivars and feral populations of oilseed rape. PLoS One 13 (9):e0204407. doi:10.1371/journal.pone.0204407 Li S, Chen L, Zhang L, Li X, Liu Y, Wu Z, Dong F, Wan L, Liu K, Hong D, Yang G (2015) BnaC9.SMG7b functions as a positive regulator of the number of seeds per silique in Brassica napus by regulating the formation of functional female gametophytes. Plant Physiology 169 (4):2744-2760. doi:10.1104/pp.15.01040 Li SY, Zhu YY, Varshney RK, Zhan JP, Zheng XX, Shi JQ, Wang XF, Liu GH, Wang HZ (2020) A systematic dissection of the mechanisms underlying the natural variation of silique number in rapeseed germplasm. Plant biotechnology journal 18:568-580. doi:10.1111/pbi.13224 Liu J, Hua W, Hu ZY, Yang HL, Zhang L, Li RJ, Deng LB, Sun XC, Wang XF, Wang HZ (2015) Natural variation in ARF18 gene simultaneously affects seed weight and silique length in polyploid rapeseed. Proceedings of the National Academy of Sciences 112 (37):E5123-5132. doi:10.1073/pnas.1502160112 Luo T, Zhang J, Khan MN, Liu J, Xu Z, Hu L (2018) Temperature variation caused by sowing dates significantly affects floral initiation and floral bud differentiation processes in rapeseed ( Brassica napus L.). Plant Science 271:40-51. doi:10.1016/j.plantsci.2018.03.004 Ma BL, Biswas DK, Herath AW, Whalen JK, Ruan SQY, Caldwell C, Earl H, Vanasse A, Scott P, Smith DL (2015) Growth, yield, and yield components of canola as affected by nitrogen, sulfur, and boron application. Journal of Plant Nutrition and Soil Science 178 (4):658-670 Makino A (2011) Photosynthesis, grain yield, and nitrogen utilization in rice and wheat. Plant Physiology 155 (1):125-129. doi:10.1104/pp.110.165076 Mena-Ali JI, Rocha OJ (2005) Effect of ovule position within the pod on the probability of seed production in Bauhinia ungulata (Fabaceae). Annals of botany 95 (3):449-455. doi:10.1093/aob/mci044 Nikolic O, Zivanovic T, Jelic M, Djalovic I (2012) Interrelationships between grain nitrogen content and other indicators of nitrogen accumulation and utilization efficiency in wheat plants. Chilean Journal of Agricultural Research 72 (1):111-116 Ohashi Y, Nakayama N, Saneoka H, Fujita K (2006) Effects of drought stress on photosynthetic gas exchange, chlorophyll fluorescence and stem diameter of soybean plants. Biologia Plantarum 50 (1):138-141 Palanivelu R, Tsukamoto T (2012) Pathfinding in angiosperm reproduction: pollen tube guidance by pistils ensures successful double fertilization. Wiley Interdisciplinary Reviews-Developmental Biology 1 (1):96-113. doi:10.1002/wdev.6 Pereira MCT, Crane JH, Montas W, Nietsche S, Vendrame WA (2014) Effects of storage length and flowering stage of pollen influence its viability, fruit set and fruit quality in ‘Red’ and ‘Lessard Thai’ sugar apple ( Annona squamosa ) and ‘Gefner’ atemoya ( A. cherimola × A. squamosa ). Scientia Horticulturae 178:55-60. doi:10.1016/j.scienta.2014.08.004 Selinski J, Scheibe R (2014) Pollen tube growth: where does the energy come from? Plant Signaling & Behavior 9 (12):e977200. doi:10.4161/15592324.2014.977200 Shi J, Zhan J, Yang Y, Ye J, Huang S, Li R, Wang X, Liu G, Wang H (2015) Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed ( Brassica napus L.). Scientific reports 5:14481. doi:10.1038/srep14481 Shi JQ, Li RY, Zou J, Long Y, Meng JL (2011) A dynamic and complex network regulates the heterosis of yield-correlated traits in rapeseed ( Brassica napus L.). PLoS One 6 (7):e21645. doi:10.1371/journal.pone.0021645 Soto G, Fox R, Ayub N, Alleva K, Guaimas F, Erijman EJ, Mazzella A, Amodeo G, Muschietti J (2010) TIP5;1 is an aquaporin specifically targeted to pollen mitochondria and is probably involved in nitrogen remobilization in Arabidopsis thaliana . The Plant Journal 64 (6):1038-1047. doi:10.1111/j.1365-313X.2010.04395.x Sun CM, Wang BQ, Yan L, Hu KN, Liu S, Zhou YM, Guan CY, Zhang ZQ, Li JN, Zhang JF, Chen S, Wen J, Ma CZ, Tu JX, Shen JX, Fu TD, Yi B (2016) Genome-wide association study provides insight into the genetic control of plant height in rapeseed ( Brassica napus L.). Frontiers in plant science 7:1102. doi:10.3389/fpls.2016.01102 Tan H, Yang XH, Zhang FX, Zheng X, Qu CM, Mu JY, Fu FY, Li JN, Guan RZ, Zhang HS, Wang GD, Zuo JR (2011) Enhanced seed oil production in canola by conditional expression of Brassica napus LEAFY COTYLEDON1 and LEC1-LIKE in developing seeds. Plant Physiology 156 (3):1577-1588. doi:10.1104/pp.111.175000 Wang CL, Hai JB, Yang JL, Tian JH, Chen WJ, Chen T, Luo HB, Wang H (2016a) Influence of leaf and silique photosynthesis on seeds yield and seeds oil quality of oilseed rape ( Brassica napus L.). European Journal of Agronomy 74:112-118. doi:10.1016/j.eja.2015.12.008 Wang JL, Tang MQ, Chen S, Zheng XF, Mo HX, Li SJ, Wang Z, Zhu KM, Ding LN, Liu SY, Li YH, Tan XL (2017) Down-regulation of BnDA1, whose gene locus is associated with the seeds weight, improves the seeds weight and organ size in Brassica napus . Plant biotechnology journal 15 (8):1024-1033. doi:10.1111/pbi.12696 Wang XD, Chen L, Wang AN, Wang H, Tian JH, Zhao XP, Chao HB, Zhao YJ, Zhao WG, Xiang J, Gan JP, Li MT (2016b) Quantitative trait loci analysis and genome-wide comparison for silique related traits in Brassica napus . BMC plant biology 16:71. doi:10.1186/s12870-016-0759-7 Wang XJ, Mathieu A, Cournède PH, de Reffye P, Zhang BG (2014) Application of a probabilistic model for analysing the abortion of seeds and pods in winter oilseed rape ( Brassica napus ). Annals of Applied Biology 165 (3):414-428. doi:10.1111/aab.12152 Wang XJ, Mathieu A, Cournede PH, Allirand JM, Jullien A, Reffye PD, Zhang BG (2011) Variability and regulation of the number of ovules, seeds and pods according to assimilate availability in winter oilseed rape ( Brassica napus L.). Field Crops Research 122 (1):60-69. doi:10.1016/j.fcr.2011.02.008 Wei B, Zhang J, Pang C, Yu H, Guo D, Jiang H, Ding M, Chen Z, Tao Q, Gu H, Qu LJ, Qin G (2015) The molecular mechanism of sporocyteless/nozzle in controlling Arabidopsis ovule development. Cell Research 25 (1):121-134. doi:10.1038/cr.2014.145 Xu JS, Song X, Cheng Y, Zou XL, Zeng L, Qiao X, Lu GY, Fu GP, Qu Z, Zhang XK (2014a) Identification of QTLs for branch number in oilseed rape ( Brassica napus L.). Journal of Genetics and Genomics 41 (10):557-559. doi:10.1016/j.jgg.2014.06.006 Xu WJ, Fu Y, Dong HL, Chen ZF, Zhang QY, Mao SS, He YJ, Qian W (2014b) Morphologic and physiologic characterization of an ovule abortion mutant in Brassica napus . Scientia Agricultura Sinica 47 (10):2062-2068. doi:10.3864/j.issn.0578-1752.2014.10.020 Yang Y, Shen YS, Li SD, Ge XH, Li ZY (2017a) High density linkage map construction and QTL detection for three silique-related traits in Orychophragmus violaceus derived Brassica napus population. Frontiers in plant science 8:1512. doi:10.3389/fpls.2017.01512 Yang Y, Shi J, Wang X, Liu G, Wang H (2016) Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis. Scientific reports 6:24124. doi:10.1038/srep24124 Yang Y, Zhu KY, Li HL, Han SQ, Meng QW, Khan SU, Fan CC, Xie KB, Zhou YM (2018) Precise editing of CLAVATA genes in Brassica napus L. regulates multilocular silique development. Plant biotechnology journal 16 (7):1322-1335. doi:10.1111/pbi.12872 Yang YH, Wang Y, Zhan JP, Shi JQ, Wang XF, Liu GH, Wang HZ (2017b) Genetic and cytological analyses of the natural variation of seed number per pod in rapeseed ( Brassica napus L.). Frontiers in plant science 8:1890. doi:10.3389/fpls.2017.01890 Zhang C, Tateishi N, Tanabe K (2010) Pollen density on the stigma affects endogenous gibberellin metabolism, seed and fruit set, and fruit quality in Pyrus pyrifolia . Journal of experimental botany 61 (15):4291-4302. doi:10.1093/jxb/erq232 Zhang LW, Li SP, Yang GS, Chen L (2012) Identification and mapping of a major dominant quantitative trait locus controlling seeds per silique as a single Mendelian factor in Brassica napus L. Theoretical and Applied Genetics 125:695-705. doi:10.1007/s00122-012-1861-3) Zheng HY, Wu HM, Pan XY, Jin WW, Li XX (2017) Aberrant meiotic modulation partially contributes to the lower germination rate of pollen grains in maize ( Zea mays L.) under low nitrogen supply. Plant and Cell Physiology 58:342-353. doi:10.1093/pcp/pcw195 Zhu XY, Zhang L, Kuang C, Guo Y, Huang CQ, Deng LB, Sun XC, Zhan GM, Hu ZY, Wang HZ, Hua W (2018) Important photosynthetic contribution of silique wall to seed yield-related traits in Arabidopsis thaliana. Photosynthesis Research 137 (3):493-501. doi:10.1007/s11120-018-0532-x Zinn KE, Tunc-Ozdemir M, Harper JF (2010) Temperature stress and plant sexual reproduction: uncovering the weakest links. Journal of Experimental Botany 61 (7):1959-1968. doi:10.1093/jxb/erq053 Tables Table 1 Means of yield and yield components of contrasting NUtE oilseed rape genotypes in the field and pot-culture experiments. Siliques per plant Seeds per silique 1000-seed weight (g) Yield per plant (g) Field experiment Pot experiment Field experiment Pot experiment Field experiment Pot experiment Field experiment Pot experiment Nt-responder 310.1 a 242.9 a 17.8 a 22.1 a 5.02 b 4.00 a 30.7 a 26.3 a Nt-nonresponder 266.2 b 207.3 b 15.4 b 17.5 b 5.81 a 4.13 a 24.7 b 21.7 b Nt-efficiency 252.2 a 210.6 a 16.8 a 19.1 a 4.71 a 3.36 b 21.8 a 21.0 a Nt-inefficiency 221.2 b 209.3 a 14.9 b 14.8 b 5.05 a 4.56 a 14.1 b 16.6 b Means of each trait between Nt-responder and Nt-nonresponder or between Nt-efficient and Nt-inefficient genotypes in each experiment followed by different letters indicate significant differences according to the ANOVA-protected LSD 0.05 test. Supplementary Files SupplementaryData.docx Supplementary data Cite Share Download PDF Status: Published Journal Publication published 21 Dec, 2021 Read the published version in Euphytica → Version 1 posted Editorial decision: Minor revisions 14 Nov, 2021 Reviews received at journal 30 Oct, 2021 Editor assigned by journal 20 Aug, 2021 First submitted to journal 20 Aug, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-830636","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":47317647,"identity":"eee2b29e-4e72-4452-9b3a-147ffa6c9ac7","order_by":0,"name":"Xiao Guo","email":"","orcid":"https://orcid.org/0000-0003-1519-8075","institution":"Northwest Agriculture and Forestry University","correspondingAuthor":false,"prefix":"","firstName":"Xiao","middleName":"","lastName":"Guo","suffix":""},{"id":47317648,"identity":"bc6a0c2e-a4fb-4284-b248-173e95f75d61","order_by":1,"name":"Bao-Luo Ma","email":"","orcid":"","institution":"Ottawa Research and Development Centre","correspondingAuthor":false,"prefix":"","firstName":"Bao-Luo","middleName":"","lastName":"Ma","suffix":""},{"id":47317649,"identity":"ac406427-14a8-4c8f-bb94-81bc3d0bc184","order_by":2,"name":"Neil B. McLaughlin","email":"","orcid":"","institution":"Ottawa Research and Development Centre","correspondingAuthor":false,"prefix":"","firstName":"Neil","middleName":"B.","lastName":"McLaughlin","suffix":""},{"id":47317650,"identity":"5d82e8e7-a1b5-4500-abff-a8f5e139a7ee","order_by":3,"name":"Yajun gao","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA8ElEQVRIiWNgGAWjYHCCBDDJDyIqQMQBYrVINgCJM0RqgQCDA8RqMTh+4Jl0wS87OePbPYYfDrYxyPHdSGD8XIBPy5mENOmZfcnGZnfOGEsAtRhL3khglp6B1z1ALbw9zInbbuRukP7YxpC44UYCGzMPPi3nH4C01CdunpG7+QfQlnrCWm4AbeH5cThxg0TuNpDDEgwIaZG88SDZmrfhuLHEjfxvFgfOSRjOPPOwWRqfFr7zOYm3ef5Uy/HPSEu+caDMRp7vePLBz/i0KBzgSWBgbIPzJYCYsQGPBgYG+Qb2AwwMf/CqGQWjYBSMgpEOAJqhU+SVEU94AAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0003-4386-3562","institution":"Northwest Agriculture and Forestry University","correspondingAuthor":true,"prefix":"","firstName":"Yajun","middleName":"","lastName":"gao","suffix":""},{"id":47317651,"identity":"af68289c-2577-467e-b760-831b40482906","order_by":4,"name":"Xiaoming Wu","email":"","orcid":"","institution":"Oil Crops Research Institute Chinese Academy of Agricultural Sciences","correspondingAuthor":false,"prefix":"","firstName":"Xiaoming","middleName":"","lastName":"Wu","suffix":""},{"id":47317652,"identity":"fc031783-5a9e-41a2-80f7-c2d597514d86","order_by":5,"name":"Biyun Chen","email":"","orcid":"","institution":"Oil Crops Research Institute Chinese Academy of Agricultural Sciences","correspondingAuthor":false,"prefix":"","firstName":"Biyun","middleName":"","lastName":"Chen","suffix":""}],"badges":[],"createdAt":"2021-08-20 14:52:23","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-830636/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-830636/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10681-021-02954-1","type":"published","date":"2021-12-21T11:54:45+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":12765920,"identity":"dd348e6e-5f40-4682-8773-0567e0103a4c","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":693075,"visible":true,"origin":"","legend":"Stamen traits under field and pot experiments. (A-D) the stamen phenotype among the high and low NUtE genotypes; anther number (E-H), pollen number (I-L) and pollen viability (M-P) of the contrasting NUtE genotypes. Each bar represents a mean plus SE of four biological replicates, and different letters above the bars within a figure indicate significant differences according to the ANOVA-protected LSD0.05 test.","description":"","filename":"floatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/e67fad731d1c7cd500600239.png"},{"id":12765912,"identity":"56189b39-2632-4c65-b675-f2a0c7c3fa36","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":911010,"visible":true,"origin":"","legend":"Initial ovule number under field and pot experiments. (A-D) the initial ovule phenotype under light microscopy; (E-H) the initial ovule numbers of contrasting NUtE genotypes. Each bar represents a mean plus SE of four biological replicates, and different letters above the bars indicate significant differences according to the ANOVA-protected LSD0.05 test.","description":"","filename":"floatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/6f94dc56c6b0941f39b32c64.png"},{"id":12765915,"identity":"737d727b-ea10-406a-b1b2-f60801062f9f","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":1296978,"visible":true,"origin":"","legend":"Silique net photosynthetic rate (A-B), silique surface area (C-D), effective ovule number (E-F) and silique biomass (G-H) among contrasting NUtE genotypes from pre-embryo, globular, heart, torpedo and maturation stages. * indicates significant difference among the contrasting NUtE genotypes (P \u003c 0.05).","description":"","filename":"floatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/7f0a764d6ce541a699114074.png"},{"id":12765916,"identity":"51f063f0-9383-4472-b133-bf78a57130b5","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":660966,"visible":true,"origin":"","legend":"BnARF18 (A-B), BnLCR (C-D), BnaC9.SMG7b (E-F) and BnDA1 (G-H) expression levels at the globular, heart and torpedo stages. The low NUtE genotypes (Nt-nonresponder and Nt-inefficient) were used as references (Nt-nonresponder for Nt-responder and Nt-inefficient for Nt-efficient, respectively). Data are expressed as the means of four biological replicates, and * indicates significant differences between the contrasting NUtE genotypes (P \u003c 0.05).","description":"","filename":"floatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/2a2bc49c43df3dbc4b0189b3.png"},{"id":12765917,"identity":"7d663715-2728-4c0a-9fb4-44ef8a2ed3d0","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":497029,"visible":true,"origin":"","legend":"Ovule development phenotypes of contrasting NUtE genotypes under field and pot experiment. (A-D) the image of ovule abortion at maturation phase; (E-H) the mean of the ovule abortion rate of contrasting NUtE oilseed rape genotypes. Each bar represents the mean of four biological replicates, and different letters above the bars indicate significant differences according to the ANOVA-protected LSD0.05 test.","description":"","filename":"floatimage5.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/601b2410e391ccae19fb53a8.png"},{"id":12765913,"identity":"0d4fc93f-b15d-4e82-acdf-a51227c4bfa1","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":264158,"visible":true,"origin":"","legend":"NUtE (A-D) and NUpE (E-H) of the contrasting NUtE genotypes under field and pot experiments. Each bar represents the mean and standard error of the mean of four biological replicates, and different letters above the bar indicate significant differences according to the ANOVA-protected LSD0.05 test. ","description":"","filename":"floatimage6.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/6ec21af109ba7ddc7196e615.png"},{"id":12766057,"identity":"bb87dd0d-7448-400a-8201-245da1d1e6ab","added_by":"auto","created_at":"2021-08-25 23:32:28","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":581521,"visible":true,"origin":"","legend":"Principal component analysis in high N (Fig.7A) and low N (Fig.7B) and the path coefficient for EON, pollen grain number, pollen activity and ovule abortion rate in high N (Fig. 7C) and low N (Fig.7D) of different oilseed rape genotypes under field and pot experiments. Nt-respn, Nt-responder; Nt-nonrs, Nt-nonresponder; Nt-effic, Nt-efficient; Nt-ineff, Nt-inefficient; EON, effective ovule number; N A, the number of anther; N P, the number of pollen; P V, pollen viability; N I, the number of initial ovule; O R, ovule abortion rate; P H, plant height; S D, stem diameter; F H, first valid branch height; F N, number of first valid branches; S A, silique surface area; 1000 W, 1000-seed weight; Y P, yield per plant; NUtE, nitrogen utilization efficiency; NUpE, nitrogen uptake efficiency.","description":"","filename":"floatimage7.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/d34a4ebe1e44647712ba64c9.png"},{"id":12765918,"identity":"4c537459-231b-4bbe-a262-978a847815cc","added_by":"auto","created_at":"2021-08-25 23:29:28","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":745784,"visible":true,"origin":"","legend":"A simplified model for the main reasons and critical period for the difference in the effective ovule numbers among contrasting NUtE oilseed rape genotypes. Compared to the low efficient genotypes, the high NUtE genotypes produced more pollen grain number, with higher pollen viability and less ovule abortion rate. After flowering, the EON of contrasting NUtE oilseed rape genotypes went through five distinct stages (proembryo, globular, heart, torpedo and maturation). There was no significant difference in EON from the proembryo to globular stages, but the high NUtE oilseed rape genotypes displayed higher EON than the low NUtE genotypes at the heart stage, and the difference in EON maintained until the maturation stage. Therefore, with higher pollen number and pollen viability, and lower ovule abortion rate, the high NUtE genotypes produced higher EON, yield and NUtE than the low NUtE genotypes.","description":"","filename":"floatimage8.png","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/613f2664c38ea8a73aa8006c.png"},{"id":16649319,"identity":"b402d049-028a-491e-ac6f-6f2b1bbf0ae6","added_by":"auto","created_at":"2021-12-21 11:54:50","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4038422,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/32098866-e79d-4448-af2b-a9baf261710c.pdf"},{"id":12766056,"identity":"70f6abd8-38ea-4a40-983a-a1a2fbd46266","added_by":"auto","created_at":"2021-08-25 23:32:28","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":539397,"visible":true,"origin":"","legend":"Supplementary data","description":"","filename":"SupplementaryData.docx","url":"https://assets-eu.researchsquare.com/files/rs-830636/v1/af4a25b3040a6cb6369ccdfb.docx"}],"financialInterests":"","formattedTitle":"\u003cp\u003ePollen Grain Number and Viability Contribute to Variation in Effective Ovule Number Among Oilseed Rape Genotypes\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003eOilseed rape (\u003cem\u003eBrassica napus\u003c/em\u003e L.) is the world\u0026rsquo;s second largest crop source of vegetable oil following soybean, with increasing significance in the international market. Yield, one of the most important traits in oilseed rape breeding, is the product of seeds per silique (SPS), siliques per plant (SPP) and thousand seed weight (TSW)(Ma et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). As with other crops, yield components are critical to yield formation and often compensate each other, depending on genotypes and interaction with environment (Yang et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2017a\u003c/span\u003e; Shi et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Additionally, the three yield components (SPS, SPP and TSW) of oilseed rape display different degree of negative correlation, which indicates that the final seed yield can be improved by increasing the individual yield components (such as SPS) (Yang et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2017b\u003c/span\u003e). SPS, which is usually negatively correlated with the other two yield components (SPP and TSW) for competition among carbohydrate sinks, has received much attention (Zhang et al. \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Wang et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2016b\u003c/span\u003e; Yang et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2017a\u003c/span\u003e; Gardner and Latta \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). On the other hand, in the study conducted by He et al. (\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2017\u003c/span\u003e), the SPS was positively correlated with nitrogen utilization efficiency (NUtE; seed dry weight / shoot N accumulation), which indicated that SPS can be used as a simple and applicable indicator for the evaluation of NUtE in oilseed rape.\u003c/p\u003e \u003cp\u003eDuring the pollination process, pollen grains have to reach and germinate on the stigma, and then form pollen tubes to accurately reach the female gametes. After that, fertilization is initialized once the male gametes meet the female gametophyte within the ovule (Palanivelu and Tsukamoto \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Thus, successful reproduction relies on the pollen grains number, grain viability and female ovule number (Gao et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Li et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Following fertilization, the developmental process from embryogenesis to final seed can be divided into five distinct stages: pre-embryo, globular, heart, torpedo and maturation (Andriotis et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Tan et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Hehenberger et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). At the same time, the fertilized ovules may fail at any one of the stages due to the internal and external environmental conditions, and the remaining ovules (effective ovule number, EON) eventually develop into viable seeds (SPS). Therefore, to some degree, the EON determines the SPS. While no significant difference in the initial ovule number was observed among the diverse oilseed rape genotypes (Li et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2015\u003c/span\u003e; Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), large variations (from 5 to 35) in the SPS have been reported (Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Based on this, we hypothesized that there is a critical period from the pre-embryo to maturation that determines the variation in EON among contrasting NUtE genotypes.\u003c/p\u003e \u003cp\u003eBiologically, the EON of oilseed rape genotypes is primarily determined by the number of flowers initiated, the number of ovules per placenta, the proportion of ovules to be fertilized (the proportion of fertile ovules \u0026times; the proportion of fertile ovules to be fertilized), the proportion of fertilized ovules to develop into seeds as well as the frequency of embryonic abortions (the proportion of fertilized ovules to develop into seeds) (Jeong et al. \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Yang et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2017b\u003c/span\u003e). Embryonic abortion, which occurs at different developmental stages within the fruit, has been shown to be non-random and based on the genotype of the embryo and environmental factors as well as their interaction (Calvino \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), such as pollen grain germination (Wang et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), silique photosynthesis (Khan et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Labra et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2017\u003c/span\u003e), nutrient accumulation (Li et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) and the number of flowers (Wang et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Brookes et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Although the developmental processes from the ovule fertilization to seed maturity is well known among oilseed rape genotypes, there is a paucity of knowledge that pinpointed the main factors that determine the EON, especially in pollen indices. Thus, establishment of the link between EON and pollen indices will provide basic information for breeders to take steps to reduce ovule abortion and increase yield of oilseed rape genotypes.\u003c/p\u003e \u003cp\u003eIt has been reported that the EON is strongly and positively correlated with NUtE, and nitrogen utilization-efficient genotypes had higher EON than the inefficient genotypes (He et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Therefore, 18 highly pure rapeseed lines with contrasting NUtE and broad genetic background were used to systematically dissect the main causes and critical period of the natural variation of EON in rapeseed under both a field and a pot experiments from 2016\u0026ndash;2017 and 2017\u0026ndash;2018. In detail, the main objectives of the current study were to: i) unravel the main factors that restrict EON, and ii) determine the critical period of EON formation. It was hoped that our results would be helpful for the researchers that work in breeding for high yielding and NUtE oilseed genotypes for sustainable agriculture in the future.\u003c/p\u003e"},{"header":"Materials And Methods","content":"\u003cp\u003ePlant materials\u003c/p\u003e \u003cp\u003eThe results from a previous study, in which 50 different oilseed rape genotypes were examined under both pot and field experiments, allowed us to categorize the tested genotypes into four groups, based on their NUtE values: Nt-responder, Nt-nonresponder, Nt-efficient, and Nt-inefficient. Nt-responder referred to genotypes with an NUtE value above the mean at high N (0.3 g N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil or 180 kg N ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e), while genotypes with a NUtE value below the mean were categorized into the Nt-nonresponder group. At the low N supply (0.1 g N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil or 0 kg N ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e), genotypes displayed a NUtE value above the mean were called Nt-efficient, and Nt-inefficient genotypes were those genotypes with a NUtE value below the mean (He et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). According to the rank of responses reported previously, we chose 18 genotypes with distinctive NUtE for this study. These included 5 Nt-responder, 5 Nt-nonresponder, 4 Nt-efficient and 4 Nt-inefficient and are described in Table S1.\u003c/p\u003e \u003cp\u003eExperimental design\u003c/p\u003e \u003cp\u003eIn this study, a field and a pot experiments were conducted at Yangling district, Shaanxi province, P.R. China (34⁰ 24\u0026rsquo; N, 108⁰ 08\u0026rsquo; E) from 2016\u0026ndash;2018. In the field trail, the soil was with a pH of 7.65, containing 13.7 g kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e organic matter, 1.19 g of total N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, 24.7 mg of available N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, 15.7 mg of Olsen-P kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e and 76.9 mg of available K kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. In the pot experiment, the soil had the following properties: 9.21 g kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e organic matter, 0.92 g of total N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, 22.0 mg of available N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, 10.7 mg of Olsen-P kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e and 73.6 mg of available K kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e, and a pH of 7.5.\u003c/p\u003e \u003cp\u003eOur previous study indicated that the Nt-responder and Nt-nonresponder genotypes were expressed only under high N supply conditions, while the Nt-efficient and Nt-inefficient were expressed under low N supply conditions. Therefore, the 10 high N genotypes were planted with high N level (150 kg N ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in the field experiment and 0.30 g of N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil in the pot experiment), while the 8 low N genotypes were planted with low N level (0 kg N ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e in field experiment and 0.10 g of N kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil in pot experiment). In the field experiment, the high or low plot of 10 m in length \u0026times; 2 m wide or 8 m length \u0026times; 2 m width, consisted of forty or thirty-two rows of oilseed rape with row spacing of 50 cm. Each experiment was arranged in a randomized complete block design with four replicates. In both experiments, sufficient phosphorus as superphosphate (P\u003csub\u003e2\u003c/sub\u003eO\u003csub\u003e5\u003c/sub\u003e 135 kg ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e and 0.20 g kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil) and potassium (K\u003csub\u003e2\u003c/sub\u003eO 150 kg ha\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e and 0.30 g kg\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e dry soil) were supplied.\u003c/p\u003e \u003cp\u003eObservations and measurements\u003c/p\u003e \u003cp\u003eIn order to monitor the dynamics of silique development, flowering and fertilization process of the genotypes were categorized into: pre-embryo (1\u0026ndash;4 days after flowering), globular (5\u0026ndash;8 days after flowering), heart (9\u0026ndash;14 days after flowering), torpedo (15\u0026ndash;22 days after flowering), and maturation (23\u0026ndash;30 days after flowering) stages (Andriotis et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Tan et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Hehenberger et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). For accurate monitoring, fully-opened flowers on the main inflorescence of each plant were marked with colored strings each day. In the pot experiment, five marked siliques of each plant at each developmental stage were chosen for the determination of net photosynthetic rate, silique length, width, surface area, EON, biomass and RNA expression levels.\u003c/p\u003e \u003cp\u003eAt the flowering stage, we removed the petals from the sampled buds with tweezers, and recorded the stamen and anther number. Then, the pollen grain number and pollen viability were measured according to Lankinen et al. (\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2018\u003c/span\u003e): Pollen grains were sprinkled on a microscope slide in Hoekstra medium with 16% sucrose. The microscope slides were placed in a dark constant temperature incubator with a temperature of 25˚C for 3 h. Pollen germination was terminated by adding 100% glycerol. Pollen grain number and germination rate (pollen viability) were determined under a light microscope (Axioplan 2; Zeiss) as the percentage of germinated pollen grains from 100 pollen grains in a randomly chosen area with the following equation: Pollen viability\u0026thinsp;=\u0026thinsp;the total number of germinated pollen grains / the total number of pollen grains \u0026times; 100%. The initiation ovule number and abortion rate were determined according to Wang et al. (\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2011\u003c/span\u003e): Ovule abortion rate\u0026thinsp;=\u0026thinsp;the aborted ovule number / the initial ovule number.\u003c/p\u003e \u003cp\u003eMeasurement of agronomic traits\u003c/p\u003e \u003cp\u003eAt maturation stage, plant height was measured from the base of the stem to the tip of the main inflorescence (Sun et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The point of measurement of stem diameter was set at 10 cm from the base of the main stem (Ohashi et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). First, valid branch height was measured as the height from the base of the stem to the effective primary branches at the bottom of the main stem. Then, the number of first valid branches was measured according to Xu et al. (\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2014a\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eMeasurement of yield and N use efficiency\u003c/p\u003e \u003cp\u003eNumber of siliques per plant was measured as the number of effective pods on the main inflorescence, branch inflorescence and the whole plant (Shi et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Seed weight of each plant was measured by weighing 500 fully developed seeds with four replications; the weight of 500 seeds was then converted to 1,000-seed weight for easy comparison with other studies (Fan et al. \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Seed yield per plant was measured as the average dry weight of seeds of the four randomly selected plants from each genotype (Shi et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Each of the ground seed and straw samples was weighed, and then digested and determined for N concentration by following the Kjeldahl procedure. Seed and shoot (seed\u0026thinsp;+\u0026thinsp;straw) N accumulations were calculated as the product of dry matter by the respective N concentration. According to He et al. (\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2017\u003c/span\u003e), N utilization efficiency (NUtE) and N uptake efficiency (NUpE) were estimated with the following equations: NUtE\u0026thinsp;=\u0026thinsp;seed yield / shoot N accumulation; NUpE\u0026thinsp;=\u0026thinsp;shoot N accumulation / N supply.\u003c/p\u003e \u003cp\u003eMeasurement of siliques indices\u003c/p\u003e \u003cp\u003eAt each stage (pre-embryo, globular, heart, torpedo, and maturation), silique net photosynthetic rate was measured on the marked siliques between 09:00 and 11:00, using a Portable Photosynthesis System (Li-6400, LI-COR, Lincoln, NE, USA). Afterwards, the marked siliques were removed from the main inflorescence. Half of each of the sampled siliques was rapidly frozen in liquid nitrogen, and stored at \u0026minus;\u0026thinsp;80\u0026deg;C for later RNA transcript analysis, and the other half was used for the determination of surface area (Hua et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2012\u003c/span\u003e) and normal ovule numbers which were counted under a light microscope. The sampled siliques were oven-dried at 105\u0026deg;C for 30 min and then at 70\u0026deg;C until a constant dry weight was reached and then were weighed.\u003c/p\u003e \u003cp\u003eQuantitative RT-PCR analysis\u003c/p\u003e \u003cp\u003eSilique samples of the 18 diverse NUtE genotypes were frozen in liquid nitrogen immediately after collection and stored at -80℃. Approximately 100 mg of tissue was ground in liquid nitrogen and total RNA was extracted using an EZNATM Plant RNA Kit (Omega Bio-Tek Inc., Norcross, GA, USA). The sample was used for cDNA synthesis using TransScript\u0026reg; First-Strand cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China) according to the manufacturer\u0026rsquo;s instructions. Quantitative PCR analyses were performed in an QuantStudio\u0026reg; Design and Analysis (QuantStudio 5, Life Technologies, CA, USA), using a TransStart Tip Green qPCR SuperMix (TransGen Biotech, China). All of the primers used in the analysis are listed in Table S2. The low NUtE genotypes (Nt-nonresponder and Nt-inefficient) were used as references (Nt-nonresponder for Nt-responder and Nt-inefficient for Nt-efficient, respectively). Data were expressed as the mean of four biological replicates\u0026thinsp;\u0026plusmn;\u0026thinsp;SD.\u003c/p\u003e \u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eStatistical Analysis\u003c/h2\u003e \u003cp\u003eAll the data from both the field and pot experiments, were subjected to analysis of variance. Estimates of correlation coefficients and tests of significance, were performed using the SPSS version 17.0. Principal component analysis was performed using the Canoco 5.0 software. When the analysis showed significance, mean comparisons were made according to the least-significant difference test at \u003cem\u003eP\u003c/em\u003e\u0026thinsp;\u003cspan type=\"Underline\" class=\"Underline\" name=\"Emphasis\"\u003e\u0026le;\u003c/span\u003e\u0026thinsp;0.05 (LSD\u003csub\u003e0.05\u003c/sub\u003e). All the figures were created using Origin 9.0 software, Microsoft Excel and PowerPoint 2016.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cp\u003eStamen and ovule characteristics\u003c/p\u003e\n\u003cp\u003eOur results showed that the stamens of the contrasting NUtE genotypes were identical and normal in the\u0026nbsp;field and pot experiments\u0026nbsp;(Fig. 1A-D), and no significant differences in anther and filament number were found (Fig. 1E-H, Supplementary Fig. 1). On average, the high NUtE (Nt-responder and Nt-efficient)\u0026nbsp;genotypes showed 44.1% (52.7% and 33.1%) greater pollen grain number and 23.5% (21.2% and 26.2%) higher pollen viability than the low NUtE (Nt-nonresponder and Nt-inefficient) genotypes (Fig. 1I-P). However, there were no significant differences in the initial ovule numbers among the high and low NUtE genotypes under both the field and pot experiments (Fig. 2A-H).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSilique physical indices\u003c/p\u003e\n\u003cp\u003eWe also observed the silique net photosynthetic rate of Nt-efficient genotypes was greater than the Nt-inefficient genotypes from the heart to maturation stages (Fig. 3B), and silique surface area for the high NUtE genotypes was higher than for the low NUtE genotypes from the heart to torpedo stages. The same trend was found in the silique length and width (Supplementary Fig. 2). Finally, higher effective ovule number and silique biomass were recorded for the high NUtE than for the low NUtE genotypes from the heart to maturation stages (Fig. 3E-H).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSilique RNA transcript levels\u003c/p\u003e\n\u003cp\u003eIn order to understand the observed differences in silique development, 18 of the contrasting NUtE genotypes were selected for measurement of the expression profile of silique development related-genes (\u003cem\u003eBnARF18\u003c/em\u003e, \u003cem\u003eBnaC9.SMG7b\u003c/em\u003e, \u003cem\u003eBnCLV3\u003c/em\u003e and \u003cem\u003eBnDA1\u003c/em\u003e). Our results found that\u0026nbsp;there was no significant difference in the expression level of the four genes among the high and low NUtE genotypes at the globular stage (Fig. 4). However, from the heart to torpedo stages, the RNA expression levels were higher by 15-40 fold in \u003cem\u003eBnARF18\u003c/em\u003e, 5-21 fold in \u003cem\u003eBnaC9.SMG7b\u003c/em\u003e, 5-15 fold in \u003cem\u003eBnCLV3\u0026nbsp;\u003c/em\u003eand 4-18 fold in \u003cem\u003eBnDA1\u003c/em\u003e for the high NUtE than for the low NUtE genotypes.\u003c/p\u003e\n\u003cp\u003eOvule abortion rate\u003c/p\u003e\n\u003cp\u003eDifferences in ovule abortion rate among the contrasting NUtE genotypes also existed in both field and pot experiments, with 39.3% (45.4% in Nt-responder and 34.3% in Nt-efficient) lower ovule abortion rate for the high NUtE than for the low NUtE genotypes (Fig. 5A-H).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAgronomic traits, yield\u0026nbsp;and NUE\u003c/p\u003e\n\u003cp\u003eNo significant difference of agronomic traits (plant height, stem diameter, first valid branch height and number of first valid branches) among the contrasting NUtE genotypes was recorded under\u0026nbsp;field or pot experiments (Table S3). However, the number of siliques per plant, seeds per silique and yield per plant were respectively 12.4%, 21.1% and 29.4% greater for the high NUtE genotypes than for the low NUtE genotypes.\u0026nbsp;In contrast,\u0026nbsp;the low NUtE genotypes displayed 14.3% higher 1000-seed weight than the high NUtE genotypes in both experiments (Table 1).Compared to the low NUtE, the high NUtE genotypes showed 14.3% and 27.7% larger\u0026nbsp;NUtE and NUpE, respectively (Fig. 6).\u003c/p\u003e\n\u003cp\u003eCorrelations between effective ovule number and the agronomic traits\u003c/p\u003e\n\u003cp\u003ePrincipal component analysis showed that the first two principal components explained about 66%-83% of the total variation (Fig. 7A-B). The first principal component was mainly spanned by traits of EON, anther number, filament number, pollen grain number,\u0026nbsp;pollen viability, initial ovule number, stem diameter, first valid branch height, 1000-seed weight, yield per plant, and NUtE on the positive side of the scale, while ovule abortion rate spanned on the negative side. The vectors indicated NUtE, the yield per plant, pollen grain number and pollen viability exhibited obviously positive correlations with EON, while there was a strong and negative correlation between ovule abortion rate and EON in both the field and the pot experiments.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eUsing the path coefficient analysis, we illustrated the direct and indirect effects of the three traits (pollen grain number, pollen viability and ovule abortion rate) on EON and their interrelationships (Fig. 7C-D). The pollen grain number and pollen viability displayed a positive direct effect on EON, while the ovule abortion rate displayed a negative direct effect on EON. A positive indirect effect of pollen grain number (or pollen viability) was evident on EON via the pollen viability (pollen grain number), however a negative direct effect occurred between EON and ovule abortion rate.\u0026nbsp;\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eOur study showed that although no significant differences in anther number and initial ovule number between the contrasting NUtE genotypes were found, the high NUtE genotypes displayed the greater pollen grain number, higher pollen viability and lower ovule abortion rate, leading in larger EON, yield and NUtE, compared to the low NUtE genotypes. After flowering, the fertilized ovules went through proembryo, globular, heart, torpedo and maturation stages. The results of silique net photosynthetic rate, surface area, biomass, and RNA expression levels revealed that there was no significant difference of EON from the proembryo to globular stages, but the high NUtE oilseed rape genotypes displayed higher EON than the low NUtE genotypes at the heart stage, until the maturation stage. We suggested that the crop physiologists and plant breeders should take effective measures to increase the pollen number and vigor, and decrease ovule abortion rate before the heart stage, so as to boost EON, seed yield and NUE. A simplified model (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e) was drawn to illustrate the mechanistic understanding of the main reasons and critical period of EON formation among the contrasting NUtE oilseed rape genotypes.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eGreater pollen grain number, higher pollen viability and lower ovule abortion rate contribute to the larger EON in the high NUtE genotypes\u003c/p\u003e \u003cp\u003eThe main biological processes that determine SPS are basically known in oilseed rape (Yang et al. \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2017b\u003c/span\u003e), however the pollen indices and ovule abortion rate responsible for its natural variation among contrasting NUtE genotypes as well as their relative contributions are poorly known. The specific aspect of present study was to report that the high NUtE genotypes produced greater seed yields and higher NUtE (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). These high NUtE genotypes displayed a larger EON, which was associated with greater pollen grain number, higher pollen viability, and lower ovule abortion rate, compared to the low NUtE genotypes. The process from the initial ovule to final seeds is affected by many factors, such as pollen sterility (Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), floral bud quantity (Luo et al. \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), inflorescence position (Wang et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), planting density (Khan et al. \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), shading (Labra et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) and silique photosynthesis (Wang et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2016a\u003c/span\u003e). In this study, we explored some factors that may have caused variation in EON among the contrasting NUtE oilseed rape genotypes. After ovule fertilization, the ovule abortion is frequently associated with seed yield (Mena-Ali and Rocha \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2005\u003c/span\u003e). In agreement with this study, we observed a negative correlation between EON and the ovule abortion rate (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eC-D), and the ovule abortion rate of the high NUtE genotypes was significantly lower than the low NUtE genotypes (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e), resulting in greater seed yield and NUtE (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). This suggested that breeding towards high NUtE genotypes can be a relatively quick trait-based strategy to reduce ovule abortion and improve seed yield and NUtE. Isoda et al. (\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2010\u003c/span\u003e) demonstrated the lower supply of photosynthates can cause more pod abortion, thereby reducing seed yields. Therefore, further increase in potential yield will likely rely on improving photosynthesis (Zhu et al. \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). According to our results (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA-B), the greater silique net photosynthetic rate of the high NUtE genotypes than that of the low NUtE genotypes leads to higher photoassimilates supply to the developing pods, and ultimately greater crop yield and NUtE (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). Additionally, the majority of reproductive abortion occurs early in development (Coello and Martinez-Barajas \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2016\u003c/span\u003e), and negative correlations between EON and the abortion rate were recorded, suggesting ovule abortion rate may have acted as the \u0026ldquo;upstream\u0026rdquo; of EON. Therefore, implementation of strategies targeted on reducing ovule abortion rate, such as sufficient supply of water (Gusmao et al. \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2012\u003c/span\u003e), appropriate temperature (Zinn et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2010\u003c/span\u003e) and light (Wang et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2011\u003c/span\u003e) during the early ovule developmental period, may increase EON, yield and NUtE in oilseed rape genotypes.\u003c/p\u003e \u003cp\u003eIt is generally believed that there were no significant differences in initial ovule number among different oilseed rape genotypes (Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Li et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). These previous studies are fully in agreement with our present results (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e), indicating the initial ovule number could not contribute to EON variation among oilseed rape genotypes. While no significant difference in the number of anthers among the high and low NUtE genotypes was found (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA-H), the high NUtE genotypes exhibited 44.1% greater pollen grain number and 23.5% higher pollen viability than the low NUtE genotypes (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eI-P). It suggested that the gametophytic competition is intensified and the chance of fertilization yielding vigorous zygotes increases when larger pollen grain number and viability are selected in breeding programs (Zhang et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Chen et al. (\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) reported that the decreased pollen viability resulted in lower grain yield in rice lines, and lower seed-set in sorghum was due to decreased pollen grain number and viability. Similar result was observed in \u003cem\u003ePyrus pyrifolia\u003c/em\u003e, where the higher number of pollen grains improved pollen germination rate, seed viability, fruit yield and fruit quality (Zhang et al. \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). The present study demonstrated significant and positive correlations of EON with both the pollen grain number and pollen viability (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003eC-D). It indicated that the greater pollen grain number and pollen viability in the high NUtE genotypes led to the higher EON, and the latter ultimately resulted in greater seed yield and NUtE (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). In oilseed rape plants, ovules fail to develop into mature seeds because of fewer viable pollens available to fertilize all the ovules, subsequently resulting in aborted ovules (Fu et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Dai et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). That means the low viable pollen grain number and low pollen viability could therefore, be the \u0026ldquo;upstream\u0026rdquo; cause of ovule abortion, and those two pollen indices ultimately regulate the variation in EON by controlling ovule abortion rate. Pollen grains are produced within the anther and released to the environment after maturation (Selinski and Scheibe \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). During pollen development, nitrogen is essential to ensure efficient growth. Inadequate N supply significantly reduced pollen grain number and pollen viability in maize (Zheng et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Soto et al. \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e2010\u003c/span\u003e). Nikolic et al. (\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2012\u003c/span\u003e) demonstrated a significant correlation between the plant N accumulation and grain yield in wheat. In our study, we found 19.2% greater shoot N accumulation at pre-flowering for the high NUtE than for the low NUtE genotypes (data not shown). Accordingly, we speculated that the pre-flowering shoot N accumulation in the high NUtE genotypes may have served as the N reserve that allowed the plant to produce larger pollen grain number with higher pollen viability, leading to less ovule abortion, greater seed yield and larger NUtE, compared to the low NUtE oilseed rape genotypes.\u003c/p\u003e \u003cp\u003eLarger effective ovule number for the high NUtE than for the low NUtE genotypes at the heart stage\u003c/p\u003e \u003cp\u003eBecause of the complex conditions, especially germplasm resources, the EON formation of diverse oilseed rape genotypes showed great variation during the reproductive stage, the most important process in the life cycle of plants (Yang et al. \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e; Pereira et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). In the present study, we first tested the anther and filament number, as well as initial ovule number, and found there was no significant difference among contrasting NUtE genotypes (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). Based on this, we surmised that the critical period for EON variation occurred after fertilization during one or more stages of the pre-embryo, globular, heart, torpedo or maturation (Andriotis et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2010\u003c/span\u003e; Tan et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2011\u003c/span\u003e; Hehenberger et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Our results clearly identified the heart stage as the critical stage at which the difference in dynamic change of EON formation became significant among the contrasting NUtE genotypes and remained until the maturation (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e\u0026ndash;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Xu et al. (\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2014b\u003c/span\u003e) and Yang et al. (\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) respectively claimed the globular and torpedo stage as the critical period for EON formation among diverse oilseed rape genotypes, which was inconsistent with our results. Our fine turning procedure with more frequent sampling in the current study allowed us to pinpoint the critical stage which may explain the inconsistency with the literature. The rapeseed silique is one of the main photosynthetic organs in oilseed rape, and it also serves as an important sink organ that stores carbohydrates for filling the developing seeds (Wang et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2016a\u003c/span\u003e; Chen et al. \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). In this study, the net photosynthetic rate and biomass of the silique were both greater for the high NUtE than for the low NUtE genotypes at the heart stage (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA-B and \u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eG-H), leading to an increase in silique biomass accumulation and sustained photosynthesis which are the major physiological determinants of EON increases (Makino \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). It was expected that plants with greater silique surface area would have higher silique photosynthesis and thus would increase seed yield (Kebede et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). Compared to the low NUtE oilseed rape genotypes, our results showed significantly larger silique surface area, silique length and width in high NUtE genotypes at heart stage, implying that the heart stage is the critical period of EON development among contrasting NUtE oilseed rape genotypes. This was consistent with the results of silique net photosynthetic ratio, biomass and EON.\u003c/p\u003e \u003cp\u003eFertilized ovules development is a complex and important process in the development of flowering plants, and it is indispensable for the next generation production (Gao et al. \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Wei et al. \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). The silique is the main organ in which fertilized ovules are formed and filled, and this is closely correlated to the growth conditions of plants and is all complex quantitative trait controlled by polygenes (Wang et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2016b\u003c/span\u003e). Recent literatures have reported the identification of genes that regulate fertilized ovules or silique development in oilseed rape plants. For example, \u003cem\u003eBnARF18\u003c/em\u003e gene simultaneously regulated silique development and seed weight (Liu et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). \u003cem\u003eBnaC9.SMG7b\u003c/em\u003e gene was found to alter fruit development and seed formation with regard to seed morphology (Li et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). The mutation of \u003cem\u003eBnCLV3\u003c/em\u003e produced more multilocular siliques with a significantly higher number of seeds per silique and a higher seed weight (Yang et al. \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), and \u003cem\u003eBnDA1\u003c/em\u003e controlled the seeds weight and organ size (Wang et al. \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Here, we selected the above four silique development genes to understand the relationship between the EON development and the expression level of genes obtained from the high and low NUtE genotypes. Our results clearly showed higher silique RNA transcript levels of the high NUtE than of the low NUtE genotypes from the heart to torpedo stages (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). It supported the conjecture that green siliques of oilseed rape functioned as a photosynthesis center, and may contribute assimilates and nutrients into the silique biomass and developing seeds, generating larger seeds (Wang et al. \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2016a\u003c/span\u003e; Hua et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Bennett et al. \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). This may be one of the reasons that the silique indices of the high NUtE genotypes were higher than the low NUtE genotypes, and further research on selection the silique development related-genes may be one of the prerequisites to enhance EON, and ultimately higher yield and NUtE (Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e and Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). Taken together, the significant differences in silique development related-genes transcript levels among the high and low NUtE genotypes were found at the heart stage, which led to the larger silique surface area, greater biomass and higher EON, accompanied with higher silique net photosynthetic rate for the high NUtE genotypes, and difference in the dynamic silique indices maintained until maturation. With these results as the supporting evidence, it was concluded that the heart stage is the critical period of determining the EON among contrasting NUtE genotypes in silique development.\u003c/p\u003e"},{"header":"Conclusion","content":"\u003cp\u003eCompared to the low NUtE oilseed rape genotypes, the high NUtE genotypes produced a larger number of pollen grains with higher pollen viability before flowering, and displayed a lower ovule abortion rate after ovule fertilization. At the heart stage, the high NUtE genotypes formed a larger EON than the low NUtE genotypes, and the dynamic difference in EON between the genotypes was maintained until the maturation. Consequently, the high NUtE genotypes produced more SPS and greater yield with larger NUtE than the low NUtE genotypes. With this evidence, we suggest that breeding programs on oilseed rape improvement to achieve higher NUtE and seed yield should be targeted at increasing pollen grain number and pollen viability, and decreasing the ovule abortion rate before the heart stage.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eAcknowledgments\u0026nbsp;\u003c/strong\u003eThis study was financially supported, in part by the grants from the National Key R\u0026amp;D Program of China (2018YFD0200907), the Special Fund for Agro-scientific Research in the Public Interest (201503124), and the Innovative Research Team Plan of the Agriculture Ministry. The senior author was sponsored by the MOE-AAFC PhD Training Program. This is a joint contribution of Northwest A \u0026amp; F University and Agriculture and Agri-Food Canada (AAFC). AAFC-ORDC contribution no. 20-xxx.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConflict of interest\u0026nbsp;\u003c/strong\u003eThe authors declare no competing financial interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u0026nbsp;\u003c/strong\u003eXG conducted the two experiments, and analyzed the data and wrote the manuscript; YX and YN participated in pot experiment; XW and BC provided the experimental materials; BM, MN and YG drafted the work and revised the manuscript critically. All authors read and approved the final manuscript.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eAndriotis VM, Pike MJ, Kular B, Rawsthorne S, Smith AM (2010) Starch turnover in developing oilseed embryos. New Phytologist 187 (3):791-804. doi:10.1111/j.1469-8137.2010.03311.x\u003c/li\u003e\n\u003cli\u003eBennett EJ, Roberts JA, Wagstaff C (2011) The role of the pod in seed development: strategies for manipulating yield. New Phytologist 190 (4):838-853. doi:10.1111/j.1469-8137.2011.03714.x\u003c/li\u003e\n\u003cli\u003eBrookes RH, Jesson LK, Burd M (2010) Reproductive investment within inflorescences of \u003cem\u003eStylidium armeria\u003c/em\u003e varies with the strength of early resource commitment. Annals of botany 105 (5):697-705. doi:10.1093/aob/mcq026\u003c/li\u003e\n\u003cli\u003eCalvino A (2014) Effects of ovule and seed abortion on brood size and fruit costs in the leguminous shrub\u003cem\u003e\u0026nbsp;Caesalpinia gilliesii\u0026nbsp;\u003c/em\u003e(Wall. ex Hook.) D. Dietr. Acta Botanica Brasilica 28 (1):59-67\u003c/li\u003e\n\u003cli\u003eChen L, Chen L, Zhang X, Liu T, Niu S, Wen J, Yi B, Ma C, Tu J, Fu T, Shen J (2018) Identification of miRNAs that regulate silique development in \u003cem\u003eBrassica napus\u003c/em\u003e. Plant Science 269:106-117. doi:10.1016/j.plantsci.2018.01.010\u003c/li\u003e\n\u003cli\u003eChen M, Chen JJ, Fang JY, Guo ZF, Lu SY (2014) Down-regulation of S-adenosylmethionine decarboxylase genes results in reduced plant length, pollen viability, and abiotic stress tolerance. Plant Cell, Tissue and Organ Culture 116:311-322. doi:10.1007/s11240-013-0405-0)\u003c/li\u003e\n\u003cli\u003eCoello P, Martinez-Barajas E (2016) Changes in nutrient distribution are part of the mechanism that promotes seed development under severe nutrient restriction. Plant Physiology and Biochemistry 99:21-26. doi:10.1016/j.plaphy.2015.11.022\u003c/li\u003e\n\u003cli\u003eDai S, Kai W, Liang B, Wang J, Jiang L, Du Y, Sun Y, Leng P (2018) The functional analysis of \u003cem\u003eSlNCED1\u0026nbsp;\u003c/em\u003ein tomato pollen development. Cellular and Molecular Life Sciences 71 (18):3457-3472. doi:10.1007/s00018-018-2809-9\u003c/li\u003e\n\u003cli\u003eFan CC, Cai GQ, Qin J, Li QY, Yang MG, Wu JZ, Fu TD, Liu KD, Zhou YM (2010) Mapping of quantitative trait loci and development of allele-specific markers for seed weight in \u003cem\u003eBrassica napus\u003c/em\u003e. Theoretical and Applied Genetics 121 (7):1289-1301. doi:10.1007/s00122-010-1388-4)\u003c/li\u003e\n\u003cli\u003eFu S, Yin L, Xu M, Li Y, Wang M, Yang J, Fu T, Wang J, Shen J, Ali A, Zou Q, Yi B, Wen J, Tao L, Kang Z, Tang R (2018) Maternal doubled haploid production in interploidy hybridization between \u003cem\u003eBrassica napus\u003c/em\u003e and \u003cem\u003eBrassica\u003c/em\u003e allooctaploids. Planta 247 (1):113-125. doi:10.1007/s00425-017-2772-y\u003c/li\u003e\n\u003cli\u003eGao H, Zhang Y, Wang W, Zhao K, Liu C, Bai L, Li R, Guo Y (2017) Two membrane-anchored aspartic proteases contribute to pollen and ovule development. Plant Physiology 173 (1):219-239. doi:10.1104/pp.16.01719\u003c/li\u003e\n\u003cli\u003eGardner KM, Latta RG (2007) Shared quantitative trait loci underlying the genetic correlation between continuous traits. Molecular Ecology 16 (20):4195-4209. doi:10.1111/j.1365-294X.2007.03499.x\u003c/li\u003e\n\u003cli\u003eGusmao M, Siddique KHM, Flower K, Nesbitt H, Veneklaas EJ (2012) Water deficit during the reproductive period of grass pea (\u003cem\u003eLathyrus sativus\u0026nbsp;\u003c/em\u003eL.) reduced grain yield but maintained seed size. Journal of Agronomy and Crop Science 198 (6):430-441. doi:10.1111/j.1439-037X.2012.00513.x\u003c/li\u003e\n\u003cli\u003eHe HY, Yang R, Li YJ, Ma A, Cao LQ, Wu XM, Chen BY, Tian H, Gao YJ (2017) Genotypic variation in nitrogen utilization efficiency of oilseed rape (\u003cem\u003eBrassica napus\u003c/em\u003e) under contrasting N supply in pot and field experiments. Frontiers in plant science 8:1825. doi:10.3389/fpls.2017.01825\u003c/li\u003e\n\u003cli\u003eHehenberger E, Kradolfer D, Kohler C (2012) Endosperm cellularization defines an important developmental transition for embryo development. Development 139 (11):2031-2039. doi:10.1242/dev.077057\u003c/li\u003e\n\u003cli\u003eHua W, Li RJ, Zhan GM, Liu J, Li J, Wang XF, Liu GH, Wang HZ (2012) Maternal control of seed oil content in\u003cem\u003e\u0026nbsp;Brassica napus\u003c/em\u003e: the role of silique wall photosynthesis. The Plant Journal 69 (3):432-444. doi:10.1111/j.1365-313X.2011.04802.x\u003c/li\u003e\n\u003cli\u003eIsoda A, Mao HX, Li ZY, Wang PW (2010) Growth of high-yielding soybeans and its relation to air temperature in Xinjiang, China. Plant Production Science 13 (2):209-217. doi:10.1626/pps.13.209\u003c/li\u003e\n\u003cli\u003eJeong N, Suh SJ, Kim MH, Lee S, Moon JK, Kim HS, Jeong SC (2012) \u003cem\u003eLn\u003c/em\u003e is a key regulator of leaflet shape and number of seeds per pod in soybean. The Plant cell 24 (12):4807-4818. doi:10.1105/tpc.112.104968\u003c/li\u003e\n\u003cli\u003eKebede B, Rahman H, Chevre AM (2014) Quantitative trait loci (QTL) mapping of silique length and petal colour in \u003cem\u003eBrassica rapa\u003c/em\u003e. Plant Breeding 133 (5):609-614. doi:10.1111/pbr.12193\u003c/li\u003e\n\u003cli\u003eKhan S, Anwar S, Kuai J, Noman A, Shahid M, Din M, Ali A, Zhou GS (2018) Alteration in yield and oil quality traits of winter rapeseed by lodging at different planting density and nitrogen rates. Scientific Reports 8 (1):634. doi:10.1038/s41598-017-18734-8\u003c/li\u003e\n\u003cli\u003eLabra MH, Struik PC, Evers JB, Calderini DF (2017) Plasticity of seed weight compensates reductions in seed number of oilseed rape in response to shading at flowering. European Journal of Agronomy 84:113-124. doi:10.1016/j.eja.2016.12.011\u003c/li\u003e\n\u003cli\u003eLankinen A, Lindstrom SAM, D'Hertefeldt T (2018) Variable pollen viability and effects of pollen load size on components of seed set in cultivars and feral populations of oilseed rape. PLoS One 13 (9):e0204407. doi:10.1371/journal.pone.0204407\u003c/li\u003e\n\u003cli\u003eLi S, Chen L, Zhang L, Li X, Liu Y, Wu Z, Dong F, Wan L, Liu K, Hong D, Yang G (2015) \u003cem\u003eBnaC9.SMG7b\u003c/em\u003e functions as a positive regulator of the number of seeds per silique in \u003cem\u003eBrassica napus\u003c/em\u003e by regulating the formation of functional female gametophytes. Plant Physiology 169 (4):2744-2760. doi:10.1104/pp.15.01040\u003c/li\u003e\n\u003cli\u003eLi SY, Zhu YY, Varshney RK, Zhan JP, Zheng XX, Shi JQ, Wang XF, Liu GH, Wang HZ (2020) A systematic dissection of the mechanisms underlying the natural variation of silique number in rapeseed germplasm. Plant biotechnology journal 18:568-580. doi:10.1111/pbi.13224\u003c/li\u003e\n\u003cli\u003eLiu J, Hua W, Hu ZY, Yang HL, Zhang L, Li RJ, Deng LB, Sun XC, Wang XF, Wang HZ (2015) Natural variation in \u003cem\u003eARF18\u003c/em\u003e gene simultaneously affects seed weight and silique length in polyploid rapeseed. Proceedings of the National Academy of Sciences 112 (37):E5123-5132. doi:10.1073/pnas.1502160112\u003c/li\u003e\n\u003cli\u003eLuo T, Zhang J, Khan MN, Liu J, Xu Z, Hu L (2018) Temperature variation caused by sowing dates significantly affects floral initiation and floral bud differentiation processes in rapeseed (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Plant Science 271:40-51. doi:10.1016/j.plantsci.2018.03.004\u003c/li\u003e\n\u003cli\u003eMa BL, Biswas DK, Herath AW, Whalen JK, Ruan SQY, Caldwell C, Earl H, Vanasse A, Scott P, Smith DL (2015) Growth, yield, and yield components of canola as affected by nitrogen, sulfur, and boron application. Journal of Plant Nutrition and Soil Science 178 (4):658-670\u003c/li\u003e\n\u003cli\u003eMakino A (2011) Photosynthesis, grain yield, and nitrogen utilization in rice and wheat. Plant Physiology 155 (1):125-129. doi:10.1104/pp.110.165076\u003c/li\u003e\n\u003cli\u003eMena-Ali JI, Rocha OJ (2005) Effect of ovule position within the pod on the probability of seed production in \u003cem\u003eBauhinia ungulata\u003c/em\u003e (Fabaceae). Annals of botany 95 (3):449-455. doi:10.1093/aob/mci044\u003c/li\u003e\n\u003cli\u003eNikolic O, Zivanovic T, Jelic M, Djalovic I (2012) Interrelationships between grain nitrogen content and other indicators of nitrogen accumulation and utilization efficiency in wheat plants. Chilean Journal of Agricultural Research 72 (1):111-116\u003c/li\u003e\n\u003cli\u003eOhashi Y, Nakayama N, Saneoka H, Fujita K (2006) Effects of drought stress on photosynthetic gas exchange, chlorophyll fluorescence and stem diameter of soybean plants. Biologia Plantarum 50 (1):138-141\u003c/li\u003e\n\u003cli\u003ePalanivelu R, Tsukamoto T (2012) Pathfinding in angiosperm reproduction: pollen tube guidance by pistils ensures successful double fertilization. Wiley Interdisciplinary Reviews-Developmental Biology 1 (1):96-113. doi:10.1002/wdev.6\u003c/li\u003e\n\u003cli\u003ePereira MCT, Crane JH, Montas W, Nietsche S, Vendrame WA (2014) Effects of storage length and flowering stage of pollen influence its viability, fruit set and fruit quality in \u0026lsquo;Red\u0026rsquo; and \u0026lsquo;Lessard Thai\u0026rsquo; sugar apple (\u003cem\u003eAnnona squamosa\u003c/em\u003e) and \u0026lsquo;Gefner\u0026rsquo; atemoya (\u003cem\u003eA. cherimola\u0026nbsp;\u003c/em\u003e\u0026times; \u003cem\u003eA. squamosa\u003c/em\u003e). Scientia Horticulturae 178:55-60. doi:10.1016/j.scienta.2014.08.004\u003c/li\u003e\n\u003cli\u003eSelinski J, Scheibe R (2014) Pollen tube growth: where does the energy come from? Plant Signaling \u0026amp; Behavior 9 (12):e977200. doi:10.4161/15592324.2014.977200\u003c/li\u003e\n\u003cli\u003eShi J, Zhan J, Yang Y, Ye J, Huang S, Li R, Wang X, Liu G, Wang H (2015) Linkage and regional association analysis reveal two new tightly-linked major-QTLs for pod number and seed number per pod in rapeseed (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Scientific reports 5:14481. doi:10.1038/srep14481\u003c/li\u003e\n\u003cli\u003eShi JQ, Li RY, Zou J, Long Y, Meng JL (2011) A dynamic and complex network regulates the heterosis of yield-correlated traits in rapeseed (\u003cem\u003eBrassica napus\u0026nbsp;\u003c/em\u003eL.). PLoS One 6 (7):e21645. doi:10.1371/journal.pone.0021645\u003c/li\u003e\n\u003cli\u003eSoto G, Fox R, Ayub N, Alleva K, Guaimas F, Erijman EJ, Mazzella A, Amodeo G, Muschietti J (2010) TIP5;1 is an aquaporin specifically targeted to pollen mitochondria and is probably involved in nitrogen remobilization in \u003cem\u003eArabidopsis thaliana\u003c/em\u003e. The Plant Journal 64 (6):1038-1047. doi:10.1111/j.1365-313X.2010.04395.x\u003c/li\u003e\n\u003cli\u003eSun CM, Wang BQ, Yan L, Hu KN, Liu S, Zhou YM, Guan CY, Zhang ZQ, Li JN, Zhang JF, Chen S, Wen J, Ma CZ, Tu JX, Shen JX, Fu TD, Yi B (2016) Genome-wide association study provides insight into the genetic control of plant height in rapeseed (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Frontiers in plant science 7:1102. doi:10.3389/fpls.2016.01102\u003c/li\u003e\n\u003cli\u003eTan H, Yang XH, Zhang FX, Zheng X, Qu CM, Mu JY, Fu FY, Li JN, Guan RZ, Zhang HS, Wang GD, Zuo JR (2011) Enhanced seed oil production in canola by conditional expression of Brassica napus \u003cem\u003eLEAFY COTYLEDON1\u003c/em\u003e and \u003cem\u003eLEC1-LIKE\u003c/em\u003e in developing seeds. Plant Physiology 156 (3):1577-1588. doi:10.1104/pp.111.175000\u003c/li\u003e\n\u003cli\u003eWang CL, Hai JB, Yang JL, Tian JH, Chen WJ, Chen T, Luo HB, Wang H (2016a) Influence of leaf and silique photosynthesis on seeds yield and seeds oil quality of oilseed rape (\u003cem\u003eBrassica napus\u0026nbsp;\u003c/em\u003eL.). European Journal of Agronomy 74:112-118. doi:10.1016/j.eja.2015.12.008\u003c/li\u003e\n\u003cli\u003eWang JL, Tang MQ, Chen S, Zheng XF, Mo HX, Li SJ, Wang Z, Zhu KM, Ding LN, Liu SY, Li YH, Tan XL (2017) Down-regulation of BnDA1, whose gene locus is associated with the seeds weight, improves the seeds weight and organ size in \u003cem\u003eBrassica napus\u003c/em\u003e. Plant biotechnology journal 15 (8):1024-1033. doi:10.1111/pbi.12696\u003c/li\u003e\n\u003cli\u003eWang XD, Chen L, Wang AN, Wang H, Tian JH, Zhao XP, Chao HB, Zhao YJ, Zhao WG, Xiang J, Gan JP, Li MT (2016b) Quantitative trait loci analysis and genome-wide comparison for silique related traits in \u003cem\u003eBrassica napus\u003c/em\u003e. BMC plant biology 16:71. doi:10.1186/s12870-016-0759-7\u003c/li\u003e\n\u003cli\u003eWang XJ, Mathieu A, Courn\u0026egrave;de PH, de Reffye P, Zhang BG (2014) Application of a probabilistic model for analysing the abortion of seeds and pods in winter oilseed rape (\u003cem\u003eBrassica napus\u003c/em\u003e). Annals of Applied Biology 165 (3):414-428. doi:10.1111/aab.12152\u003c/li\u003e\n\u003cli\u003eWang XJ, Mathieu A, Cournede PH, Allirand JM, Jullien A, Reffye PD, Zhang BG (2011) Variability and regulation of the number of ovules, seeds and pods according to assimilate availability in winter oilseed rape (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Field Crops Research 122 (1):60-69. doi:10.1016/j.fcr.2011.02.008\u003c/li\u003e\n\u003cli\u003eWei B, Zhang J, Pang C, Yu H, Guo D, Jiang H, Ding M, Chen Z, Tao Q, Gu H, Qu LJ, Qin G (2015) The molecular mechanism of sporocyteless/nozzle in controlling \u003cem\u003eArabidopsis\u003c/em\u003e ovule development. Cell Research 25 (1):121-134. doi:10.1038/cr.2014.145\u003c/li\u003e\n\u003cli\u003eXu JS, Song X, Cheng Y, Zou XL, Zeng L, Qiao X, Lu GY, Fu GP, Qu Z, Zhang XK (2014a) Identification of QTLs for branch number in oilseed rape (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Journal of Genetics and Genomics 41 (10):557-559. doi:10.1016/j.jgg.2014.06.006\u003c/li\u003e\n\u003cli\u003eXu WJ, Fu Y, Dong HL, Chen ZF, Zhang QY, Mao SS, He YJ, Qian W (2014b) Morphologic and physiologic characterization of an ovule abortion mutant in \u003cem\u003eBrassica napus\u003c/em\u003e. Scientia Agricultura Sinica 47 (10):2062-2068. doi:10.3864/j.issn.0578-1752.2014.10.020\u003c/li\u003e\n\u003cli\u003eYang Y, Shen YS, Li SD, Ge XH, Li ZY (2017a) High density linkage map construction and QTL detection for three silique-related traits in \u003cem\u003eOrychophragmus violaceus\u003c/em\u003e derived \u003cem\u003eBrassica napus\u003c/em\u003e population. Frontiers in plant science 8:1512. doi:10.3389/fpls.2017.01512\u003c/li\u003e\n\u003cli\u003eYang Y, Shi J, Wang X, Liu G, Wang H (2016) Genetic architecture and mechanism of seed number per pod in rapeseed: elucidated through linkage and near-isogenic line analysis. Scientific reports 6:24124. doi:10.1038/srep24124\u003c/li\u003e\n\u003cli\u003eYang Y, Zhu KY, Li HL, Han SQ, Meng QW, Khan SU, Fan CC, Xie KB, Zhou YM (2018) Precise editing of \u003cem\u003eCLAVATA\u003c/em\u003e genes in \u003cem\u003eBrassica napus\u003c/em\u003e L. regulates multilocular silique development. Plant biotechnology journal 16 (7):1322-1335. doi:10.1111/pbi.12872\u003c/li\u003e\n\u003cli\u003eYang YH, Wang Y, Zhan JP, Shi JQ, Wang XF, Liu GH, Wang HZ (2017b) Genetic and cytological analyses of the natural variation of seed number per pod in rapeseed (\u003cem\u003eBrassica napus\u003c/em\u003e L.). Frontiers in plant science 8:1890. doi:10.3389/fpls.2017.01890\u003c/li\u003e\n\u003cli\u003eZhang C, Tateishi N, Tanabe K (2010) Pollen density on the stigma affects endogenous gibberellin metabolism, seed and fruit set, and fruit quality in \u003cem\u003ePyrus pyrifolia\u003c/em\u003e. Journal of experimental botany 61 (15):4291-4302. doi:10.1093/jxb/erq232\u003c/li\u003e\n\u003cli\u003eZhang LW, Li SP, Yang GS, Chen L (2012) Identification and mapping of a major dominant quantitative trait locus controlling seeds per silique as a single Mendelian factor in \u003cem\u003eBrassica napus\u0026nbsp;\u003c/em\u003eL. Theoretical and Applied Genetics 125:695-705. doi:10.1007/s00122-012-1861-3)\u003c/li\u003e\n\u003cli\u003eZheng HY, Wu HM, Pan XY, Jin WW, Li XX (2017) Aberrant meiotic modulation partially contributes to the lower germination rate of pollen grains in maize (\u003cem\u003eZea mays\u003c/em\u003e L.) under low nitrogen supply. Plant and Cell Physiology 58:342-353. doi:10.1093/pcp/pcw195\u003c/li\u003e\n\u003cli\u003eZhu XY, Zhang L, Kuang C, Guo Y, Huang CQ, Deng LB, Sun XC, Zhan GM, Hu ZY, Wang HZ, Hua W (2018) Important photosynthetic contribution of silique wall to seed yield-related traits in Arabidopsis thaliana. Photosynthesis Research 137 (3):493-501. doi:10.1007/s11120-018-0532-x\u003c/li\u003e\n\u003cli\u003eZinn KE, Tunc-Ozdemir M, Harper JF (2010) Temperature stress and plant sexual reproduction: uncovering the weakest links. Journal of Experimental Botany 61 (7):1959-1968. doi:10.1093/jxb/erq053\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003e\u003cstrong\u003eTable 1 \u0026nbsp;Means of yield and yield components of contrasting NUtE oilseed rape genotypes in the field and pot-culture experiments.\u003c/strong\u003e\u003c/p\u003e\n\u003cdiv align=\"center\"\u003e\n \u003ctable border=\"1\" cellpadding=\"0\" cellspacing=\"0\" width=\"0\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd rowspan=\"2\" valign=\"top\" width=\"14.78537360890302%\"\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" width=\"21.303656597774246%\"\u003e\n \u003cp\u003eSiliques per plant\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" width=\"21.303656597774246%\"\u003e\n \u003cp\u003eSeeds per silique\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" width=\"21.303656597774246%\"\u003e\n \u003cp\u003e1000-seed weight\u0026nbsp;\u003c/p\u003e\n \u003cp\u003e(g)\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd colspan=\"2\" width=\"21.303656597774246%\"\u003e\n \u003cp\u003eYield per plant\u003c/p\u003e\n \u003cp\u003e(g)\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003eField\u0026nbsp;experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003ePot experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003eField\u0026nbsp;experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003ePot experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003eField\u0026nbsp;experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003ePot experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003eField\u0026nbsp;experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"12.5%\"\u003e\n \u003cp\u003ePot experiment\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"14.78537360890302%\"\u003e\n \u003cp\u003eNt-responder\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e310.1 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e242.9 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e17.8 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e22.1 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e5.02 b\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e4.00 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e30.7 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e26.3 a\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"14.78537360890302%\"\u003e\n \u003cp\u003eNt-nonresponder\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e266.2 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e207.3 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e15.4 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e17.5 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e5.81 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e4.13 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e24.7 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e21.7 b\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"14.78537360890302%\"\u003e\n \u003cp\u003eNt-efficiency\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e252.2 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e210.6 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e16.8 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e19.1 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e4.71 a\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e3.36 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e21.8 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e21.0 a\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"14.78537360890302%\"\u003e\n \u003cp\u003eNt-inefficiency\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e221.2 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"10.651828298887123%\"\u003e\n \u003cp\u003e209.3 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e14.9 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e14.8 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e5.05 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e4.56 a\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"bottom\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e14.1 b\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd valign=\"top\" width=\"10.651828298887123%\"\u003e\n \u003cp\u003e16.6 b\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n \u003c/table\u003e\n\u003c/div\u003e\n\u003cp\u003eMeans of each trait between Nt-responder and Nt-nonresponder or between Nt-efficient and Nt-inefficient genotypes in each experiment followed by different letters indicate significant differences according to the ANOVA-protected LSD\u003csub\u003e0.05\u003c/sub\u003e test.\u0026nbsp;\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"euphytica","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"euph","sideBox":"Learn more about [Euphytica](https://www.springer.com/journal/10681)","snPcode":"10681","submissionUrl":"https://submission.springernature.com/new-submission/10681/3","title":"Euphytica","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"effective ovule number, pollen grain number, pollen viability, ovule abortion rate, heart stage","lastPublishedDoi":"10.21203/rs.3.rs-830636/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-830636/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eWhile no significant differences in initial ovule number were found among oilseed rape genotypes, there was a large variation in effective ovule number (EON), which determines the final seeds per silique (SPS), a critical component of yield. In this study, we selected 18 oilseed rape genotypes with contrasting nitrogen utilization efficiency (NUtE) to unravel the main factors responsible for different EON and determine the critical period of EON formation under both a field and a pot experiments from 2016-2018. The high NUtE genotypes displayed significantly higher NUtE by 14.3%, along with greater yield per plant (29.4%) and SPS (21.1%) than the low NUtE genotypes. The greater productivity of the high NUtE genotypes was associated with 44.1% higher pollen grain number, 23.5% greater pollen viability, and 39.3% lower ovule abortion rate, compared to the low NUtE genotypes. In addition, the heart stage was the critical ovule development period for delineating the variability of EON among contrasting NUtE oilseed rape genotypes, when the high NUtE genotypes displayed higher silique net photosynthetic rate, surface area, biomass, and RNA expression levels. Taken together, this study indicated the pollen grain number, pollen viability and ovule abortion rate contributed to the final variation in EON and the heart stage was the critical period of determining the EON differences among contrasting NUtE genotypes. Increasing pollen grain number and pollen viability, and decreasing ovule abortion rate before heart stage should be the prerequisite for breeders to improve yield and NUtE of oilseed rape genotypes.\u003c/p\u003e","manuscriptTitle":"Pollen Grain Number and Viability Contribute to Variation in Effective Ovule Number Among Oilseed Rape Genotypes","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-08-25 23:29:26","doi":"10.21203/rs.3.rs-830636/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Minor revisions","date":"2021-11-14T11:41:31+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2021-10-30T10:36:34+00:00","index":0,"fulltext":""},{"type":"editorAssigned","content":"","date":"2021-08-21T03:56:34+00:00","index":"","fulltext":""},{"type":"submitted","content":"Euphytica","date":"2021-08-20T04:44:16+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"euphytica","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"euph","sideBox":"Learn more about [Euphytica](https://www.springer.com/journal/10681)","snPcode":"10681","submissionUrl":"https://submission.springernature.com/new-submission/10681/3","title":"Euphytica","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"stoa","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"41502882-4d82-458e-89db-61c05c002de7","owner":[],"postedDate":"August 25th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":6655877,"name":"Agronomy"},{"id":6655878,"name":"Molecular Biology"},{"id":6655879,"name":"General Biochemistry"}],"tags":[],"updatedAt":"2021-12-21T11:54:45+00:00","versionOfRecord":{"articleIdentity":"rs-830636","link":"https://doi.org/10.1007/s10681-021-02954-1","journal":{"identity":"euphytica","isVorOnly":false,"title":"Euphytica"},"publishedOn":"2021-12-21 11:54:45","publishedOnDateReadable":"December 21st, 2021"},"versionCreatedAt":"2021-08-25 23:29:26","video":"","vorDoi":"10.1007/s10681-021-02954-1","vorDoiUrl":"https://doi.org/10.1007/s10681-021-02954-1","workflowStages":[]},"version":"v1","identity":"rs-830636","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-830636","identity":"rs-830636","version":["v1"]},"buildId":"_2-kVJe1T_tPrBINL-cwx","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}