- Research article
- Open Access
Identification of key genes and its chromosome regions linked to drought responses in leaves across different crops through meta-analysis of RNA-Seq data
BMC Plant Biology volume 19, Article number: 194 (2019)
Our study is the first to provide RNA-Seq data analysis related to transcriptomic responses towards drought across different crops. The aim was to identify and map which genes play a key role in drought response on leaves across different crops. Forty-two RNA-seq samples were analyzed from 9 published studies in 7 plant species (Arabidopsis thaliana, Solanum lycopersicum, Zea mays, Vitis vinifera, Malus X domestica, Solanum tuberosum, Triticum aestivum).
Twenty-seven (16 up-regulated and 11 down-regulated) drought-regulated genes were commonly present in at least 7 of 9 studies, while 351 (147 up-regulated and 204 down-regulated) were commonly drought-regulated in 6 of 9 studies. Across all kind of leaves, the drought repressed gene-ontologies were related to the cell wall and membrane re-structuring such as wax biosynthesis, cell wall organization, fatty acid biosynthesis. On the other hand, drought-up-regulated biological processes were related to responses to osmotic stress, abscisic acid, water deprivation, abscisic-activated signalling pathway, salt stress, hydrogen peroxide treatment. A common metabolic feature linked to drought response in leaves is the repression of terpenoid pathways. There was an induction of AL1 (alfin-like), UGKYAH (trihelix), WRKY20, homeobox genes and members of the SET domain family in 6 of 9 studies. Several genes involved in detoxifying and antioxidant reactions, signalling pathways and cell protection were commonly modulated by drought across the 7 species. The chromosome (Chr) mapping of these key abiotic stress genes highlighted that Chr 4 in Arabidopsis thaliana, Chr 1 in Zea mays, Chr 2 and Chr 5 in Triticum aestivum contained a higher presence of drought-related genes compared to the other remaining chromosomes. In seedling studies, it is worth notice the up-regulation of ERF4 and ESE3 (ethylene), HVA22 (abscisic acid), TIR1 (auxin) and some transcription factors (MYB3, MYB94, MYB1, WRKY53 and WRKY20). In mature leaves, ERF1 and Alfin-like 1 were induced by drought while other transcription factors (YABBY5, ARR2, TRFL2) and genes involved phospholipid biosynthesis were repressed.
The identified and mapped genes might be potential targets of molecular breeding activities to develop cultivars with enhanced drought resistance and tolerance across different crops.
Abiotic stresses are playing a major limiting factor affecting crop growth, yield and productivity [1, 2]. Among environmental stresses, drought is one of the most serious and increasing environmental factors affecting agricultural production. Drought stress affects water uptake and plant adaptation and long-term evolution of plant species to climate change . To withstand in adverse environmental stress conditions like drought, the plant requires a substantial change in the metabolism, which includes regulation of transcription and gene expression and extensive transcriptome reprogramming upon the occurrence of the stress . Therefore, transcriptomic studies offer great insight into the mechanisms of plant stress responses. Among the small plant molecules, hormones play an important role in the modulation of the complex plant physiological and molecular responses to drought. Abscisic acid is the key hormone modulating water loss and cellular growth maintenance . However, this is only one among the many key players in the complex molecular networks underlying crop responses to environmental stresses. The outcome of the responses is regulated by complex crosstalk where small molecules (such as hormones) play a specific role of inhibition/induction of key proteins in stress signal reception, transmission and responses such as kinases, phosphatases, and transcription factors, defensive responsive genes . Some key transcription factor (TF) families such as MYB, WRKY, bZIPs have been found to be involved in a different manner depending on the type of stress. Some TFs have been object of genetic engineering to improve stress tolerance in model and crop plants [7, 8]. Transcriptomic studies are essential in gaining insight into the crop responses to drought by identifying specific genes involved in plant responses to water stresses highlighting peculiarities of each crop and identifying which genes are the base of diverse drought tolerance and resistance mechanism. Since data of each study are typically related to only one season, this may lead to reduced reliability of the conclusions driven by each study. Indeed, it is essential to find a pipeline to compare data across species in order to strengthen the meaning of every single study, validating published works across species and reducing the environmental variability that affects their reliability. This kind of works, named meta-analysis is lacking in crops, especially at the transcriptomic level. Therefore, it is highly desirable to put more efforts in developing extensive studies to systematically understand drought-stress-related mechanisms in crops, which will accelerate the development of new crop varieties with improved stress resistance to increase agricultural sustainability and food supply for a highly growing world population.
RNA-Sequencing (RNA-Seq) is a rapid technique for genome-wide gene expression analysis [9,10,11]. With the emergence of this technique, the high-throughput transcriptomic technologies have been revolutionized. This technique can be considered as an efficient way to identify genes and gene families encoding proteins involved in different metabolic pathways related to the object of the study. Next-generation sequencing methods have enabled to understand the gene expression data in both quantitative and qualitative manner  and can be used for obtaining sequences on a large scale with high sequencing depth. To assist crop breeding in developing drought-tolerant crops, it is crucial to gain insight into the complex networks of crop environmental stress responses by elucidating the molecular basis of drought-stress transduction pathways and drought tolerance mechanisms. Omic approaches have been used to validate RNA-Seq data related to environmental stress responses [13, 14]. However, transcriptomic studies present some drawbacks represented by the following: 1) high presence of false-positive results that requires validation with other platforms, 2) data are generally affected by environmental, experimental, developmental and genetic conditions, 3) experiments are typically not repeated and conducted in only one season, 4) data are highly affected by the environment, especially when performed in field conditions, 5) few replicates are usually performed due to the high costs of these analysis and the scarce integration between transcriptomic and other omic platforms. Meta-analysis improve the reproducibility of RNA-Seq studies because: 1) it filters the most meaningful information linked with the object of study, 2) eliminates data affected by environmental variability, 3) reduces false positive results, 4) increase the number of virtual replicates, 5) integrates multiple datasets. Meta-analysis studies should be integrated with statistical modelling used for each sample in order to unveil intrinsic mechanisms . In the present work, we have conducted a meta-analysis of 9 RNA-Seq studies conducted in 7 crops to deliver conserved and reliable genomic information exploitable by breeding to enhance drought resistance in crops. We analyzed (at a most comprehensive manner as possible) RNA-Seq data in crops (herbaceous, tree fruit crops, model plant) under drought using the same bioinformatics pipeline to deliver functional genomics knowledge that will guide molecular breeding to enhance drought tolerance and resistance in crops.
Drought transcriptomic responses in all kind of leaves
Based on the search criteria described in Methods, we found 22 RNA-Seq studies: 7 were performed in roots, 12 in leaves, 3 in fruits. Among leaf studies, 3 of them have no raw data available. Indeed, the analysis was performed using 9 studies [16,17,18,19,20,21,22,23,24]: 5 dealing with mature leaves, 1 in young leaves and the other 3 in seedlings. The 9 studies comprise of 2 fruit tree crops and 5 herbaceous ones.
The articles and crops selected for the study, number of up- and down-regulated genes were listed in Table 1. The analysis resulted in the identification of a total of 108,903 genes in which 53,988 were up-regulated and 54,915 were down-regulated. For each of the analysis, the total number of genes range from 9420 to 17,230. The number of genes up-regulated was in a range of 2866 to 8184 and down-regulated genes were span from 4582 to 9046. The results of the analysis along the gene ID, Arabidopsis thaliana ortholog and log2 fold change value are provided in Additional file 1: Table S1. The two Vitis vinifera studies form a cluster showing an overall transcriptomic similarity towards the analysis to drought (Fig. 1; Additional file 2: Figure S1). Although the transcriptomic responses in the two maize studies were very similar, the closeness of one maize study to wheat study was higher than between the two maize studies. The similarity in drought responses among apple, Arabidopsis thaliana and tomato were related due to the fact that they dealt with drought responses in seedling leaves.
Gene set and pathway enrichment analysis
DAVID software was used to identify the common biological processes affected by drought at transcriptomic level considering the drought-regulated genes in at least 6 of 9 studies. Metabolic pathways divided in up- or down-regulated by drought along with GO ID, its GO term, count, p-values and Benjamini values were shown (Table 2). No GO-terms related to the biological process were commonly drought-regulated in at least 7 of 9 studies. Among at least 6 of 9 articles, 38 GO-terms were down-regulated while 28 were up-regulated. Among them, it is worthy to mention some of the biological pathways that are known to be repressed by the drought stress such as wax biosynthesis and cell wall organization, fatty acid biosynthesis, protein phosphorylation. On the opposite, we identified some GO-terms that were up-regulated in response to water stress such as response to osmotic stress, response to abscisic acid, response to water deprivation, abscisic-activated signalling pathway, response to salt stress, response to hydrogen peroxide.
Abiotic stress responses
Genes mapped in the abiotic stress-related GO-terms identified by DAVID were shown in Fig. 2. Among the drought up-regulated genes involved in osmotic and salt stress, it is worth to mention the sucrose-related protein kinase and the CBL interacting protein kinase, the salt overly sensitive 1, and pyrophosphorylase 6. In the category of “response to water deprivation”, there was up-regulation of homeobox 7, lipid transfer protein 3, open stomata 1, calcineurin B-like protein. Four genes were up-regulated by drought and involved in “abscisic acid-activated signalling” while 8 of the drought up-regulated genes were involved in “response to abscisic acid”. Drought repressed three genes involved in fatty acid biosynthesis such as 3-ketoacyl-coa synthase 1 and 6 and 3-ketoacyl-coa thiolase 2.
Secondary metabolism, cellular responses, signalling
MapMan web-tool was used to identify transcriptomic effects of drought in key selected categories such as secondary metabolism, cellular responses and signalling (Additional file 3: Figure S2). Among the secondary metabolism, the drought-repressed genes were involved in terpene pathways such as terpene cyclase, phytoene synthase, farnesyl pyrophosphate synthase 2. Cellular response genes were mostly inhibited by drought. MADS transcription factors like AGL8 (agamous-like MADS-box) and SCL3 (scarecrow-like protein) were enhanced. Relating to signalling mechanisms, genes encoding for 2 leucin rich repeat genes, 3 protein kinases, a proline-rich extension like receptor kinase, a lectin protein were repressed. On the other hand, a protein kinase (AT5G56890), and two serine/threonine kinases were up-regulated.
Transcription factors and hormones
Among the drought-up regulated transcription factors, it is worthy to mention the induction of AL1 (alfin-like), UGKYAH (trihelix), WRKY20, zinc ion binding, two homeobox genes (one CDF2 and an SDG26 (SET domain)) (Fig. 3). Among the repressed ones, there were two bHLH members, a MYB factor (TKI1), two ABA-related TF (ABI3VP1), and ARR2 (cytokinin-related). Figure 4 summarized the drought-regulated genes involved in hormone-related categories. Ethylene and salicylic acid pathways were repressed by drought whereas auxin, abscisic acid, cytokinin, ethylene pathways were mostly up-regulated. Water deprivation down-regulated three genes responsive to ethylene and three responsive to salicylic acid (such as glutathione-s-transferase 2) while it up-regulated several genes responsive to auxin, abscisic acid, cytokinin and ethylene activated signalling pathway. Among the auxin-responsive genes it is worth to mention the enhancement of indole-3-butyric acid response 5 and the phytochrome associated protein 2. Relating to abscisic acid there was an up-regulation in homeobox 7, lipid transfer protein 3, shaker potassium ion channel, SNF1, potassium transport 3. The cytokinin responsive gene, heat shock protein 93, and three ethylene-related genes, ERF1, SKP1 and DREB were enhanced.
Protein-protein network analysis
The protein-protein interaction (PPI) network analysis comprises of 351 drought-related genes commonly regulated in at least 6 of 9 studies. Minimum default settings were used to reduce the number of interacting proteins and the complexity of the networks (Fig. 5). Some key genes with a high number of interactions (> 20) were highlighted. Among the up-regulated hub (highly interacting) proteins it is worthy to notice some key proteins that may play a key role in drought response such as LOS1 (Low expression of osmotically responsive genes 1), HSP90–4 (heat shock protein 90–4), SKP1B (SKP1-like protein 1B), CR88 (chlorate resistant). Interestingly, drought down-regulated highly interactive proteins such as ATL5 (ring H2 finger protein), UBQ3 (polyubiquitin 3), TTL1 (TP repeat-containing thioredoxin), ATJ20, (chaperone protein DNAJ 20), CDKA-1 (Cyclin Dependent Kinase A-1). PPI network analysis was performed for drought-regulated in common between the three seedling studies and between the five studies on mature leaves (Additional file 4: Figure S3, Additional file 5: Figure S4). In seedling, two major hub proteins WDR5A (histone methylase component) and ASHH1 (histone lysine N methyltransferase) were up-regulated (Additional file 4: Figure S3). In mature leaves, an LRR receptor-like serine/threonine protein kinase was repressed while CDKF-1 (cyclin dependent kinase F-1) was up-regulated (Additional file 5: Figure S4).
Chromosome mapping of key drought-regulated genes in crops
Key genes encoding genes related to abiotic stress responses, transcription factors, hormone metabolism (obtained from DAVID software) (respectively visualized in Figs. 2, 3, and 4) were mapped in the respective chromosomes of the 7 crops (Additional file 6: Figure S5; Additional file 7: Table S2). There were a total of 55 genes. Interestingly, we observed that in some species there was not a homogeneous distribution of these genes across chromosomes since some chromosomes contained a higher number of them. While in apple, potato and tomato there was a similar distribution of these genes in the chromosomes, whereas, in Zea mays, Triticum aestivum and Arabidopsis thaliana there was a higher presence of these genes in some of the chromosomes. In maize, a total of 29 abiotic stress-related genes were mapped to chromosome 1 implying that the chromosome 1 regions should contain more genes involved in drought resistance than the other chromosome regions. In Arabidopsis thaliana, 17 genes were mapped to chromosome 4. In Triticum aestivum chromosome 2 (2A + 2B + 2C + 2D) and 5 (5A + 5B + 5C + 5D genome) mapped respectively 15 and 12 genes. This work allowed to identify which chromosome might contain more genes involved in drought resistance and will guide the identification of new molecular markers linked with drought resistance.
The most common drought-regulated genes
The 27 genes that were drought-regulated in at least 7 of 9 studies were considered to be closely linked to drought responses. They were mapped in each plant genome (Additional file 8: Figure S6; Additional file 9: Table S3). Ethylene response factor 1, protein kinase 2B, WRKY20 and Alfin-like 1 were up-regulated in at least 8 of 9 studies. At the same time, there was an inhibition of ABA-responsive element binding protein 3 and Protein phosphatase 2A subunit A2 in at least 8 of 9 studies. Among these highly conserved genes, a zinc-finger, ocre domain-containing protein 1, Rho GTPase-activating gaco-like protein and serine carboxypeptidase-like 27 were up-regulated while SLY1 gene family was down-regulated (Additional file 9: Table S3).
Drought-regulated transcriptomic responses at different leaf stage
Attention was paid on the drought-responsive genes at different leaf developmental stages. Comparing the three studies dealing with drought transcriptomic responses in seedlings (tomato, Arabidopsis thaliana and apple), 934 commonly drought-regulated genes were identified (Fig. 6). On the other hand, 465 genes were commonly drought-regulated in at least 4 of 5 studies in mature leaf tissues. Finally, comparing the two lists of drought-regulated genes, 912 genes were specifically drought-regulated in seedlings, 443 in mature leaves and 22 in common between the two types of leaves (Additional file 10: Table S4 and Additional file 12: Table S5). These results demonstrated that transcriptome reprogramming in response to drought depends on different leaf developmental stage.
Drought-responsive genes in seedlings
Among the genes that were regulated by drought in seedlings, we paid our attention to those belonging to key categories playing an important role in drought response modulation such as hormones, transcription factors and abiotic defence responses (Additional file 10: Table S4). Relating to hormones, two ethylene-related genes (ERF4 and ethylene-responsive element binding protein (ESE3, AT5G25190)) were up-regulated in all the three seedling studies. In addition, there were other 5 up-regulated genes involved in auxin (TIR1, auxin-responsive protein (AT4G38840)), abscisic acid (HVA22), gibberellin (KAO2, GA4). On the other hand, there were 8 down-regulated hormone-related genes: oxidoreductase B2 and AIR9 (auxin-related), AREB3 (abscisic acid), BAK1 (brassinosteroid-related), CKX7(cytokinins), gibberellin-20-oxidase2 and gibberellin-2-beta-dioxygenase (gibberellins). Three genes involved thioredoxin pathways were also up-regulated: APRL5, PDIL5–1, ATY1. Unexpectedly there were also some heat stress-related genes repressed such as (HSP17.8, ARL1, GFA2, HSP98.7). Regarding with transcription factors, there were some key categories that were commonly up-regulated among crops such as MYBs (MYB3, MYB94, MYB1), bHLH, and homeobox. Relating to the WRKY family, two genes were up-regulated (WRKY53 and WRKY20) while one gene was repressed (WRKY22). Interestingly, the SET-domain family was mainly up-regulated.
Drought-responsive genes in mature leaves
We focused our attention on genes involved in the same categories that are commonly drought-regulated in at least 4 of 5 studies (cellular responses, hormones, transcription factors) (Additional file 10: Table S4). In total 4 hormone-related genes were drought-repressed such as an auxin-responsive (RRT4; O-fucosyltransferase family protein), two ABA-related genes (NCED4, HVA22A) and one salicylic acid-related (UDP-glucosyltransferase). ERF1, a key player in jasmonic acid-ethylene crosstalk was up-regulated by drought in mature leaves. Unexpectedly we observed that most of the genes encoding transcription factors were repressed by drought including YABBY5, ARR2, BLH6, TRFL2, three zinc finger proteins and other 5 genes. Alfin-like 1 was the only up-regulated transcription factor. Relating to another primary metabolism, it is worth notice that two genes involved in phospholipid biosynthesis were repressed (phosphatidylserine synthase and galactolipid galactosyltransferase).
Meta-analysis of transcriptomic responses to abiotic stresses has been performed comparing different studies that were conducted in the model plant Arabidopsis thaliana , rice , sunflower . Our study used a similar approach expanding the analysis across different crop species. Our aim was to shed light into drought response mechanisms conserved across crops instead of identifying specific responses in each crop. Our purpose was mainly to answer the two following unresolved questions:
Which genes and molecular mechanisms are conserved across species and can be considered strictly modulating drought responses in plants?
How leaf development affects crop molecular responses to drought and which genes are playing a key role in drought resistance at different developmental stages?
Common drought responses across plant species in all kind of leaves
We answered the first question by identifying drought-regulated genes in common across plant species. Twenty-seven genes were up- or down-regulated in response to drought in at least 7 of 9 studies. Some of them required particular attention considering that they have been previously linked to drought responses in single studies. They were ERF1 (involved in ethylene signalling), WRKY20 and Alfin-like 1, zinc finger ocre domain protein 1 (transcription factors), serine carboxypeptidase 27 and protein kinase 2B (involved in signalling). The involvement of these genes in drought responses is discussed below.
Among the drought-regulated 351 genes in at least 6 of 9 studies, our attention was paid on the hormone, transcription factor and stress defence categories. Among those genes related to osmotic stress, the up-regulation of sucrose nonfermenting1–related protein kinase2 (SNF1-related protein kinase; also named as SnRK2) was conserved across species. This member belongs to a family of genes that have been previously associated with osmotic stresses [28, 29].  showed that these members are induced by osmotic stress and that three of them are activated through an ABA-dependent manner. Its role is extremely important in guard cells where it is playing a key role as a central hub to mediate ABA signalling . Taken together, our meta-analysis confirmed that this gene should be an important player in sensing water deprivation in leaf tissues in different crops. This gene should be considered as a target for crop genetic engineering for the development of molecular markers associated with drought-resistance in crops.
Our work identified two genes involved in abiotic stress signalling: a calcineurin B-like (CBL) calcium sensor protein and a CBL interacting protein kinase 1. Calcineurin B-like proteins (CBLs) represent a unique family of plant calcium sensors that relay signals by interacting with a family of protein kinases, designated as CBL-interacting protein kinases (CIPKs). A previous study indicated that CIPK23 play this important role in water stress response by interaction with the calcium sensors CBL1 and CBL9 that synergistically regulates CIPK23. As suggested by , the different combination of CIPK and CBL members should be responsible for cell-specific signalling responses (osmotic stress or potassium uptake) in different organs (leaves or roots). Based on these findings, we can speculate that the simultaneous induction of calcineurin B-like calcium sensor protein (AT4G17615) and CBL interacting protein kinase 1 (AT3G17510) across 6 of 9 studies implies that the two proteins should play a major role in the activation of rapid drought sensing. This hypothesis implies that these two genes might be also considered as good targets for molecular breeding to enhance abiotic stress resistance.
Relating to GO term “response to water deprivation”, three additional genes were up-regulated such as homeobox 7, lipid transfer protein 3 (LTP3), short dehydrogenase reductase 1 (ABA2). Homeobox 7 belongs to Homeodomain-leucine zipper (HD-Zip) family proteins which are transcription factors related to environmental stress responses in plants. A member of homeobox family has been shown to confer resistance to drought in Helianthus annuus (sunflower) through over-expression . LTP3 is known to bind to lipids and its over-expression enhanced drought tolerance through the action of MYB96 that directly binds to its promoter . ABA2 is a NAD- or NADP-dependent oxidoreductases involved in ABA biosynthesis. This gene is responsive to ABA exogenous treatment . Our analysis showed that a SOS3 like calcium binding protein was commonly induced by drought in 6 of 9 crops. This gene encodes a member of the calcineurin B-like calcium sensor gene family and mediates salt tolerance by regulating ion homeostasis in Arabidopsis. We also observed an up-regulation of salt over sensitive 1 (SOS1) that is a key player of the Salt-Overly-Sensitive (SOS) pathway, essential for maintaining a normal ion ratio in the cytoplasm in salt conditions . Salinity is biphasic stress composed by an initial change of osmotic conditions followed by a subsequent stage of ionic modifications. Indeed, SOS1 plays an important role in the second phase of salinity stress. Transgenic over-expression of this gene has shown to induce drought tolerance in Arabidopsis thaliana demonstrating that improved resistance to salt stress can be obtained by limiting Na + accumulation in plants . Being drought mainly osmotic stress, the induction of this gene implies a possible role of this gene in the response to osmotic changes too. This might be explained by the fact that water deprivation has the consequence to increase the levels of soil ion concentrations which indirectly causing salt stress. Our meta-analysis highlighted the repression of fatty acid biosynthesis in response to drought in leaves. It is known that water deficit inhibits fatty acid desaturation and drought resistance has been linked with a reduction of fatty acid metabolism in cotton resulting in greater stability of the membrane system .
Relating to hormones, we found several drought-regulated genes in common between 6 of 9 studies. Three genes were involved in ABA biosynthesis and signalling, two genes in auxin response, two genes involved in ethylene-related pathways. ABA2 Arabidopsis thaliana mutants showed a reduced drought tolerance in comparison to wild type  implying that the up-regulation of this gene should be a benefit for drought resistance. Our meta-analysis showed another unexpected result: ABA3 was repressed in response to water stress in 6 of 9 crops. This gene is a basic leucine zipper (bZIP)-type ABRE-binding protein that was shown to be up-regulated by drought in vegetative tissues . Although at first glance, our results on both AREB2 and HVA22 seem to be in contrast with published findings, the repression of this gene in response to drought might be due to differences in the analyzed time points and drought intensity between studies.
Relating to ethylene biosynthesis, we found that ACS12 was constitutively repressed by drought. It is generally accepted that ethylene is involved in mediating plant responses to abiotic stress . ACS cereal mutants showed to have delayed leaf senescence in drought conditions. Mutant leaves continue to be photosynthetically active under water stress implying that leaf function is maintained . These findings showed that ethylene may serve to determine the onset of natural senescence and regulate drought-induced senescence. Based on these findings, we may speculate that the repression of ACS12 in leaves should be beneficial to inhibit ethylene biosynthesis and consequently improve drought resistance. ERF1 is known to be involved in plant disease resistance  but its role in abiotic stress responses is less clear. ERF proteins are characterized by an ERF DNA binding domain. These transcription factors bind to multiple cis-elements such as DRE/CRT and CE1 elements, involved in stress responses . In Arabidopsis thaliana, the expression of ERF1 enhanced tolerance to drought.  hypothesized that ERF1 was linked with enhanced drought resistance in rice through the induction of ABA2. Since we found that both ABA2 and ERF1 were induced in 6 of 9 crops, our findings confirm this hypothesis rendering these two genes potential targets for enhancing resistance to drought. Among the conserved drought up-regulated transcription factors, it is worth to mention WRKY20, a member of WRKYs. This finding is in agreement with published works that demonstrated an increased drought tolerance due to the over-expression of WRKY20 in Arabidopsis thaliana .
Drought-regulated genes at different leaf developmental stages
We answered the second question by identifying the drought-regulated genes commonly expressed among the three transcriptomic studies dealing with seedling responses and among the five studies performed on mature leaves.
Our findings highlighted that drought has very different transcriptomic effects on leaves depending on their developmental stage. Indeed, the identification of expression QTLs for drought resistance in leaves should clearly take in high consideration which developmental stage is considered. Among the 22 drought-regulated genes commonly expressed seedling and mature leaves it is worth noticing a transcription factor (alfin-like 1) and a heat shock protein (HSC70–7). Several members of Alfin-like TFs were up-regulated in response to different abiotic stresses in Brassica oleracea. The role Alfin-like TFs in enhancing salt stress and drought resistance is well-known when it is over-expressed in roots [46, 47]. Alfin-like 1 is a transcription co-activator  that contains a typical PHD finger binding promoter element of PRP2, a salt inducible gene . Our meta-analysis lets hypothesize that the role of this transcription factor might have a similar function in leaves.
Drought-regulated genes in seedlings
The up-regulation of two ethylene signalling genes in seedlings (ERF4 and one EREBP (ESE3)) implies that ethylene might have a promoting effect in drought response at early leaf development. These findings agree with previously published works that showed that the over-expression of ERF4 promoted adaptation to salt stress and drought . This gene is a transcriptional repressor that suppressing a repressor of defence response genes positively regulates shoot growth and water-stress tolerance in rice during early growth stages . ESE3 belongs to a family of ethylene response factor (ERF) genes that are involved in enhancing salt tolerance. Results of our work confirmed that the up-regulation of ethylene signalling should play a key role in drought resistance. HVA22 is an ABA-responsive gene regulated by environmental stresses. The up-regulation of HVA22 has been shown to be tissue-specific and in response to drought in barley . Our analysis confirmed this evidence showing an opposite trend of expression between seedlings and mature leaves. Our results challenged the hypothesis that this gene should enhance drought resistance in mature leaves.
MYB is a large family of transcription factors well-known to be involved in drought. The transgenic over-expression of MYB1 enhanced drought resistance . MYB94 activates cuticular wax biosynthesis in Arabidopsis thaliana and might be important in drought response . Our analysis confirmed the role of these two MYB factors in seedling response to drought implying that they should be considered potential targets for enhancing drought resistance. The induction of MYB factors in drought is reported in the selected articles [16, 18, 19, 22]. We also found that WRKY53 was up-regulated in seedlings in response to drought confirming previous findings that showed WRKY53 drives the inhibition of stomatal closure by reducing H2O2 content facilitating stomatal opening by promoting starch degradation  and consequently inhibiting drought tolerance. The induction of WRKY20 in response to drought should allow a positive effects on drought tolerance in crops since the over-expression of this gene improved plant yields in soybean and enhanced drought tolerance in alfalfa . Our protein-protein interaction analysis showed that WDR5A was up-regulated in all three seedling studies in response to drought. This confirmed the important role of this protein in drought responses. WDR5A is a regulating nitric oxide accumulation and NOS-like activity in guard cells to modulate stomatal closure for adaptive plant response to drought . In seedlings, this gene should drive the closure of the stomata and the survival of leaf cells under water deprivation.
Drought-regulated genes in mature leaves
In mature leaves, our meta-analysis found that ERF1 was up-regulated. The same result can be seen in . The role of this gene in drought response has been previously discussed. Considering the mature leaf datasets all together, the role of ERF1 in modulating the expression of anti-oxidant and detoxifying proteins that protect cell components in leaf mature tissues is highlighted. ERF1 should work as a ‘regulatory gene’ under different stress conditions, changing the expression of ‘functional genes’ acting as detoxification and osmotic adjustment enzymes or proteins to protect cells from damage.
BAG6 is a Calmodulin (CaM)-binding transcription activators (CAMTA), which translates calcium signatures into different biochemical, and molecular pathways  and acts as a multi-functional protein that regulates apoptotic-like processes involved in different abiotic stresses. We found that this gene was up-regulated in mature leaves across the different crops. Indeed, we speculate that this gene should be involved in signalling mechanisms in response to drought stress. Calcium (Ca2+) works as a secondary messenger in plants and it is involved in different responses to different environmental stresses . These transcription factors modulate many functional genes involved in stress tolerance in plants including drought  and regulate the expression of ERFs . Based on these findings we may hypothesize that the two genes BAG6 and ERF1 might be linked in a common signalling response to drought in crops in mature leaves. Interestingly our analysis showed that HSP70 was repressed by drought. The heat shock protein 70s (Hsp70s) and heat shock factors (Hsfs) play key roles in protecting plant cells or tissues from various abiotic stresses . It was observed that heat shock proteins play as activators or repressors, suggesting that these proteins might be modulated by both the activation and the repression mechanisms under stress condition . Indeed, the effect of the repression of HSP70 in mature leaves under drought has to be further investigated. Cyclin dependent kinases (CDKs) are signalling proteins induced by stresses such as drought . Since this gene was a highly interacting protein in drought-related gene networks, we may speculate that the induction of CDKF1 in mature leaves should play an important role in the promotion of drought resistance in crops.
Taken together all these findings, we propose a model of plant response to drought shown in Fig. 7. The first plant response should be the induction of the biosynthesis of key hormones such as ABA and ethylene driving the activation of key signalling proteins (ERF1, ABA2 and HB7). These proteins should promote the fine-tuned transcriptional modulation through the cross-talk of a complex network of transcription factors (Alfin-like 1, WRKY20, SDG26). The up-regulation of key proteins in the signal transduction (CAMTA2, KIN2, and SNF7) should provoke the induction of proteins involved in physiological defensive responses represented by stomatal closure, inhibition of fatty acid biosynthesis, an increase of osmotic potential and protection of protein folding. Molecular breeding for drought resistance should focus on these genes.
This work allowed identifying genes involved in drought responses that were conserved across crop species. The mapping of genes in the different chromosomes of the 7 analyzed crops will help in developing new molecular markers for drought resistance. This gene mapping allowed identifying which chromosomes contain more DNA regions have to be taken in higher consideration for the breeding of crops in relation to drought resistance. Our approach will be the example of similar studies dealing with the identification of key genes involved in resistance to other limiting biotic and abiotic stresses in crops.
Search strategy for selection of RNA-Seq studies
Articles published dealing in response to drought in both tree fruit crop and herbaceous species were collected. These studies were identified from Scopus and PubMed if they respect the following criteria as follows: (i) consist of RNA-seq analysis, (ii) included at least one of the following terms in title and abstract: drought, leaf, stress, abiotic stress, water stress, (iii) studies provided raw data submitted in public databases. These criteria resulted in a selection of 9 articles comprising of total 42 samples (Table 3). If different time point series were present in the same study, a single time point was selected (10 DAS or whichever is nearer). This selection was done since most of the analyzed studies were performed at 10DAS. The selected studies were grouped based on leaf developmental stage (seedling leaves, young and mature leaves). Raw data were downloaded and analyzed through a pipeline previously used in similar meta-analysis . The complete pipeline used for this study was provided in Fig. 8. The drought-modulated genes identified from each of the 9 studies were given in Additional file 11: Table S6.
Read alignment, gene differential expression and annotation
For all the 9 articles, the latest available version of the corresponding crop genome and its annotation file were downloaded from Phytozome (https://phytozome.jgi.doe.gov). The raw data files were downloaded from NCBI SRA (https://www.ncbi.nlm.nih.gov/sra) and EMBL ArrayExpress (https://www.ebi.ac.uk/arrayexpress/) according to the accession number given in the article and converted to FASTQ format using SRA toolkit version 2.3.5. Raw data underwent pre-processing by trimming low-quality bases followed by adaptor sequence removal to obtain high-quality clean reads using cutadapt version 1.8.1. The pre-processed high-quality reads were mapped to the corresponding genome with HISAT2 version 2.1.0  using the default parameters. The resulted output of HISAT2 was then used for the identification of differentially expressed genes using Cuffdiff tool in Cufflinks version 2.2.1 pipeline with default parameters. Up- and down-regulated genes with p-value < 0.05 were considered for downstream functional analysis. The DEGs selected were annotated using corresponding crop genome mapping file downloaded from the Phytozome. Custom made in-house Perl script was used for the selection of genes and mapping.
Statistical and cluster analysis
The DEGs corresponding to each study separately analyzed, with p-value < 0.05, were then taken for the statistical analysis. Using p.adjust function of R, all the statistical tests were corrected for multiple comparisons using the Benjamini-Hochberg false discovery rate . This approach can make the FDR at the desired level of α (in this study 0.05) by adjusting the P-values. R software was used for the statistical analysis. Differences among the selected studies were adjusted using the sample normalization. In order to remove systematic variation between different species, the normalization procedure served as a crucial pre-processing step to adjust for the different sample sequencing depths and other confounding technical effects. We used the geometric normalization method where FPKMs and fragment counts are scaled via the median of the geometric means of fragment counts across all libraries, as described in . The dendrogram was generated for identifying the clustering patterns of the considered studies. The grouping of the clusters for dendrogram was done using the Euclidean distance measure.
Gene set enrichment analysis
We mapped the entire differentially regulated gene IDs of each plant species to Arabidopsis thaliana and found out the corresponding best hit TAIR ID using the annotation file downloaded from Phytozome. We used MapMan  with the Arabidopsis thaliana mapping file (http://mapman.gabipd.org/) to map and visualize the metabolic overview, hormone regulation, secondary metabolism, transcription factors, and protein targeting. Firstly, we visualized the drought-regulated genes in common in at least 6 of 9 studies. Secondly, we visualized the drought-regulated genes in common between the three studies in seedlings; and finally in common between the five in mature leaves. The PageMan  analysis, plugin of MapMan, was used to visualize differences among metabolic pathways using Wilcoxon tests, no correction, and an over-representation analysis (ORA) cutoff value of 3. All the homologous TAIR IDs of the 9 studies were searched against the Database for Annotation, Visualization and Integrated Discovery (DAVID) version 6.8  Web server (https://david.ncifcrf.gov/). The gene ontology information related to the biological process was extracted from the DAVID result.
Gene mapping in crop chromosomes
The drought-regulated genes involved in abiotic stress, hormone metabolism and transcription factors were selected for the chromosome mapping. We found out the corresponding chromosome number, start and end of the drought-regulated gene IDs from the annotation file downloaded from Phytozome (https://phytozome.jgi.doe.gov) using custom made Perl script. These genes were then mapped towards the chromosome according to their chromosome number, start and end points.
Protein-protein interaction network
NetworkAnalyst , a web-based tool for network-based visual analytics of protein-protein interaction networks, was used. The list of homologous TAIR IDs from 6 of 9 studies were uploaded and mapped against the STRING interactome database with default parameters (confident score cutoff = 900 and with experimental evidence) provided in NetworkAnalyst. The networks between drought-regulated genes in seedlings and in mature leaves corresponding to the list of the visualized genes in MapMan were also obtained. To study the key connectives and to simplify the large network, we selected “Minimum Network”.
Database for Annotation, Visualization and Integrated Discovery
Differentially Expressed Genes
Boyer JS. The U.S. drought of 2012 in perspective: a call to action. Global Food Security. 2013;2:139–43.
Yu C. China’s water crisis needs more than words. Nature. 2011;470:307–7.
Lynch JP. Steep, cheap and deep: an ideotype to optimize water and N acquisition by maize root systems. Ann Bot (Lond). 2013;112:347–57.
Janiak A, Kwasniewski M, Sowa M, Gajek K, Żmuda K, Kościelniak J, et al. No time to waste: Transcriptome study reveals that drought tolerance in barley may be attributed to stressed-like expression patterns that exist before the occurrence of stress. Front Plant Sci. 2018;8.
Peleg Z, Blumwald E. Hormone balance and abiotic stress tolerance in crop plants. Plant Biol. 2011;14:290–5.
Krasensky J, Jonak C. Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J Exp Bot. 2012;63(4):1593–608.
Qin F, Shinozaki K, Yamaguchi-Shinozaki K. Achievements and challenges in understanding plant abiotic stress responses and tolerance. Plant Cell Physiol. 2011;52:1569–82.
Wang HY, Wang HL, Shao HB, Tang XL. Recent advances in utilizing transcription factors to improve plant abiotic stress tolerance by transgenic technology. Front Plant Sci. 2016;7:67.
Martinelli F, Reagan RL, Uratsu SL, Phu ML, Albrecht U, Zhao W, et al. Gene regulatory networks elucidating huanglongbing disease mechanisms. PLoS One. 2013;8:e74256. https://doi.org/10.1371/journal.pone.0074256.
Martinelli F, Ibanez A, Reagan R, Davino S, Dandekar A. Stress responses in citrus peel: comparative analysis of host responses to huanglongbing disease and puffing disorder. Sci Hortic. 2015;192:409–20. https://doi.org/10.1016/j.scienta.2015.06.037.
Zhao S, Li C, Guo Y, Sheng Q, Shyr Y. RnaSeqSampleSize: Real data based sample size estimation for RNA sequencing. BMC Bioinformatics. 2018;19:1. https://doi.org/10.1186/s12859-018-2191-5.
Martinelli F, Scalenghe R, Giovino A, Marino P, Aksenov AA, Pasamontes A, Dandekar A. Proposal of a Citrus translational genomic approach for early and infield detection of Flavescencedorée in Vitis. Plant Biosystems - An International Journal Dealing with All Aspects of Plant Biology. 2014;150(1):43–53. https://doi.org/10.1080/11263504.2014.908976.
Tosetti R, Martinelli F, Tonutti P, Barupal D. Metabolomics Approach To Studying Minimally Processed Peach (PrunusPersica) Fruit. ActaHorticulturae. 2012;934:1017–21. https://doi.org/10.17660/actahortic.2012.934.135.
Rizzini F, Bonghi C, Chkaiban L, Martinelli F, Tonutti P. Effects of postharvest partial dehydration and prolonged treatments with ethylene on transcript profiling in skins of wine grape berries. ActaHorticulturae. 2010;877:1099–104.
Zhai P, Yang L, Guo X, Wang Z, Guo J, Wang X, Zhu H. MetaComp: Comprehensive analysis software for comparative meta-omics including comparative metagenomics. BMC Bioinformatics. 2017:18(1). https://doi.org/10.1186/s12859-017-1849-8.
Clauw P, Coppens F, Beuf KD, Dhondt S, Daele TV, Maleux K, et al. Leaf responses to mild drought stress in natural variants of Arabidopsis. Plant Physiol. 2015;167:800–16.
Song K, Kim HC, Shin S, Kim K-H, Moon J-C, Kim JY, et al. Transcriptome analysis of flowering time genes under drought stress in maize leaves. Front Plant Sci. 2017;8.
Corso M, Vannozzi A, Maza E, Vitulo N, Meggio F, Pitacco A, et al. Comprehensive transcript profiling of two grapevine rootstock genotypes contrasting in drought susceptibility links the phenylpropanoid pathway to enhanced tolerance. J Exp Bot. 2015;66:5739–52.
Li P, Cao W, Fang H, Xu S, Yin S, Zhang Y, et al. Transcriptomic profiling of the maize (Zea mays L.) leaf response to abiotic stresses at the seedling stage. Frontiers in. Plant Sci. 2017;8.
Pieczynski M, Wyrzykowska A, Milanowska K, Boguszewska-Mankowska D, Zagdanska B, Karlowski W, et al. Genomewide identification of genes involved in the potato response to drought indicates functional evolutionary conservation with Arabidopsis plants. Plant Biotechnol J. 2017;16:603–14.
Orcheski B, Brown S. High-throughput sequencing reveals that pale green lethal disorder in apple (Malus) stimulates stress responses and affects senescence. Tree Genet Genomes. 2016;13.
Liu M, Yu H, Zhao G, Huang Q, Lu Y, Ouyang B. Profiling of drought-responsive microRNA and mRNA in tomato using high-throughput sequencing. BMC Genomics. 2017;18.
Salman MH, Zhang C, Pervaiz T, Zheng T, Zhang CB, Lide C, et al. Gene regulation mechanism in drought-responsive grapevine leaves as revealed by transcriptomic analysis. bioRxiv. 2016:65136.
Liu Z, Xin M, Qin J, Peng H, Ni Z, Yao Y, et al. Temporal transcriptome profiling reveals expression partitioning of homeologous genes contributing to heat and drought acclimation in wheat (Triticum aestivum L.). BMC Plant Biol. 2015;15.
Rest JS, Wilkins O, Yuan W, Purugganan MD, Gurevitch J. Meta-analysis and meta-regression of transcriptomic responses to water stress in Arabidopsis. Plant J. 2016;85(4):548–60. https://doi.org/10.1111/tpj.13124.
Muthuramalingam P, Krishnan SR, Pothiraj R, Ramesh M. Global Transcriptome Analysis of Combined Abiotic Stress Signaling Genes Unravels Key Players in Oryza sativa L.: An In silico Approach. Front Plant Sci. 2017:8. https://doi.org/10.3389/fpls.2017.00759.
Ramu VS, Paramanantham A, Ramegowda V, Mohan-Raju B, Udayakumar M, Senthil-Kumar M. Transcriptome Analysis of Sunflower Genotypes with Contrasting Oxidative Stress Tolerance Reveals Individual- and Combined- Biotic and Abiotic Stress Tolerance Mechanisms. PLoS One. 2016:11(6). https://doi.org/10.1371/journal.pone.0157522.
Boudsocq M. Osmotic signaling in plants. Multiple pathways mediated by emerging kinase families. Plant Physiol. 2005;138:1185–94.
Boudsocq M, Droillard MJ, Barbier-Brygoo H, Laurière C. Different phosphorylation mechanisms are involved in the activation of sucrose non-fermenting 1 related protein kinases 2 by osmotic stresses and abscisic acid. Plant Mol Biol. 2006;63:491–503.
Kobayashi Y, Yamamoto S, Minami H, Kagaya Y, Hattori T. Differential activation of the rice sucrose nonfermenting1-related protein kinase2 family by hyperosmotic stress and abscisic acid. Plant Cell. 2004;16(5):1163–77.
Yoshida R, Umezawa T, Mizoguchi T, Takahashi S, Takahashi F, Shinozaki K. The regulatory domain of SRK2E/OST1/SnRK2.6 interacts with ABI1 and integrates Abscisic acid (ABA) and osmotic stress signals controlling Stomatal closure inArabidopsis. J Biol Chem. 2005;281:5310–8.
Cheong YH, Pandey GK, Grant JJ, Batistic O, Li L, Kim B-G, et al. Two calcineurin B-like calcium sensors, interacting with protein kinase CIPK23, regulate leaf transpiration and root potassium uptake in Arabidopsis. Plant J. 2007;52:223–39.
Dezar CA, Fedrigo GV, Chan RL. The promoter of the sunflower HD-zip protein gene Hahb4 directs tissue-specific expression and is inducible by water stress, high salt concentrations and ABA. Plant Sci. 2005;169:447–56.
Guo L, Yang H, Zhang X, Yang S. Lipid transfer protein 3 as a target of MYB96 mediates freezing and drought stress in Arabidopsis. J Exp Bot. 2013;64(6):1755–67.
Zhou J, Wang J, Li X, et al. H2O2 mediates the crosstalk of brassinosteroid and abscisic acid in tomato responses to heat and oxidative stresses. J Exp Bot. 2014;65(15):4371–83.
Huang Z, Long X, Wang L, Kang J, Zhang Z, Zed R, et al. Growth, photosynthesis and H -ATPase activity in two Jerusalem artichoke varieties under NaCl-induced stress. Process Biochem. 2012;47:591–6.
Shi H, Lee BH, Wu SJ, Zhu JK. Overexpression of a plasma membrane Na+/H+ antiporter gene improves salt tolerance in Arabidopsis thaliana. Nat Biotechnol. 2003;21:81–5.
Pham-Thi AT, Borrel-Flood C, Vieira da Silva J, Justin AM, Mazliak P. Effects of water stress on lipid metabolism in cotton leaves. Phytochemistry. 1985;24:723–7.
Ozfidan C, Turkan I, Sekmen AH, Seckin B. Time course analysis of ABA and non-ionic osmotic stress-induced changes in water status, chlorophyll fluorescence and osmotic adjustment in Arabidopsis thaliana wild-type (Columbia) and ABA-deficient mutant (aba2). Environ Exp Bot. 2013;86:44–51.
Uno Y, Furihata T, Abe H, Yoshida R, Shinozaki K, Yamaguchi-Shinozaki K. Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc Natl Acad Sci. 2000;97:11632–7.
Young TE, Meeley RB, Gallie DR. ACC synthase expression regulates leaf performance and drought tolerance in maize. Plant J. 2004;40:813–25.
Singh K. Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002;5:430–6.
Zhang M, Liu W, Bi Y-P. Dehydration-responsive element-binding (DREB) transcription factor in plants and its role during abiotic stresses. Hereditas (Beijing). 2009;31:236–44.
Zhang G, Guo G, Hu X, et al. Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. Genome Res. 2010;20(5):646–54.
Luo X, Bai X, Sun X, Zhu D, Liu B, Ji W, et al. Expression of wild soybean WRKY20 in Arabidopsis enhances drought tolerance and regulates ABA signalling. J Exp Bot. 2013;64:2155–69.
Winicov I. Alfin1 transcription factor over expression enhances plant root growth under normal and saline conditions and improves salt tolerance in alfalfa. Planta. 2000;210:416–22.
Ramu VS, Swetha TN, Sheela SH, Babitha CK, Rohini S, Reddy MK, et al. Simultaneous expression of regulatory genes associated with specific drought-adaptive traits improves drought adaptation in peanut. Plant Biotechnol J. 2015;14:1008–20. https://doi.org/10.1111/pbi.12461.
Lee WY, Lee D, Chung WI, Kwon CS. Arabidopsis ING and Alfin1-like protein families localize to the nucleus and bind to H3K4me3/2 via plant homeodomain fingers. Plant J. 2009;58:511–24.
Bastola DR, Pethe VV, Winicov I. Alfin1, a novel zinc-finger protein in alfalfa roots that binds to promoter elements in the salt-inducible MsPRP2 gene. Plant Mol Biol. 1998;38:1123–35.
Seo Y, Park JB, Cho YJ, Jung C, Seo H, Park S-K, Nahm B, Song J. Overexpression of the ethylene-responsive factor gene BrERF4 from Brassica rapa increases tolerance to salt and drought in Arabidopsis plants. Mol Cell. 2010;30:271–7.
Joo J, Lee YH, Kim YK, Nahm BH, Song SI. Abiotic stress responsive rice ASR1 and ASR3 exhibit different tissue-dependent sugar and hormone-sensitivities. Mol Cell. 2013;35(5):421–35.
Shen Q, Chen CN, Brands A, Pan SM, Ho THD. The stress- and abscisic acid-induced barley gene HVA22: developmental regulation and homologues in diverse organisms. Plant Mol Biol. 2001;45:327–40.
Dai X, Xu Y, Ma Q, Xu W, Wang T, Xue Y, Chong K. Overexpression of a R1R2R3 MYB gene, OsMYB3R-2, increases tolerance to freezing, drought, and salt stress in transgenic Arabidopsis. Plant Physiol. 2007;143:1739–51.
Lee SB, Kim H, Kim RJ, Suh MC. Overexpression of Arabidopsis MYB96 confers drought resistance in Camelina sativa via cuticular wax accumulation. Plant Cell Rep. 2014;33:1535–46.
Sun Y, Yu D. Plant Cell Rep. 2015;34:1295 https://doi.org/10.1007/s00299-015-1787-8.
Ning W, Zhai H, Yu J, Liang S, Yang X, Xing X, et al. Overexpression of Glycine soja WRKY20 enhances drought tolerance and improves plant yields under drought stress in transgenic soybean. Mol Breed. 2017;37:19.
Liu W-C, Li Y-H, Yuan H-M, Zhang B-L, Zhai S, Lu Y-T. WD40-REPEAT 5a functions in drought stress tolerance by regulating nitric oxide accumulation in Arabidopsis. Plant Cell Environ. 2016;40:543–52.
Evans TI, Hell JW, Shea MA. Thermodynamic linkage between calmodulin domains binding calcium and contiguous sites in the C-terminal tail of Ca(V)1.2. Biophys Chem. 2011;159(1):172–87.
Reddy AS, Ali GS, Celesnik H, Day IS. Coping with stresses: roles of calcium- and calcium/Calmodulin-regulated gene expression. Plant Cell. 2011;23:2010–32.
Janiak A, Kwaśniewski M, Szarejko I. Gene expression regulation in roots under drought. J Exp Bot. 2015;67:1003–14.
Pandey N, Ranjan A, Pant P, Tripathi RK, Ateek F, Pandey HP, et al. CAMTA 1 regulates drought responses in Arabidopsis thaliana. BMC Genomics. 2013;14:216.
Wen F, Wu X, Li T, et al. Genome-wide survey of heat shock factors and heat shock protein 70s and their regulatory network under abiotic stresses in Brachypodium distachyon. PLoS One. 2017;12(7):e0180352. Published 2017 Jul 6. doi:https://doi.org/10.1371/journal.pone.0180352
Campo S, Baldrich P, Messeguer J, Lalanne E, Coca M, San SB. Overexpression of a calcium-dependent protein kinase confers salt and drought tolerance in Rice by preventing membrane lipid peroxidation. Plant Physiol. 2014;165(2):688–704.
Balan B, Caruso T, Martinelli F. Gaining insight into exclusive and common Transcriptomic features linked with biotic stress responses in Malus. Front Plant Sci. 2017;8. https://doi.org/10.3389/fpls.2017.01569.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12:357–60. https://doi.org/10.1038/nmeth.3317.
Benjamini Y, and Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. Series B (Methodological) 1995; 57(1): 289–300. Retrieved from http://www.jstor.org/stable/2346101.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11.
Thimm O, Bläsing O, Gibon Y, Nagel A, Meyer S, Krüger P, et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 2004;37:914–39. https://doi.org/10.1111/j.1365-313X.2004.02016.x.
Usadel B, Nagel A, Steinhauser D, Gibon Y, Bläsing OE, Redestig H, et al. PageMan: an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments. BMC Bioinformatics. 2006;7:535. https://doi.org/10.1186/1471-2105-7-535.
Huang DAW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. https://doi.org/10.1038/nprot.2008.211.
Xia J, Benner MJ, & Hancock RE. NetworkAnalyst - integrative approaches for protein–protein interaction network analysis and visual exploration. Nucleic Acids Research 2014; 42(W1). doi:https://doi.org/10.1093/nar/gku443.
We would like to thank Manuel Philip and Jayasankar M. J., Agrigenome Labs Private Ltd., Cochin, India for correcting the paper and generating the high quality images.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Table S1. The table consists of gene_ID, corresponding arabidopsis orthologs and log2 fold change values for significantly regulated genes (FDR < 0.05). The analysis results of each selected study has been given in seperate sheet along with the crop names. (XLSX 4134 kb)
Figure S1. Clustering heatmap showing the hierarchical relationship among the studies selected for the analysis. Resulted log2FC values of the analysis for generating the tree was indicated. (TIF 163 kb)
Figure S2. MapMan overview showing transcriptomic effects of drought in key categories selected such as secondary metabolism, cellular responses and signaling. Genes were identified as Arabidopsis orthologs of each genes of the analyzed plant species. Red means up-regulated and green means down-regulated. (TIF 226 kb)
Figure S3. Protein-protein interaction network analysis predicted for genes commonly regulated in three seedling leaf studies performed in Arabidopsis thaliana, Malus X domestica and Solanum lycopersicum based on Arabidopsis knowledgebase. (TIF 332 kb)
Figure S4. Protein-protein interaction network analysis predicted for genes commonly regulated in five mature leaf studies performed in Vitis vinifera, Solanum tuberosum, Triticum aestivum, Zea mays (study1) and Zea mays (study 2)based on Arabidopsis knowledgebase. (TIF 185 kb)
Figure S5. Key genes encoding transcription factors, hormone metabolism and abiotic stress responses obtained from DAVID software were mapped in the respective chromosomes of the 7 crops. (TIF 536 kb)
Table S2. Genes mapped in the chromosome which are involved in transcription factors, abiotic stress response and hormone metabolism. Arabidopsis ortholog, description and corresponding gene IDs were indicated. (XLSX 17 kb)
Figure S6. The 27 key genes that were drought-regulated in at least 7 of 9 studies were mapped in the respective chromosomes of the 7 crops. (TIF 536 kb)
Table S3. Genes mapped in the chromosome which are involved in atleast 7 of 9 studies. Arabidopsis ortholog, description and corresponding gene IDs were indicated. (XLSX 12 kb)
Table S4. Genes involved in the leaf developmental stage that are commonly drought-regulated were shown. Common in seedling studies and common in mature studies were given in separate sheets. GeneID, description, expression pattern and binname were indicated. (XLSX 14 kb)
Table S6. Genes involved in drought response corresponding to each of the study were shown. TAIR ID, description and log2FC were indicated. (XLSX 321 kb)
Table S5. Genes involved in drought response at different leaf developmental stages (unique for seedling, unique for mature and common for seedling and mature) were shown. TAIR ID, log2FC and expression pattern were indicated. (XLSX 56 kb)
About this article
Cite this article
Benny, J., Pisciotta, A., Caruso, T. et al. Identification of key genes and its chromosome regions linked to drought responses in leaves across different crops through meta-analysis of RNA-Seq data. BMC Plant Biol 19, 194 (2019). https://doi.org/10.1186/s12870-019-1794-y