Comparative analysis of expressed sequence tags (ESTs) between drought-tolerant and -susceptible genotypes of chickpea under terminal drought stress
© Deokar et al; licensee BioMed Central Ltd. 2011
Received: 22 December 2010
Accepted: 22 April 2011
Published: 22 April 2011
Chickpea (Cicer arietinum L.) is an important grain-legume crop that is mainly grown in rainfed areas, where terminal drought is a major constraint to its productivity. We generated expressed sequence tags (ESTs) by suppression subtraction hybridization (SSH) to identify differentially expressed genes in drought-tolerant and -susceptible genotypes in chickpea.
EST libraries were generated by SSH from root and shoot tissues of IC4958 (drought tolerant) and ICC 1882 (drought resistant) exposed to terminal drought conditions by the dry down method. SSH libraries were also constructed by using 2 sets of bulks prepared from the RNA of root tissues from selected recombinant inbred lines (RILs) (10 each) for the extreme high and low root biomass phenotype. A total of 3062 unigenes (638 contigs and 2424 singletons), 51.4% of which were novel in chickpea, were derived by cluster assembly and sequence alignment of 5949 ESTs. Only 2185 (71%) unigenes showed significant BLASTX similarity (<1E-06) in the NCBI non-redundant (nr) database. Gene ontology functional classification terms (BLASTX results and GO term), were retrieved for 2006 (92.0%) sequences, and 656 sequences were further annotated with 812 Enzyme Commission (EC) codes and were mapped to 108 different KEGG pathways. In addition, expression status of 830 unigenes in response to terminal drought stress was evaluated using macro-array (dot blots). The expression of few selected genes was validated by northern blotting and quantitative real-time PCR assay.
Our study compares not only genes that are up- and down-regulated in a drought-tolerant genotype under terminal drought stress and a drought susceptible genotype but also between the bulks of the selected RILs exhibiting extreme phenotypes. More than 50% of the genes identified have been shown to be associated with drought stress in chickpea for the first time. This study not only serves as resource for marker discovery, but can provide a better insight into the selection of candidate genes (both up- and downregulated) associated with drought tolerance. These results can be used to identify suitable targets for manipulating the drought-tolerance trait in chickpea.
Chickpea (Cicer arietinum L.), the fourth most important grain-legume crop, is grown in more than 45 countries, mostly in arid and semiarid zones. Approximately 90% of the crop is grown under rainfed conditions, wherein yield is significantly affected by abiotic stresses such as drought, heat, and cold [1–3]. Drought-related yield losses can occur in 40%-60% of the total chickpea production . Terminal drought, which occurs at the pod filling and seed-developing stage of the crop and increases in severity at the end of the season, is a major constraint to chickpea production [1, 5, 6]. The identification of differentially expressed genes between 2 genotypes differing in drought tolerance and a set of their progenies can therefore be an important indicator of drought-associated genes in chickpea.
Functional genomics approaches have been used in recent years to understand the stress-responsive mechanism in plants. Candidate genes involved in drought tolerance mechanisms have been identified, characterized, and assessed for their comparative transcriptional activity by using whole-genome sequencing or expressed sequence tag (EST) libraries. Several functional genomics studies have been performed in chickpea to identify the abiotic stress-responsive transcripts by approaches such as suppression subtractive hybridization (SSH), Super serial analysis of gene expression (SuperSAGE), microarray, and EST sequencing [7–9]. Additional file 1 summarizes results of previous studies on identifying ESTs associated with drought stress in chickpea.
SSH has been widely used to compare patterns of gene expression in tissues under different conditions. However, it has not yet been used to identify differentially expressed transcripts (both up- and downregulated) in chickpea in response to drought stress at the flowering stage of plants. In all earlier studies, except the one by Varshney et al. , water stress was imposed by either completely withdrawing water or allowing uprooted young seedlings to wilt at room temperature. However, under field conditions, water stress progresses gradually and a similar type of stress is simulated in the laboratory by the "dry down experiment," which allows comparison of different genotypes and their response toward drought . Moreover, stress response of a plant at the seedling stage can be very different from that at the reproductive stage, the latter being an important and yield-determining stage in chickpea.
In the present study, we constructed several reciprocal SSH libraries by using drought-tolerant and -susceptible genotypes as well as extreme recombinant inbred lines (RILs) for the high root biomass (HRB) and low root biomass (LRB) under terminal drought stress. This approach differs from that used in earlier studies in the following aspects: (1) use of 2 chickpea genotypes differing in their drought-tolerance capacity and their RIL progenies; (2) drought stress imposed at the flowering stage in a gradual manner by the dry down method; (3) plant samples analyzed when each plant experienced the same amount of stress, as judged by their transpiration ratio; and (4) reciprocal subtraction of transcripts by using control and stress conditions as well as susceptible and tolerant genotypes to enable a good comparison and identify both up- and downregulated genes. Thus, the EST set we used is novel and represents genes that are up and downregulated in response to terminal drought stress, and can thereby help several genes that have not been shown to be previously associated with drought stress in chickpea. The differentially expressed ESTs were analyzed using macro-array, northern blotting, and quantitative PCR.
The drought-tolerant characteristics of chickpea line ICC 4958 and drought-susceptible characteristics of ICC 1882 have been attributed to their large and prolific and small root system, respectively. An RIL mapping population (264 RILs) of ICC 4958 (large root) and ICC 1882 (small root) has been developed and phenotyped at the International Crops Research Institute for the Semi-Arid Tropics (ICRISAT), Patancheru (17° 30' N; 78° 16' E; altitude 549 m). The root phenotyping experiment was conducted in PVC cylinders with 18 cm diameter and 120 cm height, filled with soil-sand mixture in open field conditions. Plants were sampled at 35 days after sowing and different measurements were recorded as described by Kashiwagi et al. . Ten RILs for extreme phenotype of high root biomass and low root biomass were selected on the basis of phenotypic evolution [Varshney et al. unpublished] to prepare bulk cDNA SSH libraries.
Dry down procedure
Dry down, a gradual and progressive water deficit stress, was given to plants . Experiments were conducted in triplicate in a glass house receiving natural solar radiation, with air temperature regulated between 23°C and 28°C (night/day). Seeds of ICC 4958, ICC 1882, 10 RILs each for HRB and LRB, were sown in plastic pots of 8-in. diameter. Water stress (WS) treatment was initiated 35 days after the emergence of plants. All pots were saturated with water and left overnight to drain excess water. Next day, the surfaces of pots were covered with plastic beads to prevent water loss through the soil surface. Weight (in g) of individual pots was recorded daily in the morning at approximately 10.30 h. Daily loss of water through transpiration was calculated as the difference in pot weight on the current day from that on the previous day. Control plants were maintained at approximately 80% field capacity by daily compensation of water loss due to transpiration. To expose WS plants to a progressive water deficit, they were allowed to lose a maximum of 80 g of water per day; any additional loss was compensated by adding water to the pots. The transpiration of each plant was then calculated as the difference in its weight on successive days, plus water added on the previous day. Transpiration data were analyzed as described previously . Well watered (WW) pots were maintained at a normalized transpiration ratio (NTR) value of 1 and WS treatment was continued until the ratio of the transpiration of the stressed plant to the average transpiration of WW plants reached ≤0.1, that is, when the transpiration of WS plants was <1% of the WW plants, a stage defined as the endpoint for the water deficit treatment . WS plants reached this stage in 10 to13 day of initiation of stress treatment. At this stage, shoot and the root tissues from WW and WS plants were separately harvested, frozen in liquid nitrogen, and stored at -80°C for RNA extraction.
RNA and mRNA isolation
Total RNA was isolated by using the Trizol reagent (Invitrogen, Carlsbad, CA), and mRNA was further isolated by using the PolyATract mRNA Isolation System (Promega, Madison, WI). To construct bulk libraries, equal amounts of total RNA (100 μg from each RIL) isolated from 10 RILs of extreme HRB and 10 RILs of extreme LRB were pooled separately and used for mRNA isolation.
Suppression Subtractive Hybridization (SSH)
Summary of drought responsive SSH libraries and ESTs.
Name of Library
No of clones
Total no of ESTs
Type of transcripts clones
Up regulated in shoot tissue under drought stress
Down regulated in shoot tissue under drought stress
Up regulated in root tissue under drought stress
Down regulated in root tissue under drought stress
Up regulated in roots of resistant genotype (ICC 4958) under drought stress
Down regulated in roots of resistant genotype (ICC 4958) under drought stress
Up regulated in roots of extremes bulks of RILs of HRB under drought stress
Down regulated in roots of extremes bulks of RILs of HRB under drought stress
All sequences were checked for quality and then analyzed by Seqman™ II 5.08 (DNASTAR, Inc.. Lasergene Gene Corporation, Ann Arbor, MI) to detect and remove pGEMT-Easy vector sequences. A Perl script EST trimmer  was used to trim adaptors, poly A/T ends. EST sequences which were less than 100 bp long were removed. Manual sequence processing was also performed to confirm results. ESTs from individual libraries were assembled into contigs, using default parameters of CAP3 . Incorporation of ESTs into a contig required at least 95% sequence identity and a minimum 40-bp overlap. ESTs from all 8 libraries also underwent CAP3 analysis to produce a differentially expressed unigene dataset.
The NCBI BLAST program  version 2.2.6 was used to perform BLASTN and BLASTX similarity searches. BLASTN analysis was performed to determine sequence homology at the nucleotide level of this unigenes set with EST databases of Medicago truncatula, Glycine max, Lotus japonicus, and Phaseous vulgaris and also with ESTs of model plant species such as Arabidopsis thaliana, Oryza sativa, and Populus alba downloaded from NCBI. The cutoff expectation (E)-value threshold for BLASTN searches was ≤1e-5. BLASTX was performed against NCBI non-redundant (nr) database using Blast2GO with an E-value cutoff of <1e-06.
Functional categorization and GO enrichment analysis
Functional annotation was performed by using Blast2GO (version 2.2.3) , following the standard procedure of BLASTX for unigenes dataset (parameters: nr database, high scoring segment pair (HSP) cutoff length 33, report 20 hits, maximum E-value 1.0E-3), followed by mapping and annotation (parameters: E-value hit filter 1.0E-6, annotation cutoff 55, GO weight 5, HSP-hit coverage cutoff 20). GO terms were summarized according to their molecular functions, biologic processes, and cellular components. Enzyme mapping of annotated sequences was performed by using direct GO to Enzyme mapping and used to query the Kyoto Encyclopaedia of Genes and Genomes (KEGG) to define the KEGG orthologs (KOs). These KOs were then plotted into the whole metabolic atlas by using the KEGG mapping tool .
GO enrichment analysis was performed by using the Fisher exact test, as implemented in the GOSSIP module  integrated in Blast2GO package. For GO enrichment analysis, all GO terms with a cut-off threshold of pFDR(p) ≤ 0.05 were considered differentially enriched between 2 set of EST libraries. To study the genotype-specific response for ICC 4958 and ICC 1882 under drought stress, GO enrichment analysis was performed between ESTs developed from the SSH libraries AB1-1 and AB2-1, which were constructed to identify transcripts induced in response to drought in the tolerant genotype ICC 4958 and the susceptible genotype ICC 1882, respectively.
Macroarray and Northern Hybridization
To screen the differentially expressed ESTs identified in present work, two different macroarray experiments were conducted. In the first experiment, a nylon macroarray in 96-well format, using unigenes from AS1-1 and AS2-1 libraries, was constructed and total RNA from WW and WS plants of ICC 4958 were used to evaluate the differentially expressed unigenes under water stressed condition. Where as in second experiment, a nylon macroarray in a 96-well format, using unigenes from AB1-1 and AB2-1 libraries, was constructed and total RNA from water-stressed ICC 4958 and ICC 1882 were used to evaluate the genotype-specific response under water stress condition.
Equal amounts of purified PCR amplified products (100 ng) was spotted onto nylon membranes (Amersham Pharmacia Biotech, Uppsala, Sweden), using the dot-blot apparatus in 96 formats. Each blot was prepared in duplicate. PCR-amplified products of actin cDNA (GenBank: EU529707) as a housekeeping gene for normalization of the signals between the blots and neomycin phosphotransferase (NPTII) as a negative control for signal background correction were spotted on the membrane and cross-linked using UV. RNA samples were labeled during first-strand cDNA synthesis. Total RNA (5 μg ) was reverse transcribed, using SuperScript III RT enzyme (Superscript II, Life Technologies, Grand Islands, NY) in the presence of α-[32P] dCTP and used as probes. The nylon membrane were prehybridized with formamide hybridization buffer for 42°C for 6 h, the denatured probe was added, and hybridized for 24 h. Washed membranes were exposed to X-ray film (BIOMAX MR Film, Kodak) and developed after 7 days of incubation at -80°C. The image of the developed film was acquired by SYNGENE-G-Box gel documentation and analysis system (Syngene, Synoptics Ltd, Cambridge, UK) and signal intensity of each spot was calculated by the Gene tool software. Transcript levels for each unigenes were calculated as the average intensity from triplicate experiments. The intensity of each spot was normalized with respect to the intensity of actin gene. Change in level of expression was expressed as the expression ratio of normalized signal intensities of respective unigenes in control versus treatments. On the basis of macroarray results, genes exhibiting significant induction were validated by Northern blotting.
For northern blotting total RNA (20 μg) from WW and WS plants was separated by electrophoresis on a 1.2% FA agarose gel and transferred to an Immobilon™-Ny+ membrane (Millipore, USA) following the method of Sambrook et al. . PCR-amplified individual cDNA fragments (amplified with M13 forward and reverse universal sequencing primers) were purified from the agarose gel and used as probes. cDNA-amplified actin (EU529707) was the housekeeping gene control. Probes were labeled with α32P-dCTP, using the DecaLabel™ DNA labeling kit (Fermentas Life Sciences). Northern blots were scanned using a PharosFx Plus PhosphorImager (Biorad).
Quantitative real-time RT PCR
PCR primers for quantitative real-time PCR were designed with the parameters of optimum primer GC content of 50%, primer Tm > 55-65°C, primer length 18-30 nucleotides, and an expected amplicon size of 80-200 bp (see additional file 2 for primer sequences). SYBR green qPCR was performed in 96-well plates, using the Stratagene Mx3000P system and SYBR FAST qPCR Master Mix (2x) Universal (KAPA Biosystems). All qPCR reactions were run in triplicates with a no-template control to check for contaminations. PCR was conducted under the following conditions: 3 min at 95°C (enzyme activation), 40 cycles each of 3 sec at 95°C (denaturation) and 30 s at 60°C (anneal/extend). Finally, a melting curve analysis was performed from 65° to 95°C in increments of 0.5°C, each lasting 5 s, to confirm the presence of a single product and absence of primer-dimers. Two internal controls GAPDH (glyceraldehyde-3-phosphate dehydrogenase, AJ010224) and HSP90 (GR406804) were used to normalize the variations in cDNA samples . Fold changes were calculated by the 2-δδCt method .
Results and discussion
Water stress treatment
A graph of NTR values of ICC 4958, ICC 1882, and 20 RILs during the stress treatment indicates that all parental lines and RILs experienced same degree of stress (Additional file 3). The dry down procedure to impose water stress in pot experiments has been successfully employed in various plant systems, including chickpea [21–25].
Considering that terminal drought is a major constraint in achieving optimal crop yields in chickpea, all experiments were conducted at the flowering stage to identify molecular responses of chickpea under water stress. In many functional genomics studies on drought response in chickpea, drought stress has been induced by withdrawing water supply or by uprooting seedlings and allowing them to wilt at room temperature [26–28]. However, the physiologic and molecular responses to these treatments are likely to be different from those experienced by the plant during natural terminal drought conditions, wherein drought stress is gradual and allows the plant to go through various stages of adaptation. Another major limitation of all these studies is the variation in the quantum of stress experienced by different plants. Depending on their genotype as well as environmental and experimental conditions, plants experience varying degrees of stress when water is withdrawn or they are allowed to wilt for a specified duration. In our study, we sampled ICC 4958 and ICC 1882 and 20 RILs at a stage when they undergo the same degree of stress, as determined by the transpiration ratio.
cDNA SSH libraries
A total of 6432 clones were generated from the 8 SSH libraries, of which 6053 ESTs were sequenced. After a quality check, 5494 high-quality ESTs were obtained (Table 1). Four SSH libraries were constructed from resistance parent ICC 4958. In total, 2034 upregulated and 1620 downregulated ESTs were identified: 753 upregulated ESTs from library AS1-1 (shoot tissue) and 1281 from AR1-1 (root tissue), and 821 downregulated ESTs from AS2-1 (shoot tissue) and 799 from AR2-1 (root tissue). In addition, 2 reciprocal libraries were constructed using root tissues of ICC 4958 and ICC 1882: there were 503 upregulated ESTs from AB1-1 in ICC 4958 and 529 uprgulated ESTs from AB2-1 in ICC 1882. Furthermore, 400 ESTs were generated from library Bulk1-1 (constructed from the bulk of 10 extreme RILs for HRB) and 408 from library Bulk2-1 (constructed from 10 extreme RILs for LRB).
In chickpea, root growth, osmotic adjustment, and stem reserve utilization are associated with drought tolerance. Root traits such as biomass, length, density, and depth have been proposed as drought-avoidance traits under terminal drought conditions [29, 30]. Roots are considered a primary site for stress signal perception, where a signaling mechanism cascade initiates gene expression in response to drought stress. These transcriptional changes can result in successful adaptations, protecting plants against environmental stress . The differentially expressed ESTs identified in our study provide a list of gene regulated in response to terminal drought stress in root tissue of chickpea.
The SSH strategy can be used as an alternative and complementary transcript profiling tool to the GeneChip microarrays, especially to identify novel genes and transcripts present in low abundance . Thus, the SSH technology will have more utility in a system where genome sequence information and microarray chip are not available for transcript profiling.
In 2001, 47 ESTs up- or downregulated by water stress were first identified in chickpea . cDNA libraries from a drought-responsive genotype in chickpea were constructed and differentially expressed ESTs were identified using in silico approach [9, 34]. SSH libraries have been constructed from chickpea seedling after dehydration stress [27, 35] and between root tissue of 2 chickpea cultivars . Transcriptome analysis by using SuperSAGE and high-throughput 454 sequencing has generated 17,493 unique 26-bp tags (SAGE UniTags) from roots of the drought-tolerant chickpea variety ICC 588 . However, absence of a reference sequence for chickpea and the short read length of sequences (26-bp) limit the utility of this approach.
A total 5494 high-quality sequences (average length 505 bp) were generated after removing short and low-quality sequences. A total of 3062 unigenes (638 contigs and 2424 singletons) were derived from cluster assembly and sequence alignment; each contig had 2-113 ESTs with an average length of the 527 bp. The majority of contigs (84.9%) contained 5 or fewer ESTs, whereas only 2.97% contigs were made from 20 or more ESTs (Additional file 4), indicating a high degree of normalization and subtraction efficiency. All EST sequences have been deposited in the dbEST division of GenBank (HO062174-HO068058). The unigene (UG) set developed in this study is henceforth referred to as UG-TDS (unigenes responsive to terminal drought stress). CAP3 assembly analysis of our datasets with all chickpea EST sequences (34,587) deposited in NCBI dbESTs identified 1576 unigenes (51.4% of total unigenes) as singlets and are new entries to the chickpea database.
To determine the efficiency of normalization and subtraction of SSH libraries, we compared our ESTs with those generated by using non-normalized cDNA libraries. We have previously reported more than 20,000 chickpea root ESTs in response to drought and salt stress in ICC 4958 by using the same procedure to obtain tissue samples for constructing the libraries . CAP3 assembly and clustering analysis of ESTs identified 126 contigs with 1 EST from our SSH libraries and more than 5 ESTs from non-normalized libraries. Some ESTs such as HO063066 (pathogenesis-related protein), HO063205 (plasma membrane intrinsic protein), and HO067852 (Type 1 metallothionein), had single representations in SSH libraries, whereas more than 60 clones were present in non-normalized cDNA libraries. These results support the utility and efficacy of our SSH approach to reduce the redundancy and identify specific transcripts with small-scale sequencing. Dataset analysis with all chickpea EST sequences (34,587) deposited in NCBI dbESTs identified 1576 new unigenes (51.4% of the total unigenes).
Nucleotide-level diversity analysis
Functional characterization of the chickpea unigene dataset
Pathway classification of transcripts
Of the 1808 annotated sequences, 656 were annotated with 812 Enzyme Commission (EC) codes and mapped to 108 different KEGG pathways. Of the 108 pathways contained within the metabolism category (metabolic pathways), 46 were represented by 43.44% of the 656 unigenes. KEGG metabolic pathways well represented by unigenes were biosynthesis of plant hormones (44 enzymes), biosynthesis of phenylpropanoids (29 enzymes) and terpenoids and steroids (24 enzymes), biosynthesis of alkaloids derived from histidine and purine (25 enzymes) and from the shikimate pathway (24 enzymes), starch and sucrose metabolism (24 enzymes), and arginine and proline metabolism (10 enzymes). Several hormone pathways, such as of abscisic acid, ethylene, salicylic acid, and jasmonic acid, are involved in one or more environmental stresses, including drought stress and other abiotic stresses processes [38–42]. A representative KEGG map for biosynthesis of plant hormones is given in Additional file 8.
Gene ontology (GO) enrichment analysis
Differential expression analysis of unigenes under drought stress
In our study, different LEA groups of genes were found in UG-TDS: 2 unigenes encoding HVA protein (HO065000, unigene_Contig11), 5 encoding LEA proteins (HO063258, HO065296, HO0065083, UG-TDS_Contig311 and UG-TDS_Contig524), 6 encoding dehydrin (UG-TDS_Contig232, Contig320, Contig622, UG-TDS HO064933, UG-TDS HO065247 and UG-TDS HO066163), and 1 encoding ERD proteins (HO065032). Among these LEA group 2 members [LEA (HO063258) and dehydrin (HO065247)] were found highly up regulated in drought stress (Figure 8). Earlier studies in chickpea have also reported the induction of LEA proteins under drought stress [27, 34]. The expression profile of LEA genes under stress supports the role of LEA proteins as protective molecules that enable cells to survive protoplasmic water depletion . Studies on overexpression of LEA genes also support the protective role of LEA proteins by improving the stress tolerance of transgenic plants. Expression of the barley gene HVA1 in wheat and rice increases drought tolerance , and overexpression of wheat LEA genes PMA80 and PMA1959 increases dehydration tolerance in transgenic rice .
Different members of aquaporins subfamilies were found in UG-TDS: which includes, 6 unigenes encoding plasma membrane intrinsic protein (UG-TDS HO062890, HO064502, HO064741, HO064425, HO064603 and HO064612), 5 unigenes encoding tonoplast intrinsic protein (UG-TDS HO064719, HO064351, UG-TDS_Contig278, UG-TDS_Contig19 and UG-TDS_Contig156) 2 unigenes encoding NOD26-like intrinsic protein (UG-TDS HO066903 and HO062732). The maximum numbers of the unigenes encoding aquaporin were found in root libraries and downregulation of one of the member (HO062890) under drought stress was conformed in northern blot analysis (Figure 8). This is similar to downregulation of transcripts and reduction in protein levels of most the Arabidopsis aquaporin genes under drought condition , which may be an adaptive strategy for plants to minimize water flow through cell membranes and uphold leaf turgor to minimize water loss. In tobacco, NtPIP1.1 and NtPIP2.1 expression is downregulated to reduce osmotic hydraulic conductance in the roots under drought stress , supporting the role of aquaporins in drought stress maintenance.
Eleven chickpea unigenes from UG-TDS were classified as members of the AP2/ERF superfamily: 10 under the ERF family and 1 under the RAV family. Three members of this family (ERF1, ERF-2, and RAV) were analyzed by Northern blot under drought stress conditions. ERF1 was downregulated whereas ERF2 was upregulated under stress conditions. Biosynthesis of ethylene and regulation of its activation pathway are important to mediate plant developmental processes and stress responses in plants [57, 58]. The AP2/ERF family of transcription factors, especially the CBF/DREB and ERF subfamily, has been extensively studied in response to drought stress . CAP2, a member of the chickpea AP2 family, is responsive to various abiotic stress and its overexpression in tobacco increases the tolerance to dehydration and salt stress . Northern blot analysis showed that UG-TDS HO066286 coding for RAV (related to ABI3/VP1) transcription factor was downregulated under drought stress (Figure 8). Arabidopsis RAV1 is a brassinosteroid (BR) down-regulated gene. High level of BR is accompanied by a very low level of RAV1 transcripts and vice versa . The involvement of BR pathway in the enhancement of tolerance to chilling, thermo, salt, mild drought injury and pathogen attack has been confirmed in several studies [62, 63]. Therefore, the down regulation of RAV during terminal drought stress in our study may indicate the involvement of the BR pathway.
In chickpea, 3 members of the NAC gene family (CarNAC1, CarNAC3 and CarNAC5) are strongly induced by drought, salt, cold, and wounding . We have identified 8 new members of this TF family in UG-TDS with one NAC gene (HO067315) that have increased expression under drought stress validated by northern blot result (Figure 8). Expression profiling and overexpression analysis of NAC genes in several plants supports their involvement in stress tolerance [65–67, 11] K. Nakashima, L.P. Tran, D.V. Nguyen, M. Fujita, K. Maruyama, D. Todaka, Y. Ito, N. Hayashi, K. Shinozaki and K. Yamaguchi-Shinozaki, Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice, Plant J. 51 (2007), pp. 617-630. Full Text via CrossRef | View Record in Scopus | Cited By in Scopus (50).
The HDZip gene (HO062575) was among the up regulated transcription factor, as reflected from northern blot results. Two members of this gene family (HO062575 and UG-TDS_Contig226) have been identified form UG-TDS. The functional information available on plant HDZIP genes suggest that at list some of these genes are involved in response to different environmental conditions . Overexpression of sunflower HD-Zip gene Habt-4 confers drought tolerance in Arabidopsis , this is suggestive of important role of HD-Zip protein in regulation of expression of genes involved in drought tolerance.
Hypothetical proteins are genes of unknown functions predicted from the Arabidopsis or rice genome sequence. Two such genes HP-1 and HP-2 were significantly induced in WS plants. Several hypothetical genes have now been characterized by advanced bioinformatics tools by identifying similarity of conserved function domains. For example, in Arabidopsis, the family of BAG proteins initially annotated as hypothetical proteins are now annotated as bag gene family members, their function as regulators of apoptosis-like processes has also been characterized.
Functional characterization of such unknown hypothetical proteins can shed light on the mechanism of drought adaptation in chickpea.
We found transcript levels of the chlorophyll a/b-binding protein to be downregulated during stress. Most of the strongly downregulated transcripts were related to photosynthesis, photorespiration, and metabolism of amino acids and carbohydrates. In a dehydration shock treatment, the transcript level of chlorophyll a/b-binding protein remained unchanged ; similarly, in barley, chlorophyll a/b-binding protein transcript (NP_917525) levels do not change under dehydration shock treatment but are downregulated by drought stress treatment , indicating differential response of genes under dehydration and drought stress.
Comparative transcript profiles of ICC 4958 and ICC 1882 under drought stress
Genes in clusters I and II were associated with metabolic process [e.g., ethylene biosynthesis (HO062211, HO062180), flavonoid synthesis (HO062384), and amino acid biosyntheses (HO062526, HO062310 and HO062183)] and also these genes shown to be involved in drought response in several other plants [71, 72]. Upregulation of genes involved in ion binding and transport activities [e.g., ATP-binding proteins (HO063146), lipid transfer proteins (HO062394, HO062798), UDP-galactose transporters (HO062219), metal ion binding (HO062399), sulfate transporters (HO063202), tonoplast intrinsic proteins (HO062783), were also upregulated in ICC 4958. In an earlier study, we reported by in silico differential expression analysis the upregulation of the tonoplast intrinsic protein in the roots of ICC 4958, which mediates the regulation of root hydraulic conductivity in response to environmental stimuli . Several stress-related genes [e.g., pathogenesis-related proteins (HO062911, HO062939) and peroxidase (HO062698), chaperone binding (HO062569) and small heat shock protein (HO062866)] upregulated in ICC 4958 and have been shown to be induced by wounding, salt, and cold stress in other plant species [73, 74] indicates multiple stress induction of these genes. Similarly known stress-responsive transcription factors and regulators such as the AP2/ERF domain-containing transcription factor (HO062802), MYB transcription factor (HO062363), DNA repair and transcription factor XPB1 (HO062308), and transcription regulators (HO062392) were also upregulated in ICC 4958. A similar differential induction of these genes or gene categories in drought-tolerant genotypes in response to drought stress during the reproductive stage has been reported in barley .
Cluster III contained unigenes that were upregulated in ICC 1882 but not ICC 4958. One upregulated unigene encoded fructose-bisphosphate aldolase (HO063129), whose downregulation could inhibit gluconeogenesis for conserving energy in drought-stressed plants (41).
We report the sequencing, assembly, and annotation of 5494 high-quality drought-responsive EST sequences from chickpea. This dataset was generated from SSH libraries constructed using drought-tolerant and -susceptible chickpea genotypes and bulks of their progenies exhibiting HRB and LRB phenotypes. SSH libraries allowed cloning genes that are specifically up- and downregulated from the roots and shoots of chickpea in response to terminal drought. Moreover, we identified more than 1500 novel unigenes in chickpea that are associated with terminal drought stress. Besides several transcripts coding for known stress-related proteins, several novel genes with unknown functions that may have a potential role in drought tolerance in chickpea were also identified. This study also provides a comparative overview of genotype-specific expression patterns of more than 830 unigenes in root tissues of chickpea in response to drought. The up- and downregulation of some unigenes was confirmed by real-time qPCR. The EST dataset and the information about transcription of several genes can be useful for the research community and help identify potential candidate genes for drought tolerance in chickpea. Our study can also serve as an important resource for developing functional markers, full-length gene isolation, TILLING, drought-responsive promoter isolation, and in drought functional genomic studies involving overexpression, e-QTL, and manipulation of drought tolerance in chickpea.
We gratefully acknowledge the "National Fund for Basic, Strategic and Frontier Application Research in Agriculture (NFBSFARA), Indian Council of Agricultural Research (ICAR), New Delhi, India, for financial support. Authors are thankful to Dr. Pooran Gaur of ICRISAT for sharing the parental genotypes and mapping population used in the study. Thanks are also due to P Lekha and Pavana Hiremath for their help in planning of some experiments and discussions. We are grateful to Dr Vani Shankar for editing the manuscript.
- Kumar J, Abbo S: Genetics of flowering time in chickpea and its bearing on productivity in semiarid environments. Adv Agron. 2001, 72: 107-138.View ArticleGoogle Scholar
- Singh KB, Malhotra RS, Halila MH, Knights EJ, Verma MM: Current status and future strategy in breeding chickpea for resistance to biotic and abiotic stresses. Expanding the production and use of coolseason food legumes. Edited by: Muehlbauer FJ, Kaiser WJ. 1994, Dordrecht, The Netherlands: Kluwer Academic Publishers, 572-591.View ArticleGoogle Scholar
- Croser JS, Clarke HJ, Siddique KHM, Khan TN: Low-temperature stress: implications for chickpea (Cicer arietinum L.) improvement. Crit Rev Plant Sci. 2003, 22: 185-219. 10.1080/713610855.View ArticleGoogle Scholar
- Ahmad F, Gaur PM, Croser J: Chickpea (Cicer arietinum L.). Genetic Resources, Chromosome Engineering and Crop Improvement - Grain Legumes. Edited by: Singh R and Jauhar P. 2005, USA: CRC Press, 1: 185-214.Google Scholar
- Leport L, Turner NC, French RJ, Barr MD, Duda R, Davies SL, Tennant D, Siddique KHM: Physiological responses of chickpea genotypes to terminal drought in a Mediterranean-type environment. European Journal of Agronomy. 1999, 11: 279-291. 10.1016/S1161-0301(99)00039-8.View ArticleGoogle Scholar
- Siddique KHM, Brinsmead RB, Knight R, Knights EJ, Paul JG, Rose IA: Adaptation of chickpea (Cicer arietinum L.) and faba bean (Vicia faba L.) to Australia. Linking research and marketing opportunities for pulses in the 21st century. Edited by: Knight R. 2000, Dordrecht, The Netherlands: Kluwer, 289-303.View ArticleGoogle Scholar
- Molina C, Rotter B, Horres R, Udupa SM, Besser B, Bellarmino L, Baum M, Matsumura H, Terauchi R, Kahl G: SuperSAGE: the drought stress-responsive transcriptome of chickpea roots. BMC Genomics. 2008, 9: 553-10.1186/1471-2164-9-553.PubMedPubMed CentralView ArticleGoogle Scholar
- Mantri NL, Ford R, Coram TE, Pang EC: Transcriptional profiling of chickpea genes differentially regulated in response to high-salinity, cold and drought. BMC Genomics. 2007, 8: 303-10.1186/1471-2164-8-303.PubMedPubMed CentralView ArticleGoogle Scholar
- Varshney RK, Hiremath PJ, Lekha P, Kashiwagi J, Balaji J, Deokar AA, Vadez V, Xiao Y, Srinivasan R, Gaur PM, Siddique KH, Town CD, Hoisington DA: A comprehensive resource of drought- and salinityresponsive ESTs for gene discovery and marker development in chickpea (Cicer arietinum L.). BMC Genomics. 2009, 15 (10): 523-View ArticleGoogle Scholar
- Ray J, Sinclair T: : The effect of pot size on growth and transpiration of maize and soybean during water deficit stress. J Exp Bot. 1998, 49: 1381-1386. 10.1093/jexbot/49.325.1381.View ArticleGoogle Scholar
- Kashiwagi J, Krishnamurthy L, Gaur PM, Chandra S, Upadhyaya HD: Estimation of gene effect of the drought avoidance root characteristics in chickpea (C.arietinum L.). Field Crop Research. 2007, 105: 64-69.View ArticleGoogle Scholar
- EST trimming. [http://pgrc.ipk-gatersleben.de/misa/download/est_trimmer.pl]
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9 (9): 868-877. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucl Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Okuda S, Yamada T, Hamajima M, Itoh M, Katayama T, Bork P, Goto S, Kanehisa M: KEGG atlas mapping for global analysis of metabolic pathways. Nucleic Acids Res. 2008, 36: W423-426. 10.1093/nar/gkn282.PubMedPubMed CentralView ArticleGoogle Scholar
- Blüthgen N, Brand K, Cajavec B, Swat M, Herzel H, Beule D: Biological profiling of gene groups utilizing Gene Ontology. Genome Inform. 2005, 16: 106-15.PubMedGoogle Scholar
- Sambrook , Russell : Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor Laboratory Press, 3Google Scholar
- Garg R, Sahoo A, Tyagi AK, Jain M: Validation of internal control genes for quantitative gene expression studies in chickpea (Cicer arietinum L.). Biochem Biophys Res Commun. 2010, 396 (2): 283-8. 10.1016/j.bbrc.2010.04.079.PubMedView ArticleGoogle Scholar
- Livaka KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C (T)) method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticleGoogle Scholar
- Kholová J, Hash CT, Lava Kumar P, Yadav RS, Kočová M, Vadez V: Terminal drought-tolerant pearl millet [Pennisetum glaucum (L.) R. Br.] have high leaf ABA and limit transpiration at high vapour pressure deficit. J Exp Bot. 2010, 61 (5): 1431-1440. 10.1093/jxb/erq013.PubMedPubMed CentralView ArticleGoogle Scholar
- Bhatnagar-Mathur P, Vadez V, Devi MJ, Lavanya M, Vani G, Sharma KK: Genetic engineering of chickpea (Cicer arietinum L.) with the P5CSF129A gene for osmoregulation with implications on drought tolerance. Mol Breeding. 2009, 23: 591-606. 10.1007/s11032-009-9258-y.View ArticleGoogle Scholar
- Sinclair TR, Holbrook NM, Zwieniecki MA: Daily transpiration rates of woody species on drying soil. Tree Physiology. 2005, 25: 1469-1472.PubMedView ArticleGoogle Scholar
- Vadez V, Sinclair TR, Serraj R, Purcell LC: Manganese application alleviates the water deficit-induced decline of N2 fixation. Plant, Cell & Environment. 2000, 23: 497-505. 10.1046/j.1365-3040.2000.00562.x.View ArticleGoogle Scholar
- Fang X, Turner NC, Yan G, Li F, Siddique KHM: Flower numbers, pod production, pollen viability, and pistil function are reduced and flower and pod abortion increased in chickpea (Cicer arietinum L.) under terminal drought. J Exp Bot. 2010, 61 (2): 335-345. 10.1093/jxb/erp307.PubMedPubMed CentralView ArticleGoogle Scholar
- Rabbani MA, Maruyama K, Abe H, Khan MA, Katsura K, Ito Y, Yoshiwara K, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Monitoring expression profiles of rice genes under cold, drought, and high-salinity stresses and abscisic acid application using cDNA microarray and RNA gel-blot analyses. Plant Physiol. 2003, 133: 1755-1767. 10.1104/pp.103.025742.PubMedPubMed CentralView ArticleGoogle Scholar
- Boominathan P, Shukla R, Kumar A, Manna D, Negi D, Verma PK, Chattopadhyay D: Long term transcript accumulation during the development of dehydration adaptation in Cicer arietinum. Plant Physiol. 2004, 135: 1608-1620. 10.1104/pp.104.043141.PubMedPubMed CentralView ArticleGoogle Scholar
- Rensink WA, Buell CR: Microarray expression profiling resources for plant genomics. Trends in Plant Science. 2005, 10: 603-609. 10.1016/j.tplants.2005.10.003.PubMedView ArticleGoogle Scholar
- Turner NC, Wright GC, Siddique KHM: Adaptation of grain legumes (pulses) to water limited environments. Advances in Agronomy. 2001, 71: 193-231.View ArticleGoogle Scholar
- Kashiwagi J, Krishnamurthy L, Upadhyaya HD, Krishna H, Chandra S, Vadez V, Serraj R: Genetic variability of drought-avoidance root traits in the mini-core germplasm collection of chickpea (Cicer arietinum L.). Euphytica. 2005, 146: 213-222.View ArticleGoogle Scholar
- Périn C, Rebouillat J, Brasileiro ACM, Diévart A, Gantet P, Breitler JC, Johnson AAT, Courtois B, Ahmadi N, de Raissac M: Novel insights into the genomics of rice root adaptive development. Rice Genetics V. Edited by: Brar DS, Mackill DJ, Hardy B. 2007, Singapore: World Scientific, 117-141.View ArticleGoogle Scholar
- Cao W, Epstein C, Liu H, DeLoughery C, Ge N, Lin J, Diao R, Cao H, Long F, Zhang X: Comparing gene discovery from Affymetrix GeneChip microarrays and Clontech PCR-select cDNA subtraction: a case study. BMC Genomics. 2004, 5: 26-10.1186/1471-2164-5-26.PubMedPubMed CentralView ArticleGoogle Scholar
- Romo S, Labrador E, Dopico B: Water stress-regulated gene expression in Cicer arietinum seedlings and plants. Plant Physiol Biochem. 2001, 39: 1017-1026. 10.1016/S0981-9428(01)01318-3.View ArticleGoogle Scholar
- Gao WR, Wang XSH, Liu P, Chen Ch, Li JG, Zhang JS, Ma H: Comparative analysis of ESTs in response to drought stress in chickpea (Cicer arietinum L.). Biochic and Biophysic Res Commun. 2008, 376: 578-583. 10.1016/j.bbrc.2008.09.030.View ArticleGoogle Scholar
- Jain D, Chattopadhyay D: Analysis of gene expression in response to water deficit of chickpea (Cicer arietinum L.) varieties differing in drought tolerance. BMC Plant Biology. 2010, 10: 24-10.1186/1471-2229-10-24.PubMedPubMed CentralView ArticleGoogle Scholar
- Buhariwalla HK, Jayashree B, Eshwar K, Crouch JH: Development of ESTs from chickpea roots and their use in diversity analysis of the Cicer genus. BMC Plant Biology. 2005, 5: 16-10.1186/1471-2229-5-16.PubMedPubMed CentralView ArticleGoogle Scholar
- Doyle JJ, Luckow MA: The rest of the iceberg. Legume diversity and evolution in a phylogenetic context. Plant Physiol. 2003, 131: 900-910. 10.1104/pp.102.018150.PubMedPubMed CentralView ArticleGoogle Scholar
- Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y: KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008, 36: D480-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka Y, Sano T, Tamaoki M, Nakajima N, Kondo N, Hasezawa S: Ethylene inhibits abscisic acid-induced stomatal closure in Arabidopsis. Plant Physiol. 2005, 138: 2337-2343. 10.1104/pp.105.063503.PubMedPubMed CentralView ArticleGoogle Scholar
- Aprile A, Mastrangelo AM, De Leonardis AM, Galiba G, Roncaglia E, Ferrari F, Bellis LD, Turchi L, Giuliano G, Cattivelli L: Transcriptional profiling in response to terminal drought stress reveals differential responses along the wheat genome. BMC Genomics. 2009, 10: 279-10.1186/1471-2164-10-279.PubMedPubMed CentralView ArticleGoogle Scholar
- Gong P, Zhang J, Li H, Yang C, Zhang C, Zhang X, Khurram Z, Zhang Y, Wang T, Fei Z: Transcriptional profiles of drought-responsive genes in modulating transcription signal transduction, and biochemical pathways in tomato. J Exp Bot. 2010, 61: 3563-3575. 10.1093/jxb/erq167.PubMedPubMed CentralView ArticleGoogle Scholar
- Zheng J, Fu J, Gou M, Huai J, Liu Y, Jian M, Huang Q, Guo X, Dong Z, Wang H, Wang G: Genome-wide transcriptome analysis of two maize inbred lines under drought stress. Plant Mol Biol. 2010, 72 (4-5): 407-421. 10.1007/s11103-009-9579-6.PubMedView ArticleGoogle Scholar
- Abebe T, Melmaiee K, Berg V, Wise RP: Drought response in the spikes of barley: gene expression in the lemma, palea, awn, and seed. Functional & Integrative Genomics. 2009, 10 (2): 191-205.View ArticleGoogle Scholar
- Miyama M, Hanagata N: Microarray analysis of 7029 gene expression patterns in burma mangrove under high-salinity stress. Plant Sci. 2007, 172: 948-957. 10.1016/j.plantsci.2007.01.004.View ArticleGoogle Scholar
- Vernon DM, Bohnert HJ: A novel methyl trans-ferase induced by osmotic stress in the facultative halophyte Mesembryanthemum crystallinum. EMBO Journal. 1992, 11: 2077-2085.PubMedPubMed CentralGoogle Scholar
- Yoshida KT, Wada T, Koyama H, Mizobuchi-Fukuoka R, Naito S: Temporal and spatial patterns of accumulation of the transcript of myo-inositol-1-phosphate synthase and phytin containing particles during seed development in rice. Plant Physiol. 1999, 119: 65-72. 10.1104/pp.119.1.65.PubMedPubMed CentralView ArticleGoogle Scholar
- Abreu EFM, Aragão FJL: Isolation and characterization of a myo-inositol-1-phosphate synthase gene from yellow passion fruit (Passiflora edulis f. flavicarpa) expressed during seed development and environmental stress. Annals of Botany. 2007, 99: 285-292.PubMedPubMed CentralView ArticleGoogle Scholar
- Kaur H, Shukla RK, Yadav G, Chattopadhyay D, Majee M: Two divergent genes encoding L-myo-inositol 1-phosphate synthase1 (CaMIPS1) and 2 (CaMIPS2) are differentially expressed in chickpea. Plant, Cell & Environment. 2008, 31 (11): 1701-1716. 10.1111/j.1365-3040.2008.01877.x.View ArticleGoogle Scholar
- Larson SR, Raboy V: Linkage mapping of maize and barley myo-inositol 1-phosphate synthase DNA sequences: correspondence with a low phytic acid mutation. Theoretical and Applied Genetics. 1999, 99: 27-36. 10.1007/s001220051205.View ArticleGoogle Scholar
- Johnson MD, Sussex IM: 1-L-myo-inositol 1-phosphate synthase from Arabidopsis thaliana. Plant Physiol. 1995, 107: 613-619.PubMedPubMed CentralGoogle Scholar
- Delaunay AJ, Verma DPS: Proline accumulation and osmoregulation in plants. Plant J. 1993, 4: 215-223. 10.1046/j.1365-313X.1993.04020215.x.View ArticleGoogle Scholar
- Ingram J, Bartels D: The molecular basis of dehydration tolerance in plants. Annu Rev Plant Physiol Plant Mol Biol. 1996, 47: 377-403. 10.1146/annurev.arplant.47.1.377.PubMedView ArticleGoogle Scholar
- Sivamani E, Bahieldin A, Wraith JM, Al-Niemi T, Dyer WE, Ho T-HD, Qu R: Improved biomass productivity and water use efficiency under water deficit conditions in transgenic wheat constitutively expressing the barley HVA1 gene. Plant Sci. 2000, 155: 1-9. 10.1016/S0168-9452(99)00247-2.PubMedView ArticleGoogle Scholar
- Cheng Z, Targolli J, Huang X, Wu R: Wheat LEA genes, PMA80 and PMA enhance dehydration tolerance of transgenic rice (Oryza sativa L.). Mol Breeding. 1995, 10: 71-82.View ArticleGoogle Scholar
- Alexandersson E, Fraysse L, Sjövall-Larsen S, Gustavsson S, Fellert M, Karlsson M, Johanson U, Kjellbom P: Whole gene family expression and drought stress regulation of aquaporins. Plant Mol Biol. 2005, 59: 469-484. 10.1007/s11103-005-0352-1.PubMedView ArticleGoogle Scholar
- Mahdieh M, Mostajeran A, Horie T, Katsuhara M: Drought Stress Alters Water Relations and Expression of PIP-Type Aquaporin Genes in Nicotiana tabacum Plants. Plant and Cell Physiol. 2008, 49 (5): 801-813. 10.1093/pcp/pcn054.View ArticleGoogle Scholar
- Zhao XC, Schaller GE: Effect of salt and osmotic stress upon expression of the ethylene receptor ETR1 in Arabidopsis thaliana. FEBS Lett. 2004, 562: 189-192. 10.1016/S0014-5793(04)00238-8.PubMedView ArticleGoogle Scholar
- Cao Y, Song F, Goodman RM, Zheng Z: Molecular characterization of four rice genes encoding ethylene-responsive transcriptional factors and their expressions in response to biotic and abiotic stress. Journal of Plant Physiology. 2006, 163: 1167-1178. 10.1016/j.jplph.2005.11.004.PubMedView ArticleGoogle Scholar
- Agarwal PK, Agarwal P, Reddy MK, Sopory SK: Role of CBF/DREB transcription factors in abiotic and biotic stress tolerance in plants. Plant Cell Reports. 2006, 25: 1263-1274. 10.1007/s00299-006-0204-8.PubMedView ArticleGoogle Scholar
- Shukla RK, Raha S, Tripathi V, Chattopadhyay D: Expression of CAP2, an APETALA2-Family Transcription Factor from Chickpea, Enhances Growth and Tolerance to Dehydration and Salt Stress in Transgenic Tobacco. Plant Physiol. 2006, 142: 113-123. 10.1104/pp.106.081752.PubMedPubMed CentralView ArticleGoogle Scholar
- Hu YX, Wang YH, Liu XF, Li JY: Arabidopsis RAV1 is down-regulated by brassinosteroid and may act as a negative regulator during plant development. Cell Research. 2004, 14: 8-15. 10.1038/sj.cr.7290197.PubMedView ArticleGoogle Scholar
- Krishna P: Brassinosteroid-mediated stress responses. J Plant Growth Regul. 2003, 22: 289-297. 10.1007/s00344-003-0058-z.PubMedView ArticleGoogle Scholar
- Divi UK, Krishna P: Brassinosteroids confer stress tolerance. Plant stress biology: genomics goes systems biology. Edited by: Hirt H. 2009, Weinheim: Wiley-VCH, 119-135.View ArticleGoogle Scholar
- Peng H, Yu X, Cheng H, Shi Q, Zhang H, Li J, Ma H: Cloning and characterization of a novel NAC family gene CarNAC1 from chickpea (Cicer arietinum L.). Mol Biotechnol. 2010, 44 (1): 30-40. 10.1007/s12033-009-9202-8.PubMedView ArticleGoogle Scholar
- Hegedus D, Yu M, Baldwin D, Gruber M, Sharpe A, Parkin I, Whitwill S, Lydiate D: Molecular characterization of Brassica napus NAC domain transcriptional activators induced in response to biotic and abiotic stress. Plant Mol Biol. 2003, 53: 383-397.PubMedView ArticleGoogle Scholar
- Nakashima K, Tran LS, Van Nguyen D, Fujita M, Maruyama K, Todaka D, Ito Y, Hayashi N, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of a NAC-type transcription factor OsNAC6 involved in abiotic and biotic stress-responsive gene expression in rice. Plant J. 2007, 51: 617-630. 10.1111/j.1365-313X.2007.03168.x.PubMedView ArticleGoogle Scholar
- Guo YF, Gan SS: AtNAP, a NAC family transcription factor, has an important role in leaf senescence. Plant J. 2006, 47: 601-612.View ArticleGoogle Scholar
- Henriksson E, Olsson AS, Johannesson H, Johansson H, Hanson J, Engstrom P, Soderman E: Homeodomain leucine zipper class I genes in Arabidopsis: expression patterns and phylogenetic relationships. Plant Physiol. 2005, 139: 509-518. 10.1104/pp.105.063461.PubMedPubMed CentralView ArticleGoogle Scholar
- Dezer CA, Gago GM, Gonzalez DH, Chan RL: Hahb-4, a sunflower homeobox-leucine zipper gene, is a developmental regulator and confers drought tolerance to Arabidopsis thaliana plants. Transgenic research. 2005, 14: 429-440. 10.1007/s11248-005-5076-0.View ArticleGoogle Scholar
- Talame V, Ozturk NZ, Bohnert HJ, Tuberosa R: Barley transcript profiles under dehydration shock and drought stress treatments: a comparative analysis. J Exp Bot. 2007, 58: 229-240.PubMedView ArticleGoogle Scholar
- Akashi K, Miyakel C, Yokota A: A novel compatible solute in drought-tolerant wild watermelon leaves, is an efficient hydroxyl radical scavenger. FEBS Letters. 2001, 508: 438-442. 10.1016/S0014-5793(01)03123-4.PubMedView ArticleGoogle Scholar
- Capell T, Bassie L, Christou P: Modulation of the polyamine biosynthetic pathway in transgenic rice confers tolerance to drought stress. Proceedings of National Academy of Sciences, USA. 2004, 101: 9909-9914. 10.1073/pnas.0306974101.View ArticleGoogle Scholar
- Sung DY, Vierling E, Guy CL: Comprehensive expression profile analysis of the Arabidopsis Hsp70 gene family. Plant Physiol. 2001, 126: 789-800. 10.1104/pp.126.2.789.PubMedPubMed CentralView ArticleGoogle Scholar
- Radauer C, Breiteneder H: Evolutionary biology of food plant allergens. The Journal of Allergy and Clinical Immunology. 2007, 120: 518-525. 10.1016/j.jaci.2007.07.024.PubMedView ArticleGoogle Scholar
- Guo P, Baum M, Grando S, Ceccarelli S, Bai G, Li R, Von Korff M, Varshney RK, Graner A, Valkoun J: Differentially expressed genes between drought tolerant and drought sensitive barley genotypes in responsive to drought stress during the reproductive stage. J Exp Bot. 2009, 60 (12): 3531-3544. 10.1093/jxb/erp194.PubMedPubMed CentralView ArticleGoogle Scholar