Expression and Promoter Analysis of Six Heat Stress-Inducible Genes in Rice

Abstract

During the long evolutionary process, plant gradually formed a series of strategies and mechanisms to cope with stress environment such as drought, heat, cold, and high salinity. Six highly heat responsive genes were identified in rice by microarray data analysis. The qRT-PCR analysis confirmed that the expression of these six genes were highly heat inducible and moderately responded to salt stress, polyethylene glycol, and abscisic acid treatment, but little affected by cold treatment. Promoters of the three highly heat-inducible genes (OsHsfB2cp, PM19p, and Hsp90p) were used to drive GUS gene expression in rice. The results of the GUS gene expression, histochemical staining, and GUS activities in panicles and flag leaves of the transgenic rice plants confirmed high heat-induced GUS activities and moderate drought-induced activities. The three promoters exhibited similar high activity lever in rice leaf under heat, but OsHsfB2cp and PM19p showed much higher activities in panicles under heat stress. Our work confirmed that the OsHsfB2c and PM19 promoters were highly heat inducible and further characterization and reconstruction of cis-elements in their promoters could lead to the development of highly effective heat-inducible promoters for plant genetic engineering.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Rice is the most important staple food in Asia. More than 90% of the world’s rice is grown and consumed in Asia, where 60% of the world’s population lives. Rice plants are sensitive to various abiotic stresses [1]. As such, rice carries an odd portfolio of tolerances and susceptibilities to abiotic stresses as compared to other crops [2]. However, many rice-growing environments demand still greater tolerance than what has been found in most improved germplasm [3]. Global warming has been likely to induce heat stress more frequently and severely, resulting in serious yield reduction and thus threatening the security of grain production. Abiotic stresses such as high temperature or drought can induce diverse physiological and molecular responses in plants. These include extensive cellular damage and inhibition of photosynthesis, osmotic adjustment, induction of repair systems and chaperones, and changes of gene expression and metabolism [4]. Functional characterization of stress-inducible genes is important not only for understanding the molecular mechanisms of stress tolerance, but also for practical application in improving stress tolerance of crops by gene manipulation. To improve stress tolerance of crop plants, stress tolerant/resistant genes can be overexpressed in transgenic plants, using either constitutive or stress-inducible promoters. Although constitutive expression of stress resistant genes can improve resistance of transgenic plants to abiotic stresses, it may also cause stunted growth and reduction of yield in plants. Promoters that were subject to specific regulations are useful for manipulating foreign gene expression in plant cells, tissues, or organs with desirable patterns and under controlled conditions [5]. Therefore, inducible promoters including chemical inducible promoters, which drive gene expression only when plants are exposed to stresses or bioreactors to adapt the environment and maintain normal growth, are of great importance in biotechnology, molecular breeding, and agriculture [6].

In the present study, we have identified six abiotic stress responsive genes in rice from microarray data and further confirmed the expression pattern of these six genes in different organs and under abiotic stresses in rice by quantitative real-time PCR (qRT-PCR). Promoters of three highly heat-inducible genes (OsHsfB2cp, PM19p, and Hsp90p) were used to drive GUS gene expression, and the expression in transgenic rice panicles and flag leaves under heat and drought treatments was performed. These data and further work to come would be valuable for developing stress-inducible promoters in application.

2. Materials and Methods

2.1. Identification of Candidate Stress Inducible Genes from Microarray Data

In our previous study, the genome-wide gene expression analysis in response to heat shock was performed in rice panicle [7] and flag leaf [8] and the original microarray data were deposited in GEO database under accession numbers GSE38665 and GSE45259. Expression of six heat-responsive genes was derived from the above microarray data. The expression data of the six heat-responsive genes in rice seedlings under abiotic stresses were collected through the GEO accession numbers GSE14275 [9] and GSE6901 [10]. The six candidate genes in this study were named OsHsfA2a (Os03g0745000), Dehydrin (Os11g0453900), DnaJ protein putative (Os06g0195800), OsHsfB2c (Os09g0526600), AWPM-19 (Os05g0381400), and Hsp90 (Os04g0107900).

2.2. Plants Materials, Growth Condition, and Treatments

Rice seeds (Oryza sativa L. ssp. japonica) were imbibed in water in the dark at 28°C for 3 days. Germinated seeds were planted in plastic pots and grown in a climate chamber (Binder, Tuttlingen, Germany) at 25°C and 80% RH (relative humidity) with a 12 h light/12 h dark cycle and irrigated with 1/2 MS (Murashige and Skoog medium) liquid culture. After three weeks, the seedlings were subject to the abiotic stress treatments. For salt, polyethylene glycol (PEG), and abscisic acid (ABA) treatments, seedlings were irrigated with 1/2 MS liquid culture containing 100 mmol L⁻¹ NaCl, 10% PEG, and 100 μmol L⁻¹ ABA, respectively. For cold and heat shock treatments, seedlings were exposed to 4°C and 42°C, respectively, and irrigated with 1/2 MS liquid culture. After stress treatment, samples were collected separately at 0 min, 30 min, and 1, 2, 4, 8, and 12 h. For tissue specific expression analysis, roots were obtained from the three-week-old rice seedlings grown under the same condition as the controls described in the stress treatment section. Stems, leaves, sheaths, and spikes were obtained from field grown rice plants at booting stage. Rice embryos were isolated from 48 h imbibed seeds. Three replications for each experiment were performed in this study. All samples were flash-frozen in liquid nitrogen and stored at −80°C for subsequent RNA extraction.

2.3. RNA Isolation and qRT-PCR Analysis

Total RNA was extracted using TRIZOL reagent (Invitrogen, Carlsbad, CA, USA), according to the manufacturer’s protocol. The RNA samples were quantified by absorbance at 260 nm and by electrophoresis on a 1% agarose gel stained with ethidium bromide. Five micrograms of each RNA sample was incubated with 1 U DNase I (Fermentas, Vilnius, Lithuania) for 30 min at 37°C to remove residual DNA contamination. DNase I-treated RNA samples were then reverse transcribed with Oligo ${(\text{dT})}_{\mathrm{18}}$ primers using ReverTra Ace (Toyobo, Osaka, Japan) for 60 min at 42°C. The cDNA samples were then diluted threefold with buffer (10 mmol L⁻¹ Tris-HCl, pH 8.0, 1 mmol L⁻¹ EDTA).

qRT-PCR was performed on an ABI 7300 real-time PCR system (Applied Biosystems) using SYBR Premix Ex Taq II (Takara, Tokyo, Japan) under the following conditions: 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 60°C for 31 s. A rice actin gene (RAc1, accession number AB047313) expression level was used as a positive internal control (Table 1). Three biological replicates and three technical replicates for each experiment were performed.

Table 1

Primer sequences of the 6 candidate genes for qRT-PCR.

Gene name	Forward primer (5′ to 3′)	Reverse primer (5′ to 3′)
Os03g0745000	GCGACACCGAGAGCTTCTGGATC	TCGTCATCCTCCTCGTCGTTGT
Os11g0453900	TCCAGCTCCAGCTCATCCTCT	GTCGTTGTTGTGGTGTGCTCCT
Os06g0195800	TGCGGAGAGCTTCGAGGAGTT	ACCGCTTGTTCACGCCCTTG
Os09g0526600	GTCCAGCTCCAGCCAAACGAT	CCTTACGCCACTACCGCATTCC
Os05g0381400	CGCTGCTGGTGCTGAATCTGAT	AGGATGGCGAAGACGAGGAAGT
Os04g0107900	ACGCCAAGGAGAGGAAGAGGTC	ATGTTCGCCGTCCAGCCGTA
RAc1	CTTCAACACCCCTGCTATG	TCCATCAGGAAGCTCGTAG

Tissue-specific expression of the six genes were analyzed using $2^{\mathrm{\Delta }{C}_T}$ , where $\mathrm{\Delta }{C}_T=C_T$ $(\text{Target})-{\text{C}}_{\text{T}}\mathrm{\hspace{0.17em}\hspace{0.17em}}(\text{Actin})$ . The relative expression change of six genes in response to abiotic stresses was quantified using the $2^{-\mathrm{\Delta }\mathrm{\Delta }{C}_T}$ method [11]. The data were expressed as mean ± standard error.

2.4. Bioinformatics Analysis of the Three Inducible Promoters

The 1.5 kb promoter sequences upstream of the start codon ATG of the 3 stress-inducible genes were extracted from the GenBank rice genome database, and promoter cis-elements were analyzed using PLANTCARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/) and PLACE (http://www.dna.affrc.go.jp/PLACE/).

2.5. Plasmid Construction for the Three Stress-Inducible Promoters

The promoters of choice were amplified by PCR using promoter-specific primers (Table 2). The PCR was performed at 95°C for 2 min, followed by 36 cycles of 94°C for 30 s, 60°C for 30 s, and 72°C for 1 min, with a final cycle at 72°C for 8 min using 100 ng of genomic DNA. PCR products were detected by 1% agarose gel electrophoresis. The target fragments were recovered by DNA purification kit and ligated to the vector pMD19-T, and the product was transformed into Escherichia coli DH5α. Positive recombinant clones were selected by PCR and sequenced by Nanjing Jin Site Biotechnology Company. To construct the recombinant expression vectors for Agrobacterium tumefaciens, the pCAMBIA1301 vector was first digested with HindIII/BglII to remove the 35 S promoter sequence, then end-filled, and ligated as a circular vector. The target promoter fragments in the pMD19-T vectors were digested by BamHI/PstI and inserted into the BamHI/PstI sites of the modified pCAMBIA1301 vector to get the specific promoter-GUS recombinant plasmids

Table 2

Primers for cloning the 3 candidate promoters.

Promoters	Primer sequence (5′ to 3′)	Product size (bp)
OsHsfB2cp	Forward GGGTACCTATCAAAATCAATCTCAAAATTAAGCTTG	1220
	Reverse  ACTGCAGGGCGAGGACCGACCGTAGTAG
PM19p	Forward TAAGCTTGGCGCGTATCGACCATGGTT	1218
	Reverse  CGGATCCGCAACACAATTAGCACTGGTTTAGCT
Hsp90p	Forward CGGTACCTTTCTGATTGTATTTCTACGTGGACACT	1437
	Reverse  ACTGCAGAATGTTGCCTGCTGCTTGAATG

2.6. Plant Transformation

The promoter-GUS recombinant plasmids were mobilized into Agrobacterium tumefaciens LBA4404 using the freeze-thaw method [12]. Embryogeniccalli of Nipponbare were grown and transformed with Agrobacterium [13]. Transgenic calli were selected on selection medium consisting of MS salts and vitamins, 30 g/L sucrose, 2 mg/L 2,4-D, 3 g/L phytagel (P8169, Sigma-Aldrich) supplemented with 50 mg/L hygromycin B for 30 days with one transfer to fresh medium after 15 days. Transgenic plants were regenerated under hygromycin selection on regeneration medium (MS, 30 g/L sucrose, 2 mg/L kinetin, 0.5 mg/L NAA, 3 g/L phytagel). The regenerated plants were allowed to set roots for 10 days in rooting medium (one-half strength MS, 2 g/L phytagel) and followed by growing in water for 10 days before being transferred to pots in the greenhouse and finally grown to maturity. All plants were fertile with normal phenotype. The ${\text{T}}_{\mathrm{0}}$ transgenic rice plants grown in the greenhouse were confirmed for the presence of gusA. Subsequently, ${\text{T}}_{\mathrm{1}}$ and ${\text{T}}_{\mathrm{2}}$ generation of transgenic plants were raised. About 30 ${\text{T}}_{\mathrm{0}}$ transgenic plants for each of the promoter constructs were derived and grown up to ${\text{T}}_{\mathrm{2}}$ generation. At least three ${\text{T}}_{\mathrm{2}}$ transgenic lines for each promoter were considered for further characterization.

2.7. Heat and Drought Treatments on Transgenic Rice Lines

To evaluate the GUS gene expression and GUS activity in the promoter: GUS transgenic rice, heat, and drought stresses were performed at the reproductive stage. For heat stress treatment, whole plants were exposed to 40°C for 8 h. For drought treatment, the plants were withheld from water supply for 3 days until the leaves became rolled. The control plants were kept at 30°C and fully watered. The untreated plants were grown at normal conditions. The flag leaves and panicles were sampled and flash-frozen in liquid nitrogen and stored at −80°C. Three replications are performed for each experiment in this study.

2.8. GUS Gene Expression Analysis by qRT-PCR

The qRT-PCR was performed with gusA specific primers under the following conditions: 95°C for 30 s, followed by 40 cycles at 95°C for 5 s, and 60°C for 31 s using GUS5A (5′CGACGCTCACACCGATACCATC3′) and GUS3 (5′TCTCCTGCCAGGCCAGAAGTTC3′). The quantitative variation between different samples was evaluated. The amplification of rice actin gene (RAc1) using primers AF (5′CTTCAACACCCCTGCTATG3′) and AR (5′TCCATCAGGAAGCTCGTAG3′) was performed as internal control to normalize all data. The mean values for the expression levels of the target genes were calculated from three independent experiments. The relative expression level was calculated as $2^{-\mathrm{\Delta }\mathrm{\Delta }{C}_T}$ .

2.9. GUS Histochemical Staining

GUS histochemical staining was performed as described by Jefferson et al. [14]. The GUS reaction mixture consisted of the solution containing sodium phosphate (Na₄PO₄, pH 7.0), Triton-X-100, potassium ferricyanide (K₃FeCN₆), potassium ferrocyanide (K₄FeCN₆), and 5-bromo-4-chloro-3-indolyl-β-D-glucuronide (X-gluc). GUS histochemical staining on the spikelets and flag leaves from the three transgenic rice lines was conducted under heat and drought stresses. The samples were incubated in the solution at 37°C for 24 h. The reaction was stopped by adding 70% (v/v) ethanol and the pigments and chlorophylls were removed by repeated ethanol treatment.

2.10. Quantitative GUS Enzyme Activity Assay

Fluorometric assay for GUS enzyme activity was carried out according to the method of Jefferson et al. [14]. The extracted proteins were mixed with GUS assay buffer (2 mM 4-methylumbelliferyl-D-glucuronide (4-MUG), 50 mM sodium phosphate buffer with pH 7.0, 10 mM β-mercaptoethanol, 10 mM Na₂EDTA, 0.1% sodium lauroylsarcosine, and 0.1% Triton X-100). The addition of the stop buffer (0.2 M Na₂CO₃) halted the reaction. Next, 4-MUG was hydrolyzed by GUS to produce 4-methylumbelliferone fluorochrome (4-MU). GUS activity was determined in triplicate with a microplate spectrofluorometer (Megellan, Taken). The excitation wavelength was 365 nm, and the emission wavelength was 455 nm.

3. Results

3.1. Expression of Six Selected Rice Genes under Abiotic Stresses as Revealed by Microarray Data

Six genes (Os03g0745000, Os11g0453900, Os06g0195800, Os09g0526600, Os05g0381400, and Os04g0107900) were identified as highly heat-responsive genes from the microarray data of rice panicles under heat shock (Figure 1(a)). Our further microarray analysis in rice flag leaves under heat stress also confirmed that these six genes were highly heat inducible, although the expression levels decreased more quickly on time course treatment in flag leaves (Figure 1(b)). Microarray data of these six genes in rice seedling collected through the GEO accession numbers GSE14275 and GSE6901 showed that five of them were highly induced by heat, drought, and salt but less sensitive to cold stress, with the exception that Os04g0107900 was extremely heat inducible but insensitive to drought stress (Figure 2).

[figures omitted; refer to PDF]

[figure omitted; refer to PDF]

3.2. Expression Pattern of the Six Selected Genes in Different Tissues

To examine the expression pattern of the six genes in different tissues, total RNA from the roots, stems, leaves, leaf sheath, spikes, and embryos were extracted. The results of qRT-PCR analysis showed that the six genes were expressed in all the organs under normal conditions. All of the genes were expressed at higher level in all organs, with the exception that the expression of Os09g0526600 was much lower in leaf sheath (Figure 3).

[figure omitted; refer to PDF]

3.3. Expression Patterns of the Six Genes under Abiotic Stresses and ABA Treatment as Validated by qRT-PCR

To validate the stress induced expression level of the six genes obtained from the above microarray data, expression patterns of these six genes under different abiotic stresses were further analyzed by qRT-PCR. The expression of the six genes was induced by 100 mmol L⁻¹ NaCl at 30 min and reached the highest level at 2 or 4 h and then began to decrease, except for the Os06g0195800 that reached the highest level at 8 h (Figure 4(a)). The expression of six genes was early upregulated by PEG treatment (Figure 4(b)). Compared to the responses to salt and PEG stresses, the expression of the six genes was rapidly induced to much higher level by heat treatment at 30 min (Figure 4(c)) along with 8.5-fold increase for Os11g0453900 and around 20-fold increase for the other five genes. The highest levels of these genes under heat stress were between 63- and 45-fold, which are much higher than the highest levels under salt and PEG stresses (less than 7-fold). The induced expression of Os09g0526600 and Os05g0381400 remained relatively steady under heat stress. The expression of the six genes was less affected by cold treatment with the highest expression of around 2-fold (Figure 4(d)). In addition, expression of all the six genes was affected by ABA to some extent (Figure 5).

[figures omitted; refer to PDF]

[figure omitted; refer to PDF]

3.4. Bioinformatics Analysis on the Promoters of Three Stress-Inducible Genes

Three highly stress-inducible genes (Os09g0526600, Os05g0381400, and Os04g0107900) were selected for promoter element and promoter-driven GUS expression analysis. Approximately 1.2 to 1.5 kb of the upstream region of the translation initiation codon (ATG) were predicted to be the promoter region of the stress-inducible genes and named OsHsfB2c (1220 bp), PM19 (1218 bp), and HSP90 (1437 bp), respectively. The promoter analysis by PlantCARE and PLACE showed some stress response-related elements, such as G-box, GC-motif, HSE, and ABRE in the three promoters (Figure 6).

[figure omitted; refer to PDF]

3.5. Expression of GUS Gene Driven by the Three Stress-Inducible Promoters in Transgenic Rice

The promoters of the three genes were fused to the gusA reporter gene to evaluate the heat and drought stress responses of these promoters in transgenic rice plants by monitoring the expression of gusA reporter gene and GUS enzyme activities. Total RNA was extracted from flag leaf and panicle in the presence or absence of drought and heat stress for qRT-PCR. The GUS gene expression levels significantly increased in both panicles and leaves for all the three genes under heat stress. In the leaves, the expression levels were upregulated around 8- to 11-fold. In panicles, the expression levels of OsHsfB2cp and PM19p were greatly increased to more than 20-fold, while the expression of Hsp90p was only increased about 6-fold. The GUS gene expression was less induced under drought both in panicles and leaves (Figure 7). GUS histochemical staining was prominently visualized in panicles and leaves in response to heat and drought stress (Figure 8). The three promoters showed a low activity in panicle and flag leaf under drought stress. Quantification of GUS activity using fluorescence assay demonstrated that the GUS activities driven by three promoters exhibited quite similar patterns under heat and drought stresses in both panicles and leaves of the transgenic rice lines (Figure 9). These data confirmed that three promoters of choice were highly heat inducible but weakly induced by drought treatment.

[figure omitted; refer to PDF]

[figure omitted; refer to PDF]

[figure omitted; refer to PDF]

4. Discussion

Microarray has been widely used to systematically investigate the molecular reactions by which plants respond and adapt to complicated environment [15–17]. It was also applied to profile the tissue- and stress-specific gene expression [18–20], as well as transcriptome of rice developing caryopses under high temperature [21]. Based on our microarray data and rice abiotic stress microarray data available in the public GEO database, we identified six highly heat-inducible genes in rice for the promoter analysis with the aim to develop heat or abiotic stress inducible promoters for plant genetic engineering.

The transcripts of the six genes were accumulated in all the tissues we checked under normal condition. It has been elucidated that rice heat shock transcription factor (OsHsf) and heat shock proteins (OsHsp) were usually induced in response to abiotic stresses [22–29]. Os03g0745000 (OsHsfA2a) and Os09g0526600 (OsHsfA2c) were quickly and exceptionally highly induced under heat stress [22–24, 26]. OsHsps were upregulated under heat, PEG, and ABA treatments [25, 29]. Our qRT-PCR analysis confirmed that the expression of six genes was induced under heat stress, moderately induced by salt and drought treatments, and less affected by cold treatment. High expression of abiotic stress tolerance genes in these important organs under stresses is a key requirement for a stress inducible promoter. These implied that the promoters of the selected genes may be useful for developing heat inducible promoters.

The GUS gene reporter system is one of the most effective techniques employed in the study of gene regulation in plant molecular biology [27]. In this study, we further evaluated the heat and drought stress responses of three promoters in panicle and flag leaf of the transgenic rice plants by qRT-PCR, histochemical staining, and fluorescence quantification. GUS gene expression and GUS protein activity in the leaves and panicles of the OsHsfB2cp and OsPM19p transgenic plants were significantly upregulated under heat stress but less induced under drought, which is supported by the bioinformatics analysis on the promoter sequences revealing that various stress-responsive elements such as HSE, ABRE, and G-box are located in the OsHsfB2cp and OsPM19p. These results further suggested that OsHsfB2cp and PM19p may be useful for developing heat inducible promoters.

High temperature can significantly affect various activities of plant. Leaf is the main organ for photosynthesis and is easily affected by heat stress. Reproduction organs are vulnerable to heat stress. Rice is most susceptible to heat injury during flowering, as pollen viability is particularly sensitive to heat stress; even 1-2 h of high temperature at anthesis would result in high spikelet sterility [22]. It was also reported that high temperature caused rice male sterility during pollen development [30]. A joint analysis of metabolome and transcriptome was performed to comprehend the impact of high temperature on rice grain filling [21]. The inducible expression of heat tolerant genes in reproductive organs will be helpful for rice and other crops to avoid production loss under heat stress. The present study showed that OsHsfB2cp and PM19p can drive high inducible gene expression in leave and panicles and may be useful for developing heat inducible promoters with further characterization and reconstruction of their cis-elements.

5. Conclusions

Plant growth and development, especially crop yield can be greatly affected by abiotic stresses. Global warming has been likely to induce heat stress more frequently and severely, resulting in serious yield reduction and thus threatening the security of grain production. Stress-inducible promoters can be used to drive specific gene expression under stresses for increased abiotic stress tolerance. Six highly heat responsive genes were identified from rice microarray data. Our qRT-PCR analysis confirmed that these six genes were highly heat inducible and moderately responsive to salt stress, PEG, and ABA treatments. Promoters of the three highly heat-inducible genes (OsHsfB2cp, PM19p, and Hsp90p) were used to drive GUS gene expression in rice. GUS gene expression, histochemical staining, and GUS activities analysis confirmed that OsHsfB2cp and PM19p were highly heat inducible in panicles and flag leaves. Our work suggests that OsHsfB2c and PM19 could be used as highly heat inducible promoters to drive specific gene expression under heat stress for increased stress tolerance.

Authors’ Contribution

Wirat Rerksiri and Xianwen Zhang contributed equally to this work.

References

[1] J.-P. Gao, D.-Y. Chao, H.-X. Lin, "Understanding abiotic stress tolerance mechanisms: recent studies on stress response in rice," Journal of Integrative Plant Biology, vol. 49 no. 6, pp. 742-750, DOI: 10.1111/j.1744-7909.2007.00495.x, 2007.

[2] H. R. Lafitte, A. Ismail, J. Bennet, "Abiotic stress tolerance in rice for Asia: progress and the future, in new directions for a diverse planet," Plant Cell, vol. 14 no. 3, pp. 537-545, 2004.

[3] S. H. Li, X. H. Tian, Y. P. Huang, A. Y. Liu, "Preliminary analysis of high temperature harm in middle-season rice at flowering stage in Jianghan Plain in the latest 50 years," Chinese Journal of Agrometeorology, vol. 28 no. 1, 2007.

[4] S. Sarangi, J. Ghosh, A. Bora, S. Das, A. B. Mandal, "Agrobacterium-mediated genetic transformation of indica rice varieties involving Am-SOD gene," Indian Journal of Biotechnology, vol. 10 no. 1, 2011.

[5] S.-M. Yu, S.-S. Ko, C.-Y. Hong, H.-J. Sun, Y.-I. Hsing, C.-G. Tong, T.-H. D. Ho, "Global functional analyses of rice promoters by genomics approaches," Plant Molecular Biology, vol. 65 no. 4, pp. 417-425, DOI: 10.1007/s11103-007-9232-1, 2007.

[6] G. Corrado, M. Karali, "Inducible gene expression systems and plant biotechnology," Biotechnology Advances, vol. 27 no. 6, pp. 733-743, DOI: 10.1016/j.biotechadv.2009.05.006, 2009.

[7] X. W. Zhang, J. P. Li, A. L. Liu, J. Zou, X. Y. Zhou, J. H. Xiang, W. Rerksiri, Y. Peng, X. Y. Xiong, X. B. Chen, "Expression profile in rice panicle: insights into heat response mechanism at reproductive stage," PLos ONE, vol. 7 no. 11, 2012.

[8] X. W. Zhang, W. Rerksiri, A. L. Liu, X. Y. Zhou, H. R. Xiong, J. H. Xiang, X. B. Chen, X. Y. Xiong, "Transcriptome profile reveals heat response mechanism at molecular and metabolic levels in rice flag leaf," Gene, vol. 530 no. 2, pp. 185-192, 2013.

[9] W. Hu, G. Hu, B. Han, "Genome-wide survey and expression profiling of heat shock proteins and heat shock factors revealed overlapped and stress specific response under abiotic stresses in rice," Plant Science, vol. 176 no. 4, pp. 583-590, DOI: 10.1016/j.plantsci.2009.01.016, 2009.

[10] M. Jain, A. Nijhawan, R. Arora, P. Agarwal, S. Ray, P. Sharma, S. Kapoor, A. K. Tyagi, J. P. Khurana, "F-Box proteins in rice. Genome-wide analysis, classification, temporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress," Plant Physiology, vol. 143 no. 4, pp. 1467-1483, DOI: 10.1104/pp.106.091900, 2007.

[11] K. J. Livak, T. D. Schmittgen, "Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method," Methods, vol. 25 no. 4, pp. 402-408, DOI: 10.1006/meth.2001.1262, 2001.

[12] R. Hofgen, L. Willmitzer, "Storage of competent cells for Agrobacterium transformation," Nucleic Acids Research, vol. 16 no. 20, article 9877,DOI: 10.1093/nar/16.20.9877, 1988.

[13] S. Toki, N. Hara, K. Ono, H. Onodera, A. Tagiri, S. Oka, H. Tanaka, "Early infection of scutellum tissue with Agrobacterium allows high-speed transformation of rice," Plant Journal, vol. 47 no. 6, pp. 969-976, DOI: 10.1111/j.1365-313X.2006.02836.x, 2006.

[14] R. A. Jefferson, T. A. Kavanagh, M. W. Bevan, "GUS fusions: beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants," EMBO Journal, vol. 6 no. 13, pp. 3901-3907, 1987.

[15] S. Gonzali, E. Loreti, G. Novi, A. Poggi, A. Alpi, P. Perata, "The use of microarrays to study the anaerobic response in Arabidopsis," Annals of Botany, vol. 96 no. 4, pp. 661-668, DOI: 10.1093/aob/mci218, 2005.

[16] P. O. Brown, D. Botstein, "Exploring the new world of the genome with DNA microarrays," Nature Genetics, vol. 21 no. 1, pp. 33-37, DOI: 10.1038/4462, 1999.

[17] D. J. Lockhart, E. A. Winzeler, "Genomics, gene expression and DNA arrays," Nature, vol. 405 no. 6788, pp. 827-836, DOI: 10.1038/35015701, 2000.

[18] M. Fujita, Y. Horiuchi, Y. Ueda, Y. Mizuta, T. Kubo, K. Yano, S. Yamaki, K. Tsuda, T. Nagata, M. Niihama, H. Kato, S. Kikuchi, K. Hamada, T. Mochizuki, T. Ishimizu, H. Iwai, N. Tsutsumi, N. Kurata, "Rice expression atlas in reproductive development," Plant and Cell Physiology, vol. 51 no. 12, pp. 2060-2081, DOI: 10.1093/pcp/pcq165, 2010.

[19] L. Wang, W. Xie, Y. Chen, W. Tang, J. Yang, R. Ye, L. Liu, Y. Lin, C. Xu, J. Xiao, Q. Zhang, "A dynamic gene expression atlas covering the entire life cycle of rice," Plant Journal, vol. 61 no. 5, pp. 752-766, DOI: 10.1111/j.1365-313X.2009.04100.x, 2010.

[20] M. Endo, T. Tsuchiya, K. Hamada, S. Kawamura, K. Yano, M. Ohshima, A. Higashitani, M. Watanabe, M. Kawagishi-Kobayashi, "High temperatures cause male sterility in rice plants with transcriptional alterations during pollen development," Plant and Cell Physiology, vol. 50 no. 11, pp. 1911-1922, DOI: 10.1093/pcp/pcp135, 2009.

[21] H. Yamakawa, M. Hakata, "Atlas of rice grain filling-related metabolism under high temperature: joint analysis of metabolome and transcriptome demonstrated inhibition of starch accumulation and induction of amino acid accumulation," Plant and Cell Physiology, vol. 51 no. 5, pp. 795-809, DOI: 10.1093/pcp/pcq034, 2010.

[22] S. V. K. Jagadish, R. Muthurajan, R. Oane, T. R. Wheeler, S. Heuer, J. Bennett, P. Q. Craufurd, "Physiological and proteomic approaches to address heat tolerance during anthesis in rice ( Oryza sativa L.)," Journal of Experimental Botany, vol. 61 no. 1, pp. 143-156, DOI: 10.1093/jxb/erp289, 2010.

[23] H. Chauhan, N. Khurana, P. Agarwal, P. Khurana, "Heat shock factors in rice ( Oryza sativa L.): genome-wide expression analysis during reproductive development and abiotic stress," Molecular Genetics and Genomics, vol. 286 no. 2, pp. 171-187, DOI: 10.1007/s00438-011-0638-8, 2011.

[24] A. Singh, D. Mittal, D. Lavania, M. Agarwal, R. C. Mishra, A. Grover, "OsHsfA2c and OsHsfB4b are involved in the transcriptional regulation of cytoplasmic OsClpB (Hsp100) gene in rice ( Oryza sativa L.)," Cell Stress and Chaperones, vol. 17 no. 2, pp. 243-254, DOI: 10.1007/s12192-011-0303-5, 2012.

[25] S. F. Ye, S. W. Yu, L. B. Shu, J. H. Wu, A. Z. Wu, L. J. Luo, "Expression profile analysis of 9 heat shock protein genes throughout the life cycle and under abiotic stress in rice," Chinese Science Bulletin, vol. 57 no. 4, pp. 336-343, DOI: 10.1007/s11434-011-4863-7, 2012.

[26] A.-L. Liu, J. Zou, X.-W. Zhang, X.-Y. Zhou, W.-F. Wang, X.-Y. Xiong, L.-Y. Chen, X.-B. Chen, "Expression profiles of class a rice heat shock transcription factor genes under abiotic stresses," Journal of Plant Biology, vol. 53 no. 2, pp. 142-149, DOI: 10.1007/s12374-010-9099-6, 2010.

[27] S. Fior, P. D. Gerola, "Impact of ubiquitous inhibitors on the GUS gene reporter system: evidence from the model plants Arabidopsis, tobacco and rice and correction methods for quantitative assays of transgenic and endogenous GUS," Plant Methods, vol. 5 no. 1, article 19,DOI: 10.1186/1746-4811-5-19, 2009.

[28] H. Chauhan, N. Khurana, P. Agarwal, P. Khurana, "Heat shock factors in rice ( Oryza sativa L.): genome-wide expression analysis during reproductive development and abiotic stress," Molecular Genetics and Genomics, vol. 286 no. 2, pp. 171-187, DOI: 10.1007/s00438-011-0638-8, 2011.

[29] J. Zou, A. Liu, X. Chen, X. Zhou, G. Gao, W. Wang, X. Zhang, "Expression analysis of nine rice heat shock protein genes under abiotic stresses and ABA treatment," Journal of Plant Physiology, vol. 166 no. 8, pp. 851-861, DOI: 10.1016/j.jplph.2008.11.007, 2009.

[30] M. Endo, T. Tsuchiya, K. Hamada, S. Kawamura, K. Yano, M. Ohshima, A. Higashitani, M. Watanabe, M. Kawagishi-Kobayashi, "High temperatures cause male sterility in rice plants with transcriptional alterations during pollen development," Plant and Cell Physiology, vol. 50 no. 11, pp. 1911-1922, DOI: 10.1093/pcp/pcp135, 2009.