- College of Agronomy and Biotechnology, China Agricultural University, Beijing, China
The kernel setting of maize varies greatly because of the timing and intensity of water deficits. This variation can limit leaf productivity (source), the translocation of assimilated sugars (flow), and yield formation (sink). To explain the decline in kernel setting of maize under water deficits from the perspective of source-flow-sink, a 3-year experiment was conducted under a rain shelter. Five water regimes were studied. One regime included well-irrigated (CK) treatment. Four regimes involved water deficits: irrigation was withheld during the 6- to 8-leaf stage (V6−8), the 9- to 12-leaf stage (V9−12), the 13-leaf stage to tasseling stage (V13−T), and the silking stage to blister stage (R1−2). Water deficit effects on kernel setting began when the water deficit occurred at V9 and became more significant with time. Kernel weight was reduced by 12 and 11% when there were water deficits during V9−12 and V13−T, respectively. This was the result of reduced leaf area (limited source) and an altered vascular bundle in the ear peduncles (limited assimilate flow). The reduced vascular bundle number, rather than the ear peduncle cross-sectional area, significantly affected the final kernel weight when exposed to a water deficit prior to the silking stage. The water deficits prior to and close to the flowering stage significantly reduced ear kernel number; that is, 14 and 19% less during V13−T and R1−2, respectively, compared with the kernel number during the CK treatment. This reflects a smaller sink under water deficit conditions. Additionally, ovary size was reduced the most in the V13−T water deficit compared with other treatments. After rewatering, the water deficit before or during flowering stage continued to have residual effects on grain-filling in the late growth period. The grain-filling rate decreased under the V9−12 water deficit; the grain-filling duration shortened under the R1−2 water deficit; and both negative effects occurred under the V13−T water deficit. This study clearly indicated that (1) the water deficit during the vegetative organ rapid growth period both limited leaf source development and assimilate flow and slowed down kernel development, and (2) the water deficit just before and during flowering reduced kernel sink. Deficits at both times could retard grain-filling and reduce maize yield. The results of the present study might guide irrigation practices in irrigated maize or inform the management of sowing time in rainfed maize, to desynchronize the water deficit and the plant’s reactions to such deficits at different stages.
Introduction
Maize is the most important cereal crop in the world, and maize production faces great challenges in increasing or maintaining yields, based on severe abiotic stresses that are brought about by climate change (Lobell et al., 2011; Lesk et al., 2016). The most severe abiotic stress restricting maize yields is the presence of water deficits; these deficits account for approximately 50% of total maize yield (tons) loss due to abiotic stresses in China (China Statistical Yearbook, 2016).
Maize grain yield is sensitive to water deficits (Schussler and Westgate, 1991; Setter et al., 2001). Kernel number and kernel weight are the most impacted components (Borrás and Gambín, 2010). The kernel number is more flexible and is more closely correlated with grain yield than the kernel weight, especially under abiotic stress (Andrade et al., 2002). Kernel numbers can be reduced by up to 100% because of water deficits, depending on the time and intensity of the deficits (Schussler and Westgate, 1991; NeSmith and Ritchie, 1992a; Cakir, 2004). The loss of kernel numbers under a water deficit can be attributed to incompletely developed florets (Rossini et al., 2016), inhibited silk emergence (Fuad-Hassan et al., 2008; Oury et al., 2016), and zygotic abortion (Zinselmeier et al., 1999; McLaughlin and Boyer, 2004a). In contrast, kernel weight is more stable under a water deficit, and any declines usually occur during the late growth period because of the reduced leaf area (NeSmith and Ritchie, 1992b; Mansouri-Far et al., 2010). The reduction of kernel weight under a water deficit is related to impaired assimilation supply (Schussler and Westgate, 1995; Setter et al., 2001), disturbed sugar and hormonal metabolism (Ober et al., 1991; Zinselmeier et al., 1999), and decreased cell division and starch accumulation in the endorsperm (Setter and Flannigan, 2001; Yu and Setter, 2003). In other crops, kernel weight also relates to the development of parental organs prior to flowering (Yang et al., 2009; Hasan et al., 2011; Castillo et al., 2017); however, there is little relevant research on maize.
Kernel setting is associated with a source-sink relationship (Borrás and Otegui, 2001; Borrás et al., 2003b); this relationship is an important determinant of maize yield (Borrás and Gambín, 2010; Yu et al., 2015). Source activity and sink capacity can be affected by water deficits by reducing leaf area, accelerating leaf senescence (Prochazkova et al., 2001), and decreasing leaf photosynthetic rate (Pn) and its associated metabolism (Setter et al., 2001). The limited sucrose levels under water stress decrease the number of endosperm cells and starch granules, thereby, reducing sink capacity (Ober et al., 1991; Setter and Flannigan, 2001). The source and sink are not affected independently by water deficits; rather, they restrict each other. Carbohydrate shortages result in kernel abortion and fewer kernels. This can impair sink strength (Zinselmeier et al., 1999; McLaughlin and Boyer, 2004b). Consequently, the decreased kernel capacity results in sugar accumulation in the leaf and the stem, especially post-silking, which downregulates photosynthesis in the leaf (Rossi et al., 2015). The altered source-sink relationship reduces plant biomass accumulation and reduces the translocation of assimilated sugars to the kernels (Jurgens et al., 1978). When considering source-sink relationships under water deficits, the flow of assimilate from source to sink is significantly reduced under water deficits. However, more relevant evidence is needed to demonstrate these outcomes.
Therefore, this study analyzed the causal factors that reduce maize yield under water deficits in terms of source, assimilate transport, and sink, both separately and together as a system. To achieve this objective, we measured (1) grain yield and yield components; (2) the leaf area in different canopy layers and ear leaf photosynthesis at mid grain-filling; (3) vascular bundle size and number in the ear peduncle; and (4) ovary development and grain-filling dynamics. The study was conducted over 3 years under different water deficit treatments.
Materials and Methods
Experimental Site and Agronomic Management
The experiments were conducted in a rain shelter at the Shangzhuang Experimental Station at the China Agricultural University, Beijing, China (40° 08′N, 116° 10′E) during 2014–2016. The rain shelter contained 18 big ponds (2 m long, 4 m wide, and 1.8 m deep) and 18 small ponds (2 m long, 2 m wide, and 1.8 m deep). Each pond was cemented on the four sides and the bottom, and each pond was filled with calcareous fluro-aquic soil. The soil contained 9240 mg kg−1 organic matter, 519 mg kg−1 total nitrogen (TN), 68.3 mg kg−1 extractable potassium (K2O), and 16.6 mg kg−1 extractable phosphorous (P) at soil depths of 0–40 cm. The field capacity was 0.32 g g−1 at a soil depth of 0–20 cm and 0.25 g g−1 at a soil depth of 20–40 cm. The soil bulk density was 1.38 g cm−3 at a soil depth of 0–20 cm and 1.59 g cm−3 at a soil depth of 20–40 cm.
Prior to sowing, 60 kg N ha−1, 90 kg phosphorus pentoxide (P2O5) ha−1, and 150 kg K2O ha−1 were applied evenly to the soil surface. The fertilizer together with soil was then plowed under with a shovel, and an additional 120 kg N ha−1 was applied by furrowing at the 12-leaf stage (V12). The most popular maize hybrid Zhengdan 958 in China was used. Seeds were planted on June 10th in all 3 years and were harvested after the seeds reached physiological maturity. Three to four seeds per hole were sowed. The seedlings were thinned to 7.5 plants m−2 at the 3-leaf stage (V3) with a row spacing of 60 cm and plant spacing of 20 cm. Weeds, pests, and diseases were well controlled during the experimental seasons.
Experimental Design
The experiment was arranged in a completely randomized block design with three replications in 2014 and four replications in 2015 and 2016. Five treatments were conducted: (1) well-irrigated (CK) treatment with irrigation at the 6-leaf (V6), 9-leaf (V9), 13-leaf (V13), silking (R1), and blister (R2) stages; (2) introducing a water deficit during the 6–8-leaf stage (V6−8) by skipping irrigation once at V6; (3) introducing a water deficit during the 9–12-leaf stage (V9−12) by skipping irrigation once at V9; (4) introducing a water deficit during the 13-leaf stage to tasseling stage (V13−T) by skipping irrigation once at the V13 stage; and (5) introducing a water deficit during the silking stage to blister stage (R1−2) by skipping irrigation once at silking. Full irrigation was applied approximately 10 days before sowing each year. During the maize growth season, for the CK treatment, the soil water content was maintained above 50% of field capacity at a soil depth of 0–60 cm (Andrade et al., 2002). For the four water deficit treatments, the water content dropped below 50% of field capacity during the water deficit period associated with the different treatments. The relative soil water content is shown in Supplementary Figure S1.
Sampling and Measurements
Relative Soil Water Content
Soil samples were collected using an auger at three different soil layer depths (0–20, 20–40, and 40–60 cm) to determine the water content at 1 or 2 days before V6, V9, V13, R1, and R2 from 8 big ponds and 8 small bigs (2 replication per pond). The soil samples were oven dried at 105°C for 24 h. Soil water content (SWC) was calculated using the following formula:
The relative soil water content (RSWC) was calculated using the following formula:
The amount of irrigation (IW) was modified according to the method described by Tari (2016). The total irrigation amount was calculated by adding the IW for the three different soil layers (0–20, 20–40, and 40–60 cm). The IW at each layer was calculated using the following formula:
In this expression, IW is the amount of irrigation at each layer (m3); BD is the soil bulk density (kg m−3); D is the depth of each soil layer (0.2 m); A is the area of each plot (m2); FC is field capacity; and SWC is the gravimetric soil water content.
Leaf Area and Photosynthetic Rate
Three plants in 2014 and 2015, six plant in 2016 were randomly selected in each treatment at the tasseling stage to measure the length and width of each leaf. The leaf area was calculated as the leaf length × leaf width × 0.75 (Pearce et al., 1975). We recorded the leaf area at the “ear layer” (including ear leaf, one leaf above it, and one leaf below it), at the “above ear layer” (all leaves above ear layer), and at the “ear and above layer” (including leaves at both ear and above ear layers) (Supplementary Figure S2). The leaf area for each of these three layers was calculated by adding its included leaf areas (Subedi and Ma, 2005; Ning et al., 2017). The reduction of leaf area was expressed as the percentage of leaf area reduction at a specific layer for the plants exposed to the water deficit treatments compared with the plants exposed to the CK treatment.
Another three plants were randomly selected in each plot at the mid grain-filling stage [85 days after sowing (DAS)] in 2016. The net Pn of the ear leaf of each plant was measured with LI-6400 XT (Li-Cor Inc., Lincoln, NE, United States) at 10:00–12:00 am. Prior to measuring the leaf Pn, 400 μmol mol−1 CO2 and 500 μmol s−1 flow rate were set. The light flux density and temperature in the leaf chamber were consistent with ambient conditions.
Vascular Bundle and Ear Peduncle
Ear peduncles next to the base of the ears were sampled twice, one was just before tasseling (49 DAS), another is 10 days after silking (65 DAS) in 2016. The peduncles collected at 49 DAS were immediately fixed after sampling with a mixture of 10% formaldehyde, 50% ethanol, and 5% acetic acid over 24 h in a refrigerator at −4°C. The samples were then dehydrated solutions with different concentrations of ethanol: 75% ethanol for 4 h, 85% ethanol for 2 h, 90% ethanol for 2 h, 95% ethanol for 1 h, and twice in 100% ethanol for 0.5 h. The dehydrated peduncles were infiltrated with mixtures of xylene and alcohol five times (ethanol:xylene (v/v) = 2:1 for 10 min once, ethanol:xylene = 1:1 for 10 min once, ethanol:xylene = 1:2 for 10 min once, and then 100% xylene for 10 min twice). Afterward, the ear peduncles were embedded in paraffin wax and were cut into 4 μm sections using a Leica RM 2016 microtome (Leica Shanghai Instrument Co., Ltd., Shanghai, China). Finally, the peduncle sections were de-waxed and stained with 0.5% safranin and 0.5% fast green according the method described by Liu et al. (2006). Micrographs were taken using a Nikon Elipse Ci (Nikon Instruments Inc., China) and then analyzed using a CaseViewer (3DHISTECH Ltd. The Digital Pathology Company).
Cross sections of three ear peduncles were collected from big ponds (3 replicates) at 65 DAS. Razor blades were used to cut freehand slices (Xu et al., 2017), which were then stained with safranin (0.5%, w/v). Micrographs of the cross sections were taken with an Olympus SZ61 stereomicroscope (Olympus Imaging China Co., Ltd., Beijing, China). Vascular bundle size and number and ear peduncle area were measured using an Image-Pro Plus 6.0 (Media Cybernetics Inc., 2006).
Ovary Size
Fresh ovaries at the 5th kernel ring of the ear were manually cut longitudinally through the middle of the embryo (Leroux et al., 2014). The samples were then stained with 0.05% toluidine blue O for 15 min. The samples were then rinsed with distilled water for 1 min. Micrographs were taken with an Olympus SZ61 stereomicroscope.
Grain-Filling and Its Characteristics
The ears of 18 plants were randomly selected for each treatment and then bagged before silking occurred in 2015. Artificial pollination was conducted 4 days after silking according to the method described by Cárcova and Otegui (2001). Three representative ears in each treatment were harvested to measure kernel weight at 5, 15, 25, 35, 45, 60 days after pollination.
Kernels were classified as either inferior or superior based on their position in the ear. After sampling, the ears were cut into two parts at a point that was located at two-thirds of the distance between the bottom of the ear and the top. The kernels from the top one-third of the ear were defined as inferior kernels, while the rest were defined as superior kernels. Four full rows of kernels from the top one-third of the ear and two full rows of kernels from the bottom two-thirds of the ear were collected separately and dried to a constant weight at 70°C to determine kernel dry weight (Gao et al., 2017).
The grain-filling process was divided into three periods according to the study by Borrás and Gambín (2010): lag period, the active grain-filling period, and the maturation drying period. The three grain-filling periods were determined by calculating the second derivative of a logistic growth equation [W = a/(1 + be−ct)], where ‘W’ represents kernel weight (mg); ‘a’ represents the potential kernel weight; ‘b’ and ‘c’ are coefficients, which were determined by regression; and ‘t’ represents the number of days after pollination (Wang et al., 2010). The characteristic parameters of the grain-filling rate (Vlag, Vactive, and Vdrying) and the duration (Tlag, Tactive, and Tdrying) of each grain-filling period were calculated based on the derivative of the logistic growth equation (Raynaud, 2010; Gao et al., 2017). The starting time of the active grain-filling period was defined as t1 = [ln b – ln (2 + 31/2)]/c and ending time was t2 = [ln b + ln (2 + 31/2)]/c (Supplementary Figure S3). The duration of the entire grain-filling process was t3 = (ln b + 4.595)/c. Therefore, the duration of lag period (Tlag) was Tlag = t1, active grain-filling period (Tactive) was Tactive = t2 – t1, and maturation drying period (Tdrying) was Tdrying = t3 – t2. The corresponding grain-filling rates of lag period (Vlag), active grain-filling period (Vactive), and maturation drying period (Vdrying) were the daily kernel weight increase during its corresponding duration.
Grain Yield and Yield Components
All the ears from each plot were harvested, counted, and weighed once physiological maturity was reached. From these ears, 20 representative ears from big ponds and 15 ears from small ponds during 2014, 20 representative ears from all ponds during 2015, and the rest of the ears from all the ponds during 2016 were selected for each plot. The row number and kernel number of the two opposite rows of each ear were counted to determine the kernel number per ear. Kernels from these ears were then manually threshed and oven-dried at 70°C to a constant weight. Grain yield was calculated based on the dry kernel weight of each ear and ear number in each plot. The yield was adjusted to 14% moisture content. In addition, three replicates of 1000 kernels in each plot were counted, oven-dried, and weighed to determine the 1000 kernel weight.
Data Analysis
Analysis of variance (ANOVA) was conducted using the SAS 9.0 statistical package (SAS Institute, 2004). The least significant difference (LSD) was used to determine significant differences among treatments at the 0.05 probability level. Figure drawing, grain-filling curve fitting, and correlation analysis were performed using SigmaPlot 12.5 (Systat Software Inc., 2013).
Results
Grain Yield and Yield Components
The timing of water deficit significantly affected the maize yield and yield components (Table 1 and Supplementary Table S1). The maize yield significantly decreased by 8.3% (P < 0.001) as a result of the V6−8 water deficit, compared with the result of the CK treatment. There were no significant differences in the kernel number and thousand kernel weight (TKW) when comparing the V6−8 water deficit treatment and the CK treatment. Both the kernel number and TKW significantly decreased as a result of the V9–R1 water deficit; this led to a significant decline in maize yield (23.8–33.7%). In addition, the lowest TKW was observed as a result of the V9−T water deficit; the smallest kernel number was observed as a result of the V13-R2 water deficit.
TABLE 1. Average reduction of maize grain yield (%), kernel number per ear (%), and kernel weight (%) that resulted from water deficits during 2014–2016.
Leaf Area at Different Layers and Photosynthetic Rate at Mid Grain-Filling Stage
The water deficit treatments led to reductions in maize leaf area at the different layers. The reduction depended on water deficit timing (Table 2 and Supplementary Table S1). The V6−8 water deficit led to a significant decrease in leaf area in the ear layer, 7.9% (P < 0.001), compared with the leaf area observed with the CK treatment. The V9−T water deficit led to a significant decrease in leaf area for both the ear and above ear layers (P < 0.001). However, the leaf area reduction at the ear layer was greater as a result of the V9−12 water deficit compared with the V13−T water deficit (16.8% vs. 10.5%). In contrast, an opposite trend was observed at the above ear layer, with the V13−T water deficit leading to a greater reduction in leaf area (35.4%) compared with the V9−12 water deficit (20.3%). No leaf area reduction was observed as a result of the R1−2 water deficit compared with the leaf area under the CK treatment. The TKW reduction was positively correlated with the leaf area reduction, and the correlation became significant for leaf area reduction at above ear layer as well as with the ear and above ear layers (Figure 1). However, no significant correlation was observed between leaf area reduction and kernel number reduction (data not shown).
TABLE 2. Average reduction of leaf area (%) at the ear layer, the above ear layer, and the ear and above layers around the tasseling stage that resulted from the water deficits during 2014–2016.
FIGURE 1. The relationship between kernel weight and leaf area reduction at above ear layer (A), ear layer (B), and ear and above ear layer (C) in 2014 (red), 2015 (blue), and 2016 (green). CK, well-irrigated treatment; V6−8, water deficit during the 6–8-leaf stage; V9−12, water deficit during the 9–12-leaf stage; V13−T, water deficit from the 13-leaf stage to tasseling stage, R1−2, water deficit from the silking stage to blister stage. ∗indicates significance at P < 0.05.
The effects of water deficits on ear leaf Pn at the mid grain-filling stage strongly depended on water deficit timing (Figure 2). The V6−8 and V9−12 water deficits led to similar Pn values compared with the CK treatment. In contrast, water deficits after V13 significantly reduced ear leaf Pn at the mid grain-filling stage, especially for the R1−2 water deficit.
FIGURE 2. Photosynthetic rate at 85 days after sowing (mid grain-filling stage) in 2016. CK, well-irrigated treatment; V6−8, water deficit during the 6–8-leaf stage; V9−12, water deficit during the 9–12-leaf stage; V13−T, water deficit from the 13-leaf stage to tasseling stage, R1−2, water deficit from silking stage to blister stage. Values with different letters are significantly different at P < 0.05.
Vascular Bundle at Ear Peduncle
The size and number of vascular bundles of ear peduncles were reduced to some extent under water deficits; the reductions were associated with water deficit timing (Figure 3). Vascular bundle development was delayed as a result of V6−8 (Figure 3B1) and V9−12 (Figure 3C1) water deficit treatments compared with the CK treatment (Figure 3A1) at 49 DAS; however, the vascular bundle size became similar to that observed with the CK treatment (Figures 3A2–C2,F,G) at 65 DAS. The V13−T water deficit treatment postponed the vascular bundle development at 49 DAS (Figures 3A1,D1) and significantly decreased the vascular bundle size at 65 DAS. This was particularly the case for marginal vascular bundles (Figures 3A2,D2,F,G). With the R1−2 water deficit treatment, vascular bundle development was similar to that of the CK treatment at 49 DAS (Figures 3A1,E1); however, the vascular bundle size decreased significantly at 65 DAS (Figures 3A2,E2,F,G). The number of vascular bundles decreased significantly under the V9-VT water deficit treatment (Figure 3H), while the ear peduncle area decreased when the water deficit occurred after V6 (Figure 3I).
FIGURE 3. Transverse section of ear peduncle at 49 (A1–E1) and 65 (A2–E2) days after sowing (DAS); the area of the marginal vascular bundle (VB) (F) and central VB (G); and the total vascular bundle number in ear peduncle (H) and ear peduncle area (I) at 65 DAS (12 days after silking). The marginal VB, closely arranged at the edge of the ear peduncle, was consistently undeveloped and was small in size. The central VB, scattered and distributed in the center of the ear peduncle, was consistently well-developed and had a large size. CK, well-irrigated treatment; V6−8, water deficit during the 6–8-leaf stage; V9−12, water deficit during the 9–12-leaf stage; V13−T, water deficit from the 13-leaf stage to tasseling stage, R1−2, water deficit from the silking stage to blister stage. Values with different letters are significantly different at P < 0.05.
The vascular bundle number reduction was positively and significantly correlated with the reduction in TKW when exposed to water deficits (Figure 4A); however, the ear peduncle area reduction had no significant effects on TKW reduction (Figure 4B).
FIGURE 4. Correlation of kernel number reduction with vascular bundle number reduction (A) and ear peduncle area reduction (B) in 2016. CK, well-irrigated treatment; V6−8, water deficit during the 6–8-leaf stage; V9−12, water deficit during the 9–12-leaf stage; V13−T, water deficit from the 13-leaf stage to tasseling stage, R1−2, water deficit from the silking stage to blister stage. ∗indicates significance at P < 0.05.
Ovary Development and Grain-Filling Dynamics
The ovary size became smaller with delays in the water deficit prior to the silking stage (V6–VT). Ovary sizes recovered slightly when exposed to the water deficit after silking (R1−2) (Figure 5A). Ovary development was delayed as a result of the V9−T water deficit compared with ovary development during the CK treatment at 4 days after pollination (Figure 5A). Accordingly, the final kernel weight was smallest for the V9−12 and V13−T water deficit treatments (Figures 5B,C), especially for the inferior kernels. When the water deficit occurred during the early (V6−8) and late (R1−2) periods, there were no effects or slight effects on TKW.
FIGURE 5. Longitudinal-sectional micrographs (A) of the ovary at 5th kernel ring, stained with toluidine blue at 4 days after pollination, and grain-filling dynamics of inferior (B) and superior (C) kernels in 2015. CK, well-irrigated treatment; V6−8, water deficit during the 6–8-leaf stage; V9−12, water deficit during the 9–12-leaf stage; V13−T, water deficit from the 13-leaf stage to tasseling stage, R1−2, water deficit from the silking stage to blister stage. Values with different letters are significantly different at P < 0.05.
The logistic equation fit well with the grain-filling process when exposed to the water deficit based on their high correlation values, ranging from 0.989 to 0.999 (data not shown). The timing of water deficit affected grain-filling dynamics by changing the duration of grain-filling and the grain-filling rate, especially for inferior kernels as a result of the V13−T and R1−2 water deficit (Table 3). Furthermore, the effective filling period for the inferior kernels was shortened by 3.2–3.4 days under the V13–R2 water deficit, while the maturation drying period decreased by 4.0–4.3 days compared with the CK treatment. There was less than a one-day variation in the grain-filling duration under the V6–V12 (V6−8 and V9−12) water deficit treatment. In contrast, the grain-filling rate changed tremendously in both superior and inferior kernels. For superior kernels, the grain-filling rate during the lag period decreased more significantly as a result of the V9−12 water deficit (11.9%) treatment, whereas the grain-filling rate during the two late periods decreased more significantly as a result of the V13−T water deficit treatment (15.0%) compared with the CK treatment. For inferior kernels, the greatest reduction in grain-filling rate during lag period was observed as a result of the V13−T water deficit treatment (22.2%). The greatest reduction in the effective filling period (6.7%) and the maturation drying period (6.6%) occurred as a result of the V9−12 water deficit treatment.
TABLE 3. Characteristic parameters of grain-filling curve fitted by using the logistic equation, duration of lag periods (Plag), active grain-filling period (Pactive), maturation drying period (Pdrying), grain-filling rate during lag periods (Vlag), active grain-filling period (Vactive), and maturation drying period (Vdrying) that resulted from the water deficit in 2015.
Discussion
The balance between source and sink during grain-filling is of great importance for kernel setting (Jones and Simmons, 1983; Borrás et al., 2002). This balance can be influenced by water deficits, especially around the pollination period (Zinselmeier et al., 1995; Setter et al., 2001). The assimilation supply increases the kernel setting by feeding sucrose during water deficits around the pollination period (Zinselmeier et al., 1999; McLaughlin and Boyer, 2004a,b). This study found that the reduced leaf area, accelerated leaf senescence, and decreased photosynthesis resulting from water deficits reduced leaf source activity, negatively affected kernel setting, and reduced maize yield (Tables 1, 2 and Figure 2). The reduction of vascular bundle size and number at the ear peduncle when exposed to water deficit treatments also limited the transportation of sucrose (Figure 3). This represented a limited potential capacity of assimilate flow.
The study also found a slower ovary development, a shortened duration of grain-filling, and a decreased grain-filling rate as a result of water deficit treatments. This resulted in a reduction in kernel weight, especially for inferior kernels (Table 3 and Figure 5). Both the decline in source activity and sucrose flow strength hindered the development of the kernels and decreased the kernel dry weight, especially for the inferior kernels (Figure 5 and Table 3). Water deficits impacted kernel size and kernel number by altering source strength (Table 2 and Figure 2), the ability of assimilate flow (Figure 3), and sink capacity (Setter and Flannigan, 2001; Borrás and Westgate, 2006). The timing of the stress governed which component was most affected (Table 1).
Source: Leaf Area Dynamics and Photosynthetic Rate During Grain-Filling
Leaf area significantly decreased as a result of water deficits (Table 2); this decline was closely associated with water deficit timing (Cakir, 2004). In this study, the leaf area at both ear and above ear layers decreased significantly as a result of water deficits during the rapid growth period (V9-VT, Table 2). This decreased kernel weight and reduced grain yield by more than 20% (Table 1) were consistent with the findings of NeSmith and Ritchie (1992a). Previous results also indicated that water deficits occurring during the development of leaves in the ear and above ear layers could greatly reduce yields, due to the significant contribution of these leaves to grain yield (Allison and Watson, 1966; Subedi and Ma, 2005).
In addition to impacting leaf development, water deficits also accelerate leaf senescence, especially after silking (Brevedan and Egli, 2003; Sade et al., 2017). As a consequence, green leaf duration is also reduced, limiting the availability of assimilate under water deficits (Chaves et al., 2002; Borrás et al., 2003a). Boyle et al. (1991) and Zinselmeier et al. (1999) indicated that carbon starvation in maize kernel, due to limited assimilation, could result in ovary abortion, reducing the kernel number. The reduction in sink strength as a result of water deficits could result in redundant sugars in the leaf, which could trigger leaf senescence (Adams et al., 2013; Rossi et al., 2015). There was likely feedback from sink limitations on assimilation, transport, and vascular changes; however, there is a lack of relevant evidence to demonstrate this.
This study confirms that water deficits can reduce leaf Pn, with persistent effects. Water deficits occurring around the silking stage, especially during the post-silking stage, can reduce the Pn in the ear leaf at a late growth stage, even in well-watered conditions (Figure 2). This indicates that the recovery of Pn from water stress around the silking stage is difficult or incomplete. This outcome is consistent with the results of Grzesiak et al. (2006), probably because of the damaged enzymes of the photosynthetic system (Tezara et al., 1999; Chaves et al., 2002; Ali and Ashraf, 2011) and the accelerated leaf senescence (Escobar-Gutiérrez and Combe, 2012).
Potential Flow of Assimilate: Vascular Bundle Size and Number in the Ear Peduncle Limits Transportation as a Result of the Water Deficit
The long-distance transportation of sucrose in maize depends on the vascular bundle system that connects leaf vein, stem, and ear peduncle (Ruan et al., 2010; Baker et al., 2016). The vascular bundle in the ear peduncle represents the end of sucrose transportation for ear growth. Larger sizes and an increased number of the vascular bundles in the ear peduncle benefit kernel setting, especially for kernel weight (He et al., 2005, 2007). This study found that the size and/or number of vascular bundles in the ear peduncle were significantly reduced as a result of the water deficit. This was particularly true when the water deficit occurred after V9 (Figure 3). On the one hand, the smaller size and the less number of the vascular bundles likely contributed to the reduced ability of the xylem to transport water and nutrients from the soil to the plant. This could affect the assimilation transport ability and then organ development (Tanguilig et al., 1987; Otegui et al., 1995; Hu et al., 2008). On the other hand, the limited vascular bundle system could decrease daily flow of assimilate from leaves to kernels, reducing the available sugar accumulation in the kernel (Wang, 2011; Feng, 2014). Moreover, the changes in the vascular bundles likely resulted in the incomplete recovery of leaf Pn and kernel setting in the late growth period of maize (Zinselmeier et al., 1999; McLaughlin and Boyer, 2004a,b). The results of the present study suggested that the water deficit had continuous effects on assimilate transportation, even after rewatering. This was due to the changes in the size and number of the flow system components responsible for assimilation in maize. In particular, the water stress occurring during the ear growth stage between the 9-leaf stage and approximately 2 weeks after pollination appears to have had larger effects on the size and number of flow system components. This was due to the important contribution of leaves to expanded yields during this growth period (Figure 3 and Tables 1, 2).
Sink: Hampered Kernel Development and Influenced Grain-Filling
Water deficits during the silking stage greatly affected maize yield, mainly by reducing kernel number (Table 1) in association with pollen amount and activity (Uribelarrea et al., 2002; Aylor, 2004), the anthesis-silking interval (Bolaños and Edmeades, 1990, 1996), and embryo and endosperm development (Borrás and Gambín, 2010; Leroux et al., 2014). All these factors could be influenced by water deficits at the silking stage, thereby, greatly reducing the kernel number. In contrast, kernel weight was less sensitive to water deficits during R1−2. This might be because of the relatively high distribution of assimilate to individual kernels and the lower kernel number under severe water stress (Table 1). These results indicated that maintaining a high kernel number was the key to reducing or eliminating water deficit effects at the silking stage of maize production.
Several studies have examined the impact of water deficits at the silking stage of maize growth. In contrast, few studies have focused on how the number of maize kernels is affected by water deficits during the vegetative growth period. This study found that water deficits before the silking stage could significantly reduce kernel numbers; the effects became stronger over time (Table 1). The reduction in kernel number was likely to be related to the differentiation and development of reproductive organs under water stress (Fuad-Hassan et al., 2008). Setter and Flannigan (2001) showed that water deficits at an early post-pollination stage could lead to decreases in kernel weight by inhibiting cell division and endoreduplication in the endosperm of maize. Compared with the kernel number, the kernel weight of maize was more significantly affected by the water deficits occurring before silking (NeSmith and Ritchie, 1992a; Moser et al., 2006). The reduced leaf area limited the leaf source activity under water deficits prior to silking during the grain-filling period. This resulted in reduced kernel weight. The water deficit in the vegetative growth period could impact late ovary development; this effect could still be detected during grain-filling (Figure 5). This result has also been confirmed in many other crops (Yang et al., 2009; Hasan et al., 2011; Castillo et al., 2017). Therefore, ovary size is likely a reference indicator for the final kernel weight of maize under abiotic stresses. Hypothetically, ovary sensitivity to water stress is likely associated with the drought tolerance of different maize varieties. This is a topic that should be further examined in future studies.
Author Contributions
PW and YL designed the project. BZ and YL conducted the field experiment and collected data. YL made the sections and took microphotos and analyzed data. SH and YL drew the tables and figures. SH, HT, and YL drafted the manuscript.
Funding
Funding was provided by the National Science Foundation of China (Grant No. 31571592) and the National Key Research and Development Program of China (2017YFD0300410).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01326/full#supplementary-material
References
Adams, W. W., Muller, O., Cohu, C. M., and Demmig-Adams, B. (2013). May photoinhibition be a consequence, rather than a cause, of limited plant productivity? Photosynth. Res. 117, 31–44. doi: 10.1007/s11120-013-9849-7
Ali, Q., and Ashraf, M. (2011). Induction of drought tolerance in maize (Zea mays L.) due to exogenous application of trehalose: growth, photosynthesis, water relations and oxidative defence mechanism. J. Agron. Crop Sci. 197, 258–271. doi: 10.1111/j.1439-037X.2010.00463.x
Allison, J. C. S., and Watson, D. J. (1966). The production and distribution of dry matter in maize after flowering. Ann. Bot. 30, 365–381. doi: 10.1093/oxfordjournals.aob.a084082
Andrade, F. H., Echarte, L., Rizzalli, R., Della Maggiora, A., and Casanovas, M. (2002). Kernel number prediction in maize under nitrogen or water stress. Crop Sci. 42, 1173–1179. doi: 10.2135/cropsci2002.1173
Aylor, D. E. (2004). Survival of maize (Zea mays) pollen exposed in the atmosphere. Agric. For. Meteorol. 123, 125–133. doi: 10.1016/j.agrformet.2003.12.007
Baker, R. F., Leach, K. A., Boyer, N. R., Swyers, M. J., Benitez-Alfonso, Y., Skopelitis, T., et al. (2016). Sucrose transporter ZmSut1 expression and localization uncover new insights into sucrose phloem loading. Plant Physiol. 172, 1876–1898. doi: 10.1104/pp.16.00884
Bolaños, J., and Edmeades, G. O. (1990). Importance of flowering interval on breeding for drought resistance in tropical maize. Agron. Mesoam. 1, 45–50.
Bolaños, J., and Edmeades, G. O. (1996). The importance of the anthesis-silking interval in breeding for drought tolerance in tropical maize. Field Crops Res. 48, 65–80. doi: 10.1016/0378-4290(96)00036-6
Borrás, L., Curá, J. A., and Otegui, M. E. (2002). Maize kernel composition and post-flowering source-sink ratio. Crop Sci. 2002, 781–790. doi: 10.2135/cropsci2002.7810
Borrás, L., and Gambín, B. L. (2010). Trait dissection of maize kernel weight: towards integrating hierarchical scales using a plant growth approach. Field Crops Res. 118, 1–12. doi: 10.1016/j.fcr.2010.04.010
Borrás, L., Maddonni, G. A., and Otegui, M. E. (2003a). Leaf senescence in maize hybrids: plant population, row spacing and kernel set effects. Field Crops Res. 82, 13–26. doi: 10.1016/S0378-4290(03)00002-9
Borrás, L., Westgate, M. E., and Otegui, M. E. (2003b). Control of kernel weight and kernel water relations by post-flowering source–sink ratio in maize. Ann. Bot. 91, 857–867.
Borrás, L., and Otegui, M. E. (2001). Maize kernel weight response to postflowering source–sink ratio. Crop Sci. 41, 1816–1822. doi: 10.2135/cropsci2001.1816
Borrás, L., and Westgate, M. E. (2006). Predicting maize kernel sink capacity early in development. Field Crops Res. 95, 223–233. doi: 10.1016/j.fcr.2005.03.001
Boyle, M. G., Boyer, J. S., and Morgan, P. W. (1991). Stem infusion of liquid culture medium prevents reproductive failure of maize at low water potential. Crop Sci. 31, 1246–1252. doi: 10.2135/cropsci1991.0011183X003100050033x
Brevedan, R. E., and Egli, D. B. (2003). Short periods of water stress during seed filling, leaf senescence, and yield of soybean. Crop Sci. 43, 2083–2088. doi: 10.2135/cropsci2003.2083
Cakir, R. (2004). Effect of water stress at different development stages on vegetative and reproductive growth of corn. Field Crops Res. 89, 1–6. doi: 10.1016/j.fcr.2004.01.005
Cárcova, J., and Otegui, M. E. (2001). Ear temperature and pollination timing effects on maize kernel set. Crop Sci. 41, 1809–1815. doi: 10.2135/cropsci2001.1809
Castillo, F. M., Vásquez, S. C., and Calderini, D. F. (2017). Does the pre-flowering period determine the potential grain weight of sunflower? Field Crops Res. 212, 23–33. doi: 10.1016/j.fcr.2017.06.029
Chaves, M. M., Pereira, J. S., Maroco, J., Rodrigues, M. L., Ricardo, C. P., Osório, M. L., et al. (2002). How plants cope with water stress in the field? Photosynthesis and growth. Ann. Bot. 89, 907–916. doi: 10.1093/aob/mcf105
China Statistical Yearbook (2016). Available at: http://www.stats.gov.cn/tjsj/ndsj/2016/indexch.htm
Escobar-Gutiérrez, A. J., and Combe, L. (2012). Senescence in field-grown maize: from flowering to harvest. Field Crops Res. 134, 47–58. doi: 10.1016/j.fcr.2012.04.013
Feng, H. J. (2014). Effect of Plant Population on Microstructure and Function of Vascular Bundle of Summer Maize (Zea may L.). Doctoral dissertation, Shandong Agricultural University, Shandong. Available at: http://cdmd.cnki.com.cn/Article/CDMD-10434-1014347217.htm
Fuad-Hassan, A. V., Tardieu, F., and Turc, O. (2008). Drought-induced changes in anthesis-silking interval are related to silk expansion: a spatio-temporal growth analysis in maize plants subjected to soil water deficit. Plant Cell Environ. 31, 1349–1360. doi: 10.1111/j.1365-3040.2008.01839.x
Gao, Z., Liang, X. G., Zhang, L., Lin, S., Zhao, X., Zhou, L. L., et al. (2017). Spraying exogenous 6-benzyladenine and brassinolide at tasseling increases maize yield by enhancing source and sink capacity. Field Crops Res. 211, 1–9. doi: 10.1016/j.fcr.2017.05.027
Grzesiak, M. T., Grzesiak, S., and Skoczowski, A. (2006). Changes of leaf water potential and gas exchange during and after drought in triticale and maize genotypes differing in drought tolerance. Photosynthetica 44, 561–568. doi: 10.1007/s11099-006-0072-z
Hasan, A. K., Herrera, J., Lizana, C., and Calderini, D. F. (2011). Carpel weight, grain length and stabilized grain water content are physiological drivers of grain weight determination of wheat. Field Crops Res. 123, 241–247. doi: 10.1016/j.fcr.2011.05.019
He, Q. P., Dong, S. T., and Gao, R. Q. (2005). Relationship between development of spike vascular bundle and sink capacity of ear and kernel in maize (Zea mays L.). 2005. Acta Agron. Sin. 31, 995–1000.
He, Q. P., Dong, S. T., and Gao, R. Q. (2007). Comparison of ear vascular bundles in different maize cultivars. Acta Agron. Sin. 33, 1187–1196.
Hu, Y., Burucs, Z., and Schmidhalter, U. (2008). Effect of foliar fertilization application on the growth and mineral nutrient content of maize seedlings under drought and salinity. Soil Sci. Plant Nutr. 54, 133–141. doi: 10.1111/j.1747-0765.2007.00224.x
Jones, R. J., and Simmons, S. R. (1983). Effect of altered source-sink ratio on growth of maize kernels. Crop Sci. 23, 129–134. doi: 10.2135/cropsci1983.0011183X002300010038x
Jurgens, S. K., Johnson, R. R., and Boyer, J. S. (1978). Dry matter production and translocation in maize subjected to drought during grain fill. Agron. J. 70, 678–682. doi: 10.2134/agronj1978.00021962007000040036x
Leroux, B. M., Goodyke, A. J., Schumacher, K. I., Abbott, C. P., Clore, A. M., Yadegari, R., et al. (2014). Maize early endosperm growth and development: from fertilization through cell type differentiation. Am. J. Bot. 101, 1259–1274. doi: 10.3732/ajb.1400083
Lesk, C., Rowhani, P., and Ramankutty, N. (2016). Influence of extreme weather disasters on global crop production. Nature 529, 84–87. doi: 10.1038/nature16467
Liu, J. X., Liao, D. Q., Oane, R., Estenor, L., Yang, X. E., Li, Z. C., et al. (2006). Genetic variation in the sensitivity of anther dehiscence to drought stress in rice. Field Crop Res. 97, 87–100. doi: 10.1016/j.fcr.2005.08.019
Lobell, D. B., Schlenker, W., and Costa-Roberts, J. (2011). Climate trends and global crop production since 1980. Science 333, 616–620. doi: 10.1126/science.1204531
Mansouri-Far, C., Sanavy, S. A., and Saberali, S. F. (2010). Maize yield response to deficit irrigation during low-sensitive growth stages and nitrogen rate under semi-arid climatic conditions. Agric. Water Manag. 97, 12–22. doi: 10.1016/j.agwat.2009.08.003
McLaughlin, J. E., and Boyer, J. S. (2004a). Glucose localization in maize ovaries when kernel number decreases at low water potential and sucrose is fed to the stems. Ann. Bot. 94, 75–86.
McLaughlin, J. E., and Boyer, J. S. (2004b). Sugar-responsive gene expression, invertase activity, and senescence in aborting maize ovaries at low water potentials. Ann. Bot. 94, 675–689.
Moser, S. B., Feil, B., Jampatong, S., and Stamp, P. (2006). Effects of pre-anthesis drought, nitrogen fertilizer rate, and variety on grain yield, yield components, and harvest index of tropical maize. Agric. Water Manag. 81, 41–58. doi: 10.1016/j.agwat.2005.04.005
NeSmith, D. S., and Ritchie, J. T. (1992a). Short-and long-term responses of corn to a pre-anthesis soil water deficit. Agron. J. 84, 107–113. doi: 10.2134/agronj1992.00021962008400010021x
NeSmith, D. S., and Ritchie, J. T. (1992b). Maize (Zea mays L.) response to a severe soil water-deficit during grain-filling. Field Crops Res. 29, 23–35. doi: 10.1016/0378-4290(92)90073-I
Ning, P., Fritschi, F. B., and Li, C. (2017). Temporal dynamics of post-silking nitrogen fluxes and their effects on grain yield in maize under low to high nitrogen inputs. Field Crops Res. 204, 249–259. doi: 10.1016/j.fcr.2017.01.022
Ober, E. S., Setter, T. L., Madison, J. T., Thompson, J. F., and Shapiro, P. S. (1991). Influence of water deficit on maize endosperm development enzyme activities and RNA transcripts of starch and zein synthesis, abscisic acid, and cell division. Plant Physiol. 97, 154–164. doi: 10.1104/pp.97.1.154
Otegui, M. E., Andrade, F. H., and Suero, E. E. (1995). Growth, water use, and kernel abortion of maize subjected to drought at silking. Field Crops Res. 40, 87–94. doi: 10.1016/0378-4290(94)00093-R
Oury, V., Tardieu, F., and Turc, O. (2016). Ovary apical abortion under water deficit is caused by changes in sequential development of ovaries and in silk growth rate in maize. Plant Physiol. 171, 986–996. doi: 10.1104/pp.15.00268
Pearce, R. B., Mock, J. J., and Bailey, T. B. (1975). Rapid method for estimating leaf area per plant in maize. Crop Sci. 15, 691–694. doi: 10.2135/cropsci1975.0011183X001500050023x
Prochazkova, D., Sairam, R. K., Srivastava, G. C., and Singh, D. V. (2001). Oxidative stress and antioxidant activity as the basis of senescence in maize leaves. Plant Sci. 161, 765–771. doi: 10.1016/S0168-9452(01)00462-9
Raynaud, D. (2010). Why do Diffusion Data not Fit the Logistic Model? A Note on Network Discreteness, Heterogeneity and Anisotropy. Vienna: Springer, 217–218.
Rossi, M., Bermudez, L., and Carrari, F. (2015). Crop yield: challenges from a metabolic perspective. Curr. Opin. Plant Biol. 25, 79–89. doi: 10.1016/j.pbi.2015.05.004
Rossini, M. A., Maddonni, G. A., and Otegui, M. E. (2016). Multiple abiotic stresses on maize grain yield determination: additive vs multiplicative effects. Field Crops Res. 198, 280–289. doi: 10.1016/j.fcr.2016.07.004
Ruan, Y. L., Jin, Y., Yang, Y. J., Li, G. J., and Boyer, J. S. (2010). Sugar input, metabolism, and signaling mediated by invertase: roles in development, yield potential, and response to drought and heat. Mol. Plant 3, 942–955. doi: 10.1093/mp/ssq044
Sade, N., Rubio-Wilhelmi, M. D., Umnajkitikorn, K., and Blumwald, E. (2017). Stress-induced senescence and plant tolerance to abiotic stress. J. Exp. Bot. 69, 845–853. doi: 10.1093/jxb/erx235
Schussler, J. R., and Westgate, M. E. (1991). Maize kernel set at low water potential: I. Sensitivity to reduced assimilates during early kernel growth. Crop Sci. 31, 1189–1195. doi: 10.2135/cropsci1991.0011183X003100050023x
Schussler, J. R., and Westgate, M. E. (1995). Assimilate flux determines kernel set at low water potential in maize. Crop Sci. 35, 1074–1080. doi: 10.2135/cropsci1995.0011183X003500040026x
Setter, T. L., and Flannigan, B. A. (2001). Water deficit inhibits cell division and expression of transcripts involved in cell proliferation and endoreduplication in maize endosperm. J. Exp. Bot. 52, 1401–1408. doi: 10.1093/jexbot/52.360.1401
Setter, T. L., Flannigan, B. A., and Melkonian, J. (2001). Loss of kernel set due to water deficit and shade in maize: carbohydrate supplies, abscisic acid, and cytokinins. Crop Sci. 41, 1530–1540. doi: 10.2135/cropsci2001.4151530x
Subedi, K. D., and Ma, B. L. (2005). Ear position, leaf area, and contribution of individual leaves to grain yield in conventional and leafy maize hybrids. Crop Sci. 45, 2246–2257. doi: 10.2135/cropsci2004.0653
Tanguilig, V. C., Yambao, E. B., O’toole, J. C., and De Datta, S. K. (1987). Water stress effects on leaf elongation, leaf water potential, transpiration, and nutrient uptake of rice, maize, and soybean. Plant Soil 103, 155–168. doi: 10.1007/BF02370385
Tari, A. F. (2016). The effects of different deficit irrigation strategies on yield, quality, and water-use efficiencies of wheat under semi-arid conditions. Agr. Water Manage. 167, 1–10. doi: 10.1016/j.agwat.2015.12.023
Tezara, W., Mitchell, V. J., Driscoll, S. D., and Lawlor, D. W. (1999). Water stress inhibits plant photosynthesis by decreasing coupling factor and ATP. Nature 401, 914–917. doi: 10.1038/44842
Uribelarrea, M., Cárcova, J., Otegui, M. E., and Westgate, M. E. (2002). Pollen production, pollination dynamics, and kernel set in maize. Crop Sci. 42, 1910–1918. doi: 10.2135/cropsci2002.1910
Wang, N. (2011). Different Density-Tolerant Maize Varieties Vascular Bundle Characters and Its Source-Sink System. Doctoral dissertation, Shenyang Agricultural University, Shenyang. Available at: http://cdmd.cnki.com.cn/Article/CDMD-10157-1012302251.htm
Wang, Y. H., Zhang, Y., Wang, X. Y., Meng, Z. Y., and Zhang, S. L. (2010). Analysis and curve fitting on filling characteristic of summer corn in Huang-huai-hai area. Agric. Res. Arid Areas 28, 265–270.
Xu, C., Gao, Y., Tian, B., Ren, J., Meng, Q., and Wang, P. (2017). Effects of EDAH, a novel plant growth regulator, on mechanical strength, stalk vascular bundles and grain yield of summer maize at high densities. Field Crops Res. 200, 71–79. doi: 10.1016/j.fcr.2016.10.011
Yang, Z., Van Oosterom, E. J., Jordan, D. R., and Hammer, G. L. (2009). Pre-anthesis ovary development determines genotypic differences in potential kernel weight in sorghum. J. Exp. Bot. 60, 1399–1408. doi: 10.1093/jxb/erp019
Yu, L. X., and Setter, T. L. (2003). Comparative transcriptional profiling of placenta and endosperm in developing maize kernels in response to water deficit. Plant Physiol. 131, 568–582. doi: 10.1104/pp.014365
Yu, S. M., Lo, S. F., and Ho, T. H. (2015). Source–sink communication: regulated by hormone, nutrient, and stress cross-signaling. Trends Plant Sci. 20, 844–857. doi: 10.1016/j.tplants.2015.10.009
Zinselmeier, C., Jeong, B. R., and Boyer, J. S. (1999). Starch and the control of kernel number in maize at low water potentials. Plant Physiol. 121, 25–36. doi: 10.1104/pp.121.1.25
Keywords: corn, water deficit, vascular bundle, yield, leaf area
Citation: Li Y, Tao H, Zhang B, Huang S and Wang P (2018) Timing of Water Deficit Limits Maize Kernel Setting in Association With Changes in the Source-Flow-Sink Relationship. Front. Plant Sci. 9:1326. doi: 10.3389/fpls.2018.01326
Received: 25 March 2018; Accepted: 23 August 2018;
Published: 22 October 2018.
Edited by:
Weihong Luo, Nanjing Agricultural University, ChinaReviewed by:
Edmundo Acevedo, Universidad de Chile, ChileRobert John French, Department of Primary Industries and Regional Development of Western Australia (DPIRD), Australia
Copyright © 2018 Li, Tao, Zhang, Huang and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shoubing Huang, aHVhbmdzaG91YmluZ0BjYXUuZWR1LmNu Pu Wang, d2FuZ3B1QGNhdS5lZHUuY24=