Seasonal Variation in Diurnal Photosynthesis and Chlorophyll Fluorescence of Four Genotypes of Cassava ( Manihot esculenta Crantz ) under Irrigation Conditions in a Tropical Savanna Climate

Photosynthesis performance during early vegetative growth is an important physiological trait determining yield of cassava, but limited information is currently available for the tropical savanna climate of Asia. Diurnal photosynthesis and chlorophyll fluorescence of the three-month-old plants of four commercial cassava genotypes (Rayong 9, RY9; Rayong 11, RY11; Kasetsart 50, KU50 and CMR38-125-77) grown under irrigation, were investigated in three seasons i.e., rainy, cool and hot. The mean daily net photosynthetic rate (Pn) across genotypes in the rainy season (11.75μmolCO2/m/s) was significantly lower than that in the cool season (14.60 μmolCO2/m/s). Daily mean Pn in the hot season was 14.32 μmolCO2/m/s. In the rainy season, maximum photochemical quantum yield of PSII (Fv/Fm) and effective quantum yield of PSII photochemistry (ΦPSII) were significantly higher than the other seasons, while electron transfer rate (ETR) and non-photochemical quenching (NPQ) were significantly lower. Genotypic variation was observed during the hot season in which RY11 had the highest and CMR38-125-77 the lowest mean daily Pn. The prominent mechanism to avoid damages from stress during afternoon in the hot season was to reduce leaf temperature by enhancing transpiration for RY11; to close stomata early for RY9, and to increase NPQ for CMR38-125-77.


Introduction
Cassava (Manihot esculenta Crantz) is one of the most significant crops for food security and the main food source for more than 0.8 billion people in Africa and Asia [1,2].Cassava is also used to produce starch for industrial applications including paper, textile, food and beverages, plywood, glue, animal feed and ethanol [3].Thailand is ranked as the world's largest exporter of cassava products, supplying around 67% of the global market with annual production of 31 million tons in 2016 [4].The country exports cassava products in the form of dry chips, pellets, and native and modified starch [5].Most of the cassava growing areas in Thailand are in the tropical savanna climate zone with a growing period, from planting to storage root harvesting, of eight to 12 months, which covers almost all three seasons in Thailand (hot season, March to May; rainy season, June to October, and cool season, November to February) [6].
Storage root yield of cassava varies considerably depending on cultivars, climate, growing conditions, crop management and planting time [7,8].While the highest recorded experimental yield potential was 80 to 90 t ha −1 under near optimal edaphic-climatic conditions [9,10], the average yield in Thailand was recorded at 18.83-24.15t ha −1 [11].This yield gap might be due to a lack of sufficient agro-advisory information about cassava genotypes suitable for planting in different seasons and appropriate management practices for the different growing seasons [12].
Crop yield is primarily determined by photosynthesis [13,14] and net photosynthetic rate (Pn) of cassava leaves was reported to have a significant correlation with storage root yield across environments [15].Related photosynthesis parameters such as internal CO2 concentration (Ci) was also significantly correlated with dry root yield of cassava [16].In addition, changes in chlorophyll (Chl) fluorescence parameters under different environmental conditions can be used as a rapid and sensitive measure of plant health status, photosynthetic competence as well as effects of stress on plants [17,18].Information on environmental effects in different seasons on photosynthetic performance of cassava genotypes is seriously lacking for the tropical savanna climate including Thailand.If seasonal variations in environmental conditions and responsive physiological parameters are precisely measured under specific field condition, it would help determine plantbased and/or environmental factors limiting photosynthetic performance as well as providing essential input data for effective crop growth and photosynthesis modeling [19].This information will be useful for providing a set of guidelines in improving management of cultural practice of cassava growing in different seasons.
The effect of weather variations during different periods of growing seasons in Thailand on growth rate and final yield of cassava has been reported [12].However, the physiological responses, particularly photosynthesis, of cassava genotypes in different growing seasons have not been evaluated in Thailand.Therefore, the objectives of this study were to investigate the effect of weather fluctuation in different seasons on diurnal photosynthesis and Chl fluorescence of the three-monthold plants of four commercially important cassava genotypes grown under irrigation to determine appropriate genotypes for planting in different seasons based on their growth and photosynthetic performance.

Site Description
The study site was at the experimental research station, Faculty of Agriculture, Khon Kaen University, Northeast Thailand (16°28′29.7″N, 102°48′37.3″E, altitude 195 m above sea level).Soils were of the Yasothon soil series; fine loamy, siliceous, and Oxic Paleustult [20].The climate is equatorial savanna with a dry winter [21].The seasonal mean air temperatures for the Northeastern part of Thailand for 30 years (1981-2010) were 27.6, 24.2 and 28.6 °C, with the mean minimum of 24.4,18.7 and 23.2 °C, and mean maximum of 32.6, 30.6 and 35.2 °C in rainy, cool and hot seasons, respectively.The average seasonal rainfalls, during 1981-2010, were 1103 mm, 76 mm, and 224 mm in rainy, cool and hot seasons, respectively [22].
The experimental study was conducted during June 2015 to April 2016.Environmental characteristics including ambient PAR (PARA), air temperature (Tair), air relative humidity (RHair) and rainfall were recorded every five minutes by an automatic weather station (Watchdog 2000, Spectrum Technologies Inc., Lincoln, NE, USA).Air vapor pressure deficit (VPDair) was calculated from air temperature, saturated vapor pressure and RH with the aid of Smithsonian Meteorological Tables prepared by Robert J. List, Smithsonian Institution Press, City of Washington [23] (Figure 1).The weather was under the influence of monsoon winds, specifically the southwest and northeast monsoon.In 2015-2016, the three seasons were divided as follows; rainy season or southwest monsoon season from June to October, cool season or northeast monsoon from November to February, and hot or pre-monsoon from March to May.The monthly mean environmental parameters in the experimental field are shown in Figure 1.

Plant Materials
Four cassava genotypes were selected for physiological studies, all were improved cultivars or line with high starch content suitable for industrial uses.Two plant types are recognized i.e., one nonbranching plant type (cv.Rayong 9, RY9), and one branching plant type (2 cvs.Rayong 11, RY11; and Kasetsart 50, KU50, and one line CMR-38-125-77).The four cassava genotypes were planted at three different planting dates (PD), namely, PD-Jun, PD-Nov and PD-Dec.PD-Jun plants were planted on 30 June 2015, PD-Nov on 10 November 2015 and PD-Dec on 15 December 2015.Stem cuttings (20 cm long) were planted with a plant spacing of 1 m × 1 m in the plot size of 5 m wide and 7 m long.All cassava plants were fertilized and kept well-watered (soil matric potential was maintained between 0 to −30 kPa) throughout the growing and measurement period.Nitrogen fertilizer ((NH)2SO4) formula 21-0-0 was applied at one month after planting at the rate of 46.9 kg ha −1 based on soil analysis and nutrient requirements for cassava [23,24].At two months after planting the compound fertilizer N-P2O5-K2O formula 15-0-18 was applied at a rate of 312.5 kg ha −1 (Chia tai company limited, Phranakhonsiayutthaya, Krung Thep Maha Nakhon, Thailand).The plants were irrigated by a minioverhead sprinkler irrigation system.Soil matric potential was continuously monitored by Tensiometer and water mark (Watchdog 1000, Spectrum Technologies Inc., Lincoln, NE, USA) at the soil depths of 20 and 40 cm.Water applications were administered whenever soil matric potential reduced to −30 kPa at 20 cm depth.Rainy season measurements of growth and photosynthesis were performed on PD-Jun plants when they were three months old (in October 2015).Cool and hot season measurements were performed on three-month-old PD-Nov and PD-Dec plants in February and April 2016, respectively.
Plant growth rate was measured, on the planting date and after three months of planting, as the rate of increase in stem height (cm d -1 ), and rate of leaf production (leaf d −1 ), on six randomly selected plants per genotype.Leaf area index (LAI) of the three-month-old plants was measured by gap fraction analysis using LI-191R line quantum sensor (Li-Cor Inc., Lincoln, NE, USA).Canopy structure was calculated according to the equation LAI = -(1/k)ln(Qb/Qa) where k is assumed to be close to 0.5, Qb is an average below-canopy PAR and Qa is an unobstructed PAR reading [25].

Physiological Measurement
Photosynthetic measurements were made on three-month-old plants on two days in each season.The three-month-old plants were selected for photosynthetic measurements because the most active vegetative growth occurs during the period three to six months after planting, and storage roots begin to be formed at this stage [26].Rainy season measurements were made on 25 September and 2 October 2015 on plants established on 30 June 2015 (PD-Jun).Cool season measurements were made on 24 and 26 February 2016 on plants established on 10 November 2015 (PD-Nov).Hot season measurements were made on 3 and 5 April 2016 on plants established on 15 December 2015.On each day, diurnal variation in gas exchange and chlorophyll fluorescence were measured at 2 h intervals (eight time points) from 04:30 to 18:30 local time.For each time point, the measurements were performed on the central lobe of a young fully expanded leaf of two randomly selected plants of each genotype.
Gas exchange parameters were evaluated using a portable gas exchange system, infrared gas analyzer (LI-6400xt, Li-Cor Inc., Lincoln, NE, USA) equipped with the standard 2 × 3 cm leaf chamber (6400-08 clear chamber bottom).Pn, Gs, Tr and Ci were measured under ambient environmental conditions.The instantaneous light use efficiency (LUE) was calculated from Pn/absorbed photosynthetically active radiation.Water-use efficiency (WUE) was calculated as the ratio of Pn/Tr.Ratio between internal CO2 concentration and ambient CO2 (Ci/Ca) was calculated to indicate non-stomatal limitation of photosynthesis.

Data and Statistical Analysis
Comparisons among planting dates of mean ambient photosynthetically active radiation (PARA), air temperature (Tair), air relative humidity (RHair) and air vapor pressure deficit (VPDair) were done by using paired t-test.Mean comparisons of growth, environmental, photosynthesis and chlorophyll fluorescence parameters were using one-way ANOVA for comparing means among genotypes within season and among seasons within each genotype, averages of multiple comparisons were determined by a Tukey's test under Sigmaplot version 11.0 software [29].The correlation analysis between photosynthesis (Pn, Gs, Tr, ΦPSII, ETR and NPQ) and environmental parameters including PARA, Tleaf, leaf-to-air vapor pressure deficit (VPDL) and relative humidity in the plant canopy (RHC) were analyzed in different seasons.Statistical significance was taken at p < 0.05 and p < 0.01 using MSTAT-C Version 1.42 software [30].All statistical analyses followed the procedure described by Gomez and Gomez [31].

Environment and Plant Growth
The environmental conditions under which the cassava plants of the three planting dates (PD-Jun, PD-Nov and PD-Dec) were growing were recorded from the date of planting until the plants were three-month-old as shown in Table S1 (see supplementary).The values of mean daily PAR for the three planting dates were similar ranging from 412-429 µmol/m 2 /s.Among the three planting dates, the PD-Jun plants experienced the highest daily minimum (23 °C) and highest daily mean temperature (28 °C).On the other hand, the PD-Nov and PD-Dec plants were exposed to the lowest daily minimum temperature of 19 °C.The mean daily temperatures during growth of PD-Nov and PD-Dec plants were 25 and 26 °C, respectively.The highest mean daily RH was found for PD-Jun (77%) while those recorded for PD-Nov and PD-Dec were at 57% and 52%, respectively.The lowest daily minimum RH (27-31%) was recorded in the cool season for the PD-Nov and PD-Dec plants.The PD-Jun plants which were growing in the rainy season were exposed to the lowest daily maximum (2.55 kPa) and daily mean (1.02 kPa) VPDair.The daily maximum VPDair values increased for PD-Nov (3.58 kPa) and PD-Dec (4.03 kPa) plants which were growing during the cool and coolto-hot seasons, respectively.The highest total rainfall was recorded for PD-Jun (614.3 mm) and the lowest (25.6 mm) for PD-Nov plants.In this experiment, irrigation was applied by a mini-overhead sprinkler system to maintain soil matric potential level at approximately -30 kPa or higher.The highest and lowest total irrigation were applied to PD-Dec (41.6 mm) and PD-Jun (7.2 mm) plants, respectively.
Cassava plant growth rates during the first three months after planting were measured as the rate of increase in plant height (cm d −1 ) and the rate of leaf production (leaf d −1 ) as shown in Figure 2 and Table S2.The rate of increase in plant height of plants established in June (PD-Jun plants) which were growing in the rainy season, averaged across four genotypes (1.28 cm d −1 ), was significantly higher (p < 0.001) than the rates for PD-Nov (0.75 cm d −1 ) and PD-Dec (0.72 cm d −1 ) plants which were planted and growing in the cool and cool-to-hot season, respectively (Figure 2A; Table S2).Rates of leaf production across four genotypes for all three planting dates were not significantly different.For PD-Jun plants, leaf production rates for CMR38-125-77 (0.85 leaf d −1 ), RY11 (0.79 leaf d −1 ) and RY9 (0.77 leaf d −1 ) were not significantly different while that of KU50 (0.53 leaf d −1 ) was significantly lower than the others (p < 0.001).For plants established in December (PD-Dec), CMR38-125-77 showed significantly higher leaf production rate (p = 0.004) than the other genotypes (Table S2).It was noted that KU50 plants established in June and November had significantly lower (p = 0.010) leaf production rate than those in December (Table S2).Mean leaf area index (LAI) of the three-month-old plants across genotypes for PD-Jun plants (3.32) was significantly higher than (p < 0.001) those planted in the cool season (2.32 and 2.61 for PD-Nov and PD-Dec, respectively) (Figure 2C; Table S2).However, significant difference (p < 0.001) among cultivars was observed only in PD-Dec in which CMR38-125-77 had the highest mean LAI of 3.49 followed by RY11 (2.51), RY9 (2.45) and KU50 (2.00).

Diurnal Variation in Environmental Conditions during Field Measurements of Photosynthesis
Diurnal patterns of environments (PARA, Tair, RHair and VPDair) measured during 04:30-18:30 on two sunny days of photosynthesis measurements in each of the three seasons, compared with leaf/canopy parameters (PARleaf, Tleaf, RHC and VPDL) were depicted in Figure 3.For each season, the diurnal minima, maxima and means recorded during 06:30 to 16:30 across two days of investigation were summarized in Table 1.Diurnal variation in ambient PAR (PARA) compared to incident PAR (PARleaf) were displayed in Figures 3A-B, C-D and E-F for measurements performed in rainy, cool and hot seasons, respectively.Similarly, the comparison between diurnal Tair and Tleaf were shown in Figures 3G-L.Diurnal RHair compared to RHC, and VPDair compared to VPDL were demonstrated in Figures 3M-R and Figure 3S-X, respectively.The daily minima, maxima and means of the environmental and leaf/canopy parameters recorded over the two days of photosynthesis measurements were shown in Table 1.The PARleaf which were recorded on the top canopy leaves during gas exchange measurement were similar to PARA over the canopy throughout the day except for the period during 12:00-15:00 when PARleaf tended to be lower than PARA (Figures 3A-F).There were seasonal differences (p < 0.001) in the mean daily PARleaf across genotypes, with the highest mean PARleaf in the hot season (930 µmol/m 2/ s) followed by the cool (863 µmol/m 2 /s) and the rainy (518 µmol/m 2 /s) seasons (Table 1).In the rainy season, Tleaf of all genotypes were higher than Tair on the first day (Figure 3G) probably due to high PARA (Figure 3A).However, on the second day (Figure 3H) Tair tended to be than Tleaf from 11.30 onwards.On the two days of observation in the cool season, Tleaf tended to be warmer than Tair throughout the day (Figure 3I-J).In the hot season, Tair tended to be warmer than Tleaf in the late afternoon from 14.30 to 18.30 (Figure 3K-L).The highest Tleaf was noted in the hot season reaching the maximum temperatures of 41.00-44.30°C during 10:30-12:30 (Table 1).It is worth noting that both Tair and Tleaf were significantly higher (p < 0.001) in the hot season compared to the others.The mean daily Tleaf across genotypes (34.48 °C) in the hot was significantly higher (p < 0.001) than those in the rainy (31.51 °C) and cool (27.70 °C) seasons (Table 1).Nevertheless, in any season, no significant differences in Tleaf were found among genotypes.As shown in Table 1, the highest mean relative humidity in the canopy (RHC) across genotypes and air relative humidity (RHair) were recorded in rainy season with the mean daily values of 61 and 77% which differed significantly (p < 0.05).Similarly, in the cool season, RHC (29%) was also significantly lower (p < 0.05) than RHair (42%).On the contrary, in the hot season, mean daily RHC (39%) was higher than RHair (37%), although the difference was not significant.As shown in Figure 3S-X, VPDL and VPDair were low (less than 2 kPa) in the morning, increased several fold during early afternoon (particularly in the cool and hot seasons), and declined slowly in the late afternoon.As shown in Table 1, the mean daily VPDair in the hot season (3.97 kPa) was significantly higher (p < 0.05) than the other seasons, and was approximately 4.0 and 1.8 fold higher than those in the rainy (0.99 kPa) and cool (2.19 kPa) seasons.The mean daily VPDL across genotypes in the cool (2.45 kPa) and hot (2.77 kPa) seasons were significantly higher than (p < 0.001) and approximately double that in the rainy season (1.07 kPa).It is worth noted that the mean daily VPDL was similar to VPDair in the rainy and cool seasons, but in the hot season VPDL was significantly lower (p < 0.05) than VPDair.

Diurnal Chl Fluorescence of Cassava Leaves
Diurnal patterns of Chl fluorescence parameters (ΦPSII, ETR and NPQ) of four cassava genotypes measured during 04:30-18:30 on two sunny days in each of the three seasons were depicted in Figure 4.For each season, the diurnal minima, maxima and means of ΦPSII, ETR and NPQ recorded during 06:30 to 16:30 across two days of investigation were summarized in Table 2.The values for Fv/Fm were obtained from the measurements in the dark at 04:30.The Fv/Fm values were high for all genotypes and in all seasons, although some significant differences were detected.The Fv/Fm means among genotypes were significantly different (p < 0.001) only in the rainy season (Table 2) with RY9 showing the highest Fv/Fm (0.866) which was significantly higher (p < 0.001) than RY11 (0.845) but not different from that of KU50 (0.847) and CMR38-125-77 (0.858).Seasonal variation in Fv/Fm was observed i.e. the means across genotypes in rainy (0.854) and hot (0.849) seasons were significantly higher (p < 0.001) than that in the cool (0.838) season (Table 2) Diurnal patterns of ΦPSII displayed the inverted bell-shaped curves, and in the cool and hot seasons, ΦPSII values during 6:30 to 12:30 tended to decrease more rapidly than those in the rainy season (Figure 4A-F).The differences in ΦPSII means across genotypes were noted among seasons being significantly (p < 0.001) higher (0.70) in the rainy than the hot (0.58) and cool (0.56) seasons (Table 2).However, no significance differences were found among genotypes in any season.Changes in ETR over the course of the day (Figure 4G-L) were related to the intensity of sunlight (Figure 3A-F).As shown in Table 2, daily means for ETR across genotypes in the hot (159 µmol(e -)/m 2 /s) and cool (157 µmol(e -) / m 2 /s) were significantly higher (p < 0.001) than that in the rainy (90 µmol(e -)/m 2 /s) season.However, no significance differences in mean ETR were found among genotypes in any season.Diurnal patterns of NPQ were similar to those of ETR (Figure 4M-R).Cassava genotypes exhibited highest mean NPQ across genotypes (0.45) in the hot season which was significantly different (p < 0.001) from that in the rainy season (0.34), and expressed an intermediate value (0.41) in the cool season (Table 2).In any season, mean NPQ among genotypes did not differ significantly.However, it was noted that NPQ of CMR-38-125-77 tended to be higher than those of the other genotypes during 10:30 to 14:30 in the hot season (Figure 4Q, R).

Diurnal Leaf Gas Exchange
Diurnal responses of leaf gas exchange and related parameters, i.e.Pn, Gs and Tr of four cassava genotypes measured during 04:30-18:30 on the same days as Chl fluorescence measurements were displayed in Figure 5.For each season, the diurnal minima, maxima and means of R, Pn, Gs, Tr, Ci Ci/Ca, LUE and WUE which were recorded during 06:30 to 16:30 across two days of investigation were summarized in Table 2.The values for dark respiration rates (R) in Table 2 were obtained from leaf gas exchange measurements in the dark at 04:30.
Seasonal variation in R of cassava leaves was observed in which the mean across genotypes was highest in the hot (2.24 µmolCO2/m 2 /s) followed by the significantly lower rates (p = 0.006) in the cool (1.65 µmolCO2/m 2 /s) and the rainy (1.47 µmolCO2/m 2 /s) seasons (Table 2).Significantly different R among genotypes were observed in the hot and cool seasons with KU50 showing significantly lower R (p < 0.05) than CMR38-125-77 in the cool season, and RY9 in the hot season.
Diurnal changes in Pn (Figure 5A-F) were closely related to those of PARleaf (Figure 3A-F).The highest mean Pn across genotypes was found in the cool season at 14.60 µmolCO2/m 2 /s followed by that of the hot season at 14.32 µmolCO2/m 2 /s.The lowest mean Pn across genotypes was found in the rainy season (11.75 µmolCO2/m 2 /s) which was significantly lower (p = 0.044) than that in the cool (14.60 µmolCO2/m 2 /s) season.However, Pn among the four genotypes were not significantly different in any season The patterns of diurnal response in Gs (Figure 5G-L) were related to variation in Pn (Figure 5A-F).Means of Gs across genotypes significantly (p < 0.001) differed among the three seasons being highest in the rainy (0.53 molH2O/m 2 /s) followed by the hot (0.32 molH2O/m 2 /s) and the cool (0.21 molH2O/m 2 /s) seasons (Table 2).In any season Gs among the four genotypes were not significantly different.
Diurnal responses of Tr (Figure 5M-R) followed similar patterns as those of Pn (Figures 5A-F), ETR (Figures 4G-L) and NPQ (Figures 4M-R).Cassava leaves displayed the highest mean Tr across genotypes (6.58 mmolH2O/m 2 /s) in the hot season which was significantly higher (p < 0.001) than that in the cool (4.65 mmolH2O/m 2 /s) and rainy (4.52 mmolH2O/m 2 /s) seasons (Table 2).Significant differences in mean Tr among genotypes were detected only in the hot season, with RY11 showing significantly higher mean (p = 0.014) than CMR38-125-77.
Diurnal pattern of changes in Ci/Ca was depicted in Figure 6A-F showing that Ci/Ca was high in the early morning, declining to reach minimum values mostly at 12.30, then slowly increased in the afternoon.Means of diurnal Ci/Ca across genotypes were similar in the cool and hot seasons (0.70 and 0.71, respectively) which were significantly lower (p < 0.001) than that in the rainy season (0.85) (Table 2).However, in any season, mean Ci/Ca among genotypes did not differ significantly.It was apparent that cassava plants utilized light energy for photosynthesis at different efficiencies in different seasons (Figure 6G-L).The means of LUE across genotypes was highest in rainy season (0.030) and significantly higher (p < 0.001) than that in the cool (0.023) and hot (0.019) seasons (Table 2).The means of LUE among genotypes were noted in the rainy and hot seasons.In the rainy season, RY9 had lowest LUE (0.023) which was significantly different (p = 0.045) from that of RY11 (0.035) and CMR38-125-77 (0.033).In the hot season, LUE of CMR38-125-77 was lowest (0.016) which was significantly lower (p = 0.041) than that of RY11 (0.023).In general, WUE tended to be higher in the early morning then declining throughout the day (Figures 6M-R).Cassava leaves expressed highest means of diurnal WUE across genotypes in the cool (4.52 µmolCO2/mmolH2O) season which was significantly higher (p < 0.001) than that in the rainy and hot season (2.33 and 2.05 µmolCO2/mmolH2O, respectively).In any season, no significance differences in WUE among genotypes were found.

The Relationship between Cassava Leaf Photosynthesis and Environmental Field Conditions
The relationships between leaf gas exchange and environmental parameters are demonstrated by the matrix of correlation coefficient values (r) for each season in Table 3.Among the four environmental parameters, PARleaf showed the highest positive correlations with Pn across all seasons.The correlation was highest in the rainy (r = 0.93, p < 0.01) followed by the cool (r = 0.82, p < 0.01) and the hot (r = 0.73, p < 0.01) seasons.Tleaf also expressed highly significant correlations with Pn i.e., 0.76, 0.48 and 0.46 in the rainy, cool and hot seasons, respectively.RHC, on the other hand, showed highly significant negative correlations with Pn in all seasons.Pn had significant positive correlations with VPDL only in the rainy (r = 0.57, p < 0.01) and cool (r = 0.33, p < 0.01) seasons.The relationships between Tr and all four environmental parameters were highly significant and occurred in the same ways as those for Pn.It is worth noted that the relationships between Tr and PARleaf had lower correlations than between Pn and PARleaf but the opposite occurred for RHC.The relationships between Ci and all four environmental parameters were also highly significant but occurred in the opposite directions as compared to those for Pn and Tr.The absolute values of correlation coefficient were highest between Ci and Tleaf especially in the rainy and hot seasons (r = −0.90,p < 0.01).The relationships between Pn and Tr were positively correlated in all seasons being higher in the hot (r = 0.75, p < 0.01) and cool (r = 0.74, p < 0.01) and lower in the rainy season (r = 0.67, p < 0.01).The positive correlations between Pn and Gs were highly significant only during the cool and hot seasons, whereas those between Tr and Gs were highly significant in all seasons i.e. r = 0.73, 0.45 and 0.39 in the cool, hot and rainy season, respectively.Negative correlations between Pn and Ci were highest in the rainy season (r = −0.84,p < 0.01) followed by the cool (r = −0.63,p < 0.01) and the hot (r = −0.50,p < 0.01) seasons.

Discussion
The results reported here clearly demonstrated cassava resiliency under variable climatic conditions, a trait of great importance for adaptation of the crop to future climate change/global warming in the tropics/subtropics where other staple crops might fail [32,33].The seasonal climatic pattern at the study site varied considerably during the study period.The monthly mean temperature fluctuated from 24 °C in February 2016 to 32 °C in April 2016, RH from 80% in October 2015 to 40% in April 2016 and rainfall from 323 mm in August 2015 to zero in December 2015 and February 2016 (Figure 1).However, diurnal patterns of photosynthesis (leaf Pn) of field-grown cassava under full irrigation (Figures 5A-F) displayed similar response pattern in different seasons which more or less paralleled with diurnal changes in PARA (Figure 3A-F), despite differences in air temperature and RH.In the rainy season, the daily peaks of Pn occurred at 10:30 or 12:30, depending on genotypes, showing maximum Pn between 27.50 and 36.0 µmolCO2/m 2 /s (Table 2).Much higher temperature and VPDair in the hot season could be the major causes of the shift in daily peak of Pn to 8:30 or 10:30 (Figure 5E, F) depending on genotypes while maximum Pn remained high between 31.57and 33.31 µmolCO2/m 2 /s (Table 2).Rosenthal et al. [34] measured diurnal Pn of field-grown cassava in Illinois, USA during mild summer (June-August), where total rainfalls were 377 mm, mean daily minimum and maximum temperatures were 18 and 30 °C, and found that these cassava plants also exhibited the Pn peaks around noon with Pn values varying from 18-28 µmolCO2/m 2 /s.In contrast, in a seasonally dry environment, cassava had maximum Pn of 33 µmolCO2/m 2 /s (average of 10 cultivars) in the early morning (08:00) and decreased thereafter [16].Therefore, the diurnal patterns of Pn and the range of maximum Pn values in the current study were in accordance with earlier reports.
Variations in diurnal photosynthesis were apparent among genotypes in the cool and hot season.In general, the pre-noon environments were favorable for all four genotypes with small variations in Pn (Figure 5A-F).Nevertheless, during the afternoon Pn of KU50 and CMR38-125-77 tended to be slightly lower than the others particularly in the cool season, while afternoon Pn of CMR38-125-77 had a tendency to be lower than the others in both cool and hot seasons (Figures 5C-F).This could partly be due to lower light incidence on the leaf surface of KU50 and CMR38-125-77 as a result of their prominent leaf-drooping behavior under high light intensity compared to RY9 and RY11 which hardly showed any drooping (see supplementary Figure S1).Leaf movement in cassava is known as a stress avoidance mechanism in both well-watered and stressed plants [35,36].Genotypic variation in seasonal Pn of field-grown cassava was well-documented in earlier reports [37][38][39].A recent study in four African cultivars of cassava (two landraces and two improved lines) under a greenhouse condition showed almost perfect bell-shaped diurnal response with Pn peak at 12:30 and maximum Pn varying from 22 to 27 µmolCO2/m 2 /s [40].
The relationships between physical parameters of leaf including PARleaf, Tleaf, RHC, and VPDL and photosynthetic parameters (Pn, Gs, Tr, and Ci) in different seasons are expressed in a correlation matrix shown in Table 3.In each season, Pn had the strongest positive correlations with PARleaf (Table 3).Among ecological factors, PAR together with temperature and VPD have been shown to be highly correlated with Pn [41,42].Diurnal variation in PAR differed in different seasons, and Pn showed a bell-shaped response parallel with PAR during the mild spring while Pn peaked very early in the morning during hot summer [40].In this study, although maximum PARA and PARleaf occurred in the rainy followed by the hot and the cool season (Table 1), the mean daily PARleaf values were lowest in the rainy season (p < 0.001, Table 1) due to frequent cloud cover causing highly fluctuating PARA (Figures 3A,B).Lowest mean daily PARleaf in the rainy season may attribute to lower mean daily Pn (11.75 µmolCO2/m 2 /s, p = 0.044) than the cool and hot seasons (Table 2).A recent report [43] also pointed out that solar radiation is a limiting factor in rainy season based on comparison between maximum net photosynthetic rates (Pnmax) from light response curve and the predicted Pn from actual solar radiation data.It is known that Pn in fully expanded young cassava leaves developed under sunny warm climate does not reach light saturation even up to 1800 µmol/m 2 /s [41].In each season, Pn were significantly correlated with PARleaf as well as Tleaf, RHC and VPDL (Table 3).In spite of the significant differences (p < 0.001) in mean Tair and VPDair between the cool and hot seasons (Table 1), well-watered cassava in this study performed equally well showing the daily mean Pn of 14.60 and 14.32 µmolCO2/m 2 /s, respectively (Table 2).Moreover, it is worth noting that in the hot season cassava genotype RY11 could maintain maximum daily Pn at 28.9 µmolCO2/m 2 /s while leaf temperature reached 39.8 °C.Previous studies reported optimum temperature range for cassava photosynthesis in tropical environments between 30-35 °C [38,44].This indicated that these improved cassava genotypes have been well-adapted to environments in different seasons in this climatic zone.Cassava can be widely adapted to environments and usually requires a warm climate with high solar radiation for optimum photosynthesis, growth and productivity [44].
Among leaf gas exchange parameters, in the rainy season under fluctuating light intensity but high RH, higher correlation was found between Pn and Ci (-0.84, p < 0.01) than between Pn and Gs (0.59, p < 0.01) (Table 3) indicating a stronger role of photosynthetic capacity on CO2 fixation or nonstomatal regulation as compared to stomatal control of photosynthesis.Gas exchange measurements of 15 cassava cultivars mostly during high rainfall periods also found higher correlation between Pn and Ci (-0.84) than Pn and Gs (0.40) [39].Non-stomatal limitation may be attributed to mesophyll resistance to CO2 flux, carboxylation efficiency of Rubisco and RuBP regeneration [45].The extent to which Rubisco limits photosynthesis depends largely on irradiance [46].Therefore, low mean daily Pn in the rainy season could be attributed to low activity of Rubisco under fluctuating and low mean daily PAR (Figures 3A,B; Table 1).In the cool and hot seasons, Pn was influenced by both stomatal and non-stomatal controls as indicated by highly significant correlations (p < 0.01) between Pn and Gs, and Pn and Ci (Table 3).In the cool season, higher correlation between Pn and Gs (0.74, p < 0.01) than Pn and Ci (0.63, p < 0.01) may indicate stronger role of stomatal limitation inferred by lowest Gs (0.21 molH2O/m 2 /s; p < 0.001) in the cool season (Table 3).In the cool and hot season, RHair was significantly lower (p < 0.001) than that in the rainy season (Table 1).Cassava is very sensitive to low air humidity and its stomatal conductance rapidly decreases in response to dry air irrespective of soil water conditions [44,47].In the cool and hot season, maximum Gs of approximately 0.5 molH2O/m 2 /s occurred at 08:30 or 10:30 (Figures 5I-L) coinciding with high RH and low VPDL in the early morning (Figures 3O-R,U-X), thereafter Gs declined in parallel with decreasing RH and increasing VPDL resulting in decreasing Pn during the afternoon.Nevertheless, under irrigation, cassava was able to maintain relatively high mean daily Pn (14.60 and 14.32 µmolCO2/m 2 /s; Table 2) while mean Gs values were higher than 0.15 molH2O/m 2 /s (Table 2) which is the threshold value above which the plants would be considered under non-or mild water stress conditions [48].It has been suggested that photosynthesis metabolism is substantially resistant to water stress until Gs is below 0.1-0.15molH2O/m 2 /s [49].
Diurnal changes in transpiration rates (Tr) paralleled closely those of Pn (Figure 3A-F, M-R) and had highly significant correlations with Gs (Table 3).This indicated that transpiration was greatly influenced by stomatal regulation particularly in the hot season.Transpiration rates were similar in the rainy and cool seasons (4.52 and 4.65 mmolH2O/m 2 /s), and significantly higher (p < 0.001) in the hot (6.58 mmolH2O/m 2 /s) season (Figures 5M-R, Table 2) coinciding with increasing VPDair and VPDL in the latter (Figures 3S-X, Table 1).Pronounced effects of VPDL on stomatal movement and transpiration rates were classically demonstrated in cassava [47].Similar seasonal variation was reported in orange trees in which transpiration rates of well-watered orange plants in summer were about 2.5 folds higher than that in winter [50].Highest Tr in the hot season resulted in lowest WUE (2.05 µmolCO2/mmol H2O) in the hot which was significantly lower (p < 0.001) than that in the cool season (Table 2).High WUE in the cool season (4.52 µmolCO2/mmol H2O) was due to lower Tr (4.65 mmol/m 2 /s) as a result of partially closed stomata in response to dry air (Table 2).El-Sharkawy and De Tafur [39] reported similar WUE value of 4.5 µmolCO2/mmol H2O which was averaged from numerous measurements performed from upper canopy leaves of 15 cassava cultivars.Mean Tr among genotypes were significantly different (p = 0.014) only in the hot season i.e. the values appeared in the order RY11 > KU50 = RY9 > CMR38-125-7, the same order as for Gs (Table 2).High Tr played crucial role in heat dissipation, therefore maximum Tleaf in the hot season was highest in CMR38-125-77 (44.3 °C) and lowest (41.0 °C) in RY11 and KU50 (Table 1).Lower Tr in RY9 and CMR38-125-77 was related to higher WUE, although not statistically significant, than RY11 and KU50 (Table 2).When measured under natural rainfed environments, 15 cultivars of cassava exhibited large variations in WUE from 3.89 to 4.74 µmolCO2/mmol H2O [39].
The observation that maximum quantum yield of PSII photochemistry (Fv/Fm) values of all four cassava genotypes at 04:30 (predawn) were between 0.823-0.878(Table 2) across all three seasons indicated that cassava leaves were healthy, and no chronic damages occurred in PSII [51].However, significantly (p < 0.001) lower Fv/Fm was observed in the cool season compared to the others.Negative effects of low temperature (in winter and spring) on reduction of Fv/Fm have been reported in roses [52] and temperate bamboo [53].Energy utilization by a leaf is indicated by diurnal changes in the effective quantum yield of PSII photochemistry (ΦPSII) which changed in the opposite direction as PAR and showed similar patterns in all seasons.However, as shown in Figures 4A-F, the recovery of ΦPSII in the afternoon occurred much earlier in rainy season (average ΦPSII at 14:30 was 0.77) than in the cool and hot seasons (average ΦPSII at 14:30 were 0.43 and 0.53 in the cool and hot season, respectively).In addition, mean daily ΦPSII across genotypes was significantly higher (p < 0.001) than in the hot and cool seasons (Table 2).This indicated the interactive effects between light intensity and other environmental parameters such as temperature and VPD in the cool and hot seasons.Even though the soil moistures were optimized due to irrigation, stressful environments in the cool and hot season clearly posed negative effects on energy utilization [40,54].
The patterns of diurnal changes in ETR (Figures 4G-R) were similar and parallel to the curves of PAR (Figure 3A-F).The ΦPSII and the derived ETR are dependent on ambient PAR.Hence, mean daily ETR across genotypes in the rainy season was significantly lower (p < 0.001) than the others (Table 2).Theoretically, under controlled conditions ΦPSII and ETR are accurately correlated with Pn and can be used to predict CO2 assimilation and hence productivity [28].Nevertheless, under stressful environments during the afternoon in the hot season, Pn of RY9 peaked at 08:30 on the 5th April 2016 which was the hotter day (Figure 5F) while ETR continued to increase to reach maximum at 12:30 (Figure 4L).Similar results were observed in cassava cv.RY9 under both irrigated and rainfed conditions [55] and also in other plants under irrigation such as soybean [56] and peach palms [57].This indicated that after 08:30, under limited CO2 availability due to stomatal closure after 8.30 (Figures 5K,L), higher proportion of reductants generated from electron transport could be allocated to alternative electron sinks most commonly photorespiration, Mehler reactions and cyclic electron flow [58,59].These alternative pathways served to balance photosynthesis electron transfer so that light energy is optimally used for CO2 fixation and over-reduction of electron carriers and excess generation of reactive oxygen species are prevented [60].The diurnal patterns of NPQ curves were parallel to the curves of PAR because NPQ operated to dissipate excess light energy as heat to protect PSII from photodamage [61,62].Both PAR and NPQ across genotypes were highest in the hot followed by the cool and significantly lower (p < 0.001) in the rainy season (Table 1 and 2).Moreover, similar to other previously mentioned parameters (Pn, ΦPSII and ETR), genotypic variations in NPQ were more apparent in the hot season than the others.Compared with the others CMR38-125-77 had a tendency to have lowest ΦPSII and ETR but highest NPQ (Figure 4E,F,K,L,Q,R) during 10:30-14:30 in the hot season, indicating the most active photoprotective mechanisms.The most prominent component of NPQ is qE which harmlessly dissipated excess light energy as heat through functioning of the xanthophyll cycle and PsbS protein [63].It was suggested that more PsbS protein and hence increasing qE capacity might improve crop production in adverse environments [64].Despite lowest leaf-level Pn in the hot season (Table 2), CMR-38-125-77 had better growth than the others as evidenced by significantly higher leaf production rate and LAI (Figure 2).A recent report in rice showed that rice transgenic line overexpressing PsbS protein which regulated qE had higher NPQ than wild type but comparable leaf-level ETR and Pn.Higher NPQ in the transgenic rice line efficiently protected PSII from photodamage, hence displayed better growth, higher leaf area per plant, higher photosynthetic performance, greater total biomass and finally higher grain yield than the wild type [65].
Although leaf-level Pn was lowest in the rainy season cassava plants had much better growth, as indicated by significantly higher (p < 0.001) means across genotypes of rate of increase in plant height and LAI than in the cool and hot seasons (Figure 2; Table S2).In relation to photosynthesis and early vegetative growth (0-3 months), all four genotypes performed equally well in the rainy and cool season (Figure 2).However, in the hot season, CMR-38-125-77 showed significantly higher leaf production rate and LAI than the others (Figure 2).Importantly, at the age of 3 months CMR38-125-77 growing in the hot season had LAI of 3.49 m 2 m −2 , not significantly different from those growing in the rainy season (LAI = 3.57 m 2 m −2 ), whereas hot-season LAI values of the other genotypes showed 20-37% reduction from those in the rainy season crop (Figure 2, Table S2).Although CMR38-125-77 had a moderate leaf-level Pn, its high WUE as well as efficient protective mechanisms (high NPQ and leaf drooping) were beneficial for its growth performance and canopy development.Since analyses of cassava growth and yield are usually evaluated on the basis of both LAI and Pn [14,66], it can be concluded that among the studied genotypes CMR38-125-77 is most suitable for planting in the post-rainy season (in December) to obtain good vegetative growth during the first 3 months.Cock et al. [9] suggested that cassava plants should reach LAI of 3.0 m 2 m −2 as quickly as possible in order to obtain good root yield.Similar results were obtained in a parallel experiment at the same site that CMR38-125-77 planted in December had highest LAI at 4 months after planting and subsequently gave highest storage root yield [12].Moreover, Sawatraksa et al. [67], who studied the same set of genotypes planted in December under rain-fed paddy field conditions, found that CMR38-125-77 had the highest biomass during early vegetative growth.

Conclusions
Growing cassava in the tropical savanna climate under irrigation, the environmental conditions in the rainy season were the most favorable for early vegetative growth of all four cassava genotypes, based on rate of increase in plant height and LAI at 3 MAP.Among the four genotypes, CMR38-125-77 was the most suitable genotype to be planted in December and growing during the period from cool to hot season, based on its highest rate of leaf production and LAI at 3 MAP.Mean daily net photosynthesis and electron transport rates in the rainy season were slightly lower than those in the cool and hot season due to fluctuating light intensity.Cassava plants displayed several morphological and physiological mechanisms in the hot season, to protect photosynthesis machinery from being damaged under the conditions of high light intensity, temperature and VPD, by leaf drooping, early stomatal closure, enhanced transpiration, thermal dissipation by NPQ and diversion of electrons to alternative sinks, and different genotypes may employ different strategies to varying extent.

Supplementary Materials:
The following are available online at www.mdpi.com/xxx/s1, Figure S1: leaf drooping at midday of cassava cultivars; RY9 (A), RY11 (B), KU50 (C) and CMR38-125-77 (D) growing in the field under irrigated condition.Table S1: ambient photosynthetically active radiation (PARA), air temperature (Tair), relative humidity (RHair), air vapor pressure deficit (VPDair), rainfall and irrigation for the three planting dates (PD).The environmental parameters were recorded from the date of planting until the date when photosynthesis measurements were performed on the three-month-old plants.Environmental parameters for PD-Jun, PD-Nov and PD-Dec were measured during 30 June-27 September 2015, 10 November 2015-23 February 2016 and 15 December 2015-2 April 2016, respectively.The values are daily minima, maxima and means of data collected every 5 min on each day, and the data were then averaged over days.Lower case letters indicate significance differences (P < 0.05) among planting dates.Table S2: the minima, maxima and means of growth parameters including rate of increase in plant height, rate of leaf production and leaf area index (LAI) of three-month-old plants of cassava genotypes RY9, RY11, KU50 and CMR38-125-77 planted in June 2015 (PD-Jun), November 2015 (PD-Nov) and December 2015 (PD-Dec).Means which are significantly different (p < 0.05) among seasons are denoted with capital letters.Means which are significantly different (p < 0.05) among genotypes are denoted with different lower case letters.

Figure 1 .
Figure 1.Monthly mean environmental parameters at the experimental site from May, 2015 to May, 2016.Daytime ambient photosynthetically active radiation (PARA) and daytime air temperature (Tair) (A), relative humidity (RHair) and total rainfall (B), air vapor pressure deficit (VPDair) and number of rainy days (C).Rainy season ranged from June to October, cool season from November to February, and hot from March to May.Dotted vertical lines indicate the timing of photosynthetic measurements made on three-month-old cassava plants in each season.Rainy season measurements were made in October 2015 on plants established on 30 June 2015 (PD-Jun); cool season measurements were made in February 2016 on plants established on 10 November 2015 (PD-Nov); and hot season measurements were made in April 2016 on plants established on 15 December 2015 (PD-Dec).

Figure 2 .
Figure 2. Rate of increase in plant height (A), rate of leaf production (B) and leaf area index (LAI) of three-month-old plants (C) of cassava genotypes RY9, RY11, KU50 and CMR38-125-77 planted in June 2015 (PD-Jun), November 2015 (PD-Nov) and December 2015 (PD-Dec).Means which are significantly different (p < 0.05) among genotypes for each planting date are denoted by different lower-case letters.For each genotype means which are significantly different among planting dates are denoted with capital letters.Data shows mean of six replicates ± standard deviation (SD).

Figure 3 .
Figure 3.Diurnal variation in physical parameters including PARleaf (A-F), Tleaf (G-L), RHC (M -R) and VPDL (S-X) of cassava belonging to four cultivars (x, RY9; , RY11; , KU50 and , CMR38-125-77) in comparison with ambient conditions (PARA, Tair, RHair and VPDair) in rainy (A,B,G,H,M,N,S and T), cool (C,D,I,J,O,P,U and V), and hot season (E,F,K,L,Q,R,W and X).The measurement was performed on two sunny days in each season.Leaf parameters were obtained from leaf gas exchange measurements while the environmental parameters were from the weather station.Data shows mean of two replicates ± SD.

Table 1 .
Seasonal variation in environmental and leaf/canopy parameters during the time of photosynthesis measurements i.e. two days in rainy season (25 Nov and 2 Oct 2015), cool season (24 and 26 Feb 2016) and hot season (3 and 5 Apr 2016).Photosynthetically active radiation (PARA) and PARleaf, ambient PAR and over the leaf surface; Tair and Tleaf, ambient air temperature and leaf temperature; RHair and RHC, ambient air relative humidity (RH) and RH in the plant canopy; VPDair and VPDL, ambient vapor pressure deficit and VPD of the leaf.Daily minima (min), maxima (max) and means for each genotype were from 24 measurements (six time points during 6:30-16:30 h x two leaves x two days).Means which are significantly different among seasons are denoted with capital letters.Mean of leaf parameter across genotypes which is significantly different from that of the environment is denoted with *.F and P value for testing each trait among genotypes within season (the same column). 2F and P value for testing each trait among seasons of each genotype (the same row). 1

Table 2 .
Seasonal variation in Chl fluorescence and leaf gas exchange parameters of four cassava genotypes.Fv/Fm, maximal photochemical quantum yield of PSII; ΦPSII, effective quantum yield of PSII photochemistry; ETR, electron transfer rate; NPQ, non-photochemical quenching; R, respiration rate; Pn; net photosynthetic rate; Gs, stomatal conductance; Tr, transpiration rate; Ci, intercellular CO2 concentration; Ci/Ca, ratio between intercellular and ambient CO2 concentration; LUE, light use efficiency; WUE, water use efficiency.Minima (min), maxima (max) and means were from 24 measurements (six time points during 6:30-16:30 h x two leaves x two days in each season), except the values for Fv/Fm and R in which 4 measurements were taken from two leaves of each genotype at 04:30 on two days.Means which are significantly different (p < 0.05) among genotypes are denoted with different lower case letters, whereas those among seasons are represented by different capital letters.