Novel Insights into the Influence of Seed Sarcotesta Photosynthesis on Accumulation of Seed Dry Matter and Oil Content in Torreya grandis cv. “Merrillii”

Seed oil content is an important trait of nut seeds, and it is affected by the import of carbon from photosynthetic sources. Although green leaves are the main photosynthetic organs, seed sarcotesta photosynthesis also supplies assimilates to seed development. Understanding the relationship between seed photosynthesis and seed development has theoretical and practical significance in the cultivation of Torreya grandis cv. “Merrillii.” To assess the role of seed sarcotesta photosynthesis on the seed development, anatomical and physiological traits of sarcotesta were measured during two growing seasons in the field. Compared with the attached current-year leaves, the sarcotesta had higher gross photosynthetic rate at the first stage of seed development. At the late second stage of seed development, sarcotesta showed down-regulation of PSII activity, as indicated by significant decrease in the following chlorophyll fluorescence parameters: the maximum PSII efficiency (Fv/Fm), the PSII quantum yield (ΦPSII), and the photosynthetic quenching coefficient (qP). The ribulose 1, 5—bisphosphate carboxylase (Rubisco) activity, the total chlorophyll content (Chl(a+b)) and nitrogen content in the sarcotesta were also significantly decreased during that period. Treatment with DCMU [3-(3,4-dichlorophenyl)-1,1-dimethylurea] preventing seed photosynthesis decreased the seed dry weight and the oil content by 25.4 and 25.5%, respectively. We conclude that seed photosynthesis plays an important role in the dry matter accumulation at the first growth stage. Our results also suggest that down-regulation of seed photosynthesis is a plant response to re-balance the source-sink ratio at the second growth stage. These results suggest that seed photosynthesis is important for biomass accumulation and oil synthesis of the Torreya seeds. The results will facilitate achieving higher yields and oil contents in nut trees by selection for higher seed photosynthesis cultivars.

Seed oil content is an important trait of nut seeds, and it is affected by the import of carbon from photosynthetic sources. Although green leaves are the main photosynthetic organs, seed sarcotesta photosynthesis also supplies assimilates to seed development. Understanding the relationship between seed photosynthesis and seed development has theoretical and practical significance in the cultivation of Torreya grandis cv. "Merrillii." To assess the role of seed sarcotesta photosynthesis on the seed development, anatomical and physiological traits of sarcotesta were measured during two growing seasons in the field. Compared with the attached current-year leaves, the sarcotesta had higher gross photosynthetic rate at the first stage of seed development. At the late second stage of seed development, sarcotesta showed down-regulation of PSII activity, as indicated by significant decrease in the following chlorophyll fluorescence parameters: the maximum PSII efficiency (F v /F m ), the PSII quantum yield (Φ PSII ), and the photosynthetic quenching coefficient (qP). The ribulose 1, 5-bisphosphate carboxylase (Rubisco) activity, the total chlorophyll content (Chl (a+b) ) and nitrogen content in the sarcotesta were also significantly decreased during that period. Treatment with DCMU [3-(3, 4-dichlorophenyl)-1,1-dimethylurea] preventing seed photosynthesis decreased the seed dry weight and the oil content by 25.4 and 25.5%, respectively. We conclude that seed photosynthesis plays an important role in the dry matter accumulation at the first growth stage. Our results also suggest that down-regulation of seed photosynthesis is a plant response to re-balance the source-sink ratio at the second growth stage. These results suggest that seed photosynthesis is important for biomass accumulation and oil synthesis of the Torreya seeds. The results will facilitate achieving higher yields and oil contents in nut trees by selection for higher seed photosynthesis cultivars.

INTRODUCTION
Fruit and seed growth and the formation of storage reserves in them require much carbon substrate and energy through respiration. The fruits and seeds and the surrounding pericarp/seed coat have developed mechanisms to refix some of the respired carbon (Blanke and Lenz, 1989). Besides the carbon originating from leaf photosynthesis, internal recycling of respiratory CO 2 by fruits and seeds provides an important additional contribution to meet the carbon requirements of the seeds and fruits (Aschan and Pfanz, 2003).
It has been frequently shown that carbon re-fixation in reproductive organs makes in many plant species a significant photosynthetic contribution to meeting the carbon requirement of their own growth (Quebedeaux and Chollet, 1975;Whiley et al., 1992;King et al., 1998;Proietti et al., 1999;Furbank et al., 2004;Imai and Ogawa, 2009;Lytovchenko et al., 2011;Hu et al., 2012;Hua et al., 2012;Xu et al., 2016). Seed oil synthesis is dependent upon the supply of photosynthate from photosynthetic organs of the seeds (Baud and Lepiniec, 2010). It has been also reported that the photosynthesis taking place in green embryos, or in the external green non-foliar organs (e.g., pericarp, pod, silique, peel), plays an important role in the regulation of crop oil content (Singal et al., 1995;Asokanthan et al., 1997;Ruuska et al., 2004;Schwender et al., 2004;Goffman et al., 2005;Allen et al., 2009;Hua et al., 2012;Wang et al., 2016). Therefore, enormous research efforts have been devoted to studying the traits and mechanisms of photosynthesis in reproductive organs when discussing their growth and development.
Fruits maintain a relatively high gross photosynthetic rate at the early stage of development; the rate then progressively decreases until fruit maturity in evergreen species (Proietti et al., 1999;Hieke et al., 2002;Imai and Ogawa, 2009). It is known that similarly as the chloroplasts in leaves, the chloroplasts in the outer layer of fruits and seeds contain well-organized grana thylakoids and starch grains (Platt-Aloia and Thompson, 1981;Bondada and Oosterhuis, 2003;Li et al., 2006). The number of the granular thylakoids decreases and the size of the plastoglobules increases during fruit development (Prebeg et al., 2008). Some enzymes and pigments involved in photosynthesis have been found in the external green non-foliar organs of the reproductive organs, for example, in the silique wall of rapeseed (Hua et al., 2012), the pod of soybean (Quebedeaux and Chollet, 1975), in the ear and awns of wheat (Martinez et al., 2003;Li et al., 2004), and in the peel of tomato and apple (Xu et al., 1997;Chen and Cheng, 2007). The photosynthetic rate of fruits and seeds has been estimated by CO 2 gas exchange measurements (Wullschleger and Oosterhuis, 1990;Hua et al., 2012;Xu et al., 2016) and by oxygen electrode measurements (Chen and Cheng, 2007;Hu et al., 2012;Hiratsuka et al., 2015). The contribution of fruits and seeds photosynthesis to the yield has been also extensively studied by estimating the carbon budget of CO 2 gas exchange or O 2 evolution measurements (Dejong and Walton, 1989;Hieke et al., 2002;Imai and Ogawa, 2009;Hu et al., 2012), by shading (Maydup et al., 2010;Hu et al., 2012) and by using photosynthetic inhibitors, such as 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU) or paraquat (Maydup et al., 2010;Hua et al., 2012).
Tree nuts are rich in monounsaturated (MUFA) and polyunsaturated fatty acids (PUFA), which are recommended as an important resource of healthy diet in human populations throughout the world (O'Neil et al., 2012). Torreya grandis cv. "Merrillii" (later also "Torreya, " or the "Torreya cultivar"), a tree species with significant economic value, has drupe-like fruits that provide an example of seedoil nuts, with an oil content of 55%; about which 80% are unsaturated fatty acids Wang and Xiu (2005). Yu et al. (1986) reported that active chloroplasts are observed in the fresh outer layer (sarcotesta) of seeds of Torreya. The sarcotesta remain green for a long time during seed development. However, little information is available on the underlying physiological and biochemical mechanisms of photosynthesis in Torreya seeds and the relative contribution of seed photosynthesis to seed dry matter accumulation and oil synthesis during the seed development.
The aim of the present study is to (1) examine the potential photosynthetic mechanisms underlying the variable photosynthetic rate at the different seed development stages, and (2) to evaluate the role of seed photosynthesis in dry matter accumulation and oil synthesis in the seeds. Our understanding of the photosynthetic physiology and of the role of seed photosynthesis at the different stages of seed development will help to provide a theoretical basis for achieving higher nut yields and higher oil content by effective utilization of potential photosynthesis of seed sarcotesta.

Plant Materials and Growth Conditions
Nine Torreya grandis cv. "Merrillii" grafting trees planted at Yuqian Town, Lin'an City, Zhejiang Province, China (30 • 10 ′ N, 119 • 22 ′ E) were used in this study. The trees were established by using 1 year old Torreya grandis cv. "Merrillii" scion on 2 year old Torreya grandis Fort. Ex Lindl. root stock in 2001, and they began bearing seed since 2009. The trees were spaced in a 5 m × 5 m orientation and they received 62.4 kg N ha −1 , 48 kg P ha −1 , and 43.2 kg K ha −1 in March; and 31.2 kg N ha −1 , 24 kg P ha −1 , and kg 21.6 kg K ha −1 after seed harvest. Throughout the study, the trees were maintained with standard fertilization, irrigation, and pest control practices recommended for the Torreya cultivar (Dai et al., 2008;Zhang et al., 2017).
The field measurements were carried out in two consecutive growing seasons during years 2015-2016. The development of the Torreya flowers and seeds lasted approximately 17 months, spanning 2 years. Accordingly, the seeds measured in 2015 and 2016 were based on blooming in 2014 and 2015, respectively. The blooming season lasted from middle April to early May both in 2014 and 2015. The seeds used in the study were produced by hand pollination on 12th April 2014 and 15th April 2015. The fertilization process is usually finished in early September (Liu et al., 2017). In present study, the viable seeds began to stick out from the seed scale on 15th April 2015 and 20th April 2016. Each of the seeds to be measured was labeled with a plastic tag on the day when the seed scale was breached and that day was documented as the day of protrusion for the seed (Figure 1). Seeds of approximately the same size were selected for the measurements.

Design of the Field Study
With the exception of a separate experiment (see below), four Torreya trees were sampled for the measurements. From each sampled tree, seeds and leaves for the measurements were sampled randomly from the outer part of the crown at a height of approximately 1.0-1.5 m. The number of trees and seeds/leaves sampled are indicated below for each measurement category. Sampling took place on eight dates between 18th May and 9th September in 2015 and on six dates between 12th May and 11th September in 2016. However, some of the measurements were carried out only on the sampling dates in 2015 (see below), and sampling for the microscopic examination of the chloroplast ultrastructure was carried out only during one date in 2015. For each sampling date, the time of the measurement was determined as the number of days after the seeds emerged from the seed scale (days after seeds protrusion of, DASP).

Measurements of the Dimensions and Weight of the Seeds
For measurements of the dimensions and weight, three seeds from each of the four sampled trees were sampled on each sampling date in 2015 and 2016. Length (L), width (W) and thickness (T) of the seeds were measured with electronic digital calipers. Making the simplifying assumption that the seed shape is spherical, the surface area of the seeds was calculated by Equation (1) (Baryeh, 2001): The seeds were weighed by an electronic balance. The seed fresh weight was determined immediately after sampling. After that, the seeds were dried at 80 • C for more than 24 h until a constant weight (=dry weight) was attained. In order to examine FIGURE 1 | A microscopic photograph of a Torreya grandis cv. "Merrillii" seed before (A) and after (B) the protrusion of the seed from the seed scale.
qualitatively the color of the seeds, they were also photographed at each sampling occasion in 2015.

Measurement of Chlorophyll Fluorescence in Seeds
For measurements of chlorophyll fluorescence, three seeds from each of the four sampled trees were sampled in 2015 on four sampling dates at 27, 91, 106, 138, and 149 DASP. Chlorophyll fluorescence was measured in detached seeds with a pulsemodulated fluorometer PAM 2500 (Walz, Effeltrich, Germany). Seeds were wrapped in a wet paper towel after collecting them from the trees. We made the measurements in the middle part of the seeds where the surface is flat. Care was also taken to keep the distance between the fiber probe and the seed surface constant for all measurements. Prior to measurement, seeds were darkadapted for 30 min. The minimal fluorescence (F o ) was recorded after dark adaptation and the maximum fluorescence (F m ) after a strong saturating pulse with the intensity and the width of 10,000 µmol m −2 s −1 and 600 ms, respectively. Subsequently, each seed was illuminated with an actinic light of 1,161 µmol m −2 s −1 for 5 min, which allowed the fluorescent signal to reach a steady-state level (F s ) before a saturating pulse was imposed to obtain the maximal fluorescence under light exposure (F m '). Shortly after the saturating pulse, the actinic light was turned off and a far-red light was switched on for 3 s to rapidly oxidize the primary quinone acceptor in PSII by drawing electrons from PSII to PSI, so that the minimal fluorescence (F o ') under light exposure could be determined. The maximum PSII efficiency in dark-adapted seeds was calculated as F v /F m = (F m -F o )/F m (Kitajima and Butler, 1975). The PSII quantum yield (Φ PSII ) in the light was calculated as (F m ' -F s )/F m ' (Genty et al., 1989), and the photochemical quenching coefficient (qP) as (F m ' -F s )/(F m ' -F o ') (Krause and Weis, 1991). Photochemical quantum yield of open PSII centers in the light-adapted state (F v '/F m ') was calculated as (F m ' -F o ')/F m ' (Schreiber et al., 1994).

Measurement of Leaf and Seed CO 2 Exchange
For measurements of CO 2 exchange, three seeds and their associated leaves from each of the four sampled trees were sampled on each sampling date. CO 2 exchange was measured in the leaves and seeds using a portable photosynthesis measuring instrument (LI6400, LI-COR, Lincoln, NE, USA) according to Hu et al. (2013). Both for leaves and seeds, a cluster leaf chamber (6400-22, LI-COR) with an 18-RGB light source (LI-COR, Lincoln) was used.
In order to examine closely the seed photosynthetic parameters, seed CO 2 exchange was measured at several PAR levels on all sampling dates in 2015. Subsequently, equations were fitted to the data by estimating the values of their coefficients representing the basic photosynthetic parameters.
Light response curves of seeds were determined as follows: seeds enclosed in the chamber were first kept at 1,200 µmol photon m −2 s −1 for at least 20 min; thereafter, PAR was decreased in a stepwise manner to measure the net photosynthetic rate (P n ) at ten levels of photosynthetically active radiation (PAR): 0, 50, 100, 150, 200, 400, 600, 800, 1,000, and 1,200 µmol photons m −2 s −1 . At each PAR level P n was recorded when it reached a stable level, usually this took approximately 5 min. To avoid midday depression of photosynthesis in natural high-light (Larcher, 2003), the measurements were carried out in the field between 0800-1,100 and 1,400-1,600 h. The block temperature was maintained at 27-30 • C for all measurements.
To calculate the P (g,max) on a surface area, a dry weight or an individual seed basis, we used Equation (2) for calculating the net photosynthetic rate (P n ) (Charles -Edwards, 1981): Coefficient α is the reciprocal of PAR at which 50% of the asymptotic rate of P n is attained; coefficient b is the gradient of the saturation curve of Equation (2) at the origin, that is, the maximum quantum yield. The dark respiration rate (r) was estimated as the CO 2 exchange rate in darkness with PAR = 0 µmol photon m −2 s −1 . The gross photosynthetic rate (P g ) of seeds was calculated by summing the net photosynthetic rate (P n ) and the respiration rate (r). Finally, the maximum P g at light saturation (P (g,max) ) was calculated as b/α (Imai and Ogawa, 2009). For leaves, we used the standard method of expression P n and r per one-sided leaf area, so that the leaf area was taken as one-half of the total two-sided surface area of the leaf. Therefore, for seeds P g was divided by half of the surface area of the seed when P g was expressed as per surface area. The P (g,max) and rvalues were calculated both on surface area (µmol CO 2 m −2 s −1 ), seed dry weight (µmol CO 2 g −1 DW s −1 ), and individual seed (µmol CO 2 seed −1 ) basis. The apparent quantum use efficiency (AQE) was calculated by a linear regression on P g on the PAR in the range PAR <200 µmol photon m −2 s −1 . The light saturation point (LSP) was determined using an AQ response curve analysis software (Version 1.0, LI-COR, 2/2008, Ranjan et al., 2012). In order to compare seed and leaf gross photosynthesis at the same PAR, leaf gas exchange was measured at PAR = 800 µmol photons m −2 s −1 in 2015 and at PAR = 1,300 µmol photon m −2 s −1 in 2016. The net CO 2 assimilation rates of leaves were recorded after equilibration to obtain the steady-state photosynthetic rate at the applied PAR level. In order to calculate the gross photosynthetic rate, the dark respiration rate of the leaves was also measured. For the seeds, we obtained the value of the net photosynthetic rate from their light response curves measured in 2015 at the same PAR level as used in the measurements with leaves (PAR = 800 µmol photons m −2 s −1 and PAR = 1,300 µmol photons m −2 s −1 in 2015 and 2016, respectively).

Measurements of Chlorophyll Content, Rubisco Activity and Nitrogen Content of Seed Sarcotesta
For measurements of chlorophyll content, Rubisco (ribulose 1, 5 -bisphosphate carboxylase) activity, and nitrogen content, three seeds from each of four sampled trees were sampled on each sampling date (Rubisco and nitrogen content only in 2015). The sarcotesta (each about 25 mm 2 ) were extracted with 5 ml 95% ethyl alcohol for 24 h in darkness. Absorbance of the supernatant was measured with a spectrophotometer (UV-2550, Shimadzu, Kyoto, Japan) at 649, 664, and 470 nm after centrifugation. The contents of total chlorophylls (Chl a+b ), chlorophyll a (Chl a ) and chlorophyll b (Chl b ) were then calculated according to Lichtenthaler (1987).
Sarcotesta samples were dried and ground, then the total nitrogen content was determined with an azotometer (Kjeltec-2300 FOSS, Sweden) according to the micro-Kjeldahl method (Schuman et al., 1972).

Microscopy of the Chloroplast Ultrastructure in Leaves and Sarcotesta
Tissue samples (sarcotesta and leaf slices) were taken from seeds and the associated leaves at 81 DASP (5 July 2015) for the analysis of chloroplast ultrastructure. The samples were immediately fixed in 2.5% (v/v) glutaraldehyde (0.1 M phosphate buffer, pH 7.0) for at least 4 h once cut from the plants. The samples were then immersed in 1% (v/v) osmium tetroxide for postfixation. The specimens were dehydrated using a graded series of ethanol and embedded in epoxy resin for ultrathin sectioning and transmission electron microscopy (H7650, HITACHI, Tokyo, Japan).

Determining the Effects of DCMU Application on Seed Dry Weight and Seed Quality
In order to further examine the role of seed photosynthesis in fulfilling the carbon requirements of the seeds, an experiment with DCMU [3-(3,4-dichlorophenyl)-1,1-dimethylurea] was carried out in 2015. DCMU is a chemical which is commonly used to inhibit Photosystem II in photosynthesis by binding to the D1 protein in the reaction center (Allen and Holmes, 1986;Chow et al., 1990). The seeds were treated with DCMU application according to Maydup et al. (2010) with some modifications. For this experiment, the labeled seeds (approximately 80 seeds from each tree) were divided into two parts. One part was treated with DCMU application, and the other part was used as a control so that it was treated otherwise similarly but without DCMU application. In detail, starting at 52 DASP (6th June), the experimental seeds were treated with 100 µmol m −3 of DCMU mixed with 0.1% (v/v) Tween-20 as a surfactant. Control seeds were treated with Tween-20 alone at the same concentration. The intact seed were immersed with DCMU solution in a fixed plastic valve bag for 30 min. During the immersing, the leaves were not allowed to be in contact with any DCMU solution. The DCMU application needed to decrease the photosynthetic electron transport rate (ETR) to approximately 30% was assessed by a pulse-modulated fluorometer PAM 2500 (Walz, Effeltrich, Germany) 1 day after application. Subsequently, the DCMU was applied at midday on sunny days (any rain may dilute the solution) once a week from 52 DASP until the seeds had reached maturity. The concentration of 100 µmol m −3 was chosen according to our preliminary experiment. Fan et al. (2007) also used this concentration in the solution with which spinach leaf discs were vacuum infiltrated with the aim of inhibiting Photosystem II and disrupting the redox poise needed for cyclic electron transport. Their study also shows that the DCMU application does not have any other effects beyond the inhibition of photosynthesis on plants. See supplementary Material as Table S1 for further details.
On 11th September 2015 the seeds were harvested and the length, the width, and the fresh weight (FW) of the intact seeds were measured immediately. Then the sarcotesta were removed from the seeds, and the fresh weight of the remaining hard seed was measured. After drying to a constant weight the dry weight (DW) of both the sarcotesta and the hard seed were measured. The sum of these two was taken as the dry weight of the intact seeds.
The relative contribution of the sarcotesta photosynthesis to the dry weight of the seeds was calculated as the dry weight difference between the control seeds and the seeds treated with DCMU, as a percentage of the dry weight of the control seeds. This was carried out separately both to the intact seeds and to the seeds with sarcotesta removed.
Finally, in order to measure the oil content, the kernel of both control and DCMU treated seeds was ground with liquid nitrogen. About 8 g (per extraction thimble made of thick filter paper) of full fat flour was defatted in a Soxhlet apparatus with petroleum ether solvent (boiling point range 38.7-54.8 • C) for 8 h. The defatted samples were dried overnight (10-12 h) in a fume hood to remove residual petroleum ether and then weighed to calculate the oil content. The fatty acids composition of the oil was analyzed by gas chromatography (Agilent 7890A, Agilent Technologies, USA) according to the Fatty acid methyl esters (FAME) method of Maxwell and Marmer (1983). The standard gas chromatogram of main fatty acids is shown in Figure S1.

Statistical Analyses
Data were subjected to analysis of variance (ANOVA) using SPSS statistical software (16.0, IBM, New York, USA). The data are presented as the mean ± standard deviation (SD). The differences in the average values among different sampling occasions were tested by the least significant difference (LSD). Differences at P ≤ 0.05 were considered significant. The relationships between Rubisco activity, nitrogen content and P (g,max) were determined using correlation analysis.

Seed Color, Dimensions, and Weight
The color of the seeds at 27 DASP was deep green, turned to yellowish green at the middle stage of seed development (about at 69 DASP), and then become again much greener with further seed development (Figure 2A). Similar patterns of the seed length, width, and area were found both in 2015 (Figures 2B,C) and 2016 (Figures 2D,E): They all showed first a somewhat linear increase at the early stage, but later the increase leveled off before the maximum values were attained.
The seed fresh and dry weights exhibited a significant increasing trend during the seed development. The fresh weight of the seeds increased rapidly from 27 to 45 DASP in 2015 with a mean of 0.3 g day −1 (Figure 3A) and from 28 to 46 DASP in 2016 with a mean of 0.2 g day −1 (Figure 3D), followed in both years by a slower increase before reaching a constant weight. The seed dry weight significantly increased by 471 and 126% from 27 to 69 DASP and 70 to 149 DASP, respectively ( Figure 3B). It also showed a more linear increasing trend from 28 to 140 DASP in 2016 ( Figure 3E). The mean seed dry/fresh weight ratio was in 2015 the lowest at 45 DASP, followed first a plateau for approximately 24 days, then by an increase until 138 DASP, and remained finally another plateau thereafter (Figure 3C). Similarly, the mean seed dry/fresh weight ratio was the lowest at 46 DASP in 2016 ( Figure 3F).

CO 2 Exchange in the Seeds
The photosynthetic light response curves of seeds during their developmental stages are shown in Supplementary Data Figure S2. From 27 to 91 DASP the seeds had a high light saturation point (LSP) of 1,110∼1,183 µmol m −2 s −1 . Then, the LSP significantly decreased with seed age greater than 91 DASP ( Figure 4A). The apparent quantum use efficiency (AQE) of seeds showed an increasing trend from 27 to 69 DASP, followed by a decreasing trend as seed development progressed ( Figure 4B). The dark respiration R d and the maximum rate of gross photosynthesis P (g,max) per seed surface area rapidly decreased from 27 to 45 DASP, then stabilized for about 24 days, followed by another decrease until seed-maturation ( Figure 4C). The R d and P (g,max) per seed exhibited a dramatic increase from 27 to 91 DASP, followed by a pronounced decrease during the development of the seed (Figure 4D). The R d and P (g,max) per dry weight indicated a rapid decreasing trend from 27 to 45 DASP followed by a relatively stable value for about 24 days, and then another decrease as seed development proceeded (Figure 4E).

Chlorophyll Fluorescence in Sarcotesta of the Seeds
The maximum quantum yield of PSII (F v /F m ) in the sarcotesta of the seeds was maintained at 0.74∼0.76 from 27 to 138 DASP, after that the (F v /F m ) significantly decreased with further seed development ( Table 1). The actual photochemical efficiency of PSII (Φ PSII ), the photochemical quenching coefficient (qP) and the photochemical quantum yield of open PSII centers in the light-adapted state (F v '/F m ') maintained a relatively constant value from 27 to 91 DASP. However, the Φ PSII and qP significantly decreased from 91 to 106 DASP, while the F v '/F m ' was constant. The F v /F m in the sarcoteata of the seed was significantly lower at 138 DASP than that at 106 DASP. All of the parameters Φ PSII , qP and F v '/F m ' significantly decreased from 138 to 149 DASP (Table 1).

Chlorophyll Content, Rubisco Activity and Its Relationship with P (g, max) in the Sarcotesta
The total chlorophyll content (Chl a+b ) per area of sarcotesta of the seeds showed a rapid decrease from 27 to 69 DAPS in 2015 ( Figure 5). This initial drop was about 63% as compared with the first measured value (Figure 5). After that, the total chlorophyll content showed an increase from 70 to 106 DASP, and then a significant decrease as seed development proceeded. The Chl a/b ratio of sarcotesta of the seeds remained relatively low for a long period, then reached its peak value at 138 DASP.
The Rubisco activity of sarcotesta of seed showed a decreasing trend during the seed development ( Figure 6A). The reduction in Rubisco activity was 53.9% from 27 to 69 DASP and 43.8% from 70 to 138 DASP, respectively. The Rubisco activity of sarcotesta was significantly positively correlated with its P( g, max ) ( Figure 6B).

Nitrogen (N) Content in Sarcotesta and Its Relationship with P (g, max)
The N content of sarcotesta showed a dramatic decrease from 27 to 69 DASP ( Figure 7A). After that the N content level was constant until 106 DASP, after which it rapidly decreased as seed development proceeded (Figure 7A). The N content in sarcotesta was linearly and positively correlated with P (g,max) per dry weight of the seeds (Figure 7B, P = 0.0013).

Seed and Leaf Photosynthesis and the Contribution of Seed Photosynthesis to Seed Carbon Requirements
Seed photosynthesis and leaf photosynthesis were compared by measuring both of them under the same PAR conditions. FIGURE 3 | Fresh weight (A,D), dry weight (B,E) and dry weight/fresh weight ratio (C,F) of seeds of Torreya grandis cv. "Merrillii" during seed development in 2015 and 2016. For each quantity the mean ± SD is indicated. Different letters denote significant differences at P ≤ 0.05 level. n = 4 trees. DASP, days after seed protrusion. Similar patterns were found both in 2015 ( Figure 8A) and 2016 ( Figure 8B): In the seeds, the rate of gross photosynthesis P g was very high at the beginning of the seed development. After that it decreased rapidly. In the leaves, a contrasting pattern from low rates of P g at the beginning of the growing season to higher ones later was found. Accordingly, the rate of P g was at early stages much higher in seeds than in leaves, whereas later this was reversed (Figure 8).

Chloroplast Ultrastructure in Leaves and Sarcotesta
The chloroplast ultrastructure was examined on one occasion at 81 DASP in 2015. Chloroplasts were spherical both in the leaves and in the sarcotesta (Figure 9). There was no significant difference in the chloroplast length, width, or grana number between leaves and sarcotesta ( Table 2); however, the granal lamellae number was significantly greater in sarcotesta than in leaves (P < 0.05; Figure 9 and Table 2). Leaves at 81 DASP possessed well-differentiated chloroplasts that contained numerous layers of grana and well-developed stromal lamellae (Figures 9A,B). The system of granal and intergranal thylakoids was fully developed also in seeds, accompanied by many plastoglobuli and starch grains in the chloroplasts of sarcotesta (Table 2, Figures 9C,D). However, no plastoglobuli or starch grains were observed in the chloroplasts in leaves (Figures 9A,B).

Effects of DCMU Incubation on the Seed Dry Weight and Seed Oil Quality
The photosynthetic contribution of seeds to the seed weight was assessed by measuring the reduction in their weight after DCMU application ( Table 3). This was carried out separately with intact seeds, and with seeds where the sarcotesta were removed before measuring the weight. The fresh weight of the intact seeds and seeds with sarcotesta removed was reduced by 22.0 and 5.4% from that of the control seeds, respectively. The contribution of seed photosynthesis to the dry weight of the seeds, as estimated by FIGURE 4 | CO 2 exchange parameters in seeds of Torreya grandis cv. "Merrillii" during seed development in 2015. The light saturation point LSP (A), the apparent quantum use efficiency AQE (B), the rate of CO 2 exchange on surface area (C), single seed (D) and dry weight (E) basis. In the rate of CO 2 exchange (C-E), both the maximal rate of gross photosynthesis at light saturation (P (g,max), filled circles) and respiration rate (R d , open squares) are indicated. For each quantity the mean ± SD is indicated. Different letters denote significant differences at P ≤ 0.05 level. n = 4 trees. DASP, days after seed protrusion.
the reduction in the dry weight caused by the DCMU application, was in the intact seeds and seeds with sarcotesta removed about 25.4 and 40.0%, respectively ( Table 3). The DCMU application significantly reduced also the length and width of the seeds. In addition to reducing the size and weight of the seeds, the DCMU application significantly reduced the oil content of the seeds by about 25.5% as compared with the control seeds ( Table 3).

Seed Development in Torreya
The dry/fresh weight ratio of fruits can be used to assess the different phases of fruit growth (Imai and Ogawa, 2009). In the present study, seed development was divided into two phases according to the seed growth pattern. The first phase was the rapid seed growth stage, which was characterized by a rapid increase in seed size. This occurred from 1 to 69 DASP in 2015 and from 1 to 46 in 2016 (Figure 2). During this period, the seed dry weight/fresh weight ratio decreased to its lowest value (Figures 3C,F), indicating that seed water content increased. In the second phase (the seed filling stage), the increase of seed  length and width slowed as seed dry weight continued to increase. The main oil synthesis was obviously during this period, as suggested by the result of Tian et al. (1989). They found that the oil content in Torreya seeds rapidly increased from 21.4 to 54.5% during the period from 5th July to 5th September.

Photosynthetic Characteristics in Torreya Seed Sarcotesta
The seeds had in the present study relatively high light saturation point LSP the period from 27 to 91 DASP; however, the LSP decreased rapidly from 92 to 106 DASP ( Figure 4A). The site where the Torreya cultivar was studied is in Southeastern China, having both high light intensity and high air temperature during the summer months from July to August, corresponding to the period from 75 to 138 DASP in this study. Thus, during the second stage of seed development, the seeds may have been exposed to an excess light which cannot be used for photosynthesis, so that and their PSII were then susceptible to photoinhibition owing to the decreased LSP (Galmés et al., 2007). The photoinhibition of PSII can be detected in vivo by a significant decrease in the F v /F m (Krause and Weis, 1991). However, no change in F v /F m was observed in seed sarcotesta of the Torreya cultivar between 91 and 106 DASP ( Table 1), indicating that the photoinhibition did not occur in the sarcotesta during this period. However, the Φ PSII was significantly lower at 106 DASP than that at 91 DASP. This is consistent with the findings of Zhang et al. (2011) who reported that soybean leaves dissipate the excess excitation energy thermally by the down-regulation of PSII activity to protect their photosynthetic apparatus from photodamage under drought condition. The parameter qP indicates the proportion of PSII reaction centers that are open (Maxwell and Johnson, 2000). Since Φ PSII = qP×F v '/F m ', the significant reduction of Φ PSII found in the present study from 91 to 106 DASP resulted primarily from an increased proportion of closed PSII reaction centers, as indicated by the lower values of qP; while values of F v '/F m ' were essentially the same during that period ( Table 1). These results suggest that large proportion of the PSII reaction centers were closed at 106 DASP so that they may have operated as quenchers to help PSII from photoinhibition (Huner et al., 2006). There is another possible reason for the significant reduction of LSP during 91 to 106 DASP, while with little change in F v /F m of the seed sarcotesta. Plant may decrease photosynthetic rate to re-balance the source-sink ratio when carbohydrate levels reached a threshold value in leaves (Paul and Foyer, 2001). Accordingly, we suggest that the size growth of the seeds (storage organ) was close to completion during that period from 91 to 106 DASP (Figure 2): so that more carbohydrate accumulation in seed sarcotesta may have taken place. That may have caused the reduction of photosynthetic rate and LSP. To test this hypothesis, further studies are necessary to investigate the trend of carbohydrate content in sarcotesta during the seed development.
Furthermore, the apparent quantum use efficiency (AQE) was between 0.04 and 0.06 mol CO 2 (mol photons) −1 in the first stage of seed development (from 1 to 69 DAPS, Figure 4B). These values of AQE are higher than those observed earlier the Torreya leaves associated to the seeds (0.01-0.02 mol CO 2 mol −1 , Huang, 2015). Similar results of higher AQE in capsule wall compared with its associated leaves has also been observed in cotton (Hu, 2013). These results suggest that the seed sarcotesta can effectively use low light for photosynthesis, especially in the early morning or late afternoon, which is important for yield formation.
Respiration provide the energy for the synthesis of new phytomass (Amthor, 2000). In the present study, the rate of dark respiration R d calculated per seed dry weight was very high at the initial stage of seed development, then R d rapidly decreased. These changes in the rate of dark respiration are consistent with the rate of increase in the seed dry mass (Figures 3B,E), supporting the notion of respiratory cost of synthesis of new biomass.
The P (g,max) expressed per unit surface area significantly decreased by 18.7% between 27 and 91 DASP ( Figure 4C); however, the P (g,max) per seed dramatically increased during the same time ( Figure 4D). This may be due to a rapid increase in the surface area of the seeds (Figure 2C). The seed area significantly increased by 18.2% also later from 91 to 138 DASP; however, the P (g,max) per surface area significantly decreased during the same time by 49%. Thus, a significant decrease in the P (g,max) per FIGURE 6 | Ribulose-1,5-biphosphate carboxylase (Rubisco) activity (A) and the correlation between the maximal rate of gross photosynthesis at light saturation (P (g,max) ) per area, and Rubisco activity (B) in sarcotesta of seeds of Torreya grandis cv. "Merrillii" during seed development in 2015. In (A) the mean ± SD is indicated. Different letters denote significant differences at P ≤ 0.05 level. n = 4 trees. DASP, days after seed protrusion. The solid line represent the best-fit linear regressions between the maximal rate of gross photosynthesis at light saturation (P (g,max) ) per area, and Rubisco activity: *P ≤ 0.05; **P ≤ 0.01; *** P ≤ 0.001; ns, not significant.
FIGURE 7 | Nitrogen content (A) and the correlation between the maximal rate of gross photosynthesis at light saturation P (g,max) per dry weight, and nitrogen content (B) in sarcotesta of seeds of Torreya grandis cv. "Merrillii" during seed development in 2015. In (A) the mean ± SD is indicated. Different letters denote significant differences at P ≤ 0.05 level. n = 4 trees. DASP, days after seed protrusion. The solid line represent the best-fit linear regressions between the maximal rate of gross photosynthesis at light saturation P (g,max) per dry weight, and nitrogen content: *P ≤ 0.05; **P ≤ 0.01; ***P ≤ 0.001; ns, not significant. seed during the 91 to 138 DASP stage can be attributed to the decreased photosynthetic rate per seed surface area.
Chlorophyll plays an important role in the light absorption and energy transduction, which are the basis of photosynthesis. In the present study, the initial rapid decrease in chlorophyll content per area from 27 to 69 DASP (Figure 5) could have been caused to a large extent by a dilution effect due to the fast increase in seed volume. There was some increase in the chlorophyll content from 91 to 106 DASP, when the P (g,max) per surface area decreased by 22.7%. These findings indicate that the decrease of seed photosynthesis during the period from 91 to 106 DASP was not caused by low chlorophyll content. This notion is in agreement with the general observation that light-saturated photosynthesis of leaves is largely unrelated to their chlorophyll content (Björkman, 1981).
However, the Chl a+b content significantly decreased by 23% accompanied by significantly decreased F v /F m form 106 to 138 DASP (Table 1 and Figure 5). The Chl a/b ratio was significantly higher at 138 DASP than that at 106 DASP (Figure 5), indicating that degradation of photosynthetic pigments, especially Chl b took place during that period. Therefore, the photosynthetic apparatus seemed to begin to be damaged at the late second development stage of the seeds.
Nitrogen may be a key regulator of photosynthesis in plants. In leaves, more than half of the nitrogen is allocated to photosynthetic proteins, such as Rubisco (Evans, 1989), and the amount, or activity, of Rubsico correlates strongly with photosynthetic rate (Evans, 1989;Reich and Walters, 1994;Priwitzer et al., 1998;Hrstka et al., 2005). Thus, nitrogen content has been used as a qualitative measure of the Rubisco content   (Makino, 2003). In the present study, the Rubisco activity and nitrogen content of sarcotesta both showed decreasing trend during seed development (Figures 6A, 7A). Moreover, a positive correlation was found between Rubisco activity and P (g,max) per surface area ( Figure 6A). The nitrogen content in sarcotesta significantly decreased by 25% from 106 to 138  Data are mean ± SD. Different letters denote significant differences at P ≤ 0.05 level. n = 5 trees.
DASP ( Figure 7A). We also observed a significant positive correlation between the P (g,max) calculated per dry weight and nitrogen content in the sarcotesta ( Figure 7B). Therefore, the decreased Rubisco activity caused by the decreased nitrogen content may have caused the damage of the photosynthetic apparatus. Liu et al. (2014) reported that the nitrogen content in Torreya kernels significantly increased by 50% from 19 July to 19 August, corresponding to 95 ∼ 126 DASP in present study. It has been reported that nitrogen in foliage and pericarp of almond is rapidly imported into the embryo (Weinbaum and Muraoka, 1986). Similarly, the decrease in percent nitrogen in sarcotesta in the present study suggests that the amount of nigtrogen-containing compounds, such as proteins or amino acids, decreased, as they were rapidly imported into the developing seed.

The Relative Contribution of Seed Photosynthesis to Seed Development
The photosynthetic rate per seed area significantly decreased from 32 to 69 DASP, which corresponds to the first stage of seed development between 17th May and 23rd June, but was then significantly higher than the rate in the associated current-year leaves which began to sprout 30th April and reached the maximum rate of photosynthesis 20th July. Thus, we suggest that the photosynthesis in seeds plays an important role in seed growth during the first stage of seed development. However, at the second stage of seed development, when seed size growth is close to completion, the constant photosynthetic The data indicate the contents of the various fatty acids as percentages in oil content. Data are mean ± SD. Different letters denote significant differences at P ≤ 0.05 level. n = 3 trees.
rate in the associated leaves, accompanied by the significantly decreased photosynthetic rate in seed sarcotesta during that time (Figure 8), suggest a less important role of seed photosynthesis on dry matter accumulation during that period. The chloroplasts of leaves possessed many well-organized grana thylakoids with no plastoglobuli or starch grains, while those in sarcotesta contained a large number of granas accompanied by many plastoglobuli and starch grains (Figure 9). Plastoglobuli are usually formed in senescing chloroplasts (Lichtenthaler, 2013). Our observation of these characteristics indicates that the chloroplasts of sarcotesta began to senesce in the second stage of seed development. In addition, the accumulated starch grains can dislocate the grana (Bondada and Oosterhuis, 2003) and push them toward the periphery of the chloroplast, which might impede photosynthetic activity. The import of carbon from photosynthetically active source tissues is needed for seed development (Gallardo et al., 2008), and it also affects the synthesis of fatty acids (Xiong et al., 2010). Rolletscheck et al. (2005) have shown that illumination of developing soybean seeds causes photosynthetic release of significant amounts of O 2 . This seed photosynthesis plays a key role in the local energy state, storage metabolism and flux of carbon to lipid biosynthesis. Although green leaves are the main sources of photosynthate production, sarcotesta photosynthesis was also supplies assimilates for oil synthesis. Indeed, our results shows that the reduction in dry weight and oil content of the intact seed induced by DCMU application was 25.4 and 25.5%, respectively ( Table 3). This indicates that the reduced supply of seed sarcotesta photosynthates caused by the DCMU application inhibited oil synthesis. Results from earlier studies suggests that the contents of polyunsaturated fatty acids are significantly increased as a response to abiotic stresses, such as freezing stress, salt stress, oxidative stress etc. (Watanabe et al., 2004;Zhou et al., 2010;Zhang et al., 2012). Considering that DCMU application inhibits the photosynthetic electron transport, the increased polyunsaturated fatty acids content caused by DCMU application can be explained as an adaptive response to the deficiency of the photosynthate supply from seed photosynthesis. This is in line with Wang et al. (2016) who found that green siliques photosynthesis greatly influences seed oil content and oil composition of oilseed rape.

CONCLUSIONS
Our results show that seed photosynthesis is important for the biomass accumulation and oil synthesis in Torreya grandis cv. "Merrillii" seeds. The seed photosynthesis is especially active during the first stage of seed development, when the photosynthetic production in the attached current-year leaves is much lower. It appears that the down-regulation of PSII activity, of Rubisco activity, and of Chl (a+b) , and nitrogen contents in the sarcotesta are well coordinated in order to rebalance the source-sink ratio at the late second stage of seed development. Simultaneously the photosynthetic activity of the attached current year leaves increases, so that the photosynthetic activity is higher in the leaves than in the seeds during the second stage of seed development. Our results also show that seed photosynthesis inhibition affects the fatty acid composition of the seeds. These results have implications for cultivation practices in the future, as high oil content require efficient photosynthesis during the seed development stage of the Torreya cultivar. Therefore, the high photosynthetic activity in sarcotesta may be a useful trait for the selection of cultivars with high yield and high oil contents.

ACKNOWLEDGMENTS
We thank LvJia Lu and Han Yao for their support in plant culture and technical assistance and Yin Xu for the determination of chloroplast ultrastructure. This study was supported by the National Science Foundation of China (grant 31570616), the open project funds for forestry discipline in Zhejiang province (grant KF201312), the key technical integration and demonstration promotion of the efficient ecological cultivation of Torreya young forest (grant [2015] No. TS03), selective breeding of new cultivars in Torreya. grandis (grant 2016C02052-12).

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2017. 02179/full#supplementary-material Figure S1 | The standard gas chromatogram of the main fatty acids in seeds of Torreya grandis cv. "Merrillii." Figure S2 | Light responses of net photosynthesis P n in seeds of Torreya grandis cv. "Merrillii" during seed development; expressed on surface area (A), dry weight (B), and single seed (C) basis. For the P n values mean ± SD is indicated. Different letters denote significant differences at P ≤ 0.05 level. n = 4. The light responses were measured on eight dates indicated by the symbols and DASP values in the upper right-hand corner of the uppermost panel. DASP, days after seeds protrusion. Note: The values of P n are negative because of high rates of respiration. The rates were high because the measurements were carried out with intact seed containing also the non-photosynthetic embryo in addition to the photosynthetic sarcotesta. The absolute value of the negative P n indicates the re-fixation by the sarcotesta photosynthesis of the CO 2 released from respiration. See introduction of the main text for further explanation.
Table S1 | The inhibitive effect of different concentrations of DCMU on the ETR of seed in comparison with the control.