Rhythms produced by the endogenous circadian clock play a critical role in allowing plants to respond and adapt to the environment. While there is a well-established regulatory link between the circadian clock and responses to abiotic stress in model plants, little is known of the circadian system in crop species like soybean. This study examines how drought impacts diurnal oscillation of both drought responsive and circadian clock genes in soybean. Drought stress induced marked changes in gene expression of several circadian clock-like components, such as LCL1-, GmELF4- and PRR-like genes, which had reduced expression in stressed plants. The same conditions produced a phase advance of expression for the GmTOC1-like, GmLUX-like and GmPRR7-like genes. Similarly, the rhythmic expression pattern of the soybean drought-responsive genes DREB-, bZIP-, GOLS-, RAB18- and Remorin-like changed significantly after plant exposure to drought. In silico analysis of promoter regions of these genes revealed the presence of cis-elements associated both with stress and circadian clock regulation. Furthermore, some soybean genes with upstream ABRE elements were responsive to abscisic acid treatment. Our results indicate that some connection between the drought response and the circadian clock may exist in soybean since (i) drought stress affects gene expression of circadian clock components and (ii) several stress responsive genes display diurnal oscillation in soybeans.

The soybean transcriptome displays strong variation along the day in optimal growth conditions and also in response to adverse circumstances, like drought stress. However, no study conducted to date has presented suitable reference genes, with stable expression along the day, for relative gene expression quantification in combined studies on drought stress and diurnal oscillations. Recently, water deficit responses have been associated with circadian clock oscillations at the transcription level, revealing the existence of hitherto unknown processes and increasing the demand for studies on plant responses to drought stress and its oscillation during the day. We performed data mining from a transcriptome-wide background using microarrays and RNA-seq databases to select an unpublished set of candidate reference genes, specifically chosen for the normalization of gene expression in studies on soybean under both drought stress and diurnal oscillations. Experimental validation and stability analysis in soybean plants submitted to drought stress and sampled during a 24 h timecourse showed that four of these newer reference genes (FYVE, NUDIX, Golgin-84 and CYST) indeed exhibited greater expression stability than the conventionally used housekeeping genes (ELF1-β and β-actin) under these conditions. We also demonstrated the effect of using reference candidate genes with different stability values to normalize the relative expression data from a drought-inducible soybean gene (DREB5) evaluated in different periods of the day.

Since drought can seriously affect plant growth and development and little is known about how the oscillations of gene expression during the drought stress-acclimation response in soybean is affected, we applied Illumina technology to sequence 36 cDNA libraries synthesized from control and drought-stressed soybean plants to verify the dynamic changes in gene expression during a 24-h time course. Cycling variables were measured from the expression data to determine the putative circadian rhythm regulation of gene expression. We identified 4866 genes differentially expressed in soybean plants in response to water deficit. Of these genes, 3715 were differentially expressed during the light period, from which approximately 9.55 % were observed in both light and darkness. We found 887 genes that were either up- or down-regulated in different periods of the day. Of 54,175 predicted soybean genes, 35.52 % exhibited expression oscillations in a 24 h period. This number increased to 39.23 % when plants were submitted to water deficit. Major differences in gene expression were observed in the control plants from late day (ZT16) until predawn (ZT20) periods, indicating that gene expression oscillates during the course of 24 h in normal development. Under water deficit, dissimilarity increased in all time-periods, indicating that the applied stress influenced gene expression. Such differences in plants under stress were primarily observed in ZT0 (early morning) to ZT8 (late day) and also from ZT4 to ZT12. Stress-related pathways were triggered in response to water deficit primarily during midday, when more genes were up-regulated compared to early morning. Additionally, genes known to be involved in secondary metabolism and hormone signaling were also expressed in the dark period. Gene expression networks can be dynamically shaped to acclimate plant metabolism under environmental stressful conditions. We have identified putative cycling genes that are expressed in soybean leaves under normal developmental conditions and genes whose expression oscillates under conditions of water deficit. These results suggest that time of day, as well as light and temperature oscillations that occur considerably affect the regulation of water deficit stress response in soybean plants.

Water deficit is an important climatic problem that can impair agriculture yield and economy. Genetically modified soybean plants containing the AtNCED3 gene were obtained aiming drought-tolerance improvement. The NCED3 gene encodes a 9-cis-epoxycarotenoid dioxygenase (NCED, EC 1.13.11.51), an important enzyme in abscisic acid biosynthesis. ABA activates the expression of drought-responsive genes, in water-deficit conditions, targeting defense mechanisms and enabling plants to survive under low water availability. Results from greenhouse experiments showed that the transgene AtNCED3 and the endogenous genes GmAREB1, GmPP2C, GmSnRK2 and GmAAO3 presented higher expression under water deficit (WD) in the event 2Ha11 than in WT-plants. No significant correlation was observed between the plant materials and WD conditions for growth parameters; however, gas exchange measurements decreased in the GM event, which also showed 80% higher intrinsic water use when compared to WT plants. In crop season 2015/16, event 2Ha11 showed higher total number of pods, higher number of pods with seeds and yield than WT plants. ABA concentration was also higher in GM plants under WD. These results obtained in field screenings suggest that AtNCED3 soybean plants might outperform under drought, reducing economic and yield losses, thus being a good candidate line to be incorporated in the soybean-breeding program to develop drought-tolerant cultivars.