Original Research ARTICLE
Regulatory Cross-Talks and Cascades in Rice Hormone Biosynthesis Pathways Contribute to Stress Signaling
- 1Department of Biophysics, Molecular Biology and Bioinformatics, University of Calcutta, Kolkata, India
- 2Computational Systems Biology Group, Center of Excellence in Systems Biology and Biomedical Engineering, University of Calcutta, Kolkata, India
Crosstalk among different hormone signaling pathways play an important role in modulating plant response to both biotic and abiotic stress. Hormone activity is controlled by its bio-availability, which is again influenced by its biosynthesis. Thus, independent hormone biosynthesis pathways must be regulated and co-ordinated to mount an integrated response. One of the possibilities is to use cis-regulatory elements to orchestrate expression of hormone biosynthesis genes. Analysis of CREs, associated with differentially expressed hormone biosynthesis related genes in rice leaf under Magnaporthe oryzae attack and drought stress enabled us to obtain insights about cross-talk among hormone biosynthesis pathways at the transcriptional level. We identified some master transcription regulators that co-ordinate different hormone biosynthesis pathways under stress. We found that Abscisic acid and Brassinosteroid regulate Cytokinin conjugation; conversely Brassinosteroid biosynthesis is affected by both Abscisic acid and Cytokinin. Jasmonic acid and Ethylene biosynthesis may be modulated by Abscisic acid through DREB transcription factors. Jasmonic acid or Salicylic acid biosynthesis pathways are co-regulated but they are unlikely to influence each others production directly. Thus, multiple hormones may modulate hormone biosynthesis pathways through a complex regulatory network, where biosynthesis of one hormone is affected by several other contributing hormones.
Plants being sessile creatures are forced to adapt to various adverse conditions to successfully complete their life-cycles. A repertoire of complex signaling systems was acquired by plants to respond to physiological and environmental cues. Phytohormones are one of the key mediators that afford plants the ability to rapidly respond to external cues by adjusting their metabolism. Previous research efforts have revealed importance of individual hormones in plant development and stress response (Creelman and Mullet, 1997; Yamaguchi, 2008; Zhao, 2010). Biological activity of any hormone depends on its availability which is controlled by level of its synthesis, transport and conjugation or degradation as well as efficiency of signal perception and transduction. Modulations in any of these have direct impact on downstream hormone responses (Figure 1A). Moreover, activities of various hormones have been shown to be overlapping; and cross-talk among different hormonal response pathways is a well-established phenomenon (Gazzarrini and McCourt, 2003; Robert-Seilaniantz et al., 2011; Vanstraelen and Benková, 2012). Hormonal interactions provide an additional level of complexity in regulation of plant response to internal and external cues as well as render a feedback mechanism to balance system robustness and dynamicity. Thus, the final response of a plant is not determined by activity of a single hormone but rather by a complex network where activity of a hormone is also modulated by other contributing hormones.
Figure 1. Overview of factors that may regulate hormone signaling and biosynthesis, which leads to plant stress and defense response. (A) Hormone signaling brought about by biological active hormone depends on their availability which is controlled by level of their synthesis, transport and conjugation or degradation as well as efficiency of signal perception and transduction. Modulations in any of these have direct impact on downstream stress and defense responses. (B) Transcription factors bind to cis-regulatory elements (CRE) in the promoters of hormone biosynthesis genes and regulate hormone biosynthesis and consequently hormone signaling and response. Conversely promoters of transcription factor genes may harbor hormone responsive CREs that control transcription of these genes in response to hormone signaling.
In rice both disease resistance and abiotic stress responses appear to be controlled by a complicated hormone signaling network. Salicylic acid (SA) and jasmonic acid (JA) promote resistance against pathogens (Figure 2). Ethylene (ET) may have negative or positive impact on disease resistance depending upon pathogen lifestyle (De Vleesschauwer et al., 2013). Abscisic acid (ABA), JA and ET interact positively or negatively depending upon stress conditions to promote stress responses. ABA is known to be antagonistic to SA and suppresses disease resistance but promotes abiotic stress tolerance (Kohli et al., 2013). The interactions of growth hormones are equally complex e.g., gibberellic acid (GA) is known to suppress stress responses through its antagonistic interaction with both ABA and JA (De Vleesschauwer et al., 2013). Auxin is also known to suppress defense against pathogen but augments abiotic stress tolerance (Kohli et al., 2013). Brassinosteroid (BR) through its interaction with stress hormones as well as growth hormones plays important role in biotic and abiotic stress response. Role of cytokinin (CK) in stress management is complicated and condition specific. CK interacts negatively with ABA, JA, and BR and promotes stress tolerance by delaying stress response (Figure 2). It is depleted from tissue exposed to prolonged stress, but accumulates during stress of short duration. Despite this cross-talk the metabolic pathways controlling individual hormone production are specific and apparently non-redundant. This poses the questions about regulation of independent hormone biosynthesis pathways and the possible ways in which multiple hormone production pathways may be modulated resulting in an integrated response.
Figure 2. Model showing probable hormone defense (Magnaporthe oryzae infection) and stress (drought) networking in Oryza sativa. Salicylic acid (SA), jasmonic acid (JA) ethylene (ET), brassinosteroid (BR) and abscisic acid (ABA) interact positively or negatively depending upon stress conditions. ABA is antagonistic to SA and suppresses disease resistance but promotes abiotic stress tolerance. Gibberellic acid (GA) suppress stress responses. Auxin suppress defense against pathogen but augments abiotic stress tolerance. Cytokinin (CK) interacts negatively with ABA, JA and BR. Positive and negative regulatory actions are indicated by arrows and lines with bars, respectively.
Regulation of hormone biosynthetic pathways at both transcriptional and post transcriptional level has been well documented (Xiong and Zhu, 2003; Frigerio et al., 2006; Argueso et al., 2007), however modulations of these pathways in an integrated way remains unclear. One of the possibilities may be to use common transcription regulators for several hormones thereby controlling the transcription of genes related to more than one hormone biosynthesis at the same time. Transcriptional control of genes is brought about by cis-regulatory elements (CRE), located in promoters of target genes, and their corresponding transcription factors. As the promoters may harbor multiple CREs, gene expression depends upon binding of corresponding transcription factors and the specific combinatorial logic among them (Figure 1B). CRE mediated regulation has the potential to produce gene expression patterns which may be tissue and/or condition specific. Several computational approaches intended to understand the CRE mediated combinatorial regulation at the promoter regions of genes (Pilpel et al., 2001; Sudarsanam et al., 2002; Kato et al., 2004; Yu et al., 2006; Murakami et al., 2008; Deb and Kundu, 2015) are available. In present work, we applied the existing computational approaches to study possible CRE mediated regulation of hormone biosynthetic pathways in rice leaf and the dynamics of these regulation under Magnaporthe oryzae attack (causal organism of blast disease) and drought stress. Results of this work allowed us to map the CRE-mediated regulation of differentially up-regulated rice hormone biosynthesis (HB) genes under blast and drought conditions. Through analysis of CREs associated with differentially expressed HB genes we were able to develop insights about effect of different hormones and stress conditions on regulation of all HB pathways. We were also able to identify key transcription factors for CRE elements that integrate and modulate these pathways.
2. Materials and Methods
2.1. Promoter Sequence Data
Promoter sequence data (1 kb upstream) of rice loci is collected from Rice Genome Annotation Project, version 6.1 (Ouyang et al., 2007). RepeatMasker, version open-4.0.2 (Smit etal., 2013–2015) is used to mask interspersed and simple repeats in the promoter sequences. This filtering restricts unwanted discovery of CREs into the stretches of repeated DNAs present in the promoters.
2.2. Cis-Regulatory Element Data
The information of CREs is collected from PLACE (Plant cis-acting regulatory element) database (Higo et al., 1999). This initial dataset is screened to fetch only the unique entries. Signal Scan (Prestridge, 1991) tool is used to fetch the occurrences and positions of the CREs in the promoters.
2.3. Mapping Plant Hormone Biosynthesis (HB) Pathway Genes
The information of rice metabolic pathways and their associated loci is collected from the RiceCyc database (v.3.3) (Jaiswal et al., 2006). The initial data is further screened to get the hormone biosynthesis pathways. All together information for 8 different plant hormone biosynthesis pathways (Auxin, ABA, BR ,CK, ET, GA, JA, and SA) was obtained.
2.4. Gene Expression Data
Gene expression data were collected from Gene Expression Omnibus (GEO) (Edgar et al., 2002) for Rice under two different conditions, blast infection (GSE7256) and drought (GSE26280). RMA (Robust Multiarray Average) normalization was performed using R and Bioconductor (Gentleman et al., 2004) and differentially up-regulated genes (log fold change >2; p < 0.05) from each experiment (4 dpi of GSE7256 and panicle stage leaf of GSE26280) were considered for further analysis.
2.5. Estimation of Over-Represented CREs in Individual Rice Promoters
The over-represented CREs in the promoters can be obtained by different statistical tests. We used cumulative hypergeometric statistics (Hughes et al., 2000; Sudarsanam et al., 2002) (Equation 1) to determine the over-represented CREs in individual promoters of the respective genes.
Where, N is the total number of occurrences of all CREs in rice genome, m is the total number of occurrences of a considered CRE in the genome, x is the number of occurrences of that CRE in an individual promoter of a gene, n is the total number of occurrences of all CREs in that promoter and t is the summation index. Thus, PCRE is the probability of occurrence of a considered CRE in a particular promoter. A particular CRE was considered as over-represented in a specific promoter if the respective PCRE is < 0.05. This statistics was executed for all CREs in each HB locus promoter.
2.6. Estimation of CRE Co-occurrence
In our previous work (Deb and Kundu, 2015) a unique methodology has been introduced to identify the significantly co-occurring CRE pairs in the rice promoters and to develop different level of networks using them; the same methodology is applied in the present study. We computed the co-occurrence tendency in terms of COR value for all possible CRE pairs present in the promoters of differentially up-regulated HB genes for each condition (blast and drought). A series of filtering steps (Deb and Kundu, 2015) were used to extract out only the significantly co-occurring CRE pairs.
2.7. Construction of Networks: Locus-Locus Co-enriched CRE Network and CRE Co-occurrence Network
We constructed locus-locus co-enriched CRE networks and CRE co-occurrence networks for each condition (blast and drought) separately. The locus-locus co-enriched CRE network was developed using the up-regulated HB genes and the over-represented CREs at their promoters. Each node of the network was specified by individual locus and two loci were interconnected if they have at least one common over-represented CRE.
Individual CRE co-occurrence network was constructed using the significantly co-occurring CRE pairs where nodes represent individual CRE and two CREs were connected if they significantly co-occur with each other (Deb and Kundu, 2015). We intended to find the co-occurring CRE combinations in the form of cliques. Clique(k) is interpreted as a complete (fully connected) subgraph having k number of nodes. In CRE co-occurrence network, a clique represents a combination of co-occurring CREs having statistically significant co-occurrence relationship with every other CREs of the same clique. CFinder (Adamcsek et al., 2006) was used to identify the cliques in each network.
3.1. CREs in Up-Regulated Hormone Biosynthesis Genes under Blast and Drought
We found 37 HB loci in blast infected leaf and 26 HB loci in leaf under drought stress to be up-regulated. In case of blast dataset (GSE7256), genes associated with Gibberellin biosynthesis were absent. Similarly Salicylic acid biosynthesis associated genes were absent from drought dataset (GSE26280) (Table S1). This indicates that SA pathways are unlikely to have any role in case of drought stress whereas JA pathway has important implications.
We estimated over-represented CREs in individual promoters of these two sets of loci. On an average 8 CREs (8.2 in blast and 7.8 in drought) were over-represented in individual promoters. Promoter of ABA biosynthesis locus Os07g05940.1 was associated with highest number of CREs, sixteen in number. This locus was common to both blast and drought datasets. Lowest number of CREs was associated with JA biosynthesis locus Os06g11210.1, it had a single CRE element (Table S1).
3.2. CRE Co-occurrence Networks of Hormone Biosynthesis Genes under Blast and Drought
The co-occurrence tendency for all possible CRE pairs present in the promoters of differentially up-regulated HB genes for each condition (37 for blast and 26 for drought) were calculated and significantly co-occurring pairs based on the filtering steps (Deb and Kundu, 2015) were extracted and used to generate two separate CRE co-occurrence networks. Network for blast consisted of 137 nodes (CREs) and 296 edges (co-occurring relationship) whereas, that for drought consisted of 122 nodes and 232 edges. In case of blast condition, the co-occurrence network shows clustering coefficient 0.16, network density 0.03 and characteristics path length 3.11. On the other hand, the drought related network has clustering coefficient 0.06, network density 0.03 and characteristics path length 3.59.
When we merged these two networks, a relatively large number of nodes were found to be common (102), whereas, number of common interacting edges were relatively small in number (24) (Figure 3A). This indicates that though common CREs are present in the promoters of both blast and drought up-regulated HB genes, their co-occurring interaction patterns differ from condition to condition. This further indicates a possibility of presence of unique combinations of transcription factors that regulate the biosynthetic genes at two different stresses.
Figure 3. Merged CRE co-occurrence network and locus-locus co-enriched CRE networks of hormone biosynthesis genes up-regulated under Magnaporthe oryzae infection and Drought. (A) Merged CRE co-occurrence network of hormone biosynthesis genes up-regulated under Magnaporthe oryzae infection—Blast and Drought. Nodes (individual CREs) explicitly found in the Blast are highlighted (red) in the left, while those found in the Drought are highlighted (green) in the right hand side of the network. Nodes in the middle (cyan) are found in both Blast and Drought. Pink and bottle green color of edges imply that a significant co-occurrence is found between the two connected CREs explicitly in Blast and Drought respectively. A blue edge implies that the two connected CREs significantly co-occur in both the conditions. (B,C) Locus-locus co-enriched CRE networks for up-regulated hormone biosynthesis (HB) genes. (B) Magnaportheoryzae infection and (C) during drought stress. Each node of the network represents individual locus and two loci are interconnected if they have at least one common over-represented CREs. The edge weight represents the number of common over-represented CREs among them.
3.3. Locus-Locus Co-enriched CRE Network
The locus-locus co-enriched CRE network was developed using the up-regulated HB genes having over-represented CREs at their promoters. Each node of the network represents individual locus and two loci are interconnected if they have at least one common over-represented CREs. Thus, the interacting edges represent the regulatory similarities among the promoters (loci) where the edge weight shows the number of common over-represented CREs among them. Interestingly we found that loci of same HB pathways as well as of different HB pathways share abundant edges. This indicates that different HB pathways may have common regulators (Figures 3B,C).
3.4. Cliques of CRE Co-occurrence Networks
Clique of CREs is a combination of co-occurring CREs i.e., each CRE has significant co-occurrence relationship with other members of the same clique. Clique (k = 3–8) analysis was performed on CRE co-occurrence network for each condition. In case of blast and drought 76 and 29 cliques were observed, respectively. Twenty one cliques in case of blast and a single clique in case of drought were found to encompass all the HB pathways. Whereas 55 cliques in blast and 28 cliques in drought were found to be present only in loci of particular HB pathways. It was interesting to note that in case of Blast we did not find any clique which are exclusive to a single HB pathway; in case of drought we found two (d0:[CACTFTPPCA1, GTGANTG10, PROXBBNNAPA] and d14:[CPBCSPOR, CACTFTPPCA1, PROXBBNNAPA]) exclusive to cytokinin biosynthesis genes. Under blast infection 7 cliques were associated with all 37 up-regulated HB loci, whereas under drought condition a single clique was associated with 25 (out of 26) up-regulated HB loci and another was associated with 22 up-regulated HB loci. These 9 cliques were composed of 10 CREs of which 6 were common to both blast and drought loci (Table 1). The common CREs are binding sites for GATA, GT1, RAV1, ARR1, DOF, and WRKY71 transcription factors of these GT1, RAV1, and WRKY71 are reported to be associated with stress response. GT1 influences SA inducible expression of pathogenesis related—PR proteins (Buchel et al., 1999). RAV1 is a positive regulator of leaf senescence (Woo et al., 2010) and WRKY71 is related to both biotic and abiotic stress response (Xie et al., 2006; Liu et al., 2007). These observations indicate that transcription factors that bind to these CREs may act as master regulators of defense and stress response (Table 1).
Table 1. List of CREs and corresponding transcription factors (TF) associated with up-regulated hormone biosynthesis loci under blast (Magnaporthe oryzae infection) and drought stress.
In blast infected leaves BR HB loci were most abundant in cliques (94.7%) followed by CK (93.7%), JA (85.5%) and SA (75.0%) (Figure S1). In case of drought SA up-regulated loci were not observed in expression data so were completely absent but BR (82.8%), CK (79.3%) and JA (86.2%) showed similar trend (Figure S1). Loci of ethylene biosynthesis pathway were found to be associated with 50 cliques of which 48 were associated with JA and 46 were associated with both JA and SA. This indicates the existence of some master transcription factors which regulate both ET and JA biosynthesis genes and all three (ET, JA, SA, HB genes).
Expression analysis of HB pathway associated genes under Magnaporthe oryzae infection and drought stress, demonstrates that some genes are differentially up-regulated under both conditions whereas some are exclusively up-regulated in either condition. These findings indicated that differential regulation of hormone biosynthesis may contribute to stress response. Analysis of CREs in promoters of these up-regulated HB genes allowed us to extract insights about regulation of HB pathways and CREs (and/or transcription factors) that may integrate and modulate hormone biosynthesis.
4.1. Cytokinin Conjugation is Regulated by ABA
Role of cytokinins in regulating plant growth and development by promoting cell division and stimulating sink strength is well established. Recently evidence has accumulated demonstrating importance of CK in environmental stress response as well. Endogenous CK levels were reported to alter stress tolerance of plants and conversely stress conditions were found to change CK content in plants (Kudoyarova et al., 2007; Rivero et al., 2007; Albacete et al., 2008; Ghanem et al., 2008; Havlová et al., 2008; Werner et al., 2010; Nishiyama et al., 2011). Studies of stress response in CK-deficient plants showed that CKs negatively regulate stress signaling (Werner et al., 2010; Nishiyama et al., 2011). Down-regulation of CK active content is associated with prolonged drought and salt stress (Kudoyarova et al., 2007; Albacete et al., 2008; Ghanem et al., 2008; Nishiyama et al., 2011). In agreement with above observations we found that all CK HB genes up-regulated in present study (16 in blast and 11 in drought) were glucosyltransferases, which are known to be involved in conversion of active CKs to inactive conjugates (Table S1).
We estimated and analyzed over-represented CREs in the promoters of all up-regulated CK HB loci. We found that 3 (out of 11) and 8 (out of 16) CK HB loci in drought and blast infected leaves were also associated with various ABA response CREs respectively (Table S2). As all the up-regulated CK HB loci are involved in CK conjugation, this indicates that ABA may modulate CK conjugation associated loci to reduce active CK. This is in agreement with previous experimental studies which show ABA regulates key genes in CK biosynthesis (Takei et al., 2004; Sakakibara, 2006). Incidentally we found 9-cis-epoxycarotenoid dioxygenase 1 (Os07g05940.1), a key gene of ABA biosynthesis (Xiong and Zhu, 2003) to be up-regulated in present study as well. 9-cis-epoxycarotenoid dioxygenase converts 9-cis-neoxanthin to xanthoxin in chloroplast (Schwartz et al., 1997). Xanthoxin is then transported to cytoplasm and converted to active ABA hormone through a two step reaction (Xiong and Zhu, 2003). Thus, under both blast and drought CK conjugation and ABA production loci are up-regulated simultaneously. To dissect ABA-CK interaction further we looked for over-represented CREs in promoter of 9-cis-epoxycarotenoid dioxygenase and found a CK-response element (ARR1AT) together with five ABA response elements to be enriched. These findings suggest that there may be a crosstalk between CK and ABA HB pathway regulation (Figure 4). This corresponds well with analysis of CK mutants and regulation of CK metabolic genes, which revealed its antagonism and cross-talk with ABA (Ha et al., 2012). They also observed that CK deficient plants show increased sensitivity to ABA and up-regulation of stress and/or ABA responsive genes, whereas these effects are negated with elevation in CK levels. Hence active ABA modulates CK conjugation to reduce active CK, which in turn may lead to hypersensitivity to ABA and increased stress tolerance. Conversely active CK as well as active ABA may modulate ABA production and downstream ABA responses (Figure 4).
Figure 4. Schematic representation of antagonistic interaction between cytokinin (CK) and abscisic acid (ABA) biosynthesis pathways. CK response element (CKRE) and ABA response element (ABARE) are present in promoter of 9-cis-epoxycarotenoid dioxygenase, a key enzyme of ABA biosynthesis pathway. These response elements may regulate ABA biosynthesis in response to CK and ABA signals. Conversely ABARE present in promoter of Glucosyltransferases, enzymes responsible for CK-conjugation, may modulate CK-conjugation in response to ABA signaling. CK conjugation leads to CK inactivation. Solid lines depict biosynthesis of hormones and responsible enzymes. Dashed lines depict regulation of hormone biosynthesis genes through CRE elements.
4.2. DREB2A Co-ordinates Ethylene and Jasmonic Acid Biosynthesis under Drought and Blast Attack
ET and JA interactions and their response play important role in abiotic and biotic stress management (Wang et al., 2002). In case of blast we found two ET HB genes and six JA HB genes to be up-regulated; similarly in case of drought we found two ET and four JA HB genes to be up-regulated. We compared cliques associated with these JA and ET HB loci and noted that these 7 cliques (b35:SURECOREATSULTR11 DRECRTCOREAT RAV1AAT CGCGBOXAT, b71:SURECOREATSULTR11 RAV1AAT CBFHV,d2:CBFHV LTRE1HVBLT49 ROOTMOTIFTAPOX1, d9:DRECRTCOREAT CACTFTPPCA1 GATABOX, d10:GT1GMSCAM4 CBFHV LTRE1HVBLT49, d11:POLLEN1LELAT52 CBFHV LTRE1HVBLT49, d23:CELLCYCLESC POLLEN1LELAT52 CBFHV) contain DREB protein binding elements (DRECRTCOREAT and CBFHV) under both conditions. None of dehydration-responsive elements (DREs) were found in cliques associated with up-regulated ABA loci. Moreover the transcription factor, DREB2A, that binds to these DREs (DRECRTCOREAT and CBFHV) was also found to be up-regulated under both conditions.
DREB proteins are one of the best characterized transcription factors involved in abiotic stress response (Lata and Prasad, 2011). In rice DREB2A transcription factor (Os01g07120.2) functions are considered important under osmotic stress (Nakashima et al., 2000; Sakuma et al., 2006). DREB2A expression is induced strongly by high salt stress (Lee et al., 2010) and dehydration but weakly by exogenous application of ABA (Dubouzet et al., 2003). C2H2 zinc finger family protein (ZPF) transcription factors are known to modulate DREB expression (Figueiredo et al., 2012). ZPF179 (Os01g62190.1) and ZPF252 (Os12g39400.1) expression is known to be induced by salt and drought stress (Xu et al., 2008; Sun et al., 2010). Transgenic rice constitutively expressing ZPF179 exhibits up-regulation of DREB2A. In present study under both conditions DREB2A (3.07 fold in blast and 5.35 fold in drought) as well as ZPF179 (6.32 fold in blast and 7.72 fold in drought) and ZPF252 (2.44 fold in blast and 5.77 fold in drought) were up-regulated. However, the expressions of all three DREB2A, ZPF179, and ZPF252 were higher in case of drought than in case of blast. We found that ABA response elements were enriched in ZFP179 (ABRELATERD1, ABREOSRAB21) and ZFP252 (PROXBBNNAPA, RYREPEATBNNAPA) but not in DREB2A. DREB2A promoter contained several DREB binding CREs (DRECRTCOREA, CBFHV, DRE2COREZMRAB17) indicating that DREB2A expression may be regulated by DREB and that ABA may not modulating DREB2A directly but it may do so through ZPFs (Figure 5).
Figure 5. Schematic representation of regulation of ethylene (ET) and jasmonic acid (JA) biosynthesis by DREB transcription factors. DREB expression is regulated by C2H2 zinc finger transcription factors. C2H2 zinc finger transcription factors promoter contains abscisic acid response CRE (ABARE) and DREB binding CRE (DRE). Solid lines indicate transcriptional regulation of genes from previous literature and CRE occurrences. Dashed lines indicate transcriptional regulation of genes from CRE occurrences in present study.
4.3. Cytokinin and Brassinosteroids Play a Role in Stress Response and May Interact Directly
BR and Cytokinin are both growth promoting hormones that are known to have a role in stress management (Peleg et al., 2011; Jiang et al., 2013; Nahar et al., 2013). ABA-CK and ABA-BR are known to interact antagonistically (Tran et al., 2007; Zhang et al., 2009; Nishiyama et al., 2011). In case of blast we found six BR HB genes, sixteen CK HB genes and 9-cis-epoxycarotenoid dioxygenase (an enzyme responsible for the first committed reaction in ABA biosynthesis) to be up-regulated. In drought we found five BR HB genes, eleven CK HB genes and ABA HB gene 9-cis-epoxycarotenoid dioxygenase to be up-regulated. These results indicate that CK, BR and ABA may be important for stress response under both conditions.
BR regulates cell expansion along longitudinal axis while CK regulates cell expansion along transverse axis (Kitanaga et al., 2006). BR is known to interact with different phytohormones (Choudhary et al., 2012; Gruszka, 2013) and to promote tolerance in plants to a wide range of stresses, including heat, cold, drought and salinity, (Dhaubhadel et al., 1999, 2002; Kagale et al., 2007; Koh et al., 2007). In Arabidopsis seedlings exogenous application of BR enhanced expression of 1-aminocyclopropane-1-carboxylatesynthase (ACS), a key gene of ethylene biosynthesis pathway (Muday et al., 2012). BR also prevents 26S mediated ubiquitination and degradation of ACS protein (Hansen et al., 2009). In rice BR enhances JA level as well as promotes expression of thionin genes encoding antimicrobial peptides (Kitanaga et al., 2006). Role of CK in stress management has been recently revealed and is only partially understood. In case of rice blast disease CK levels elevated in leaf around infection site and was found to act synergistically with SA to activate PR genes (Jiang et al., 2013). However, SA levels did not alter with accumulation of CK during blast infection (Jiang et al., 2013). Rice plants over-expressing isopentyltransferase, (IPT), a key CK biosynthesis gene resulted in delayed stress response and drought tolerance. Moreover the elevated level of CKs under water stress corresponded well with decreased expression of JA biosynthesis gene OPR2 and increased expression of BR biosynthesis genes Delta14-sterol reductase (DWF5) and C-8 sterol isomerase (HYD1) (Peleg et al., 2011). ABA is known to interact antagonistically with both CK (Tran et al., 2007; Nishiyama et al., 2011) and BR (Zhang et al., 2009). The primary signaling outputs of BR are regulated by ABA signaling (Zhang et al., 2009), so whether effect of CK on BR biosynthesis is direct or mediated by ABA is not clear (Peleg et al., 2011).
CK and BR HB loci promoters were inspected for known ABA (CBFHV, DRE, ABRE, DPBF, LTRE) and CK (CPBCSPOR, ARR1AT) response elements. As BR-response elements are not present in PLACE database we scanned all the CK and BR HB up-regulated loci promoters for the presence of BR response element (CGTG[T/C]G) (Reinhold et al., 2011). BR response element was found it to be present in 3 BR loci and 2 CK loci in drought condition; and in most of CK and BR loci under blast attack. In promoters of up-regulated BR HB genes Os05g29990.1 (NAD dependent epimerase) and Os08g17500.1 (cinnamoyl-CoA reductase), CK (CPBCSPOR, ARR1AT) response elements were present. These loci are associated with initial steps of BR biosynthesis. Os08g17500.1 and another loci with similar annotation Os02g08420.1 were also associated with ABA responsive element DRE2COREZMRAB17. Therefore, our results indicate that in case of blast CK and ABA both may interact directly with BR HB loci likewise BR may also interact directly with CK HB loci (Figure 6).
Figure 6. Model depicting interaction among cytokinin (CK), abscisic acid (ABA) and brassinosteroid (BR). Lines with bars represent antagonistic interactions. Solid lines depict signal network interaction and dashed lines depict regulation of hormone biosynthesis genes through CRE elements. CRE elements that may be responsible for the interactions are indicated above the arrows. Abscisic acid (ABA) interacts antagonistically with both cytokinin (CK) and brassinosteroid (BR). CK response elements (CPBCSPOR and ARR1AT) and ABA response element (DRE2COREZMRAB17) are enriched in promoters of BR biosynthesis associated genes. BR response element (CGTG[T/C]G) is present in promoters of CK biosynthesis associated genes.
4.4. Salicylic Acid and Jasmonic Acid are Unlikely to Effect Each-Others Biosynthesis Directly
SA and JA are stress hormone, SA is related to plant defense against pathogen attack whereas JA response is associated with both abiotic and biotic stress. Promoters of SA and JA HB loci were found to contain several common CREs (ASF1MOTIFCAMV, CGCGBOXAT, SURECOREATSULTR11, WRKY71OS, ACGTCBOX, CTRMCAMV35S). These common CRE were associated with 72.4 % of cliques (combination of CREs) for blast infection. These observations indicate that JA and SA HB loci may be co-regulated under fungal attack. Previous research has shown that interactions between SA and JA are complex and condition specific (Thaler et al., 2012), however antagonistic interactions seem to prevail (Pieterse et al., 2009; Tsuda et al., 2009). Cross-talk between SA and JA is proposed to play a role in fine-tuning of plant defense response (Pieterse et al., 2012). In a study expression of thirteen OPR (12-oxo-phytodienoic acid reductase—enzyme that converts linolenic acid to jasmonic acid) isogenes was monitored under different stress conditions. It was found that OPR1, OPR2, OPR10, and OPR13 in rice were induced by SA (Jang et al., 2009). SA may regulate JA biosynthesis through intermediaries i.e., NPR1 or directly.
NPR1 is known to play an important role in SA-JA antagonism, NPR1 is required for functional salicylate signaling in Arabidopsis-Pseudomonas syringae pathosystem. When NPR1 is impaired through a mutation, it resulted in enhanced expression of LOX2, an important enzyme in JA biosynthesis (Spoel et al., 2003). We also found NPR1 to be up-regulated (fold change = 2.949) under blast infection. NPR1 is known to regulate expression of genes by enhancing DNA binding of TGA transcription factors (Després et al., 2003). We did find TGA binding CRE (ASF1MOTIFCAMV) in promoters of both SA and JA HB loci but TGA expression was not up-regulated under present condition. Thus, SA through NPR1 facilitated TGA binding to ASF1MOTIFCAMV may regulate SA and JA biosynthesis. However, either this does not require up-regulation of TGA or the time point, experimental conditions and data analysis strategy used for this study were inadequate to capture expression of TGA.
JA HB genes LOX and four OPR were found to be up-regulated but none of them had any SA-responsive CRE in their promoters. Similarly out of three only one SA HB gene Os04g43800.1 contained one JA responsive CRE (T/GBOXATPIN2). Os04g43800.1 is a phenylalanine ammonia-lyase which is also associated with secondary metabolite phenylpropanoid derivative synthesis. Analysis of CRE in up-regulated JA and SA HB genes under both drought and blast revealed no CRE (except one JA responsive CRE in SA loci) that can be tied to SA or JA response directly, so it is highly unlikely that these hormones directly affect each-other's biosynthesis.
Our investigations into CRE-mediated regulation of differentially up-regulated hormone biosynthesis associated loci under blast infection and drought stress revealed that hormone biosynthesis pathways are subjected to hormone cross-talk. Plant response to blast and drought is associated with transcriptional modulation of biosynthesis genes of several hormones and these hormones modulate bio-availability of each other through a complex regulatory network. Abscisic acid and Brassinosteroid regulates Cytokinin by modulating its conjugation. Conversely Brassinosteroid biosynthesis is effected by both Abscisic acid and Cytokinin. Through our analysis we found that DREB transcription factors may be important in modulation of both Jasmonic acid and Ethylene biosynthesis. Moreover we found that though DREB response are known to be ABA independent but at least DREB2A response may be regulated by ABA signaling through C2H2 zinc finger protein transcription factors. Interestingly we found no indication of direct interaction of Jasmonic acid or Salicylic acid with each other's biosynthesis pathways. Thus, multiple hormones may modulate several hormone biosynthesis pathways through a complex regulatory networks, where biosynthesis of one hormone is effected by several other contributing hormones.
AD, RG, and SK conceived the idea and designed the study. AD and RG collected the data. AD developed and executed computational pipelines and took part in result analysis. RG performed literature survey and result analysis. AD and RG wrote the manuscript. SK edited the manuscript and supervised the study. AD and RG contributed equally to this study. All authors read and approved the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
AD acknowledges UGC (RFSMS), Government of India, for financial assistance. RG acknowledges Center of Excellence in Systems Biology and Biomedical Engineering for financial assistance. This work is supported by Center of Excellence in Systems Biology and Biomedical Engineering, University of Calcutta.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01303
Table S1. List of hormone biosynthesis associated loci with fold change >2 and adjusted p-value < 0.05 for Magnaporthe oryzae infection and drought stress in rice leaves.
Table S2. List of cytokinin biosynthesis associated loci, their known annotations, CREs found to be enriched in their promoters, ABA responsive CREs if any present among the enriched CREs and corresponding transcription factors.
Figure S1. Bar graph representation of different cliques and associated loci. (A) Cliques of CREs associated with hormone biosynthesis (HB) genes up-regulated under Magnaporthe oryzae infection. (B) Cliques of CREs associated with hormone biosynthesis genes up-regulated under drought stress. Horizontal axis represents the cliques and vertical axis represents number of loci a particular clique is associated with. Loci belonging to different HB pathways are represented by different colors. Hormone abbreviations: CK, cytokinin; SA, salicylic acid; JA, jasmonic acid; ET, ethylene; ABA, abscisic acid ; GA, gibberellic; BR, brassinosteroid.
Adamcsek, B., Palla, G., Farkas, I. J., Derényi, I., and Vicsek, T. (2006). Cfinder: locating cliques and overlapping modules in biological networks. Bioinformatics 22, 1021–1023. doi: 10.1093/bioinformatics/btl039
Albacete, A., Ghanem, M. E., Martínez-Andújar, C., Acosta, M., Sánchez-Bravo, J., Martínez, V., et al. (2008). Hormonal changes in relation to biomass partitioning and shoot growth impairment in salinized tomato (Solanum lycopersicum L.) plants. J. Exp. Bot. 59, 4119–4131. doi: 10.1093/jxb/ern251
Buchel, A. S., Brederode, F. T., Bol, J. F., and Linthorst, H. J. (1999). Mutation of GT-1 binding sites in the Pr-1A promoter influences the level of inducible gene expression in vivo. Plant Mol. Biol. 40, 387–396. doi: 10.1023/A:1006144505121
Choudhary, S. P., Yu, J.-Q., Yamaguchi-Shinozaki, K., Shinozaki, K., and Tran, L.-S. P. (2012). Benefits of brassinosteroid crosstalk. Trends Plant Sci. 17, 594–605. doi: 10.1016/j.tplants.2012.05.012
Després, C., Chubak, C., Rochon, A., Clark, R., Bethune, T., Desveaux, D., et al. (2003). The arabidopsis NPR1 disease resistance protein is a novel cofactor that confers redox regulation of dna binding activity to the basic domain/leucine zipper transcription factor TGA1. Plant Cell 15, 2181–2191. doi: 10.1105/tpc.012849
Dhaubhadel, S., Browning, K. S., Gallie, D. R., and Krishna, P. (2002). Brassinosteroid functions to protect the translational machinery and heat-shock protein synthesis following thermal stress. Plant J. 29, 681–691. doi: 10.1046/j.1365-313X.2002.01257.x
Dhaubhadel, S., Chaudhary, S., Dobinson, K. F., and Krishna, P. (1999). Treatment with 24-epibrassinolide, a brassinosteroid, increases the basic thermotolerance of brassica napus and tomato seedlings. Plant Mol. Biol. 40, 333–342. doi: 10.1023/A:1006283015582
Dubouzet, J. G., Sakuma, Y., Ito, Y., Kasuga, M., Dubouzet, E. G., Miura, S., et al. (2003). Osdreb genes in rice, Oryza sativa L., encode transcription activators that function in drought-, high-salt-and cold-responsive gene expression. Plant J. 33, 751–763. doi: 10.1046/j.1365-313X.2003.01661.x
Figueiredo, D. D., Barros, P. M., Cordeiro, A. M., Serra, T. S., Lourenço, T., Chander, S., et al. (2012). Seven zinc-finger transcription factors are novel regulators of the stress responsive gene OsDREB1B. J. Exp. Bot. 63, 3643–3656. doi: 10.1093/jxb/ers035
Frigerio, M., Alabadí, D., Pérez-Gómez, J., García-Cárcel, L., Phillips, A. L., Hedden, P., et al. (2006). Transcriptional regulation of gibberellin metabolism genes by auxin signaling in arabidopsis. Plant Physiol. 142, 553–563. doi: 10.1104/pp.106.084871
Gentleman, R. C., Carey, V. J., Bates, D. M., Bolstad, B., Dettling, M., Dudoit, S., et al. (2004). Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 5:R80. doi: 10.1186/gb-2004-5-10-r80
Ghanem, M. E., Albacete, A., Martínez-Andújar, C., Acosta, M., Romero-Aranda, R., Dodd, I. C., et al. (2008). Hormonal changes during salinity-induced leaf senescence in tomato (Solanum lycopersicum L.). J. Exp. Bot. 59, 3039–3050. doi: 10.1093/jxb/ern153
Gruszka, D. (2013). The brassinosteroid signaling pathwaynew key players and interconnections with other signaling networks crucial for plant development and stress tolerance. Int. J. Mol. Sci. 14, 8740–8774. doi: 10.3390/ijms14058740
Ha, S., Vankova, R., Yamaguchi-Shinozaki, K., Shinozaki, K., and Tran, L.-S. P. (2012). Cytokinins: metabolism and function in plant adaptation to environmental stresses. Trends Plant Sci. 17, 172–179. doi: 10.1016/j.tplants.2011.12.005
Havlová, M., Dobrev, P. I., Motyka, V., Štorchová, H., Libus, J., Dobrá, J., et al. (2008). The role of cytokinins in responses to water deficit in tobacco plants over-expressing trans-zeatin O-glucosyltransferase gene under 35S or SAG12 promoters. Plant Cell Environ. 31, 341–353. doi: 10.1111/j.1365-3040.2007.01766.x
Hughes, J. D., Estep, P. W., Tavazoie, S., and Church, G. M. (2000). Computational identification of cis-regulatory elements associated with groups of functionally related genes in saccharomyces cerevisiae. J. Mol. Biol. 296, 1205–1214. doi: 10.1006/jmbi.2000.3519
Jang, S., Cho, K., Shibato, J., Han, O., Iwahashi, H., Tamogami, S., et al. (2009). Rice osoprs: transcriptional profiling responses to diverse environmental stimuli and biochemical analysis of osopr1. J. Plant Biol. 52, 229–243. doi: 10.1007/s12374-009-9022-1
Jiang, C.-J., Shimono, M., Sugano, S., Kojima, M., Liu, X., Inoue, H., et al. (2013). Cytokinins act synergistically with salicylic acid to activate defense gene expression in rice. Mol. Plant-Microbe Interact. 26, 287–296. doi: 10.1094/MPMI-06-12-0152-R
Kagale, S., Divi, U. K., Krochko, J. E., Keller, W. A., and Krishna, P. (2007). Brassinosteroid confers tolerance in Arabidopsis thaliana and Brassica napus to a range of abiotic stresses. Planta 225, 353–364. doi: 10.1007/s00425-006-0361-6
Kato, M., Hata, N., Banerjee, N., Futcher, B., and Zhang, M. Q. (2004). Identifying combinatorial regulation of transcription factors and binding motifs. Genome Biol. 5:R56. doi: 10.1186/gb-2004-5-8-r56
Kitanaga, Y., Jian, C., Hasegawa, M., Yazaki, J., Kishimoto, N., Kikuchi, S., et al. (2006). Sequential regulation of gibberellin, brassinosteroid, and jasmonic acid biosynthesis occurs in rice coleoptiles to control the transcript levels of anti-microbial thionin genes. Biosci. Biotechnol. Biochem. 70, 2410–2419. doi: 10.1271/bbb.60145
Koh, S., Lee, S.-C., Kim, M.-K., Koh, J. H., Lee, S., An, G., et al. (2007). T-DNA tagged knockout mutation of rice OsGSK1, an orthologue of arabidopsis BIN2, with enhanced tolerance to various abiotic stresses. Plant Mol. Biol. 65, 453–466. doi: 10.1007/s11103-007-9213-4
Kohli, A., Sreenivasulu, N., Lakshmanan, P., and Kumar, P. P. (2013). The phytohormone crosstalk paradigm takes center stage in understanding how plants respond to abiotic stresses. Plant Cell Rep. 32, 945–957. doi: 10.1007/s00299-013-1461-y
Kudoyarova, G. R., Vysotskaya, L. B., Cherkozyanova, A., and Dodd, I. C. (2007). Effect of partial rootzone drying on the concentration of zeatin-type cytokinins in tomato (Solanum lycopersicum L.) xylem sap and leaves. J. Exp. Bot. 58, 161–168. doi: 10.1093/jxb/erl116
Lee, S. J., Kang, J. Y., Park, H. J., Kim, M. D., Bae, M. S., Choi, H. I., et al. (2010). DREB2C interacts with ABF2, a bZIP protein regulating abscisic acid-responsive gene expression, and its overexpression affects abscisic acid sensitivity. Plant Physiol. 153, 716–727. doi: 10.1104/pp.110.154617
Murakami, K., Imanishi, T., Gojobori, T., and Nakai, K. (2008). Two different classes of co-occurring motif pairs found by a novel visualization method in human promoter regions. BMC Genomics 9:112. doi: 10.1186/1471-2164-9-112
Nahar, K., Kyndt, T., Hause, B., Höfte, M., and Gheysen, G. (2013). Brassinosteroids suppress rice defense against root-knot nematodes through antagonism with the jasmonate pathway. Mol. Plant-Microbe Interact. 26, 106–115. doi: 10.1094/MPMI-05-12-0108-FI
Nakashima, K., Shinwari, Z. K., Sakuma, Y., Seki, M., Miura, S., Shinozaki, K., et al. (2000). Organization and expression of two arabidopsis DREB2 genes encoding dre-binding proteins involved in dehydration-and high-salinity-responsive gene expression. Plant Mol. Biol. 42, 657–665. doi: 10.1023/A:1006321900483
Nishiyama, R., Watanabe, Y., Fujita, Y., Le, D. T., Kojima, M., Werner, T., et al. (2011). Analysis of cytokinin mutants and regulation of cytokinin metabolic genes reveals important regulatory roles of cytokinins in drought, salt and abscisic acid responses, and abscisic acid biosynthesis. Plant Cell 23, 2169–2183. doi: 10.1105/tpc.111.087395
Ouyang, S., Zhu, W., Hamilton, J., Lin, H., Campbell, M., Childs, K., et al. (2007). The TIGR rice genome annotation resource: improvements and new features. Nucleic Acids Res. 35(Suppl. 1):D883–D887. doi: 10.1093/nar/gkl976
Peleg, Z., Reguera, M., Tumimbang, E., Walia, H., and Blumwald, E. (2011). Cytokinin-mediated source/sink modifications improve drought tolerance and increase grain yield in rice under water-stress. Plant Biotechnol. J. 9, 747–758. doi: 10.1111/j.1467-7652.2010.00584.x
Pieterse, C. M., Van der Does, D., Zamioudis, C., Leon-Reyes, A., and Van Wees, S. C. (2012). Hormonal modulation of plant immunity. Ann. Rev. Cell Dev. Biol. 28, 489–521. doi: 10.1146/annurev-cellbio-092910-154055
Reinhold, H., Soyk, S., Šimková, K., Hostettler, C., Marafino, J., Mainiero, S., et al. (2011). β-amylase–like proteins function as transcription factors in arabidopsis, controlling shoot growth and development. Plant Cell 23, 1391–1403. doi: 10.1105/tpc.110.081950
Rivero, R. M., Kojima, M., Gepstein, A., Sakakibara, H., Mittler, R., Gepstein, S., et al. (2007). Delayed leaf senescence induces extreme drought tolerance in a flowering plant. Proc. Natl. Acad. Sci. U.S.A. 104, 19631–19636. doi: 10.1073/pnas.0709453104
Robert-Seilaniantz, A., Grant, M., and Jones, J. D. (2011). Hormone crosstalk in plant disease and defense: more than just jasmonate-salicylate antagonism. Ann. Rev. Phytopathol. 49, 317–343. doi: 10.1146/annurev-phyto-073009-114447
Sakuma, Y., Maruyama, K., Osakabe, Y., Qin, F., Seki, M., Shinozaki, K., et al. (2006). Functional analysis of an arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell 18, 1292–1309. doi: 10.1105/tpc.105.035881
Schwartz, S. H., Tan, B. C., Gage, D. A., Zeevaart, J. A., and McCarty, D. R. (1997). Specific oxidative cleavage of carotenoids by VP14 of maize. Science 276, 1872–1874. doi: 10.1126/science.276.5320.1872
Smit, A. F. A., Hubley, R., and Green, P. (2013–2015). RepeatMasker Open-4.0. Available online at: http://www.repeatmasker.org
Spoel, S. H., Koornneef, A., Claessens, S. M., Korzelius, J. P., Van Pelt, J. A., Mueller, M. J., et al. (2003). NPR1 modulates cross-talk between salicylate-and jasmonate-dependent defense pathways through a novel function in the cytosol. Plant Cell 15, 760–770. doi: 10.1105/tpc.009159
Sudarsanam, P., Pilpel, Y., and Church, G. M. (2002). Genome-wide co-occurrence of promoter elements reveals a cis-regulatory cassette of rrna transcription motifs in saccharomyces cerevisiae. Genome Res. 12, 1723–1731. doi: 10.1101/gr.301202
Sun, S.-J., Guo, S.-Q., Yang, X., Bao, Y.-M., Tang, H.-J., Sun, H., et al. (2010). Functional analysis of a novel Cys2/His2-type zinc finger protein involved in salt tolerance in rice. J. Exp. Bot. 61, 2807–2818. doi: 10.1093/jxb/erq120
Takei, K., Yamaya, T., and Sakakibara, H. (2004). Arabidopsis CYP735A1 and CYP735A2 encode cytokinin hydroxylases that catalyze the biosynthesis of trans-zeatin. J. Biol. Chem. 279, 41866–41872. doi: 10.1074/jbc.M406337200
Tran, L.-S. P., Urao, T., Qin, F., Maruyama, K., Kakimoto, T., Shinozaki, K., et al. (2007). Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 104, 20623–20628. doi: 10.1073/pnas.0706547105
Werner, T., Nehnevajova, E., Köllmer, I., Novák, O., Strnad, M., Krämer, U., et al. (2010). Root-specific reduction of cytokinin causes enhanced root growth, drought tolerance, and leaf mineral enrichment in arabidopsis and tobacco. Plant Cell 22, 3905–3920. doi: 10.1105/tpc.109.072694
Woo, H. R., Kim, J. H., Kim, J., Kim, J., Lee, U., Song, I. J., et al. (2010). The RAV1 transcription factor positively regulates leaf senescence in arabidopsis. J. Exp. Bot. 61, 3947–3957. doi: 10.1093/jxb/erq206
Xie, Z., Zhang, Z.-L., Zou, X., Yang, G., Komatsu, S., and Shen, Q. J. (2006). Interactions of two abscisic-acid induced wrky genes in repressing gibberellin signaling in aleurone cells. Plant J. 46, 231–242. doi: 10.1111/j.1365-313X.2006.02694.x
Xu, D.-Q., Huang, J., Guo, S.-Q., Yang, X., Bao, Y.-M., Tang, H.-J., et al. (2008). Overexpression of a tfiiia-type zinc finger protein gene ZFP252 enhances drought and salt tolerance in rice (Oryza sativa L.). FEBS Lett. 582, 1037–1043. doi: 10.1016/j.febslet.2008.02.052
Yu, X., Lin, J., Masuda, T., Esumi, N., Zack, D. J., and Qian, J. (2006). Genome-wide prediction and characterization of interactions between transcription factors in Saccharomyces cerevisiae. Nucleic Acids Res. 34, 917–927. doi: 10.1093/nar/gkj487
Zhang, S., Cai, Z., and Wang, X. (2009). The primary signaling outputs of brassinosteroids are regulated by abscisic acid signaling. Proc. Natl. Acad. Sci. U.S.A. 106, 4543–4548. doi: 10.1073/pnas.0900349106
Keywords: hormone biosynthesis, transcriptional regulation, cis-regulatory element, cross-talk, Magnaporthe, drought
Citation: Deb A, Grewal RK and Kundu S (2016) Regulatory Cross-Talks and Cascades in Rice Hormone Biosynthesis Pathways Contribute to Stress Signaling. Front. Plant Sci. 7:1303. doi: 10.3389/fpls.2016.01303
Received: 05 July 2016; Accepted: 15 August 2016;
Published: 26 August 2016.
Edited by:Alexandre V. Morozov, Rutgers University, USA
Reviewed by:Victor P. Andreev, Arbor Research Collaborative for Health, USA
Kira Vyatkina, Saint Petersburg Academic University, Russia
Copyright © 2016 Deb, Grewal and Kundu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sudip Kundu, email@example.com