Abstract

During sugarcane growth, the Early Elongation stage is critical to cane yield formation. In this study, parameters of 17 sugarcane varieties were determined at the Early Elongation stage using CI-301 photosynthesis measuring system and CI-100 digital plant canopy imager. The data analysis showed highly significant differences in leaf area index (LAI), mean foliage inclination angle (MFIA), transmission coefficient for diffused light penetration (TD), transmission coefficient for solar beam radiation penetration (TR), leaf distribution (LD), net photosynthetic rate (PN), transpiration rate ( ), and stomatal conductance (GS) among sugarcane varieties. Based on the photosynthetic or canopy parameters, the 17 sugarcane varieties were classified into four categories. Through the factor analysis, nine parameters were represented by three principal factors, of which the cumulative rate of variance contributions reached 85.77%. A regression for sugarcane yield, with relative error of yield fitting less than 0.05, was successfully established: sugarcane yield = −27.19 − 1.69 × PN + 0.17 ×   + 90.43 × LAI − 408.81 × LD + 0.0015 × NSH + 101.38 ×   ( ). This study helps provide a theoretical basis and technical guidance for the screening of new sugarcane varieties with high net photosynthetic rate and ideal canopy structure.

1. Introduction

Photosynthesis is the basis of yield formation. It is generally recognized that about 90%–95% of crop yield is formed from the assimilated carbons. Therefore, how to improve crops’ photosynthetic capacity is one of the most important targets in the studies of genetics and breeding [14]. The canopy structure of a crop should be optimized so that an increased proportion of light can reach the leaves at the base of the plants [5, 6]. Studies on photosynthetic characteristics of crops such as Zea mays [79], Oryza sativa [1012], and Glycine max [13, 14] have been extensively performed. Also, inheritance of photosynthetic characteristics [1517], diurnal variation of net photosynthetic rate [18, 19], and seasonal variation [20, 21] have been reported in the separated sugarcane seedling populations. Indices such as shape, size, quantity, and spatial distribution of leaves are directly related to the light environment and light utilization rate of the population. They are the important factors affecting the light distribution in crop population and photosynthesis [22].

Optimizing a reasonable population structure, improving the light distribution in the population, and increasing the rate of light utilization are all effective ways in obtaining high crop yields. Leaf area index (LAI) determines canopy light interception, which affects the photosynthetic rate of a crop population. When LAI is at the optimal level, canopy light interception reaches the highest level for improved photosynthetic capacity and potential yield increase [23]. The relations between sugarcane leaf morphology and cane yield and between sucrose content and transpiration rate of leaves have been studied [2427]. However, studies on the canopy characteristics of sugarcane were rarely reported [28, 29], although there are many yield-related traits and indices controlling photosynthesis and canopy structure of sugarcane. The interactions between these indices are complicated due to the influences of genetic characteristics, environmental factors, and sample selections. The factor analysis method allows the categorization of a large number of correlated agronomic traits into several groups based on principal factors. It has been widely used in studies of crop germplasm resources. In sugarcane, a correlation analysis of seasonal variation of the canopy structure parameters to yield-related traits was conducted [21]; however, there has been no report on the application of factor analysis to the classification of photosynthetic and canopy parameters and their effects on crop yield in sugarcane [30, 31].

During sugarcane growth, Early Elongation stage is a critical period for cane yield formation. During this period, an appropriate number of stalks per unit area are controlled through timely cultivation to achieve high yield [25]. In this study, we aimed to investigate the effects of photosynthetic and canopy characteristics of 17 sugarcane varieties at the Early Elongation stage on cane yield. The results of this study would provide the scientific basis for breeding of sugarcane varieties with high photosynthetic efficiency and for improving cultivation technique. Hierarchical cluster analysis was performed by the numerical classification method on the photosynthetic and canopy parameters of the involved sugarcane varieties. The relationship and the stepwise regression analysis between canopy and photosynthetic characteristics and cane yield were also conducted, in order to screen new sugarcane varieties with high net photosynthetic rate (PN) and ideal canopy structure.

2. Materials and Methods

2.1. Sugarcane Varieties

Seventeen sugarcane varieties were tested, namely, FN94-0403, FN94-0744, FN95-1726, FN96-0907, FN98-10100, ROC10, GT94-116, GT95-118, GT96-211, GT96-44, GT97-18, MT70-611, YT92-1287, YT96-107, YT96-794, YT96-835, and YT96-86.

2.2. Field Experiment Design and Cane Yield Traits

Field experiment was conducted on the Experimental Farm of the Key Laboratory of Sugarcane Biology and Genetic Breeding, Ministry of Agriculture/Fujian Agriculture and Forestry University located at Jianxin, Cangshan District, Fuzhou, Fujian (longitude: 119.23E, latitude: 26.08N). A randomized block design with three replications was adopted. Plot area was 33 m2, with three rows at 1.1 m row spacing. Planting density was 45,000 two-bud setts ha−1. The preceding crop was sugarcane, and the soil type was sandy loam soil. The nutrients of the arable layer before sowing included organic matter of 10.5 g/kg, total nitrogen of 0.91 g/kg, alkaline hydrolysis nitrogen of 90.01 mg/kg, available phosphorus of 110.4 mg/kg, and rapidly available potassium of 369.5 mg/kg. The rows were covered by plastic films after planting. Planting was performed in a scheme of double buds and double rows. Herbicides were applied before planting and plastic film mulching. Weed control was repeatedly performed as needed. The basal fertilizer Calcium Super Phosphate at 750 kg ha−1 was applied during planting. During growing season, 975 kg ha−1 and 600 kg ha−1 KCl were applied for two times as topdressing. Field management level was only slightly more intensive than that of the local standard operations on cultivation, fertilization, irrigation, and pest control. Each application was completed covering all the experimental plots on the same day.

In Fuzhou, sugarcane elongation period lasts from late June to early November, while Early Elongation stage commences from late June and ends in early July. Measurements of photosynthetic characteristics and canopy parameters were conducted in early July when sugarcane seedlings reached an average height of 80.9 cm. Data on yield-related traits, including plant height ( ), stalk diameter ( ), single stalk weight (SSW), and the number of effective stalks per hectare (NSH), were collected before harvesting. All the stalks in the middle row of each plot were cut and weighed. The sampled sugarcane areas were measured. The number of millable stalks within the sampling area was also counted. Single stalk weight and cane yield were calculated by the following formulae [25]:Single stalk weight = (plant height × stalk diameter2 × 0.785)/1000; Cane yield = single stalk weight × the number of effective stalks per hectare.

2.3. Determination of Photosynthetic Parameters

Photosynthetic parameter data were collected between 8:30–11:30 AM on a sunny morning in early July, using a CI-301 photosynthesis system (CID Co., Ltd, Vancouver, WA, USA) under natural light conditions. The data included net photosynthetic rate (PN, μmolm−2 s−1), transpiration rate ( , mmolm−2 s−1), and stomatal conductivity (GS, mmolm−2 s−1). The PN equals the rate of photosynthetic CO2 fixation minus the rate of CO2 loss during respiration. The is the amount of evaporation per unit time from a leaf surface. The GS is the rate of CO2 entering stomata [17, 22]. The measurements were taken three times on each variety following the protocol described earlier with minor modifications [17, 22, 31]. The first youngest fully expanded (+1) leaves at the top canopy were measured reciprocally at the middle to upper section excluding the midrib. The direction of leaf chamber was adjusted towards sunlight to ensure that measurements were done under a uniform light intensity. A total of 18 plants were measured for each variety under the following natural conditions: light intensity (1,781.17 ± 103.21 μmolm−2 s−1), ambient temperature ( ) (28.07 ± 2.44°C), leaf temperature (30.14 ± 2.26°C), ambient relative humidity (RH) (47.71 ± 4.88%), and ambient CO2 concentration ( ) (331.02 ± 25.35 μL L−1).

Several parameters were recorded simultaneously, including effective photosynthetic radiation (PAR), relative humidity (RH), ambient temperature ( ), ambient CO2 concentration ( ), stomatal conductivity (GS), evaporation rate ( ), and between cell CO2 concentrations ( ). Because all determinations were done under the natural sunlight, only those measurements were kept when natural light intensity is at 1,600 μmol/M2 Sec, while measurements taken under 1,600 μmol/M2 Sec light intensity were discarded. CO2 concentrations were collected from the upper level of canopy after passing through the dual-bottle apparatus and were slightly lower than those collected above the ground level. This is due to the fact that sugarcane is a highly efficient C4 crop. Only mean values and standard deviations were shown in the text.

2.4. Determination of Canopy Parameters

A fisheye lens of the CI-100 digital plant canopy imager (CID Co., Ltd, Vancouver, WA, USA) was attached to an observation rod and placed in the center of the rows in a cloudless twilight afternoon with a little sunshine in early July. The rod was adjusted to the horizontal direction for taking photos when there was no shadow or any other external influences. Fifteen testing points were selected for each variety and five images were taken from each plot. LD was represented by distribution frequency of the leaf within each azimuth. The canopy indicators such as leaf area index (LAI), mean foliage inclination angle (MFIA), transmission coefficient for diffuse light penetration (TD), transmission coefficient for solar beam radiation penetration (TR), and leaf distribution (LD) were calculated by Plant Canopy Analysis software provided by CID Company (CID Co., Ltd, Vancouver, WA, USA) [21].

2.5. Statistical Methods

Variance analysis, cluster analysis, factor analysis, and regression analysis were performed using DPS software (Zhejiang University, Hangzhou, China) [30]. After standardization of original data, the distance coefficient was determined as the chi-square distance and cluster analysis by Ward’s method [30]. Varimax rotation was used to change the coordinates used in principal component analysis and factor analysis [21, 30]. A stepwise regression analysis was also conducted for the photosynthetic and canopy parameters and their correlation with cane yield [21, 30].

3. Results

3.1. Estimated Cane Yield Trait Data

Estimated cane yield trait data are shown in Table 1.

3.2. Differences in Photosynthetic Parameters among 17 Sugarcane Varieties

There were extremely significant differences in the mean PN, , and GS among the 17 sugarcane varieties (Table 2). Based on PN, the 17 sugarcane varieties were classified into four categories. Category I (high PN) included five varieties, namely, FN98-10100, GT94-116, GT96-211, YT96-794, and YT96-86, with a mean PN value of 37.09 μmolm−2 s−1. Category II included four varieties (FN94-0403, FN95-1726, FN96-0907, and YT96-107) with a mean PN value of 35.68 μmolm−2 s−1. Category III included six varieties (GT95-118, GT97-18, MT70-611, ROC10, YT92-1287, and YT96-835) with a mean PN value of 34.24 μmolm−2 s−1. Category IV included two varieties, FN94-0744 and GT96-44, with a mean PN value of 34.24 μmolm−2 s−1. The 17 sugarcane varieties could be classified into three groups based on the values. Group I (high ) included 12 varieties, namely, MT70-611, FN94-0744, GT96-211, ROC10, YT96-107, YT96-835, FN94-0403, FN95-1726, FN96-0907, FN98-10100, YT96-794, and YT96-86, with a mean value of 4.11 mmolm−2 s−1. Group II (mid- ) included four middle varieties (GT94-116, GT95-118, GT97-18, and YT92-1287) with a mean value of 3.65 mmolm−2 s−1. The variety GT96-44 alone represented Group III (low ) with a value of only 3.10 mmolm−2 s−1. The 17 sugarcane varieties were classifiable into four categories based on GS. Category I (high GS) included five varieties (GT96-211, ROC10, YT96-107, YT96-794, and YT96-86) with a mean GS value of 154.85 mmolm−2 s−1. Category II (mid-GS) included six varieties (FN94-0403, FN95-1726, FN96-0907, FN98-10100, MT70-611, and YT96-835) with a mean GS value of 142.4 mmolm−2 s−1. Category III (low GS) included five varieties (FN94-0744, GT94-116, GT95-118, GT97-18, and YT92-1287), with a mean GS value of 128.46 mmolm−2 s−1. The variety GT96-44 alone formed Category IV with the lowest GS value of 109.7 mmolm−2 s−1.

The result of cluster analysis showed that these 17 varieties could be grouped into four clusters (Table 3, Figure 1). Cluster I included seven varieties, namely, FN94-0403, FN95-1726, FN96-0907, FN98-10100, GT96-211, YT96-794, and YT96-86, which shared relatively high PN, , and GS values. Cluster II included five varieties, that is, FN94-0744, GT94-116, GT95-118, GT97-18, and YT92-1287. The PN values of these varieties were at middle levels, while all had high and GS values. Cluster III included four varieties (MT70-611, ROC10, YT96-107, and YT96-835) with middle PN, low , and low GS values. Cluster IV only had one variety, GT96-44, with low values of PN, , and GS.

3.3. Differences in the Canopy Parameters among 17 Sugarcane Varieties

There were highly significant differences in LAI, MFIA, TD, TR, and LD among the 17 sugarcane varieties at the Early Elongation stage (Table 4). According to LAI, 17 sugarcane varieties were classified into three categories. Category I included six large LAI varieties FN94-0403, GT94-116, FN96-0907, FN98-10100, GT96-44, and MT70-611, with a mean value 1.35. Category II included 4 middle LAI varieties FN94-0744, GT96-211, GT97-18, and YT92-1287, with a mean value 1.155. Category III included seven low LAI varieties FN95-1726, GT95-118, ROC10, YT96-107, YT96-794, YT96-835, and YT96-86, with a mean value 1.033.

Hierarchical cluster analysis was performed by merging all canopy parameters’ data, that is, LAI, MFIA, TD, TR, and LD, and the output was shown in Figure 2. The 17 varieties again were grouped into four clusters (Table 5, Figure 2). Cluster I included two varieties, FN94-0403 and FN98-10100, with large mean LAI and relatively small mean MFIA and TD values. Cluster II included four varieties, namely, FN96-0907, GT94-116, GT96-44, and MT70-611. Cluster III contained five varieties, that is, FN94-0744, GT96-211, GT97-18, YT92-1287, and YT96-86, with medium LAI, MFIA, TD, TR, and LD values. Cluster IV had six varieties, namely, FN95-1726, ROC10, YT96-107, YT96-835, GT95-118, and YT96-794, with small LAI values but relatively large MFIA and TD values.

3.4. Relationships between Photosynthetic and Canopy Parameters and Cane Yield

Cane yield-related trait values including photosynthetic parameters (PN, , and GS), canopy parameters (LAI, MFIA, TD, TR, and LD), and other traits (NSH, SSW) corresponding to the first three principal factors, joint degree (JD), and special variances (SV) are listed in Table 6. The cumulative contribution rates of the characteristic values of the first three principal factors reached 85.77%. Table 6 also indicates that the essential information of the first three principal factors can reflect most of the information of these photosynthetic gas exchange parameters, canopy parameters, and yield-related traits. The nine traits had much in common with a high degree of freedom. It also revealed that these three principal factors represented the nine traits well.

The factor loading matrices after varimax rotation in three principal factors were shown in Table 7. By comparison between the factor loading matrices of postvarimax rotation and the ones of prevarimax rotation, the load values of important variables in the principal factors were increased significantly. The biological significance of the principal factors after the varimax rotation also became more significant. As can be seen from Table 7, the large load values for canopy parameters, LAI, MFIA, TD, TR, and LD, indicated that these values played a major role in the first principal factor. The canopy parameters also had a close relationship with cane yield. On the other hand, the load values of the photosynthetic parameters PN, , and GS were relatively large in the second principal factor, while the load values of SSW and NSH were relatively large in the third principal factor. Collectively, these factors were entitled as the yield component factors.

The correlation analyses data listed in Table 8 indicated the following: (1) PN had a positive correlation with and GS; (2) LAI had a negative correlation with MFIA, TD, and TR, but a significantly positive correlation with NSH; (3) cane yield had a significantly negative correlation with MFIA, TD, and TR, but a significantly positive correlation with LD; and (4) SSW had significantly negative correlation with NSH, but a significantly positive correlation with   . Based on regression analysis, cane yield was predictable based on photosynthetic and canopy parameters. The regression equation is as follows: Cane yield = − 27.19 − 1.69 × PN + 0.17 ×   + 90.43 × LAI − 408.81 × LD + 0.0015 × NSH + 101.38 ×   ( ). The relative error of cane yield fitting was less than 5% for all the 17 sugarcane varieties. Parameters NSH and play a decisive role in cane yield formation with highly significant correlation coefficients of 0.902 and 0.94. Parameters LAI and LD also play a major role with significant partial correlation coefficients of 0.682 and 0.624. However, parameters PN and only influence the formation of cane yield slightly with nonsignificant partial correlation coefficients of −0.331 and 0.298. Based on these observations, we conclude that the Early Elongation (growth) stage is critical to the life cycle of sugarcane. It is during this elongation stage that NSH and start to form a canopy structure, which in turn determines the yield level of photosynthetic products. As was discussed in our previous report [21], cane yield formation could be influenced by many factors, which need to be adjusted according to canopy structure to ensure final cane yield formation.

4. Discussion and Conclusions

Regarding the genotypic differences in photosynthetic capacities among sugarcane varieties, the PN and GS could be affected greatly by the specific combining ability among sugarcane varieties [21]. However, the correlation coefficient between F1 hybrids and their parents reached a significant level in terms of general combining ability. Therefore, parental varieties with high photosynthetic efficiency should be selected in sugarcane breeding [22]. Sugarcane breeding practice demonstrates that PN and GS (water use efficiency) can be combined. Although the photosynthetic parameters PN, , and GS had high degrees of separations in F1 seedling populations [17], the high photosynthetic efficiency, once generated upon cross-recombination, can be immobilized by asexual reproduction, which is a unique feature of sugarcane production. This is why all the advanced varieties possess good PN, , and GS parameters [31].

In the present study, there were significant differences in the photosynthetic parameters among the 17 sugarcane varieties at the Early Elongation stage, including PN, , and GS, which are in accordance with the previous reports [17, 22, 31]. The sugarcane varieties were classified according to the three photosynthetic traits, respectively, and the comprehensive traits and natural types were well reflected. The clustering patterns were relatively stable, although the masking of some traits by others and some ambiguous intercategory differences did exist. Therefore, single and integrated traits should be combined in photosynthetic efficiency breeding to obtain a better product [8].

Leaf area determines the canopy light interception. The distribution of leaf area within canopy is an important index of crop canopy structure [32]. Previous studies showed that there were significant differences in canopy parameters LAI, MFIA, and TD at the seedling stage among different genotypes. Varieties with large LAI and small TD at the seedling stage had a higher number of leaves per unit space. A positive relationship exists between even leaf spatial distribution and light interception efficiency of the population and hence higher cane yield [33]. Leaf canopy structure directly affects the canopy light interception and light distribution at different layers. Researches on morphology and anatomy revealed that sugarcane populations with transparent canopy of short, narrow, thick, and straight leaves had a larger yield potential. In addition, LAI has a very close relationship with dry matter accumulation. Therefore, a proper increase of the population LAI can increase dry matter accumulation [3236]. In the present study, MFIA affects the receiving of solar radiation and light distribution in canopy. The ideal leaf population structure has the most effective leaf area with the characteristics of continuously changing the inclination angle distribution.

Canopy structure and light distribution of sugarcane at the Early Elongation stage play a crucial role in the accumulation and distribution of photosynthetic products, in the growth and development of population, and in the formation of final products [2730]. This has been clearly explained by the correlations between the main canopy parameters and cane yield in this study. Zhang et al. [35] showed that canopy MFIA of winter wheat decreased with the increase of nitrogen content of soil. The light distributions of two crop populations with the same LAI are different, if the spatial distributions of MFIA and leaves are different. They can lead to different photosynthetic rates of canopy and different amounts of dry matters were produced by populations per unit time [35]. However, further optimization is still required when the yield-related traits are indirectly selected based on the main canopy parameters since the yield of sugarcane has a negative correlation with quality-related traits in this study. What should also be stressed is that, for those varieties with large LAI at the Early Elongation stage, the LAI should be adjusted to be within a certain range at the middle to late stages. Mutual shading among leaves should be avoided, which weakens light intensity, reduces light interception, and decreases photosynthetic efficiency. Therefore, rapid and accurate diagnosis of canopy structure of sugarcane and a timely, reasonable regulation of population size play important roles in the cultivation of high-quality and high-yield sugarcane varieties [32, 33].

During the Early Elongation stage, the photosynthetic and canopy parameters play a crucial role in cane yield formation [21, 25]. In this study, highly significant differences in LAI, MFIA, TD, TR, LD, PN, , and GS were observed among 17 sugarcane varieties at the Early Elongation stage. However, during the factor analysis process, nine parameters were expressed by three principal factors, of which the cumulative variance contribution rate reached 85.77%. Both correlation analysis and factor analysis showed that, at the Early Elongation stage, there were close relationships between photosynthetic and canopy parameters and cane yield. Based on regression analysis for sugarcane yield, with relative error of yield fitting less than 5%, a cane yield prediction formula was established: Cane yield = − 27.19 − 1.69 × PN + 0.17 ×   + 90.43 × LAI − 408.81 × LD + 0.0015 × NSH + 101.38 ×   ( ), in which photosynthetic parameters, such as PN and , and canopy parameters, such as LAI and LD, are used as the indirect indices for selecting better yield-related traits at the Early Elongation stage. From all the above, this study should provide the theoretical foundation and technical guidance for identifying new sugarcane varieties with high photosynthetic rates and ideal canopy structures.

Conflict of Interests

The authors declare no conflict of interests.

Acknowledgments

This study was supported by the 948 Program on the Introduction of International Advanced Agricultural Science and Technique of Department of Agriculture (2014-S18), the Public Welfare Industrial Research Project of China (201003009-2), and the earmarked fund for the Modern Agriculture Technology of China (CARS-20). The authors especially thank Luguang Wu at School of Agriculture and Food Sciences, University of Queensland, Brisbane, Australia, for valuable editorial comments during paper preparation.