Genome-wide transcriptome analysis of soybean primary root under varying water-deficit conditions
© Song et al. 2016
Received: 10 September 2015
Accepted: 6 January 2016
Published: 15 January 2016
Soybean is a major crop that provides an important source of protein and oil to humans and animals, but its production can be dramatically decreased by the occurrence of drought stress. Soybeans can survive drought stress if there is a robust and deep root system at the early vegetative growth stage. However, little is known about the genome-wide molecular mechanisms contributing to soybean root system architecture. This study was performed to gain knowledge on transcriptome changes and related molecular mechanisms contributing to soybean root development under water limited conditions.
The soybean Williams 82 genotype was subjected to very mild stress (VMS), mild stress (MS) and severe stress (SS) conditions, as well as recovery from the severe stress after re-watering (SR). In total, 6,609 genes in the roots showed differential expression patterns in response to different water-deficit stress levels. Genes involved in hormone (Auxin/Ethylene), carbohydrate, and cell wall-related metabolism (XTH/lipid/flavonoids/lignin) pathways were differentially regulated in the soybean root system. Several transcription factors (TFs) regulating root growth and responses under varying water-deficit conditions were identified and the expression patterns of six TFs were found to be common across the stress levels. Further analysis on the whole plant level led to the finding of tissue-specific or water-deficit levels specific regulation of transcription factors. Analysis of the over-represented motif of different gene groups revealed several new cis-elements associated with different levels of water deficit. The expression patterns of 18 genes were confirmed byquantitative reverse transcription polymerase chain reaction method and demonstrated the accuracy and effectiveness of RNA-Seq.
The primary root specific transcriptome in soybean can enable a better understanding of the root response to water deficit conditions. The genes detected in root tissues that were associated with key hormones, carbohydrates, and cell wall-related metabolism could play a vital role in achieving drought tolerance and could be promising candidates for future functional characterization. TFs involved in the soybean root and at the whole plant level could be used for future network analysis between TFs and cis-elements. All of these findings will be helpful in elucidating the molecular mechanisms associated with water stress responses in soybean roots.
KeywordsSoybean Gene expression profile Water-deficit RNA-Seq Hormone interplay Metabolism Transcription factor Root
Soybean is the third most cultivated crop worldwide and the most important vegetable oil and protein source for humans. Soybean production is limited by environmental constraints, particularly drought [1, 2]. Water availability to the plant is the most important factor, as soybeans use approximately 450–700 mm of water during the growing season .
Deeper root growth is important for maintaining crop yields, especially under drought stress conditions. It has been reported that low water potential at soybean vegetative stages will decrease or stop its shoot growth, while the root continues to grow . This shoot-root response in soybean under water-deficit conditions allows the plant to search for additional soil water while maintaining higher water use efficiency . However, the limited vegetative plant growth and reduced levels of photosynthesis under drought stress negatively affect the soybean yield potential. In soybeans, short-term and moderate drought stress during vegetative growth stages generally does not impact the soybean yield in comparison with the longer-term and severe drought stress. Soybean genotypes with the ability to develop a larger root system before the occurrence of drought stress at the flowering stage were reported to have higher yield performance . Root growth was less affected when drought was imposed at the reproductive R4 stage and ceased at the R5 stage . The taproot elongation rate is the major factor influencing the depth of the soybean rooting system under water-deficit conditions . These results indicated that the lack of soil moisture at critical stages of growth will profoundly impact the productivity and that developing a large root system at the early vegetative (V) growth stage will place the soybeans in an excellent position to maintain turgor under drought conditions.
Therefore, understanding the soybean taproot response to drought is very critical for effective management of abiotic stress. Consequently, understanding the genetic regulation of the taproot response to drought will help to identify specific genes and metabolic pathways for either gene-based marker selection or genetic engineering to develop soybeans with better root-related traits. Accordingly, characterization of water-deficit induced changes in transcripts within the growth zone of the primary root is important to understand the mechanisms that control the response of root growth to water deficits. Until now, the molecular events associated with drought stress in soybean taproots are not well known and few studies have analyzed differences in transcriptional responses to different drought conditions in the same tissue or whole plant. Investigation and understanding of the responsible genomic mechanisms for drought stress will be useful to improve varieties with yield protection in soybean breeding programs. The novel genes involved in root development under drought conditions through expression pattern analysis can also be used for further functional genomic studies.
With the advent of next-generation sequencing (NGS) methods, Illumina/Solexa’s sequencing technology has been used for understanding the complexity of gene expression and regulation networks in several plant species responding to abiotic stresses, including cotton , Chinese cabbage , chickpea , maize , Ammopiptanthus mongolicus , and soybean . The research in soybean response to drought was investigated in soybean slow-wilting (PI 416937, PI 471938 and PI 567690) and fast-wilting lines (Benning, Hutcheson and Pana) [15–17]. Both genotypic and non-genotypic differential genes response to water deficit were identified, which extends the current understanding of plant hydraulic conductivity. Several high potential candidate genes were also identified to elucidate the mechanism of slow-wilting. However, all of these experiments used leaf tissue to profile the gene expression pattern under water-deficit or high vapor pressure deficit (VPD) conditions. Recently, the root transcriptomes of DT2008 and W82 soybean seedlings under short-term (2 h or 10 h) dehydration conditions were analyzed by use 66 K Affymetrix Soybean Array GeneChip . The differential expression of genes (osmoprotectant biosynthesis, detoxification or cell wall-related proteins, kinases, transcription factors and phosphatase 2C proteins) may cause higher drought tolerability of DT2008 vs. W82. This research enable us to identify early responsive genes and understand the upstream regulation mechanism of soybean root response to dehydration. However, no research has focused on how the soybean taproot responses to different drought level treatments until now.
To provide novel insights into the molecular basis of drought stress in the soybean root, genome-wide transcriptome profiling using Illumina/Solexa’s sequencing approach has been used at different drought conditions in the present work. First, the comparative analysis in the root between well-watered and drought stress conditions showed several key metabolic and hormone pathways involved in regulating the root response under varying water deficit levels. Secondly, the TF analysis led to the identification of key molecular regulators to improve drought tolerance/avoidance. These findings will benefit the elucidation of mechanisms for temporal and spatial regulation of genes in soybean under drought conditions.
Phenotypic response of the soybean shoot and leaf to varying water deficit treatments
Identification of differentially expressed genes (DEGs) under different water-deficit stresses in the soybean taproot
Transcripts associated with metabolism were down-regulated in soybean roots under mild and severe water-deficit treatments
All of these results together demonstrated that water stress has negative effects on the metabolic processes of soybean roots and the down-regulated metabolism may slow down root growth in mild and severe drought treatments.
Carbohydrate metabolism in soybean roots under water-deficit stress
It is known that carbohydrates accumulate in plants as a response to drought , but the alteration in metabolites and pathways specific to taproot responses still remain unknown. Differential expression of genes involved in sucrose and starch metabolism after imposition of water-deficit stress were found in the current study.
The hormone regulation network under mild and severe water-deficit treatments in soybean roots
Hormones play a central role in integrating the environmental and intrinsic cues during plants adapt and grow under suboptimal treatments . To elucidate the key hormones associated with root development under water-deficits, different hormone regulations of transcripts were investigated in this study. In total, 37 transcripts associated with auxin synthesis and signaling pathway genes were found to be regulated by drought. Among them, genes related to the auxin biosynthesis pathway (IAA-leucine resistant-like gene 66, indole-3-acetic acid-amido synthetase3.1 and 3.6, IAA-amino acid hydrolase IAR3) were found to be down-regulated upon water-deficit and showed up-regulation after re-watering. Concurrently, two auxin transporter genes (PIN2 and PIN5) that were down-regulated during drought treatments, exhibited an increase in the transcript levels during recovery. Several auxin responsive genes were detected to be regulated differently at varying levels of drought stress (Additional file 4). These results identified several key auxin signal and transporter genes that were regulated by water-deficit in the taproot.
It has been reported that the effect of water-deficiency on ethylene synthesis was linked with the plant’s rate of stressed [21, 22]. Rapid induction of water-deficit stress induces ethylene production and a slower imposition of water-deficit stress inhibits ethylene biosynthesis [21, 22]. In our results, 102 transcripts associated with the ethylene pathway were regulated, and most of the ethylene signal transduction and synthesis related genes were down-regulated in MS and SS treatments, but resumed at the recovery phase. They were mainly grouped into subgroup1, subgroup2 and subgroup3 (genes up-regulated or constantly under VMS and SR, down-regulated or constantly under MS and SS, Additional file 5). The ACC oxidases, ACO1 and ACO4, which catalyze the last step of the ethylene biosynthesis pathway, were found to be inhibited by water-deficit stress. Several ethylene response factors (ERFs) and ACC synthase (ACS) genes in ethylene pathway were also found to be down-regulated (Additional file 5). These results suggested the stress treatment in this study was a slow progression of drought induction. Ethylene synthesis and signal transduction played an important role in the soybean primary root growth under the water stress treatments.
Genes potentially involved in other hormonal pathways were also identified in the present study: abscisic acid (ABA, 22 genes), gibberellin (GA, 25 genes), jasmonic acid (JA, nine genes), salicylic acid (SA, nine genes) and cytokinin (16 genes) (Additional files 6 and 7). No brassinosteroid related genes were found. Although ABA is the major hormone known to be involved in drought stress response, there were only a few ABA related genes that changed in expression level under water deficient treatments. Among them, the transcripts of ABA synthesis genes,NCED (9-cis-epoxycarotenoid dioxygenase) and CCD8 (carotenoid cleavage dioxygenases), increased with respect to the level of drought imposed. This trend was found to be reversed after water recovery. Several genes involved in the crosstalk between the two hormones were also found here. Two GA synthesis related genes, GA2OX1 and GA2OX2 (gibberellin 2-oxidase), were found to be down-regulated, which were also involved in the ethylene pathway. CCD7 and CCD8 (carotenoid cleavage dioxygenases), which function in root branching by involving ABA and the strigolactones synthesis pathway, were shown to change expression patterns under water-deficit treatments. Another gene, GmJAR1 (JASMONATE RESISTANT 1), which is involved in both the auxin and JA pathways, was also inhibited during the drought stress treatment and recovered after re-watering. All of the above results indicated that the ethylene and auxin synthesis and signal transduction play a more important role than other hormones in the soybean taproot response to water deficient stress. On the other hand, hormone crosstalk probably plays a significant role in soybean root responses to drought.
Transcription factor (TF) modulation network in the soybean root
The TFs specifically regulated in both MS and SS treatment were selected for further gene family analysis (Additional file 9). Genes belonging to the bHLH (basic helix-loop-helix), AP2 (activating Protein 2), MYB (myeloblastosis) and WRKY (is a class of DNA-binding proteins that recognize the TTGAC(C/T) W-box elements) families represent most of the differentially expressed TFs. MYB represents the highest number of transcripts under drought followed by WRKY and bHLH. Six TFs were found to change under every stress level and they probably play the conserved or fundamental roles in taproot responses to stress. One of them, a heat shock related TF (Glyma04g05500), was up-regulated in all stress levels. Its homologue was found to regulate the heat stress response in Arabidopsis . Another ERF (Ethylene response factor, containing the AP2 domain) TF was down-regulated. Two MYB15 TFs and one U-box TF were found to be up-regulated under VMS and SR treatments and were down-regulated under MS and SS treatments. In Arabidopsis, the AtMYB15 gene was found to be involved in the ABA, ethylene, and jasmonic acid mediated signaling pathways. Overexpression of AtMYB15 showed an improved survival rate and reduced water loss under water deficit treatments .
Regulation of TFs in soybean whole plant responses under drought stress
The whole plant acts as a system to respond differently to stress treatments with spatio-temporal transcript regulation. Transcriptional control of stress-response gene expression is a crucial component of plant responses to a range of varying levels of environmental stresses. Understanding the regulation of TFs at a whole plant level during water deficit treatments will enable identification of key master regulatory elements involved in plant responses. The regulatory mechanisms of drought responses were specifically identified in the root water-deficit response and their regulation patterns were also evaluated in leaf and shoot tissue transcriptomes.
Conserved motif analysis in the promoters of regulated genes in the soybean root
The promoter motifs, which are conserved in clusters of co-expressed and functionally related genes, may be involved in mediating coordinated gene activity [25, 26]. Over-represented sequences analyses were conducted on the promoter regions (2,000 base pairs upstream of the start codon) and all of the enriched motifs from the different drought response groups are listed in Additional file 11.
In this analysis, two motifs (AAAAAAAA and TATATATA) over-represented in each group were identified. The second motif was normally found approximately 25–35 base pairs upstream of the start site. Several known motifs involved in the known drought response pathway or BR pathway were also identified, such as CRTGCAY (ABI3), GCCRCS (ERF1), CACRC and ACGHGK (BES1), CTCYCYC (ABI4), and GMCACGB (ABF1, BZR1, and BES1). This result indicated that those enriched motifs were probably involved in the drought regulations network induced at different stress levels.
Because groups 2 and 6 (group2: genes up-regulated under VMS and SR treatments, but down-regulated under MS and SS treatments; group6: genes up-regulated under MS and SS treatments, but down-regulated under VMS and SR treatments;) stand for oppositely regulated mechanisms in response to drought in the taproot tissues, motif analysis on genes from these two groups in the root tissues were further assessed. Fifty-seven over-represented motifs were identified in the MS and SS up-regulated drought response genes group (Group 2) and 49 over-represented motifs were isolated in the MS and SS down-regulated drought response genes group (Group 6). Only one motif (CACGHG) overlapping in these two groups. The CACGHG motif was a “G-box”, which is the binding site of the GBF4 TF in Arabidopsis. It was also conserved in the PIF, BES ((BRI1-EMS-SUPPRESSOR 1), and MYC TF target element. Enriched sequences regulated by ARR10, ANT, and Dof3 TFs were found only in Group 2 and not in Group 6. The ARR10 TF is involve in the cytokinin signaling pathway and often shows a down-regulated response to abiotic stress . ANT encodes an AP2-like ethylene-responsive TF and its function is most closely related to floral development . It has been reported that the expression pattern of the DOF TF responds to abiotic stress in several species [29, 30]. Two enriched cis-elements (CWCCACS and CAYSCAC) were found in Group 6, which were regulated by TF GAMYB. It has been reported that GmGBP1, which is one soybean GAMYB binding protein, regulates flowering and improves drought tolerance in Arabidopsis . These results indicated that different enriched motifs were involved in very specific regulation mechanism responses to different drought levels.
qRT-PCR (Real-time quantitative PCR) validation of differentially expressed transcripts from RNA-Seq
To validate the RNA-Seq expression data and its reliability, 18 DEGs were selected for qRT-PCR analysis. To compare these two different methods, the relative expression measurement from the qRT-PCR was transformed into fold change by base 2 to match with the RNA-Seq fold change value.
After comparing the results between qRT-PCR and RNA-Seq, the expression patterns from the qRT-PCR were highly consistent with the RNA-Seq results. Furthermore, a correlation between RNA-Seq and qRT-PCR was evaluated using log2 expression levels. As shown in Additional file 12, the qRT-PCR measurements were highly correlated with the RNA-Seq results (Y = 0.8476x + 0.2348, R2 = 0.8346). The comparison of data between different methods showed the accuracy of the RNA-Seq method and showed that our RNA-Seq data were accurate and efficient. The RNA-seq results can be used for gene expression profiling during the soybean root response to water-deficit stress. Notably, a few genes showed different expression patterns between the RNA-seq and qRT-PCR data. Such inconsistency was probably caused by the relatively low expression of these genes and by no significant changes being found in the RNA-seq results.
The right root system architecture will allow plants to survive from a water deficit environment with an increasingly deeper, more vertically orientated root system to access the deep soil water. Spollen et al. reported that the elongation rates of soybean primary roots were less sensitive than the shoot and the growth of the shoot had more inhibition than the root under the same water potential . Our phenotypic data from the shoot indicated that there is a good correlation between root water uptake capacity and shoot dry weight under water-limited conditions. This study unraveled the signaling networks, including hormone signaling, metabolism, and TFs in soybean primary root responses to the varying water-deficit conditions. Exploring the complex crosstalk between different regulatory levels involved in the soybean stress response will promote understanding of the regulation network of drought stress.
Metabolic changes in soybean root under water deficit conditions
It has been known that a variety of osmolytes which includes sugars (such as raffinose, trehalose, sucrose and sorbitol), sugar alcohols (such as mannitol), amino acids (such as proline), and secondary metabolites (such as isoflavonoids and saponins) will accumulate in the plants when it responses to the abiotic stress, [33–35]. These metabolism products are able to stabilize proteins and cellular structures to maintain cell turgor by osmotic adjustment. Redox metabolism removes excess levels of ROS (Reactive oxygen species) and re-establishes the cellular redox balance [36, 37]. The type of metabolism and related genes involved in the adjustment of soybean roots under water deficit conditions are still unclear. Galactinol and raffinose act as osmoprotectants and provide an adaptation to the water stress conditions . Arabidopsis plants over-expressing GolS1 and GolS2 (galactinol synthase) accumulated high levels of galactinol and raffinose and were tolerant to drought and salinity stress [38, 39]. Several GOLS1 and GOLS2 genes that were up-regulated under the VMS, MS, and SS conditions were inhibited at the water-recovery phase in the present study. These candidate genes will be very promising and can be used for soybean transgenic development to improve soybean tolerance to water deficit conditions.
Starch serves as a primary carbohydrate source and can be mobilized to soluble sugars if the plants experience an environmental challenge. It is reported that the starch content will be depleted by drought and the soluble sugars will accumulate in the leaves by enhanced beta-amylase enzyme activity [40, 41]. Here, we found the reverse adjustment for the starch metabolism in the soybean root, which indicates that metabolites involved with energy production and growth were increased or shifted from shoots to roots. The total metabolism is accelerated to adapt to the water-deficit environment when water is a limiting factor. It is reported that the shoots and roots respond to drought in an opposite metabolic way, while shoots decrease their growth metabolism (lower concentrations of sugars, amino acids, nucleosides, nitrogen, phosphorus and potassium) and roots increase metabolism in a mirrored response . BMY8 (β-amylase) antisense Arabidopsis plants accumulated high starch levels during cold shock (4 °C for 6 h) . Three BMY8 genes were also found down-regulated under the MS and SS conditions in the soybean root, which indicates the GmBMY8 gene probably plays a role under drought stress by modulating starch metabolism.
As important membrane components, the changes in the composition of lipids may help to maintain membrane integrity and preserve cell compartmentation under water stress conditions. In response to drought, the total leaf lipid content decreased progressively. The time-course of the decrease in Arabidopsis leaf lipid content correlated well with the expression of genes involved in lipid degradation . Most of the lipid biosynthesis related genes were down-regulated in the root in this study. In transgenic Arabidopsis with suppression of the inositol polyphosphate 1-phosphatase (IPP) signal by heterologous expression of the human type I inositol polyphosphate 5-phosphatase exhibited an increased drought tolerance phenotype, which was mediated in part via a DREB2A-dependent pathway . The content and composition may change when plants are exposed to different abiotic and biotic stresses . Little is known about the effects of water-deficit on lignin synthesis in roots until now. Research has shown that there was an increased expression of genes involved in lignin biosynthesis during the intermediate and final stages of water stress (from 48 h to 72 h) in rice roots (Oryza sativa L.), such as those coding for PAL, C3H, 4-coumarate: coenzyme A ligase (4CL), caffeoyl coenzyme A O-methyltransferase (CCoAOMT), cin-namyl alcohol dehydrogenase (CAD), and peroxidase . Similarly, lignin synthesis-related proteins were induced mostly in the later stage of drought stress in the roots of wild watermelon (Citrullus lanatus sp.), which may function in the enhancement of physical desiccation tolerance and drought adaptation . These results indicated that growth reduction and tolerance to water-deficit were associated with more lignin production in the roots. However, most of the lignin biosynthesis related genes in the soybean root showed a down-regulated pattern under the MS and SS conditions. The reason for this phenomenon may be caused by different regions of root and stress levels or degrees.
The analysis of metabolism related genes contributes to the understanding of soybean stress biology. Genetic engineering using genes encoding components of drought stress-related metabolic pathways has shown the potential to enhance drought resistance in soybeans [49–51]. Integration of proteomics and metabolomics analysis will provide more evidence to elucidate the nature of soybean water deficit stress responses in the future.
Hormone interaction involved in the primary root development under water-deficit
Three hormones,auxin, ethylene and ABA, were identified to be involved in drought responsive pathways in soybean . In this study, ethylene and auxin synthesis and signal transduction had more important roles than other hormones in the soybean taproot response to water deficient stress. At present, hormone interactions involved in abiotic stress were mainly derived from Arabidopsis thaliana. Until now, it was unknown how the hormone crosstalk regulated primary root elongation under water-deficit stress in soybean roots. Several genes (i.e., GA2OX1, GA2OX2, ACS, CCD7, CCD8, and JAR1) involved in different hormone crosstalk and related TFs were found in this study (Additional files 4, 5, 6 and 7), which indicated that the interactions between different hormones play a crucial role in soybean root development.
The expression of genes associated with the auxin transporter has been shown to be regulated by ethylene . Conversely, auxin was found to affect the synthesis of ethylene . Here, the expression of PIN2 and PIN5 were found to be regulated by the water-deficit stress, and several ACS genes (ACS1, ACS6, ACS7, ACS8, and ACS10) were regulated under MS and SS conditions. CCD7 and CCD8 genes were involved in the ABA and strigolactones (SL) synthesis pathway. SL, as plant hormone, functions as a positive regulator of primary root elongation and as a negative regulator of adventitious root formation due to modulation of auxin flux in the root [55–57]. Wang et al. reported that ABA treatment increased soybean CCD7 and CCD8 transcription . Based on this study, hormone balance between SL and ABA might play an important role in primary root development under MS and SS treatments. Under drought stress, higher levels of ABA may up-regulate the activities of downstream TFs (MYBs), which would modify the auxin biosynthesis and trigger auxin responsive cascades, controlling root development. The increase in ABA content and signal transduction promotes the synthesis of ethylene, then inhibits the elongation of root. This complex network of hormone crosstalk allows cell differentiation to balance cell division and the formation of the root system architecture.
All of the above results give clues about transcripts associated with hormones that modulate root growth under drought stress. While progress has been made in identifying how water deficit affects hormone biosynthesis, response, and transport, the correlation between the content of hormones in the soybean taproot and the expression levels of genes are not well understood. It will be interesting and useful to investigate changes in hormone content in soybean taproots under water-deficit conditions in the future.
Water-deficits induced TFs in soybean root
TFs are important classes of genes that regulate expression of downstream drought-responsive genes. In this study, several families, which include MYB, AP2, bHLH, WRKY, and Cys2His2 zinc-finger, were identified to be involved in soybean water-deficit responses either from tissue specific or from stress specific TFs. Among these TFs, MYB, MYC, AP2, and HD-ZIP families play a central role in drought tolerance [59–62]. These results predicated that soybean probably shares a conserved regulation mechanism with other species.
TFs DREB1A and DREB1D were found to be highly up-regulated by MS and SS conditions in the root. In transgenic Arabidopsis, overexpression of the DREB1/CBF family genes showed increasing tolerance to freezing, drought, and high salinity [62–64]. Overexpression of DREB/CBF TFs has been reported to enhance drought tolerance in soybeans [65, 66]. Drought is often accompanied by elevated air and leaf temperatures in the field. Heat shock TF (HSF) genes have been implicated not only in thermotolerance but also in plant growth and development and response to various stresses . Six HSF transcription factors were found highly induced under the MS and SS conditions in soybean roots. After water recovery, the expression level decreased compared with SS treatment (Additional file 9). Several auxin, salt, and cold pathway related TFs were found in the DEG group in the root tissue. ICE1 encodes an upstream bHLH TF that regulates the transcription of CBF genes in the cold . One soybean ICE TF (Glyma06g17330.1) was induced by the MS and SS treatments. Two auxin response TFs were affected; IAA26 (Glyma07g01800.1) was up-regulated and IAA19 (Glyma13g43780.1) was down-regulated by the stress. The induced IAA 26 in rice leaves could dynamically regulate the phytohormone IAA level to respond to drought stress . Two PIF3 (phytochrome-interacting factor) genes (Glyma19g40980.1 and Glyma10g28290.1) were induced by the MS and SS treatments in the root. It was reported that the rice PIF-like protein Os-PIL1/PIL13 was a key regulator of reduced internode elongation in rice under drought conditions . Overexpressing the STO (salt tolerance) gene in Arabidopsis conferred a higher salt tolerance phenotype than the wildtype through binding to a MYB TF . One soybean STO (Glyma11g13570.1) gene was found induced under the MS and SS conditions, but was down-regulated after water recovery. Two GA pathway TFs; SCARECROW (SCR) and RGA-LIKE 2 della protein (RGL2), were also found to be down-regulated in all levels of drought stress. SCR is known to be a key regulator of stem cell renewal and radial patterning in the root of Arabidopsis [72, 73].
These results indicate that the water-deficit response in soybean roots may be involved in many different abiotic and hormone signaling pathways that modify the root growth. Further molecular work needs to be conducted to evaluate the importance of these TFs in the water-deficit response and to determine the role of individual genes. Investigating the transcriptional regulatory network of DEGs involved in different tissues at different water deficit stress levels would provide more information for future functional analysis.
In conclusion, this study provides a large-scale transcriptome analysis of soybean taproot, leaf, and shoot in response to different water-deficit levels and demonstrated the usefulness of the DEG approach for identifying key genes across tissues and stress levels. The DEGs involved in metabolism, hormone, and regulation networks were identified and proved to be associated with the development of soybean roots under drought. Our gene expression analysis results not only provided detailed and accurate information for characterizing the drought stress response in soybeans but also provided candidates for genetic engineering for stress resistance in soybeans. Genetic modifications of these candidates could play a vital role in achieving drought tolerance and yield performance.
Plant growth conditions and stress treatment
Williams 82 soybean seeds were planted in 5-gallon pots filled with soil media (turface and sand mixed in 2:1 ratio) in the growth chamber at the Sears Plant Growth Facility at the University of Missouri. The controlled conditions of the growth chamber was 27/21 °C day/night temperature, photoperiod of 16/8 day/night, 60 % relative humidity. Three seeds were sown per pot and thinned to one plant per pot when the plants had two sets of unfolded trifoliate leaves (Vegetative (V) 2 stage). All pots were kept well-watered until V3 (three unfolded trifoliate leaves) stage . For very mild stress (VMS), drought stress was imposed by withholding water for 5 days. For mild stress (MS), drought stress was imposed by withholding water for 12 days. For severe stress (SS), drought stress was imposed by withholding water for 19 days. For water recovery after severe stress (SR), plants were re-watered for 2 days after withholding water for 19 days. Each treatment had corresponding control plants and all control plants were kept well-watered until sampling. The soil moisture of each pot was measured on a daily basis to monitor the treatment condition by using the PR2 soil moisture profile probe (Delta-T Devices Ltd., UK). The value recorded at the sampling day was used for Fig. 1b. Plant height, water potential, stomatal conductance, canopy temperature, and leaf relative water content (RWC) were recorded or determined at the sampling day in all control and stressed plants. These physiological traits were collected from an attached leaflet of the fourth trifoliate leaf from the main-stem apex at midday. Full leaves were used in the determination of RWC. Water potential was measured using a pressure chamber (PMS Instrument Co., Albany, CA, USA). A LI-COR 6400XT portable photosynthesis system (LI-COR, Lincoln, NE) was used to measure net stomatal conductance. Canopy temperatures were collected with an infrared meter (FLIP Ex severe infra-red meter, Everett, WA). RWC was determined by the ratio of tissue fresh weight to tissue turgid weight. This is termed relative tissue weight.
Five-centimeter long root tissues from the root tip, the first fully expanded trifoliate, and the entire shoot bearing third trifoliate leaf were collected along with their respective controls. Three individual plants were used for tissue sampling, measurement, and RNA isolation. Three independent biological replicates were used to reduce technical and biological variability. The collected tissues were rinsed with RNase-free water immediately to remove sand and turface, and then blotted to remove excess water. Each samples were wrapped with aluminum foil, and flash frozen in liquid nitrogen and stored in a -80 °C freezer for further RNA isolation.
RNA isolation, library construction, and RNA sequencing
Total RNA was isolated from the roots with the RNeasy Plant mini kit (Qiagen, Cat#:74904) following the manufacturer’s protocols. On-column DNase digestion with the RNase-Free DNase set were used to remove the DNA contamination (Qiagen, Valencia, CA, Cat#:79254). The RNA quantity was checked with a NanoDrop Spectrometer (ND-1000 Spectrophotometer, Thermo Scientific). The RNA integrity was analyzed using the Agilent 2100 Bioanalyzer (RNA Nano Chip, Agilent, Santa Clara, CA, USA). The RNA integrity number value ≥7 was used to select the best RNA for library construction. Total RNA, 2 μg from each sample, was used for library construction. The RNA TruSeq Stranded mRNA LT sample prep kit (Illumina, San Diego, CA, RS-122-2101) was used to prepare the RNA-Seq library following the manufacturer’s protocols. A Qubit quantitation assay and the NGS Fragment Analysis (Agilent 2100 Bioanalyzer, Santa Clara, CA) were conducted on each library. Single end reads were generated by the Illumina HiSeq 2000 (read length 1 × 100 base; Illumina, Inc. San Diego, CA).
Mapping of RNA-Seq reads
The initial base calling and quality filtering of the reads generated with the Illumina analysis pipeline (Fastq format) were performed using a custom Perl script and the default parameters of the Illumina pipeline (http://www.illumina.com). Additional filtering for poor-quality bases was performed using the FASTX-toolkit available in the FastQC software package (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). High-quality mRNA-Seq reads were aligned to the Glycine max reference genome (Gmax1.1 version) and Phytozome v9.0 gene model release using Tophat (version 140 1.4.1) . The genome indexes for Tophat were built using the Bowtie build command of Bowtie (version 0.12.7) with the reference genome file as input . Tophat was then run with default parameters and a reference GTF file using the –G option, and replicates of each treatment/sample were mapped independently to improve alignment sensitivity and accuracy for further analysis.
Sequence data analysis and differential counting
Only uniquely mapped reads were used in the analysis. The gene expression (FPKM) levels were estimated using Cufflinks software (version v2.1.1)  while differential gene expression analysis was performed using Cuffdiff (version v2.1.1) among the different sample comparisons. Only the genes with more than a 2-fold change and a p-value less than 5*10−5 were considered as significant differentially expressed genes. We used the R value and the absolute value of log2 ratio > 2 as the threshold to judge the significance of gene expression differences. The analyzed RNA-Seq gene expressions and the differential gene expression results are available for browsing in Soybean Knowledge Base (SoyKB)  at http://soykb.org. The datasets (SoyKB/Information/Dataset 17) can be browsed using the differential expression analysis suite tools, which provide access to the gene lists, Venn diagrams, Volcano plots, and function for pathway analysis.
Pathway identification and TF analysis
The different levels of water-deficit regulated genes were clustered using dChip software and the response pathways were plotted by MapMan . Multiple biological or metabolic pathways were plotted together with the mapped gene intensity of fold change by a green and red schema. The statistical cutoff of the fold change was two.
Promoter sequences (2,000 bp upstream) of the translation site of DEGs were extracted from the Phytozome database (http://www.phytozome.net/). The presence and abundance of the known cis-elements were analyzed using the DREAM tool  in MEME (based on E value <0.05). Then, similar motif searches were made from the JASPAR CORE (2014) plants database in Tomtom tool  (another MEME suite program).
Quantitative real-time PCR analysis
Quantitative real-time PCR was performed on an ABI7900HT detection system (Life Technologies, Foster City, CA) using the Maxima SYBR Green/ROX qPCR Master Mix (2×) (Cat# K0223, Thermo, USA) following the manufacturer’s protocol. PCR amplification was achieved by the following program: 50 °C for 2 min, 95 °C for 10 min, and then 40 cycles of 95 °C for 15 s, 60 °C for 1 min. The comparative Ct method for quantification was used to quantify the relative expression of specific genes . Actin (Glyma18g52780) was selected as an internal control to normalize gene expression. All primers were designed using Primer3 web-interface (http://frodo.wi.mit.edu/primer3/input.htm; ). The primer sequences are listed in Additional file 13. The reactions were performed with three biological replicates and repeated once for technical replicate.
Availability of supporting data
The data sets supporting the results of this article are included within the article (and its additional files). Illumina sequences were deposited to the NCBI Sequence Read Archive under the accessions SRP067593 (http://www.ncbi.nlm.nih.gov/sra/SRP067593).
activating Protein 2
carotenoid cleavage dioxygenases 7
carotenoid cleavage dioxygenases 8
differentially expressed genes
ethylene response factors
gibberellin 2-oxidase 1
JASMONATE RESISTANT 1
Heat shock TF
leaf water potential
Missouri Soybean Merchandising Council
RGA-LIKE 2 della protein
relative water content
very mild stress
high vapor pressure deficit
a class of DNA-binding proteins that recognize the TTGAC(C/T) W-box elements
We thank Theresa Musket for proof reading. Funding for this research was provided by the United Soybean Board and the MSMC (Missouri Soybean Merchandising Council).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Stacey G, Vodkin L, Parrott WA, Shoemaker RC. National science foundation-sponsored workshop report. Draft plan for soybean genomics. Plant Physiol. 2004;135:59.PubMedPubMed CentralView ArticleGoogle Scholar
- Sprent JI. Nodulation in legumes. Kew: Royal Botanic Gardens; 2001.Google Scholar
- Dogan E, Kirnak H, Copur O. Deficit irrigations during soybean reproductive stages and CROPGRO-soybean simulations under semi-arid climatic treatments. Field Crop Res. 2007;103:154.View ArticleGoogle Scholar
- Yamaguchi M, Sharp RE. Complexity and coordination of root growth at low water potentials: recent advances from transcriptomic and proteomic analyses. Plant Cell Environ. 2010;33:590.PubMedView ArticleGoogle Scholar
- Lenssen A. Soybean response to drought. Integrated Crop Management News Online. Iowa State University Extension. 2012. http://crops.extension.iastate.edu/cropnews/2012/06/soybean-response-drought. Accessed 22 Jun 2012.
- Hirasawa T, Tanaka K, Miyamoto D, Takei M, Ishihara K. Effects of pre-flowering moisture deficits on dry matter production and ecophysiological characteristics in soybean plants under drought treatments during grain filling. Jpn J Crop Sci. 1994;63:721.View ArticleGoogle Scholar
- Hoogenboom G, Peterson CM, Huck MG. Shoot growth rate of soybeans as affected by drought stress. Agron J. 1987;79:598.View ArticleGoogle Scholar
- Taylor HM, Burnett E, Booth GD. Taproot elongation rates of soybeans. Z. Acker Pflanzenbau Bd. 1978;146:33.
- Wang G, Zhu Q, Meng Q, Wu C. Transcript profiling during salt stress of young cotton (Gossypium hirsutum) seedlings via Solexa sequencing. Acta Physiol Plant. 2012;34:107.View ArticleGoogle Scholar
- Yu S, Zhang F, Yu Y, Zhang D, Zhao X, Wang WH. Transcriptome profiling of dehydration stress in the Chinese cabbage (Brassica rapa L. ssp. pekinensis) by tag sequencing. Plant Mol Biol Rep. 2012;30:17.View ArticleGoogle Scholar
- Wang X, Liu Y, Jia Y, Gu H, Ma H, Yu T, et al. Transcriptional responses to drought stress in root and leaf of chickpea seedling. Mol Biol Rep. 2012;39:8147.PubMedView ArticleGoogle Scholar
- Shen Y, Zhang Y, Chen J, Lin H, Zhao M, Peng H, et al. Genome expression profile analysis reveals important transcripts in maize roots responding to the stress of heavy metal Pb. Physiol Plant. 2013;147:270.PubMedView ArticleGoogle Scholar
- Zhou Y, Gao F, Liu R, Feng J, Li H. De novo sequencing and analysis of root transcriptome using 454 pyrosequencing to discover putative genes associated with drought tolerance in Ammopiptanthus mongolicus. BMC Genomics. 2012;13:266.PubMedPubMed CentralView ArticleGoogle Scholar
- Fan XD, Wang JQ, Yang N, Dong YY, Liu L, Wang FW, et al. Gene expression profiling of soybean leaves and roots under salt, saline-alkali and drought stress by high-throughput Illumina sequencing. Gene. 2013;512:392.PubMedView ArticleGoogle Scholar
- Shin JH, Vaughn JN, Abdel-Haleem H, Chavarro C, Abernathy B, Kim KD, et al. Transcriptomic changes due to water deficit define a general soybean response and accession-specific pathways for drought avoidance. BMC Plant Biol. 2015;15:26. doi:10.1186/s12870-015-0422-8.PubMedPubMed CentralView ArticleGoogle Scholar
- Prince SJ, Joshi T, Mutava RN, Syed N, Joao Vitor Mdos S, Patil G, et al. Comparative analysis of the drought-responsive transcriptome in soybean lines contrasting for canopy wilting. Plant Sci. 2015;240:65.PubMedView ArticleGoogle Scholar
- Devi MJ, Sinclair TR, Taliercio E. Comparisons of the effects of elevated vapor pressure deficit on gene expression in leaves among two fast-wilting and a slow-wilting soybean. PLoS One. 2015;10(10):e0139134.PubMedPubMed CentralView ArticleGoogle Scholar
- Ha CV, Watanabe Y, Tran UT, Le DT, Tanaka M, Nguyen KH, et al. Comparative analysis of root transcriptomes from two contrasting drought-responsive Williams 82 and DT2008 soybean cultivars under normal and dehydration conditions. Front Plant Sci. 2015;6:551.PubMedPubMed CentralView ArticleGoogle Scholar
- Seki M, Umezawa T, Urano K, Shinozaki K. Regulatory metabolic networks in drought stress responses. Curr Opin Plant Biol. 2007;10:296.PubMedView ArticleGoogle Scholar
- Peleg Z, Blumwald E. Hormone balance and abiotic stress tolerance in crop plants. Curr Opin Plant Biol. 2011;14:290.PubMedView ArticleGoogle Scholar
- Xu CC, Zou Q. Effect of drought on lipoxygenase activity, ethylene, and ethane formation in leaves of soybean plants. Acta Bot Sin. 1993;35:31.Google Scholar
- Morgan PW, Drew MC. Ethylene and plant responses to stress. Physiol Plant. 1997;100:620.View ArticleGoogle Scholar
- Schramm F, Larkindale J, Kiehlmann E, Ganguli A, Englich G, Vierling E, et al. A cascade of transcription factor DREB2A and heat stress transcription factor HsfA3 regulates the heat stress response of Arabidopsis. Plant J. 2008;53:264.PubMedView ArticleGoogle Scholar
- Ding Z, Li S, An X, Liu X, Qin H, Wang D. Transgenic expression of MYB15 confers enhanced sensitivity to abscisic acid and improved drought tolerance in Arabidopsis thaliana. J Genet Genomics. 2009;36:17.PubMedView ArticleGoogle Scholar
- Do JH, Choi DK. Clustering approaches to identifying gene expression patterns from DNA microarray data. Mol Cells. 2008;25:279.PubMedGoogle Scholar
- Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM. Systematic determination of genetic network architecture. Nat Genet. 1999;22:281.PubMedView ArticleGoogle Scholar
- Ramireddy E, Chang L, Schmülling T. Cytokinin as a mediator for regulating root system architecture in response to environmental cues. Plant Signal Behav. 2014;9:e27771.PubMedPubMed CentralView ArticleGoogle Scholar
- Klucher KM, Chow H, Reiser L, Fischer RL. The AINTEGUMENTA gene of Arabidopsis required for ovule and female gametophyte development is related to the floral homeotic gene APETALA2. Plant Cell. 1996;8:137.PubMedPubMed CentralView ArticleGoogle Scholar
- Ma J, Li MY, Wang F, Tang J, Xiong AS. Genome-wide analysis of Dof family transcription factors and their responses to abiotic stresses in Chinese cabbage. BMC Genomics. 2015;16(1):33.PubMedPubMed CentralView ArticleGoogle Scholar
- He L, Su C, Wang Y, Wei Z. ATDOF5.8 protein is the upstream regulator of ANAC069 and is responsive to abiotic stress. Biochimie. 2015. doi: 10.1016/j.biochi.2014.12.017.
- Zhang Y, Zhao L, Li H, Gao Y, Li Y, Wu X, et al. GmGBP1, a homolog of human ski interacting protein in soybean, regulates flowering and stress tolerance in Arabidopsis. BMC Plant Biol. 2013;6:13.Google Scholar
- Spollen WG, Sharp RE, Saab IN, Wu Y. Regulation of cell expansion in roots and shoots at low water potentials. In: Smith JAC, Griffiths H, editors. Water deficits, plant responses from cell to community. Oxford: Bios Scientific; 1993. p. 37–52.Google Scholar
- Urano K, Kurihara Y, Seki M, Shinozaki K. ‘Omics’ analyses of regulatory networks in plant abiotic stress responses. Curr Opin Plant Biol. 2010;13:132.PubMedView ArticleGoogle Scholar
- Wu W, Zhang Q, Zhu Y, Lam HM, Cai Z, Guo D. Comparative metabolic profiling reveals secondary metabolites correlated with soybean salt tolerance. J Agric Food Chem. 2008;56:11132.PubMedView ArticleGoogle Scholar
- Yamaguchi M, Valliyodan B, Zhang J, Lenoble ME, Yu O, Rogers EE, et al. Regulation of growth response to water stress in the soybean primary root. I. Proteomic analysis reveals region-specific regulation of phenylpropanoid metabolism and control of free iron in the elongation zone. Plant Cell Environ. 2010;33:223.PubMedView ArticleGoogle Scholar
- Gábor K, Robert L, Gabriella S, Virág S, Lívia SS, Gábor G, et al. Genetic manipulation of proline levels affects antioxidants in soybean subjected to simultaneous drought and heat stresses. Physiol Plant. 2005;124:227.View ArticleGoogle Scholar
- Valliyodan B, Nguyen HT. Understanding regulatory networks and engineering for enhanced drought tolerance in plants. Curr Opin Plant Biol. 2006;9:189.PubMedView ArticleGoogle Scholar
- Taji T, Ohsumi C, Iuchi S, Seki M, Kasuga M, Kobayashi M, et al. Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J. 2002;29:417.PubMedView ArticleGoogle Scholar
- Nishizawa A, Yabuta Y, Shigeoka S. Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 2008;147:1251.PubMedPubMed CentralView ArticleGoogle Scholar
- Todaka D, Matsushima H, Morohashi Y. Water stress enhances beta-amylase activity in cucumber cotyledons. J Exp Bot. 2000;51:739.PubMedView ArticleGoogle Scholar
- Kaplan F, Guy CL. beta-Amylase induction and the protective role of maltose during temperature shock. Plant Physiol. 2004;135:1674.PubMedPubMed CentralView ArticleGoogle Scholar
- Gargallo-Garriga A, Sardans J, Pérez-Trujillo M, Rivas-Ubach A, Oravec M, Vecerova K, et al. Opposite metabolic responses of shoots and roots to drought. Sci Rep. 2014;4:6829.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaplan F, Guy CL. RNA interference of Arabidopsis beta-amylase8 prevents maltose accumulation upon cold shock and increases sensitivity of PSII photochemical efficiency to freezing stress. Plant J. 2005;44:730.PubMedView ArticleGoogle Scholar
- Gigon A, Matos AR, Laffray D, Zuily-Fodil Y, Pham-Thi AT. Effect of drought stress on lipid metabolism in the leaves of Arabidopsis thaliana (ecotype Columbia). Ann Bot. 2004;94:345.PubMedPubMed CentralView ArticleGoogle Scholar
- Perera IY, Hung CY, Moore CD, Stevenson-Paulik J, Boss WF. Transgenic Arabidopsis plants expressing the type 1 inositol 5-phosphatase exhibit increased drought tolerance and altered abscisic acid signaling. Plant Cell. 2008;20:2876.PubMedPubMed CentralView ArticleGoogle Scholar
- Moura JC, Bonine CA, de Oliveira Fernandes Viana J, Dornelas MC, Mazzafera P. Abiotic and biotic stresses and changes in the lignin content and composition in plants. J Integr Plant Biol. 2010;52:360.PubMedView ArticleGoogle Scholar
- Yang L, Wang CC, Guo WD, Li XB, Lu M, Yu CL. Differential expression of cell wall related genes in the elongation zone of riceroots under water deficit. Russ J Plant Physiol. 2006;53:390.View ArticleGoogle Scholar
- Yoshimura K, Masuda A, Kuwano M, Yokota A, Akashi K. Programmed proteome response for drought avoidance/tolerancein the root of a C-3 xerophyte (wild watermelon) under water deficits. Plant Cell Physiol. 2008;49:226.PubMedView ArticleGoogle Scholar
- de Ronde JA, Cress WA, Krüger GH, Strasser RJ, Van Staden J. Photosynthetic response of transgenic soybean plants, containing an Arabidopsis P5CR gene, during heat and drought stress. J Plant Physiol. 2004;161:1211.PubMedView ArticleGoogle Scholar
- de Ronde JA, Laurie RN, Caetano T, Greyling MM, Kerepesi I. Comparative study between transgenic and non-transgenic soybean lines proved transgenic lines to be more drought tolerant. Euphytica. 2004;138:123.View ArticleGoogle Scholar
- Simon-Sarkadi L, Kocsy G, Várhegyi A, Galiba G, de Ronde JA. Genetic manipulation of proline accumulation influences the concentrations of other amino acids in soybean subjected to simultaneous drought and heat stress. J Agric Food Chem. 2005;53:7512.PubMedView ArticleGoogle Scholar
- Le DT, Nishiyama R, Watanabe Y, Tanaka M, Seki M, Ham le H, et al. Differential gene expression in soybean leaf tissues at late developmental stages under drought stress revealed by genome-wide transcriptome analysis. PLoS One. 2012. doi: 10.1371/journal.pone.0049522.
- Růzicka K, Ljung K, Vanneste S, Podhorská R, Beeckman T, Friml J, et al. Ethylene regulates root growth through effects on auxin biosynthesis and transport-dependent auxin distribution. Plant Cell. 2007;19:2197.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsuchisaka A, Theologis A. Unique and overlapping expression patterns among the Arabidopsis 1-amino-cyclopropane-1-carboxylate synthase gene family members. Plant Physiol. 2004;136:2982.PubMedPubMed CentralView ArticleGoogle Scholar
- Ruyter-Spira C, Kohlen W, Charnikhova T, van Zeijl A, van Bezouwen L, de Ruijter N, et al. Physiological effects of the synthetic strigolactone analog GR24 on root system architecture in Arabidopsis: another belowground role for strigolactones? Plant Physiol. 2011;155:721.PubMedPubMed CentralView ArticleGoogle Scholar
- Rasmussen A, Mason MG, De Cuyper C. Strigolactones suppress adventitious rooting in Arabidopsis and pea. Plant Physiol. 2012;158:1976.PubMedPubMed CentralView ArticleGoogle Scholar
- Koltai H. Strigolactones are regulators of root development. New Phytol. 2011;190:545.PubMedView ArticleGoogle Scholar
- Wang RK, Wang CE, Fei YY, Gai JY, Zhao TJ. Genome-wide identification and transcription analysis of soybean carotenoid oxygenase genes during abiotic stress treatments. Mol Biol Rep. 2013;40:4737.PubMedView ArticleGoogle Scholar
- Shin DJ, Moon SJ, Han S, Kim BG, Park SR, Lee SK, et al. Expression of StMYB1R-1, a novel potato single MYB-like domain transcription factor increases drought tolerance. Plant Physiol. 2011;155:421.PubMedPubMed CentralView ArticleGoogle Scholar
- Sugano S, Kaminaka H, Rybka Z, Catala R, Salinas J, Matsui K, et al. Stress-responsive zinc finger gene ZPT2-3 plays a role in drought tolerance in petunia. Plant J. 2003;36:830.PubMedView ArticleGoogle Scholar
- Yu H, Chen X, Hong YY, Wang Y, Xu P, Ke SD, et al. Activated expression of an Arabidopsis HD-START protein confers drought tolerance with improved root system and reduced stomatal density. Plant Cell. 2008;20:1134.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu Q, Kasuga M, Sakuma Y, Abe H, Miura S, Yamaguchi-Shino-zaki K, et al. Two transcription factors, DREB1 and DREB2, with an EREBP/AP2 DNA binding domain separate two cellular signal transduction pathways in drought- and low-temperature-respon-sive gene expression, respectively, in Arabidopsis. Plant Cell. 1998;10:1391.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaglo-Ottosen KR. Arabidopsis CBF1 overexpression induces COR genes and enhances freezing tolerance. Science. 1998;280:104.PubMedView ArticleGoogle Scholar
- Kasuga M, Liu Q, Miura S, Yamaguchi-Shinozaki K, Shinozaki K. Improving plant drought, salt, and freezing tolerance by gene transfer of a single stress-inducible transcription factor. Nat Biotechnol. 1999;17:287.PubMedView ArticleGoogle Scholar
- Polizel AM, Medri ME, Nakashima K, Yamanaka N, Farias JR, de Oliveira MC, et al. Molecular, anatomical and physiological properties of a genetically modified soybean line transformed with rd29A:AtDREB1A for the improvement of drought tolerance. Genet Mol Res. 2011;10:3641.PubMedView ArticleGoogle Scholar
- de Paiva Rolla AA, de Fátima Corrêa Carvalho J, Fuganti-Pagliarini R, Engels C, do Rio A, Marin SR, et al. Phenotyping soybean plants transformed with rd29A:AtDREB1A for drought tolerance in the greenhouse and field. Transgenic Res. 2014;23:75.PubMedView ArticleGoogle Scholar
- Nishizawa A, Yabuta Y, Yoshida E, Maruta T, Yoshimura K, Shigeoka S. Arabidopsis heat shock transcription factor A2 as a key regulator in response to several types of environmental stress. Plant J. 2006;48:535.PubMedView ArticleGoogle Scholar
- Chinnusamy V, Ohta M, Kanrar S, Lee BH, Hong X, Agarwal M, et al. ICE1: a regulator of cold-induced transcriptome and freezing tolerance in Arabidopsis. Genes Dev. 2003;17:1043.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang D, Pan YJ, Zhao XQ, Zhu LH, Fu BY, Li ZK. Genome-wide temporal-spatial gene expression profiling of drought responsiveness in rice. BMC Genomics. 2011;12:149. doi:10.1186/1471-2164-12-149.PubMedPubMed CentralView ArticleGoogle Scholar
- Todaka D, Nakashima K, Maruyama K, Kidokoro S, Osakabe Y, Ito Y, et al. Rice phytochrome-interacting factor-like protein OsPIL1 functions as a key regulator of internode elongation and induces a morphological response to drought stress. Proc Natl Acad Sci U S A. 2012;109:15947.PubMedPubMed CentralView ArticleGoogle Scholar
- Nagaoka S, Takano T. Salt tolerance‐related protein STO binds to a Myb transcription factor homologue and confers salt tolerance in Arabidopsis. J Exp Bot. 2003;54:2231.PubMedView ArticleGoogle Scholar
- Di Laurenzio L, Wysocka-Diller J, Malamy JE, Pysh L, Helariutta Y, Freshour G, et al. The SCARECROW gene regulates an asymmetric cell division that is essential for generating the radial organization of the Arabidopsis root. Cell. 1996;86:423.PubMedView ArticleGoogle Scholar
- Helariutta Y, Fukaki H, Wysocka-Diller J, Nakajima K, Jung J, Sena G, et al. The SHORT-ROOT gene controls radial patterning of the Arabidopsis root through radial signaling. Cell. 2000;101:555.PubMedView ArticleGoogle Scholar
- Fehr WR, Caviness CE, Burmood DT, Pennington JS. Stage of development descriptions for soybeans, Glycine max (L.) Merr. Crop Sci. 1971;11:929.View ArticleGoogle Scholar
- Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009;25:1105.PubMedPubMed CentralView ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357.PubMedPubMed CentralView ArticleGoogle Scholar
- Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31:46. doi:10.1038/nbt.2450.PubMedView ArticleGoogle Scholar
- Joshi T, Fitzpatrick MR, Chen S, Liu Y, Zhang H, Endacott RZ, et al. Soybean knowledge base (SoyKB): a web resource for integration of soybean translational genomics and molecular breeding. Nucleic Acids Res. 2014;42:1245.View ArticleGoogle Scholar
- Usadel B, Poree F, Nagel A, Lohse M, Czedik-Eysenberg A, Stitt M. A guide to using MapMan to visualize and compare Omics data in plants: a case study in the crop species, Maize. Plant Cell Environ. 2009;32:1211.PubMedView ArticleGoogle Scholar
- Bailey TL. DREME: Motif discovery in transcription factor ChIP-seq data. Bioinformatics. 2011;27:1653.PubMedPubMed CentralView ArticleGoogle Scholar
- Gupta S, Stamatoyannopoulos JA, Bailey TL, Noble WS. Quantifying similarity between motifs. Genome Biol. 2007;8:R24.PubMedPubMed CentralView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2-[Delta][Delta]CT method. Methods. 2001;25(4):402–8.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000;132:365.PubMedGoogle Scholar