Photochemical efficiency correlated with candidate gene expression promote coffee drought tolerance

The aim of this study was to identify the correlation between photochemical efficiency and candidate genes expression to elucidate the drought tolerance mechanisms in coffee progenies (Icatu Vermelho IAC 3851-2 × Catimor UFV 1602-215) previously identified as tolerant in field conditions. Four progenies (2, 5, 12 and 15) were evaluated under water-deficit conditions (water deficit imposed 8 months after transplanting seedlings to the pots) and under irrigated system. Evaluations of physiological parameters and expression of candidate genes for drought tolerance were performed. Progeny 5 showed capacity to maintain water potential, which contributed to lower qP variation between irrigated and deficit conditions. However, the increases of qN and NPQ in response to stress indicate that this progeny is photochemically responsive to small variations of Ψam protecting the photosystem and maintaining qP. Data obtained for progeny 12 indicated a lower water status maintenance capacity, but with increased qN and NPQ providing maintenance of the ɸPSII and ETR parameters. A PCA analysis revealed that the genes coding regulatory proteins, ABA-synthesis, cellular protectors, isoforms of ascorbate peroxidase clearly displayed a major response to drought stress and discriminated the progenies 5 and 12 which showed a better photochemical response. The genes CaMYB1, CaERF017, CaEDR2, CaNCED, CaAPX1, CaAPX5, CaGolS3, CaDHN1 and CaPYL8a were up-regulated in the arabica coffee progenies with greater photochemical efficiency under deficit and therefore contributing to efficiency of the photosynthesis in drought tolerant progenies.


Results
For this study, three progenies potentially drought tolerant (progenies 5, 12 and 2) and one sensitive (progeny 15) were previously selected in the field, where the average productivity of the most productive progenies (progenies 5, 12 and 2) was 70.00 bag ha −1 and the vigor score was 10. The progeny 15, on the other hand, showed productivity of 34.00 bag ha −1 and the vigor score was 5.
Leaf Ψpd with values below − 1.0 MPa confirmed the unstressed condition of irrigated plants (Table 1). In this condition, no differences were observed between the progenies in relation to the values of net CO 2 assimilation rates (A), stomatal conductance (gs), transpiration (E), vapor pressure deficit (VPD) and leaf temperature (Tleaf), quantum yield of PSII electron transport (ΦPSII), ETR photochemical quenching, non-photochemical quenching (qN) and the NPQ. For variables A, gs and E there was an increase in their values in the second evaluation followed by a slight decrease in the last, while for the variables VPD and Tleaf, highest values were observed at 35 days after irrigation withdraw. The progeny 15 presented lower leaf water potential compared to the other progenies, reaching values below − 5.0 MPa at 28 and 35 days after irrigation withdraw. It is worth to mention that at 35 days, progeny 12 reached the water potential of − 3.35 MPa, while progeny 2 and 5 maintained higher values of water potential, being 1.30 and 1. 63 MPa, respectively (Table 1). The values of net CO 2 assimilation rates (A), stomatal conductance (gs), and transpiration (E) were lower under water deficit conditions than irrigated controls for all progenies at 28  In relation to the parameters ΔF and ΔFmʹ, positive values were found indicating an increase of F and Fmʹ under water deficit (Fig. 1A,B). From the flux density of photosynthesizing photons of 600 μmol m −2 s −2 , progeny 2 presented higher ΔF while the lowest value was found in progeny 15. As for ΔFmʹ, the progeny 12 presented the highest value when submitted to 162 to 230 μmol m −2 s −2 , as compared to the other progenies.
Considering the ΔETR, the curves presented negative values, which shows that under water deficit, the ETR of the progenies were lower than when irrigated (Fig. 1C). We should draw attention to progeny 15, which presented a more negative ΔETR, and the progeny 12 with ΔETR values close to zero. Of these, progeny 15 presented a more expressive reduction of ETR and progeny 12, a small magnitude of reduction. For progenies 5 and 2, the reduction was intermediate.
As for ΔɸPSII, corresponding to the variation of the effective quantum yield, progeny 12 presented a variation close to zero, indicating that the deficit did not reduce the ɸPSII values, as compared to the irrigated one (Fig. 1D). In the other progenies, the values of the curves were negative and, therefore, the effective quantum yield under water deficit was lower in these progenies. With the increase in the density of photosynthetic photon flux, the differences between the progenies were less pronounced, with no differences occurring at 1100 μmol m −2 s −2 .
In general, with increasing photosynthesizing photon flow density, progenies 5 and 12 presented positive values of ΔqN and ΔNPQ, indicating higher qN and NPQ, when submitted to water deficit (Fig. 1e,f). In contrast, progenies 15 and 2 presented negative values of ΔqN and ΔNPQ, therefore with reduction of these parameters under water deficit. The values of ΔqP were negative indicating a reduction of qP under deficit, but the progeny Table 1. Pre-dawn water potential (Ψpd), photosynthetic rate (A), stomatal conductance (gs), transpiration (E), vapor pressure deficit (VPD) and leaf temperature (Tleaf) of 4 irrigated (I) and non-irrigated (NI) progenies at 0, 28 and 25 days of evaluation. Means followed by the same lowercase letter in the column and averages followed by the same capital letter in the row in each irrigation condition do not differ statistically from each other at the 5% probability level by the Duncan test. Means followed by asterisks in the line, for the same time and even progeny, are statistically different at the 5% probability level by the Duncan test.  www.nature.com/scientificreports/ 5, which presented higher values of delta, was observed, indicating that the qP was less affected in this progeny as compared to the others (Fig. 1g).
Regarding the maximum photochemical efficiency of PSII (evaluated by the Fv/Fm ratio), similar values were observed under irrigated and non-irrigated conditions for progenies 5 and 12 ( Fig. 2A). However, the Fv/Fm ratio was lower in the non-irrigated condition compared to the irrigated condition in progenies 15 and 2, with lower values in relation to progenies 5 and 12. In the drought-stressed progenies, the ETR/A ratio enhanced remarkably. However, there were no significant differences between them (Fig. 2B).
In order to verify the correlation between physiological characteristics and gene expression, a statistical analysis was performed, using PCA, at 35 days, when plants effectively reached a water stress (Fig. 3). The analyses considered the two PCs corresponding to the two eigenvalues with the highest values, which explained more than 65% of the variance. Figure 3A shows a greater variance explained by the PC1 for positively associated variances with CaMYB1, CaAPX1, CaGolS3, CaERF017, CaDHN1, CaNCED, CaAPX5 and CaEDR2 and negatively with the physiological characteristics. The PC1 was able to differentiate irrigation regimes, with higher values of gene expression for non-irrigated treatments. On the other hand, in PC2 the variables of greatest contribution were CaMYB1, CaAPX1, CaGolS3, CaERF017, CaDHN1, CaNCED, CaAPX5, CaEDR2 and CaPYL8a, together with the fluorescence variables. In this case, progenies 12 and 5 under water deficit conditions were discriminated together with the progenies under irrigated conditions (Fig. 3B). According to the PCA, gene expression analyses will be focused and detailed for the genes that contributed to the formation of PC2, since they were decisive in distinguishing progenies 5 and 12 that showed the greatest photochemical efficiency.
A first univariate analysis for the expression data of candidate genes, all the treatments (time and condition) were grouped by progeny and analyzed statistically. The only significant difference among the expression patterns of the different progenies was found for the CaMYB1 (Fig. 4) In this case, the pattern of relative expression, among all treatments, displayed by the progenies 5 and 12, was statistically different from progenies 2 and 15 ( Fig. 4).
The second analysis was performed only with the gene-expression data of the 35 days condition as it was done for the PCA analysis. The obtained results with the 35 day's data analysis among the treatments and progenies showed that statistical differences were observed for the expression data of all genes studied. The results of the analyses of the expression profiles of some transcription factors such as CaERF017, CaMYB1 and CaEDR2 are presented on Figs. 5A-C. For CaERF017 an expressive expression was found at 35 days, with higher values found under non-irrigated progenies 2, 5 and 12 and lower values under both, irrigated and non-irrigated progeny 15 (Fig. 5A). For the CaMYB1 gene-expression levels, progenies 5 and 12 presented the highest expression levels under the non-irrigated condition (Fig. 5B). However, as shown on Fig. 3, high levels of CaMYB1 expression were also observed for progenies 5 and 12, at the time point of 28 days. For the CaEDR2 gene, at 35 days, progenies 2 and 5 presented higher expression levels under both conditions, as compared to the other progenies, with significant higher values under the non-irrigated condition as compared to the irrigated (Fig. 5C).
For CaAPX1, expression was up-regulated in progenies 2, 5 and 12 under the non-irrigated condition (Fig. 5D). Furthermore, for CaAPX5 transcript levels, higher expression values were found for non-irrigated progenies 5 and 12 (Fig. 5E). Higher values of relative expression of CaGolS3 were observed in progeny 5 followed by 12, both under the non-irrigated condition, as compared to the other treatments studied (Fig. 5F).
Results of the CaDHN1, shows that the highest value of relative expression was observed for progeny 5 in the non-irrigated condition, whereas for progeny 12 there was no differential expression between irrigated and non-irrigated condition (Fig. 5G). The highest value of relative expression for the CaNCED gene was presented by progeny 5 non-irrigated at 35 days (Fig. 4H). The CaPYL8a gene displayed the higher values for progenies 5 and 12 for the non-irrigated condition as compared to the other treatments (Fig. 5I).
Noteworthy that for all genes tested, as shown on Fig  Ranked results were analyzed using the non-parametric Kruskal-Wallis test followed by a pairwise comparison with the Mann-Whitney test. Significance level were established as p-values < 0.05, for both tests. Means followed by the same lowercase letter and averages followed by the same capital letter in each irrigation condition do not differ statistically from each other. Means followed by asterisks for even progeny are statistically different.

Discussion
In Brazil, the damaging effects of water deficit are enhanced by the high irradiances usually observed in coffeegrowing regions 4 . The coffee tree has a decreased stomatal conductance (gs) under water stress, which leads to a huge oxidative pressure 26,27 . Therefore, progenies which can undertake adjustments in photochemical and molecular processes might better endure drought stress conditions. In this study, we sought to understand about the strategies these different progenies use in order to avoid photochemical damaging.
The results presented here indicated that coffee progenies 2, 5 and 12 that presented high productivity under field conditions in a year of severe water deficit, considered potentially tolerant, showed a different pattern of photochemical and molecular responses to the drought conditions. Progeny 5 showed capacity to maintain water potential, which contributed to lower qP variation and gas exchange between irrigated and deficit conditions, proving the drought tolerance condition observed in the field. However, qN and NPQ increases in response to the water stress indicate that this progeny is photochemically responsive to small variations of Ψpd protecting the photosystem and maintaining qP. As qP represents the portion of absorbed energy that is used in the photochemical step by photosystem II reaction centers 17 , it may also be indicative of the proportion of reaction centers that are open 13 . The highest values of the non-photochemical quenchings qN and NPQ protect the reaction centers of photosystem II by dissipating excess energy in the form of heat 28,29 . This was also evident from values of Fv/ www.nature.com/scientificreports/ Fm that remained unchanged in stressed plants, indicating that there was no photoinhibition in this progeny and that light reactions were in a fully active state 26 . Among the tolerants, progeny 12 presented lower water status maintenance capacity, as verified in the lower values of Ψpd. It should be noted, however, that the reduction of gas exchange under water deficit showed similar levels to those progenies that maintained the leaf-water potential. In spite of having reduced qP under water deficit, the increase of qN and NPQ provided maintenance of the photochemical parameters as ɸPSII and ETR. This shows that this progeny has a certain tolerance of the photosynthetic apparatus to a severe water deficit (Ψpd = − 3.35 MPa). The stress condition may lead to a slight increase in excitation pressure at the reaction center in the PSII 30 , as indicated by the decrease in qP, however, not being sufficient to decrease the ETR, which in turn reflects activities of photosystem II 28 and the effective photochemical efficiency of PSII (ɸPSII). The maintenance of Fv/Fm values in this progeny under drought also highlights that there were no PSII inhibitions 26 . Similar results were found by Li et al. 28 that observed under drought an increase of qN and decrease of qP in sugar beet leaves, suggesting that excess energy was dissipated as heat, which prevented a photochemical inhibition, due to a decreased degree of reaction-centers opening.
On the other hand, progeny 2, even with higher Ψam, proved to be photochemically sensitive, since there was a greater reduction in qP, ETR, qN, NPQ and Fv/Fm under water deficit. These results indicate that the photosynthetic apparatus may have been damaged by photoinhibition and that it was unable to protect itself against excess energy 26,28,31 .
The progeny 15 exhibited sensitivity to drought, as previously noted under field condition as this progeny had the lowest productivity in the year of occurrence of a severe drought. It presented the lowest capacity of water potential (Ψpd) maintenance after 28 days of water stress, due to the higher decay rate of this variable among the progenies studied. Progeny 15 showed photochemical sensitivity, with reduction of all fluorescence parameters under water deficit. There is a decrease in CO 2 availability for ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) due to stomatal closure induced by drought. This may causes an excess of energy in the chloroplasts that can lead to an increase in ROS levels, increasing oxidative stress 7 . Thus, drought stress effects on the photosynthetic apparatus of this progeny has typically started with stomatal effects, and ended in metabolic changes caused by severe stress 32 . It is emphasized that the plants were not exposed to different levels of oxidative pressure since there were no differences among the progenies in ETR/A ratio within the same system. However, a high ETR/A ratio was observed in all progenies under drought stress, which are likely associated with an increase in www.nature.com/scientificreports/ reducing power that might lead to oxidative stress 27 . In progeny 15 drought-stressed, a decrease in Fv/Fm was found in relation to well-watered plants, suggesting that effective oxidative damage occurred in these plants. It was also observed that in South African coffee cultivars lower values of Fv/Fm were found in drought-sensitive cultivars 26 . Also, the agronomic vigor analysis showed that progenies that didn't present effective photochemical adjustments had lowest agronomic vigor, marked depletion symptoms and yellowing (leaf scald). Overall, progenies 5 and 12 showed high photochemical efficiency and energy dissipation capacity, analyzed by the parameters of the chlorophyll a fluorescence curves. Tounekti et al. 26 studying South African coffee cultivars confirmed that photoprotection is an important factor that affects photosynthetic productivity, and that it varies between the cultivars. Interestingly, the PCA analysis revealed that the expression levels of CaMYB1, CaAPX1, CaGolS3, CaERF017, CaDHN1, CaNCED, CaAPX5, CaEDR2 and CaPYL8a genes clearly displayed the major response to drought stress and discriminated these progenies (5 and 12) under the water-deficit conditions. Accordingly, the gene expression analysis showed that for CaMYB1, CaERF017and CaEDR2, transcript levels coding for these regulatory proteins were higher in the tolerant progenies 5 and 12 as compared to the others. In C. arabica CaERF017 gene was identified as C canephora DREB-like gene 25 . The DREB (Dehydration Responsive Element Binding) gene family corresponds to key transcription factors involved in responses to various abiotic stresses, which regulate the expression of various responsive genes, acting downstream of ABA-dependent and -independent signal transduction pathways 33,34 . The CaERF017 gene was the most expressed DREB-like gene in leaves under low humidity in the drought tolerant clones 25 . The MYB transcription factors have been reported by regulating the biosynthesis of secondary metabolites, regulation of lateral root growth and of the stomatal Abundance of transcripts was normalized using CaUBQ10 gene expression as an endogenous control. The results were expressed using 15 I as the reference sample (Relative Expression = 1). Values of three replicates were presented as means ± SE (bars). Ranked results of the relative expression data of the progenies were analyzed using the non-parametric Kruskal-Wallis test followed by a pairwise comparison with the Mann-Whitney test. Significance level were established as p-values < 0.05, for both tests. Means followed by the same lowercase letter and averages followed by the same capital letter in each irrigation condition do not differ statistically from each other. Means followed by asterisks for even progeny are statistically different.  35,36 . The EDR2 gene encodes a protein that acts as a negative regulator of cell death caused by the attack of pathogens and mediated by salicylic acid 37 . Among ABA-related genes, CaNCED expression highlighted the progeny 5 without irrigation. NCED catalyzes the cleavage of 9-cis-xanthophylls to xanthoxin, which is the key regulation setp of ABA biosynthesis in plants 38 . Accumulation of this hormone can trigger the expression of many stress-responsive genes increase the tolerance of the plant 39 . In C. canephora the tolerant clone 14 showed expression 4 times larger than the other clones in response to drought. NCED gene expressions were strongly induced in response to water stress in tolerant progenies of several common bean progenies 40,41 . The ABA synthesized is perceived by ABA receptors as the PYRABACTIN RESISTANCE1 (PYR1)/PYR1-LIKE (PYL)/REGULATORY COMPONENTS OF ABA RECEPTORS (RCAR) family of proteins 42,43 . Greater expression of CcPYL3 and CcPYL7 was detected in tolerant clones of C. canephora, but with a differential expression between them 21 as observed in this study for C. arabica with greater expression in progenies 5 and 12 under drought.
The high levels of gene expression encoding isoforms of the antioxidant enzyme ascorbate peroxidase in progenies 5 and 12, mainly under water deficit condition, may be related to the maintenance of the gas exchanges and greater photochemical efficiency in these progenies under non-irrigated condition. The excess excitation energy can lead to the production of reactive oxygen species (ROS) that can oxidize lipids, proteins and nucleic acids in water deficit condition 44 . The peroxidase enzyme ascorbate (APX) plays a key role in the removal of ROS, catalyzing the conversion of H 2 O 2 to H 2 O, using ascorbate with the electron donor 45 . Thus, increased expression of enzymes antioxidants may provide greater protection of photosystem II and other membranes to oxidative stress, contributing to efficiency of the photochemical phase of photosynthesis in these genotypes that showed the highest values of Fv/Fm. Similarly, the CaGolS3 expression was detected with higher levels in the progenies 5 and 12, corroborating with the literature that shows that in C. arabica the GolS3 isoform is highly expressed under moderate and severe drought conditions 46 . In C. canephora a fourfold higher expression was also observed in the tolerant clone 14 under drought condition 18 . CaGolS3 gene encodes an isoform of the enzyme galactinol synthase that is involved in the first step of oligosaccharide synthesis in the raffinose (RFOs) Family. RFOs are involved both in osmoprotection as in ROS scavenging, thus promoting increased protection the photosynthetic apparatus against oxidative damages caused by various stresses such as drought 47 . On the other hand, the lower levels of expression of these genes in progeny 2 may have contributed to the decrease of all fluorescence values, perhaps because of a less effective protective and antioxidant action.
Noteworthy is also that the CaDHN1 gene expression was found in non-irrigated plants mainly in progenies 5 and 12. Dehydrins are known as a group II of LEA proteins (Late embryogenesis abundant), accumulated in the later stages of development in response to the most diverse environmental stresses, participating significantly in the stabilization of membranes, enzymes and nucleotides 48 . The expression of dehydrins can be strongly induced by abiotic stresses 49,50 , being more expressed in tolerant clones 18,22 .
In our study, the better photochemical efficiency and greater capacity to dissipate unused energy in photosynthesis, was clearly associated with higher gene expression levels coding for regulatory proteins, ABA-synthesis, cellular protectors and isoforms of the antioxidant enzyme ascorbate peroxidase, indicating that these genes can contribute to maintaining the integrity of the photosynthetic apparatus in progenies tolerant to water deficit. On the other hand, progeny 15 showed photochemical sensitivity, as well as the lowest expression levels of almost all genes tested. In addition these results are also in agreement with the lowest productivity of progeny 15, observed under field conditions in a year of a severe drought period.
Progenies with better photochemical efficiency and greater capacity to dissipate unused energy in photosynthesis may present better condition to maintain its function under low leaf water status and to recovery after a period of severe drought, under a large scale of radiation intensity. crossing, based on data on productivity and agronomic vigor observed in 2014/2015 crop year (a period with severe water deficits in this region). According to the EPAMIG genetic improvement program database, three more productive genotypes with maximum agronomic vigor (progenies 5, 12 and 2) and one less productive and with lowest agronomic vigor (progeny 15) were selected. Plants were scored according to a ten point arbitrary scale in which score 1 corresponds to plants with the poorest vegetative vigor and marked depletion symptoms. The score 10 indicates plants with excellent vigor, more leafy and with marked vegetative growth of the productive branches 51,52 . Productivity is measured by weighing the fruit immediately after harvest. Then a sample of about 3 L from each plot is weighed and placed for drying in the sun. After drying, the coffee is weighed, processed and weighed again to calculate productivity in processed bags per hectare (bag ha −1 ).

Material and methods
The use of coffee progenies in the present study has the appropriate permissions and/or licenses of the EPAMIG coffee breeding program. The use of coffee plants parts complies with international, Brazilian and institutional guidelines. www.nature.com/scientificreports/ Climatological Station of Lavras). The minimum, average and maximum temperatures and relative humidity and photosynthetically active radiation of the greenhouse during the experimental period were 17 °C, 26 °C, 35 °C, 39%, 63%, 87% and 300 µmol m −2 s −1 , respectively. The seedlings were grown in 20 L polyethylene pots, with a substrate consisting of a mixture of subsoil soil, sand and bovine manure (3: 1: 1, v/v/v). The fertilization was carried out according to the substrate analysis, following the recommendations of Guimarães et al. 53 . Phytosanitary treatments were carried out preventively to control the main pests and diseases of higher incidence in the region. After being transferred to the pots, the seedlings were kept in a greenhouse, with free exchange of air during a period of eight months.
Until the day before the application of the water treatments all plots were irrigated in order to maintain the substrate of the vase in the field capacity (FC). In order to control the amount of water applied in each vessel, the humidity sensor (ML2X THETA PROBE, Delta-T Devices) was used and calibrated according to soil water moisture, from soil saturation to field capacity.The relative soil moisture at 100% field capacity was 0.46 m −3 m −3 . After the eight-month period, the plants were separated into two groups. In the first group (i.e., irrigated plants (I), control), plants remained irrigated, maintaining 100% of the field capacity. The plants from the second group (i.e., drought-stressed plants (NI)) were not watered until the soil water content reached approximately 35% FC (i.e., 0.16 m −3 m −3 ; ca. 28 days). They were then maintained at 35% FC for additional 7 days, (i.e., 0. Statistical analyses. Physiological traits were evaluated by analysis of variance using the F test, with the Duncan's New Multiple Range Test (DNMRT) at 5% probability for to perform the mean differences using the R software 54 and agricolae package 55 . The physiological traits and gene expression profiles were also analyzed by Principal Component Analysis (PCA), using the R software 54 and the FactoMineR package 56 .
For the RT-qPCR results, Fv/Fm and ETR/A, ranked results were analyzed using the non-parametric kruskal-Wallis test followed by a pairwise comparison with the Mann-Whitney test. Significance levels were established as p-values < 0.05, for both tests. The spreadsheet software Excel was used to perform the analyses.
Physiological traits. All samplings and measurements were made using leaves from the third or fourth pair of apex plagiotropic branches.
Leaf samples for pre-dawn water potential (Ψpd) measurements were collected between 3 and 5 am. The collected leaves were placed in plastic bags, properly closed and identified according to each treatment. Then the samples were placed in a thermal box with ice, to avoid excessive water loss, and taken to the place where Ψpd measurements were made using the Scholander-type pressure chamber (PMS Instruments Plant Moisture-Model 1000).
Chlorophyll fluorescence was measured with a pulse amplitude modulation fluorometer (MINI PAM-1999, Heinz Walz GmbH, Effeltrich, Germany). The analyzes were carried out following the methodology described by Schreiber et al. 57 and Genty et al. 58 and the terminologies and calculations proposed by Snel and Van Kooten 59 . Leafclips were used for measurements of the minimum (F0) and maximum (Fm) fluorescence yield in darkadapted (30 min) leaf tissues. The minimum fluorescence (F0) was measured using a weak modulated measuring beam (ca. 0.03 μmol m −2 s −1 ). Subsequently, the maximal fluorescence (Fm) was measured by applying a saturating actinic light pulse (> 6000 μmol m −2 s −1 ) for 0.8 s. Thus, it was possible to calculate the maximum photochemical efficiency of photosystem II (PSII), Fv/Fm = (Fm − Fo)/Fm.
In light-adapted leaf tissues, instant light-response curves were obtained using the light-curve program of the Mini-PAM. The photosynthetic photon flux, ranging from 0 to 1620 µmols m −2 s −1 for a period of 4 min in eight steps following each other within 30 s. At the end of each light level, a saturating pulse was applied to determine fluorescence parameters. These Chl fluorescence parameters were used to effective photochemical efficiency of PSII (ΦFSII = ΔF/Fm′), the apparent electron transport rate (ETR = (ΔF/Fm′ × PPFD × 0.5 × 0.84), the photochemical quenching (qP = (Fm′-F)/(Fm′ − F0′)) and the non-photochemical quenching [NPQ = (Fm − Fm′)/ Fm′; qN = Fm − Fm′/Fv)] 60 . For ETR calculation, it was assumed that quanta were evenly distributed between PSII and PSI (0.5), and leaf light absorption was considered to be 0.84 61 . At 35 days of stress imposition the light curves were presented as delta values (Δ), i.e., the difference between the values of the progenies under water deficit in relation to the irrigated control. This evaluation also included Fv/Fm and electron transport-to-net photosynthesis rate ratio (ETR/A). Samples  qPCR real-time assay. For real-time quantitative PCR, the first strand cDNA was synthesized using 1 μg of total RNA, the ImProm-II Reverse Transcription System and oligo(dT15) according to the manufacturer's recommendations (Promega, Madison, WI, USA). Real-time quantitative PCR (RT-qPCR) experiments were carried out using the protocol recommended for the use of a 7500 Fast Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) as previously described by Marraccini et al. 21 . Primers (Table 2) were designed using Primer Express software (Applied Biosystems) and were preliminarily tested for their specificity and efficiency against a cDNA mixture from leaves. Data were analysed using 7500 Fast Software v2.0.6 (Applied Biosystems) to determine cycle threshold (Ct) values. The specificity of the PCR products generated for each set of primers was verified by analysing the Tm (dissociation) of amplified products. PCR efficiency (E) was estimated using absolute fluorescence data captured during the exponential phase of amplification of each reaction with the equation (1 + E) = 10(− 1/slope) 62 . Expression levels were calculated by applying the formula (1 + E )−ΔΔCt where ΔCt target = Ct targetgene − Ct CaUBQ and ΔΔCt = ΔCt target − ΔCt reference sample . Gene expression levels were normalized to the expression level of the CaUBQ10 63 reference gene. Genotype 15 irrigated at zero time was used as internal calibrator with relative expression equal to 1.

Data availability
All data generated or analysed during this study are included in this published article. www.nature.com/scientificreports/ Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.