- 1College of Life Sciences, Nankai University, Tianjin, China
- 2School of Life Sciences and Technology, University of Electronic Science and Technology of China, Chengdu, China
- 3Center for Informational Biology, University of Electronic Science and Technology of China, Chengdu, China
Dehydration responsive element binding proteins are transcription factors of the plant-specific AP2 family, many of which contribute to abiotic stress responses in several plant species. We investigated the possibility of increasing drought tolerance in the traditional Chinese medicinal herb, Salvia miltiorrhiza, through modulating the transcriptional regulation of AtDREB1C in transgenic plants under the control of a constitutive (35S) or drought-inducible (RD29A) promoter. AtDREB1C transgenic S. miltiorrhiza plants showed increased survival under severe drought conditions compared to the non-transgenic wild-type (WT) control. However, transgenic plants with constitutive overexpression of AtDREB1C showed considerable dwarfing relative to WT. Physiological tests suggested that the higher chlorophyll content, photosynthetic capacity, and superoxide dismutase, peroxidase, and catalase activity in the transgenic plants enhanced plant drought stress resistance compared to WT. Transcriptome analysis of S. miltiorrhiza following drought stress identified a number of differentially expressed genes (DEGs) between the AtDREB1C transgenic lines and WT. These DEGs are involved in photosynthesis, plant hormone signal transduction, phenylpropanoid biosynthesis, ribosome, starch and sucrose metabolism, and other metabolic pathways. The modified pathways involved in plant hormone signaling are thought to be one of the main causes of the increased drought tolerance of AtDREB1C transgenic S. miltiorrhiza plants.
Introduction
Drought is one of the primary abiotic stresses limiting global agricultural production (Zeppel et al., 2014; Shavrukov et al., 2015; Wang et al., 2016). Plant genetic engineering provides an important approach for improving drought resistance in crops via the introduction of transgene(s) to mitigate physiological and agronomic penalties associated with drought (Bihmidine et al., 2013; Patel et al., 2015; Mendiondo et al., 2016). One strategy for increasing drought tolerance is to genetically engineer plants with transcription factors (TF) that control the expression of several abiotic stress response genes (Cabello et al., 2014; de Paiva Rolla et al., 2014; Rehman and Mahmood, 2015). TFs, which are important regulators of gene transcription in various organisms, are classified into a number of families based on the sequences and structures of their DNA-binding domains (Lata and Prasad, 2011; Yamasaki et al., 2013). For example, the plant-specific APETALA2 (AP2)/ethylene-responsive element–binding protein (EREBP) family comprises five subfamilies: AP2, RAV, DREB, ERF, and others (Lata and Prasad, 2011; Licausi et al., 2013; Vahala et al., 2013). DREB TFs, some of the earliest TFs found to be involved in abiotic stress, comprise 56 members divided into six groups (A-1 to A-6) (Hichri et al., 2015). All DREB TFs, which contain a single AP2 domain, bind to the consensus sequence A/GCCGAC in the C repeat or dehydration response element (DRE) in the promoters of genes that are induced under low temperatures and/or water stress (Kovalchuk et al., 2013; Licausi et al., 2013; Dou et al., 2015; Park et al., 2015; Jia et al., 2016). Overexpressing DREB TFs in transgenic plants markedly improves plant survival under both drought and low-temperature conditions (Shan et al., 2007; Byun et al., 2015; Kidokoro et al., 2015; Shavrukov et al., 2015).
When designing genetic engineering strategies to combat abiotic stress in plants, the timing and developmental regulation of the targeted transgene(s) must be coordinated with the occurrence of the corresponding stress (Bihmidine et al., 2013). Misexpression of the targeted transgene(s) in the absence of stress may negatively affect crop performance (Kovalchuk et al., 2013; Shavrukov et al., 2015). For instance, ectopic expression of CBF/DREB genes via fusion to the constitutive 35S CaMV promoter in Arabidopsis (Yamasaki and Randall, 2016), tobacco (Zhou et al., 2014), potato (Storani et al., 2015), and tomato (Li et al., 2012) improved plant performance under abiotic stress but resulted in developmental abnormalities that, under non-stress environments, penalized plant performance. Therefore, precise regulation of the timing and expression level of the transgene is an important consideration when designing constructs to be introduced into a plant genome. For example, when stress-inducible promoters were used in transgenic cassettes carrying transgenes thought to improve abiotic stress tolerance in rice (Ravikumar et al., 2014), barley (Kovalchuk et al., 2013), wheat (Morran et al., 2011), and tomato (Rai et al., 2013), phenotypic development was normal in the transgenic lines under non-stress conditions, yet plant performance was still improved under stress conditions.
Salvia miltiorrhiza, commonly known as Danshen, is a widely used traditional Chinese medicinal plant; its roots are used to treat hyperlipidemia as well as cardiovascular and cerebrovascular diseases (Cui et al., 2015; Luo et al., 2015; Xu et al., 2015; Zhou et al., 2016). The recent decline in environmental conditions, especially the increased drought conditions in soil due to higher global temperatures, has resulted in a sharp reduction in the planting area and yield of S. miltiorrhiza (Wu et al., 2014; Zhang et al., 2016). Therefore, it is important to improve the stress resistance and adaptability of this plant. However, to date, efforts to improve this plant by genetic engineering have mainly focused on increasing the contents of active ingredients such as tanshinone and phenolic acids (Zhang Y. et al., 2014; Zhao et al., 2015). By contrast, few studies have focused on improving the stress tolerance of S. miltiorrhiza using transgenic methods (Wang et al., 2012; Wu et al., 2014).
Arabidopsis contains three DREB1/CBF genes in tandem on chromosome 4 in the order DREB1B/CBF1, DREB1A/CBF3, and DREB1C/CBF2 (Park et al., 2015; Jia et al., 2016). The high level of similarity among these genes, as well as their tandem organization and the observation that they have the same transcriptional orientation, suggest that share a common origin (Kang et al., 2013). Overexpression of AtDREB1A, AtDREB1B, and AtDREB1C in Arabidopsis helped confirm that they have matching functional activities; the phenotypes of these lines mimicked the multiple biochemical changes associated with cold acclimation (Gilmour et al., 2004). We previously demonstrated that the drought resistance of both AtDREB1A and AtDREB1B transgenic S. miltiorrhiza plants was significantly improved under severe drought stress (Wei et al., 2016a,b). However, the effects of AtDREB1C overexpression in transgenic S. miltiorrhiza has not previously been examined, and whether overexpressing AtDREB1C can improve the drought tolerance of transgenic plants remains to be determined.
In this study, we demonstrated that ectopic expression of Arabidopsis DREB1C led to substantial improvements in the capacity of S. miltiorrhiza to survive under severe drought stress conditions. We assessed the effects of ectopic expression of AtDREB1C on drought resistance in S. miltiorrhiza by investigating shoot growth, the photosynthetic system (stomatal conductance rate, transpiration rate, and net photosynthesis rate), antioxidative metabolism (both antioxidant enzymes and lipid peroxidation), and transcriptomic changes in transgenic plants. We found that expressing AtDREB1C under the control of the stress-inducible RD29A promoter greatly improved drought tolerance and minimized the negative effects of this gene on plant growth, as were observed using the strong constitutive CaMV35S promoter.
Materials and Methods
Vector Construction
The full-length open reading frame of AtDREB1C (GenBank accession number: AB013817) and the Arabidopsis thaliana RD29A promoter region (GenBank accession number: AY973635) were PCR-amplified from genomic DNA using primers containing the restriction sites NcoI/BstEII and BamHI/NcoI, respectively (Supplementary Table S1). The resulting PCR products were cloned into the pEASY-T1 vector (TransGen Biotech, Beijing, China). Following digestion, the AtDREB1C fragment was inserted into binary plant vector pCAMBIA3301, containing the CaMV35S constitutive promoter (Figure 1). To construct the pRD29A::AtDREB1C plasmid, the 35S promoter from p35S::AtDREB1C was replaced with the RD29A promoter using BamHI and NcoI (Figure 1). The vectors were individually transformed into Agrobacterium tumefaciens strain LBA4404 using standard heat-shock methods prior to use in the plant transformation experiments.
FIGURE 1. T-DNA regions of the binary vectors employed for Agrobacterium tumefaciens-mediated transformation. (A) Diagram of binary vector pCAMBIA3301. (B) Binary vector p35S::AtDREB1C, containing AtDREB1C driven by the 35S promoter. (C) Binary vector pRD29A::AtDREB1C, containing AtDREB1C driven by the RD29A promoter. 35S-P, cauliflower mosaic virus-CaMV35S RNA promoter; RD29A-P, Arabidopsis thaliana RD29A promoter; 35S-T, CaMV35S poly A; NOS-T, 3′ terminator region of the nopaline synthase gene; Bar, phosphinothricin (R); RB, right border; LB, left border.
Plant Materials, Genetic Transformation, and Molecular Characterization
Salvia miltiorrhiza seeds were surface sterilized with 0.1% HgCl2 and cultured on hormone-free MS agar medium containing 30 g L-1 sucrose and 8 g L-1 agar. The cultures were incubated in a growth chamber at 25 ± 2°C under a 16 h light/8 h dark cycle provided by a cool white fluorescent lamp, providing a PPFD level of 26 μmol m-2 s-1. The S. miltiorrhiza plants were transformed via the Agrobacterium-mediated leaf-disk method (Wei et al., 2016b), followed by weekly subculturing of the leaf disks. After 35 days, the shoots were cut off of the plants and transferred into selective rooting medium, leading to the production of roots within 2 weeks. The rooted plantlets were cut into segments from their internodes and propagated on a 1/2-strength MS basal medium; 1.5-month-old plantlets were used for PCR screening and for gene expression evaluations. Positive transformants were transferred to soil for further observation and stress treatments.
The cetyltrimethylammonium bromide (CTAB) method was used to isolate genomic DNA from young leaves (Wei et al., 2016b). Positive pCAMBIA3301 transgenic lines were detected using GUS-F/R. Primer sets used for the p35S::AtDREB1C transgenic lines were designed based on the sequences of the 35S promoter of pCAMBIA3301 (for the forward primer, 35S-F) and AtDREB1C (for the reverse primer, AtDREB1C-R). Similarly, RD29A-F/AtDREB1C-R was used to detect pRD29A::AtDREB1C transgenic plants (Supplementary Table S1). Primers 18S-F/R were designed based on the conserved regions of the control gene, 18S rRNA.
A modified CTAB method was used to extract total RNA from WT and transgenic lines after drought stress treatment as described previously (Wei et al., 2016b). Semi-quantitative RT-PCR was performed to investigate GUS or AtDREB1C expression in transgenic S. miltiorrhiza lines using specific primers (Supplementary Table S1). The actin gene was used as an internal control.
Stress Treatment
Transgenic S. miltiorrhiza seedlings were grown in a growth chamber (4°C) for cold treatment. For abscisic acid (ABA), drought, and salt stress treatments, the seedlings were carefully rinsed with water and transferred to 1/4 MS solutions containing 100 μM ABA, 20% PEG6000, or 250 mM NaCl, respectively. Treatment with 1/4 MS solution lacking supplements was used for the control. After 6 h, control and stress-treated seedlings were harvested, immediately frozen in liquid nitrogen, and stored at -80°C for further analysis (Wei et al., 2016a).
To analyze drought tolerance in transgenic S. miltiorrhiza plants, plantlets in sterile tubes in solid MS medium were transferred to sterile soil containing nutrients and cultured in a growth chamber at 25 ± 2°C under a 16-h/8-h photoperiod (light intensity, 50 μmol m-2 s-1). Prior to drought stress treatment, the pots for both the well-watered and drought-stress treatments were saturated with water and allowed to drain overnight. Plants at uniform developmental stages were chosen for stress treatments. On day 54, the plants were stressed by withholding water for 9 days, the relative water content (RWC) and other physiological parameters were measured in the leaves using the third-youngest leaf of three randomly selected plants per treatment. After 11 days of drought stress treatment, since all plants showed symptoms of severe wilting, the plants were watered again.
Determining Physiological Indices
O2- and H2O2 levels were measured as previously described (Wei et al., 2016a). The RWC and chlorophyll content of both WT and transgenic plants were determined as described (Wei et al., 2016b). Leaf gas exchange parameters, including stomatal conductance (gs), transpiration rate (E), and net photosynthesis rate (Pn) were measured using a portable gas exchange system (CIRAS-2, PP-System, Amesbury, MA, USA) (Wei et al., 2016a). Malondialdehyde (MDA) content and superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activity were evaluated as described (Wei et al., 2016a).
Transcriptome Analysis
Total RNA from the mixed sample (three individual plants) was isolated using the modified CTAB method as described above. Following RNase-free DNase I treatment (New England BioLabs) for 30 min at 37°C to remove residual DNA, the samples were sequenced by Biomarker Technologies Co., Ltd. (Beijing, China), and cDNA library construction and Illumina HiSeq2500 sequencing were performed as described previously (Wei et al., 2016b).
Validation of RNA-seq Data by qRT-PCR
A total of 30 pairs of gene-specific primers (Supplementary Table S1) were designed to produce amplicons for validating the RNA-seq data. Quantitative reverse-transcription PCR (qRT-PCR) was performed on an iQ5.0 instrument (Bio-Rad, USA) using SYBR Green qPCR kits (Roche) according to the manufacturer’s protocol. Relative gene expression levels were calculated using the 2-ΔΔCt method; expression levels were quantified by normalization against Actin and GAPDH. All assays for each gene were performed in triplicate synchronously under identical conditions.
Statistical Analysis
All data are presented as the means ± standard error (SE) of at least three replicates. Statistical analysis (analysis of variance; ANVOA) was performed using SAS software version 9.1 (SAS Institute Inc., Cary, NC, USA) to test for significant differences between WT and transgenic lines. The mean values of each treatment group were compared using Duncan’s test at P < 0.05.
Results
Generation of Transgenic Lines
To examine the role of AtDREB1C/CBF2 in S. miltiorrhiza, we constructed transgenic plants harboring p35S::AtDREB1C and pRD29A::AtDREB1C (Figure 1), expressing AtDREB1C under the control of the constitutive CaMV35S and the stress-inducible RD29A promoter, respectively. We selected putative transgenic lines harboring these constructs on medium containing Basta and confirmed them using PCR with gene-specific primers (Figures 2A,B). We performed RT-PCR analysis using cDNA from the transgenic plants (p35S::AtDREB1C and pRD29A::AtDREB1C), which amplified the 918 bp AtDREB1C fragment, under stress conditions (Figure 2C). A 1081 bp PCR product corresponding to the region harboring the GUS-specific primer from the pCAMBIA3301 vector control plants was detected (data not shown). We selected two independent p35S::AtDREB1C transgenic lines (p35S::AtDREB1C-01 and 04) and two pRD29A::AtDREB1C transgenic lines (pRD29A::AtDREB1C-02 and 07), which produced high levels of AtDREB1C after drought stress, for further analysis.
FIGURE 2. Identification of transgenic AtDREB1C Salvia miltiorrhiza plants via PCR. (A) PCR analysis of primary transformants using specific primers for 35S::AtDREB1C. Lanes 1 and 16, size markers; lanes 2–13, DNA from putative transformants; lane 14, untransformed control; lane 15, p35S::AtDREB1C. The 18S rRNA gene was used as a control. (B) PCR analysis of primary transformants using specific primers for RD29A::AtDREB1C. Lanes 1 and 16, size marker; lanes 2–13, DNA from putative transformants; lane 14, untransformed control; lane 15, pRD29A::AtDREB1C. The 18S rRNA gene was used as a control. (C) RT-PCR analysis of AtDREB1C expression in transgenic plants using under drought conditions. Lanes 1 and 16, size markers; lanes 2–7, DNA from PCR-positive p35S::AtDREB1C transformants; lanes 8–13, DNA from PCR-positive pRD29A::AtDREB1C transformants; lane 14, untransformed control; lane 15, p35S::AtDREB1C and pRD29A::AtDREB1C. The actin gene served as the internal control.
AtDREB1C Expression Patterns in Transgenic S. miltiorrhiza Plants
We performed RT-PCR to investigate the effects of exogenous ABA treatment and various stress conditions on AtDREB1C expression. We detected a relatively high level of AtDREB1C expression in the p35S::AtDREB1C transgenic lines (p35S::AtDREB1C-01 and 04) under both untreated and various stress-treated conditions (Figure 3). However, AtDREB1C was expressed at very low levels in the pRD29A::AtDREB1C transgenic lines (pRD29A::AtDREB1C-02 and 07) without stress, but strongly induced after exogenous ABA application and various stress treatments (Figure 3).
FIGURE 3. Induction of AtDREB1C expression under various conditions. ABA- and stress-induced AtDREB1C expression determined by RT-PCR. RNA was isolated from seedlings subjected to stress treatments as follows: 100 μM ABA (6 h), 250 mM NaCl (6 h), or 20% PEG6000 (6 h) in 1/4 Murashige and Skoog (MS) solution; cold (4°C, 6 h). Untreated indicates untreated control plants, and 1/4 MS indicates seedlings treated with 1/4 MS solution lacking supplements. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on.
Morphological Changes and Drought Tolerance in Transgenic S. miltiorrhiza Plants
We observed the growth of the AtDREB1C transgenic S. miltiorrhiza plants under normal growing conditions and compared them with WT plants and pCAMBIA3301 vector control plants. Both p35S::AtDREB1C transgenic lines showed significant growth retardation compared to the control plants, whereas the growth of the pRD29A::AtDREB1C transgenic lines was unaffected (Figure 4A). The average plant height of the p35S::AtDREB1C transgenic lines was less than 1/3 that of WT plants (Figure 4C). In addition, the leaves of p35S::AtDREB1C transgenic plants were darker and shorter than those of WT (Figure 4B). The leaf size (length and width) was also smaller in the p35S::AtDREB1C plants than in WT (Figure 4D).
FIGURE 4. Growth of AtDREB1C transgenic S. miltiorrhiza plants. (A) Wild-type and transgenic seedlings. (B) Wild-type and transgenic leaves. (C) Plant height (cm) in wild-type and transgenic lines. (D) Leaf size (cm) in wild-type and transgenic lines. The labels under the graphs correspond to the plant materials shown in (A) and (B). WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control. Bars represent the mean + SE of three independent experiments. ∗∗Significantly different at the P < 0.01 level compared to WT.
To investigate whether AtDREB1C expression driven by the different promoters was associated with increased stress tolerance, we subjected WT and transgenic plantlets to 11 days of drought stress. At the early stage of treatment (3 days), all plants showed normal growth (Figure 5). At day 6 of treatment, the leaves of WT and the control line pCAMBIA3301-18 rapidly withered, whereas the growth of p35S::AtDREB1C and pRD29A::AtDREB1C plants was unaffected by drought stress. After 9 days of drought stress, both WT and pCAMBIA3301-18 plants exhibited severe wilting. By contrast, the growth of p35S::AtDREB1C and pRD29A::AtDREB1C transgenic plants was almost normal under drought stress treatment. However, on the last day of drought treatment (day 11), all plants, including p35S::AtDREB1C and pRD29A::AtDREB1C transgenic plants, were severely wilted (Figure 5). After re-watering, both p35S::AtDREB1C and pRD29A::AtDREB1C transgenic plants recovered quickly, whereas the growth of drought-stressed WT and pCAMBIA3301-18 plants was severely inhibited; the plants failed to recover and ultimately died (Figure 5). All lines, including WT and pCAMBIA3301-18, grew well under well-watered conditions.
FIGURE 5. Morphological changes in wild-type and transgenic plants under drought stress. D, drought; R, re-watered; W, well-watered; WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control.
Ectopic Expression of AtDREB1C Reduces ROS Accumulation under Drought Stress
Since drought stress increases ROS production, we investigated whether transgenic modification would reduce ROS accumulation after drought stress by examining O2- and H2O2 accumulation in the plants using NBT and DAB staining, respectively. Under normal conditions, little NBT staining was detected in the WT and pCAMBIA3301-18, whereas after drought stress, clear NBT staining was detected in both lines (Figure 6A). By contrast, little NBT staining was detected in the AtDREB1C transgenic lines, even under drought stress. In addition, lower levels of H2O2 were detected in the transgenic lines than in the WT and vector control lines, as revealed by DAB staining (Figure 6B). Quantitative analysis of O2- and H2O2 levels indicated that the transgenic lines produced significantly lower levels of ROS than WT (Figures 6C,D).
FIGURE 6. Changes in O2- (A,C) and H2O2 (B,D) levels in wild-type (WT) and transgenic leaves from plants subjected to drought stress. Drought-stressed leaves were incubated in nitroblue tetrazolium (NBT) or diaminobenzidine (DAB) solution. Blue staining indicates the location and levels of O2- (A). Brown staining indicates H2O2 accumulation (B). O2- production rate (C) and H2O2 content (D) were measured after drought treatment. The labels under the graphs correspond to the plant materials shown in (A) and (B). Control, plants growing under normal conditions; Drought, plants growing after drought stress treatment (6 days). WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control. Bars represent the mean + SE of three independent experiments. ∗∗Significantly different at the P < 0.01 level compared to WT.
AtDREB1C Expression Affects MDA Content, as Well as SOD, POD, and CAT Activity
We measured MDA content, an indicator of lipid peroxidation, in the leaves of WT and transgenic lines. After 9 days of drought treatment, we detected significantly lower MDA levels (P < 0.01) in the transgenic lines than in WT and control line pCAMBIA3301-18 (Figure 7A). We further analyzed drought stress responses in the plants by monitoring SOD, POD, and CAT activity; these enzymes scavenge harmful ROS that accumulate under stress. The trends in the activities of these antioxidant enzymes were similar in all lines after 9 days of drought stress (Figures 7B–D). The activities of all antioxidant enzymes were higher in the AtDREB1C transgenic lines than in WT and pCAMBIA3301-18, with a slight difference observed between p35S::AtDREB1C and pRD29A::AtDREB1C transgenic plants. These findings suggest that the enhanced drought tolerance in the AtDREB1C transgenic lines (both the 35S and RD29A promoter lines) is related to their reduced MDA content and increased antioxidant enzyme activity.
FIGURE 7. Malondialdehyde levels (A) and SOD (B), POD (C), and CAT (D) activity in WT and transgenic S. miltiorrhiza plants after withholding watering for 9 days. WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control. Bars represents the mean + SE of three independent experiments. ∗∗Significantly different at the P < 0.01 level compared to WT.
Overexpressing AtDREB1C Increases RWC and Chlorophyll Content in Transgenic S. miltiorrhiza Plants
We calculated RWC and chlorophyll content in the plants after 9 days of drought stress. As shown in Figure 8A, the RWC was higher in all AtDREB1C transgenic lines than in the WT and vector control line (pCAMBIA3301-18), whereas there was no significant difference in these values between the p35S::AtDREB1C and pRD29A::AtDREB1C transgenic lines. The transgenic lines had significantly higher chlorophyll contents than WT and the vector control (pCAMBIA3301-18; P < 0.01), with p35S::AtDREB1C-04 plants showing the highest value (Figure 8B).
FIGURE 8. Effects of drought stress on relative water content (RWC) (A) and chlorophyll content (B) in WT and transgenic S. miltiorrhiza plants after withholding watering for 9 days. WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control. Bars represent the mean + SE of three independent experiments. ∗Significantly different at the P < 0.05 level compared to WT. ∗∗Significantly different at the P < 0.01 level compared to WT.
Increased Photosynthetic Capacity in AtDREB1C Transgenic Plants under Drought Stress
To assess the photosynthetic capacity of the transgenic lines, we used a portable photosynthesis system (CIRAS-2) to measure Pn, E, and gs-values over the course of drought treatment in soil. Under well-watered conditions, these values varied little among S. miltiorrhiza lines throughout the study period (data not shown). However, when subjected to drought stress, the WT and control line (pCAMBIA3301-18) exhibited a significant decline in Pn, E, and gs, whereas the AtDREB1C lines were less strongly affected by this treatment. The greatest difference was observed on day 9 of drought treatment, with no significant differences detected between the p35S::AtDREB1C and pRD29A::AtDREB1C lines (Figure 9).
FIGURE 9. Effects of drought stress on leaf gas exchange parameters in WT and transgenic S. miltiorrhiza plants after withholding watering for 9 days. (A) Net photosynthesis rate (Pn). (B) Transpiration rate (E). (C) Stomatal conductance (gs). WT, wild type. The transgenic lines are indicated by the construct names followed by the line number: p35S::AtDREB1C-01 refers to p35S::AtDREB1C line 01 and so on; pCAMBIA3301-18 was used for the positive control. Bars represent the mean + SE of three independent experiments. ∗∗Significantly different at the P < 0.01 level compared to WT.
Transcriptome Sequencing of AtDREB1C-Overexpressing S. miltiorrhiza under Drought Stress
To further investigate the molecular mechanisms underlying the enhanced drought tolerance in the AtDREB1C overexpression lines, we compared the global expression profiles of lines p35S::AtDREB1C-04 and pRD29A::AtDREB1C-07 versus WT following drought treatment. Mixed RNA samples from three separate 51-day-old plants (after 6 days of drought) were used to generate a cDNA library, yielding a total of 18,810,158, 17,445,426, and 18,966,357 clean reads for the WT, p35S::AtDREB1C-04, and pRD29A::AtDREB1C-07 cDNA libraries, respectively. The guanine–cytosine (GC) contents in these libraries were 50.16, 48.50, and 49.77%, respectively, and all Q30 percentages were higher than 85% (Supplementary Table S2). As no reference genome was available, we constructed a de novo transcriptome assembly by combining all clean reads from different cDNA libraries. High-throughout sequencing identified 13,153,489 contigs, 161,844 transcripts, and 73,748 unigenes (Supplementary Table S3). To obtain functional annotations, we assigned the unigenes to categories in the COG, GO, KEGG, KOG, Pfam, Swiss-Prot, and NR databases (Supplementary Table S4). A total of 43,428 unigenes were identified, 14,182 of which were enriched in the COG database, while 26,570 were analyzed in the GO database, and 19,338 were mapped to KEGG pathways.
We analyzed the expression of the unigenes in the WT and two transgenic lines using Bowtie (Langmead et al., 2009) and RSEM software (Li and Dewey, 2011) and normalized the values by FPKM (Fragments Per Kilobase of transcript per Million mapped reads) (Trapnell et al., 2010). A total of 1,970 differentially expressed genes (DEGs) were identified via clustering analysis (Figure 10), 704 of which were same between two comparisons (WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07; Figure 11). Compared with WT, 629 unigenes were upregulated and 351 were downregulated in p35S::AtDREB1C-04, while 1,148 unigenes were upregulated and 378 were downregulated in pRD29A::AtDREB1C-07 (Supplementary Table S5). Additionally, compared with p35S::AtDREB1C-04, 342 unigenes were upregulated and 101 were downregulated in pRD29A::AtDREB1C-07 transgenic S. miltiorrhiza (Supplementary Table S5).
FIGURE 10. Clustering analysis of DEGs in WT and AtDREB1C transgenic S. miltiorrhiza plants based on their expression profiles obtained by RNA-Seq. WT, wild type; p35S::AtDREB1C-04, p35S::AtDREB1C transgenic line 04; pRD29A::AtDREB1C-07, pRD29A::AtDREB1C transgenic line 07. The color scale corresponds to the log2 (FPKM) values of genes in various samples.
FIGURE 11. Venn diagram analysis of DEGs between two comparisons (WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07). WT, wild type; p35S::AtDREB1C-04, p35S::AtDREB1C transgenic line 04; pRD29A::AtDREB1C-07, pRD29A::AtDREB1C transgenic line 07.
To validate the results of expression profiling, we selected 40 genes that were significantly upregulated and 20 that were significantly downregulated in the AtDREB1C transgenic plants for further analysis (Supplementary Table S6). Quantitative RT-PCR was performed on a set of 30 genes (18 significantly upregulated genes and 12 significantly downregulated genes) selected at random from among the significant DEGs in the transgenic lines compared to WT plants following drought treatment (Supplementary Table S6). The relative expression levels of these genes were similar to the expression profiles determined from the respective RNA-seq data (Supplementary Figure S1). A high correlation (R2 > 0.98) was found between the qRT-PCR and RNA-seq results (Figure 12), confirming the accuracy of the RNA-seq data.
FIGURE 12. Correlation of fold changes determined by RNA-Seq (x-axis) and qRT-PCR data (y-axis). WT, wild type; p35S::AtDREB1C-04, p35S::AtDREB1C transgenic line 04; pRD29A::AtDREB1C-07, pRD29A::AtDREB1C transgenic line 07.
COG Enrichment and KEGG Pathway Analysis of DEGs
Annotations of DEGs against the COG database showed that more than 17% of DEGs in both comparisons (WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07) could not be annotated accurately and were therefore classified in the cluster ‘general function prediction only.’ Based on the number of annotated genes in each category, ‘transcription’ and ‘signal transduction mechanisms’ were two of the top three categories in the two DEG sets (WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07) according to the COG database (Figures 13A,C), which suggests that DEGs in these groups play important roles in the drought stress response of AtDREB1C transgenic S. miltiorrhiza.
FIGURE 13. COG and KEGG classification of DEGs in two comparisons [WT versus p35S::AtDREB1C-04 (A,B) and WT versus pRD29A::AtDREB1C-07 (C,D)]. WT, wild type; p35S::AtDREB1C-04, p35S::AtDREB1C transgenic line 04; pRD29A::AtDREB1C-07, pRD29A::AtDREB1C transgenic line 07.
Furthermore, we analyzed the regulatory pathways of the DEGs using KEGG analysis to obtain gene functional and genomic information. The top 50 enriched KEGG pathways of DEGs among the WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07 comparisons are shown in Figures 13B,D, respectively. Many DEGs from both comparisons were enriched in the KEGG pathways of photosynthesis, plant hormone signal transduction, ribosome, and phenylpropanoid biosynthesis. Plant hormone regulatory pathways, including auxin, cytokinin, gibberellin, ABA, ethylene, brassinosteroid, jasmonic acid, and salicylic acid pathways (Figure 14), ranked second and fourth among the top 50 KEGG pathways of DEGs identified in the WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07 comparisons, respectively (Figures 13B,D). Except for the cytokinin pathway, other hormone regulatory pathways appear to be involved in AtDREB1C-mediated enhanced drought tolerance in transgenic S. miltiorrhiza (Figure 14).
FIGURE 14. Analysis of DEGs related to plant hormone signal transduction in two comparisons [WT versus p35S::AtDREB1C-04 (A) and WT versus pRD29A::AtDREB1C-07 (B)]. Red/green indicate up-/downregulated DEGs, respectively. Blue indicates DEGs with mixed patterns of regulation. WT, wild type; p35S::AtDREB1C-04, p35S::AtDREB1C transgenic line 04; pRD29A::AtDREB1C-07, pRD29A::AtDREB1C transgenic line 07.
Analysis of auxin, gibberellin, brassinosteroid, jasmonic acid, and salicylic acid KEGG pathways showed that all DEGs were upregulated in both p35S::AtDREB1C and pRD29A::AtDREB1C transgenic plants (Figure 14). The genes encoding ERF1/2, a pivotal regulator in the ethylene pathway, exhibited mixed expression patterns in both comparisons (WT versus p35S::AtDREB1C-04 and WT versus pRD29A::AtDREB1C-07). The genes for the core components for ABA signaling, PYR/PYL, PP2C, and SnRK2, which function as ABA receptors, repressors, and positive regulators, exhibited various patterns of expression; PP2C was downregulated and SnRK2 was upregulated in the p35S::AtDREB1C transgenic lines, while PP2C exhibited mixed patterns of expression, and PYR/PYL and SnRK2 were upregulated in pRD29A::AtDREB1C transgenic S. miltiorrhiza (Figure 14).
In addition, we analyzed the DEGs between p35S::AtDREB1C-04 and pRD29A::AtDREB1C-07 transgenic lines after drought stress. According to the COG database, more than 15% of DEGs were classified in ‘secondary metabolites biosynthesis, transport and catabolism,’ which was the largest category with accurate annotation (Supplementary Figure S2A). Furthermore, ‘phenylalanine metabolism’ and ‘phenylpropanoid biosynthesis’ were the top two KEGG pathways of DEGs identified in the comparison of p35S::AtDREB1C-04 versus pRD29A::AtDREB1C -07 (Supplementary Figure S2B). Compared with p35S::AtDREB1C-04, most of the significant DEGs were upregulated in pRD29A::AtDREB1C-07. Meanwhile, almost all the significant DEGs were expressed at very low levels in WT (Supplementary Table S7).
Discussion
Water deficit, a major abiotic stress factor, limits crop yields, thereby affecting the incomes of millions of farmers in arid and semi-arid regions (Ferdous et al., 2015; Todaka et al., 2015; Zhu et al., 2016). Such losses in yields vary greatly depending on the timing, duration, and intensity of the stress. This variability, combined with other location-specific environmental stressors, such as high light and temperature, hamper efforts aimed at producing drought-tolerant plants through conventional breeding (Fang and Xiong, 2015; Thomas, 2015). Compared to traditional breeding and marker-assisted selection, directly introducing stress tolerance-related genes via genetic engineering represents a more effective, rapid approach for increasing abiotic stress tolerance in plants (Patel et al., 2015; Todaka et al., 2015). One important strategy for genetic engineering is to induce a TF gene that regulates several genes involved in abiotic stress (Rehman and Mahmood, 2015). Indeed, overexpressing stress-inducible DREB TF genes activates the expression of many target genes containing DRE elements in their promoters, leading to improved stress tolerance (Morran et al., 2011; Tillett et al., 2012; Guttikonda et al., 2014; Dou et al., 2015; Kidokoro et al., 2015). However, constitutively expressing DREB genes often has various pleiotropic effects on the growth and development of plants, such as growth retardation, dwarfism, and delayed flowering (Sharabi-Schwager et al., 2010; Shavrukov et al., 2015; Yamasaki and Randall, 2016; Wu et al., 2017). Several attempts have been made to overcome the problem of severe growth retardation by reducing the duration of DREB overexpression through the use of stress-inducible promoters (Morran et al., 2011; Datta et al., 2012; Kovalchuk et al., 2013; Ravikumar et al., 2014).
In the present study, we functionally validated the role of AtDREB1C in drought tolerance in S. miltiorrhiza. Like the TFs AtDREB1A and AtDREB1B (Wei et al., 2016a,b), ectopic expression of AtDREB1C led to a substantial improvement in the capacity of S. miltiorrhiza to survive under severe drought stress conditions, suggesting that these three TFs play similar roles in the plant stress response. When considering the effects of TFs on the growth of transgenic plants, we observed diverse phenotypes among transgenic lines expressing AtDREB1A, AtDREB1B, and AtDREB1C driven by the constitutive CaMV35S promoter: under normal conditions, p35S::AtDREB1A transgenic S. miltiorrhiza showed slightly stunted growth (Wei et al., 2016b), p35S::AtDREB1B transgenic lines exhibited normal growth (Wei et al., 2016a), and p35S::AtDREB1C plants exhibited severe dwarfism. Interestingly, no growth inhibition or phenotypic aberrations were observed in transgenic plants harboring any of these genes driven by the stress-inducible RD29A promoter. One possible explanation for these results is that while the results provide no clear evidence for specialized functions for AtDREB1A, AtDREB1B, and AtDREB1C, quantitative differences may exist, with the effects of a TF on a given trait varying among transgenic lines (Gilmour et al., 2004). Meanwhile, the RD29A promoter::DREB1 system has proven to be highly effective for overexpressing DREB1 proteins during exposure to stress (Ravikumar et al., 2014; Ma et al., 2015).
We performed physiological and transcriptome analysis to elucidate the mechanisms involved in the drought resistance of the transgenic plants. Physiological tests suggested that the higher chlorophyll content (Figure 8B), photosynthetic capacity (Figures 9A–C), and SOD, POD, and CAT activity (Figures 7B–D) enhanced the drought stress resistance of transgenic plants compared to WT. Similar effects were observed in other plant species, including tomato (Lee et al., 2003), sugarcane (Reis et al., 2014), and Lolium perenne (Li et al., 2011). Transcriptome analysis has been used to elucidate the gene regulatory networks involved in abiotic stress responses in several plant species, such as A. thaliana (Mao et al., 2016), rice (Oono et al., 2014), potato (Zhang N. et al., 2014), Brassica juncea (Bhardwaj et al., 2015), and Robinia pseudoacacia L. ‘Idaho’ (Xiu et al., 2015). Transcriptome analysis of S. miltiorrhiza identified a number of DEGs between the AtDREB1C transgenic lines and WT, suggesting that their corresponding proteins play specific roles in the drought stress response. COG analysis indicated that DEGs involved in ‘transcription,’ ‘signal transduction mechanisms,’ ‘carbohydrate transport and metabolism,’ and ‘replication, recombination and repair’ play important roles in drought tolerance (Figures 13A,C). The main pathways regulated by AtDREB1C, according to KEGG enrichment analysis, include photosynthesis, plant hormone signal transduction, phenylpropanoid biosynthesis, ribosome, starch and sucrose metabolism, and plant–pathogen interactions (Figures 13B,D).
In our previous study, global expression profiling was compared for WT and pRD29A::AtDREB1A-31 transgenic line after drought stress (Wei et al., 2016b). Similarly, numerous genes involved in signal transduction, transcriptional activation, photosynthesis, protein protection, carbohydrate transport and metabolism were significantly upregulated in pRD29A::AtDREB1A-31 plants. Further analysis showed that, the top 10 categories of COG enrichment of DEGs were almost the same between AtDREB1A and AtDREB1C transcriptome. Besides, photosynthesis, plant hormone signal transduction, ribosome, and plant–pathogen interactions were also the main pathways regulated by AtDREB1A TF. Furthermore, qRT-PCR analysis revealed that the expression of 30 upregulated genes, which selected from AtDREB1A and AtDREB1C transcriptome, were increased significantly in AtDREB1B transgenic plants (Wei et al., 2016a; Wei et al., unpublished), indicating that AtDREB1-type TFs may regulate the same downstream genes under drought stress (Gilmour et al., 2004).
Phytohormones are major regulators of plant development, growth, and physiological responses, and they mediate environmental stress responses by integrating environmental stimuli and regulatory networks (Jaillais and Chory, 2010; Depuydt and Hardtke, 2011; Xiu et al., 2015). Manipulating plant hormone metabolism and signaling processes represents a promising strategy for increasing plant stress tolerance (Cabello et al., 2014). ABA is a key regulator of plant responses to environmental stress, particularly osmotic stress (Seiler et al., 2011; Jones, 2016). The essential components of ABA signaling have recently been identified, and their mode of action has been clarified. In the current model, ABA signaling involves three core components: receptors (PYR/PYL/RCAR), protein kinases (SnRK2/OST1), and protein phosphatases (PP2C), (Cramer et al., 2011; Osakabe et al., 2014; Kong et al., 2015; Chen et al., 2016). Overexpressing FpDREB2A in R. pseudoacacia increases drought stress resistance. Analysis of the root transcriptome of FpDREB2A transgenic line revealed that PYL, PP2C, and SnRK2 expressed in an interval of down- and up-regulation, which suggested that the transgenic plants regulated ABA signal transduction and supposed to be involve in drought stress responses (Xiu et al., 2015). In this study, we found that ABA signaling pathways were altered in transgenic S. miltiorrhiza compared to WT (Figure 14). The up- or down-regulation of the core components of ABA signaling in the transgenic plants may give the plants the flexibility needed to adapt to drought stress, although the exact mechanisms remain unclear.
Auxin also plays a critical role in plant drought stress responses (Depuydt and Hardtke, 2011). The expression levels of numerous auxin-related genes are altered under dehydration stress (Jain and Khurana, 2009). The induction of YUCCA7 (YUC7), encoding a flavin monooxygenase that functions in the tryptophan-dependent auxin biosynthetic pathway, increases endogenous auxin levels and enhances drought resistance in Arabidopsis (Du et al., 2012; Lee and Luan, 2012). Auxin alters the expression of numerous genes, including primary auxin-responsive genes of the Aux/IAA, GH3, and SAUR gene families. Aux/IAA gene family members are implicated in the light-induced regulation of auxin responses. Although to date, no SAUR genes have been functionally characterized, SAUR proteins bind to calcium/calmodulin, suggesting that they play a role in auxin signaling via calcium ions (Park et al., 2007). GH3s, which are members of the acyl adenylate-forming firefly luciferase superfamily, catalyze the adenylation of specific substrates (Wang et al., 2008). GH3-mediated auxin homeostasis is an essential component of the complex auxin regulatory network underlying stress adaptation in Arabidopsis (Park et al., 2007). In the present study, we identified upregulated DEGs in the auxin pathway, including GH3 and SAUR family members (Figure 14), suggesting that these two groups of auxin-responsive genes play specific roles in AtDREB1C-mediated drought tolerance in S. miltiorrhiza.
Other hormones, including gibberellin, ethylene, brassinosteroid, jasmonic acid, and salicylic acid, also play direct or indirect roles in the abiotic stress response (Peleg and Blumwald, 2011; Cabello et al., 2014; Zhou et al., 2014; Chen et al., 2016). For example, ethylene may modulate the expression of genes considered to be effectors of ethylene signaling (Shi et al., 2012). Environmental stress induces ethylene accumulation, which increases plant survival under adverse conditions. Recent evidence indicates that phytohormones function in multiple processes. Crosstalk between various plant hormone regulatory pathways results in synergetic or antagonistic interactions, which play crucial roles in abiotic stress responses (Depuydt and Hardtke, 2011; Peleg and Blumwald, 2011; Song et al., 2014; Shi et al., 2015; Chen et al., 2016). Plant growth retardation often occurs when plant hormone homeostasis is altered. Excess ABA biosynthesis may deplete the pool of precursors needed for chlorophyll biosynthesis, thereby leading to growth retardation (Sreenivasulu et al., 2012). GH3-8 is an IAA–amido synthetase that alters auxin homeostasis; overexpressing GH3-8 leads to abnormal plant morphology, as well as retarded growth and development, in transgenic rice (Ding et al., 2008). Therefore, we speculate that constitutive expression of AtDREB1C alters plant hormone homeostasis, especially at the vegetative stage, which might be the main cause of the dwarf phenotypes of the p35S::AtDREB1C transgenic lines.
Taken together, our findings indicate that AtDREB1C plays important roles in the drought resistance of transgenic S. miltiorrhiza. The improved drought tolerance of these plants was associated with improved photosynthesis and antioxidant defense systems. Analysis of the transcriptome of AtDREB1C transgenic plants showed changes in the expression of genes involved in pathways including photosynthesis, plant hormone signal transduction, phenylpropanoid biosynthesis, ribosome, starch and sucrose metabolism, and plant–pathogen interactions. In particular, the differential expression of genes involved in plant hormone signaling are thought to be one of the main causes of increased drought tolerance in the AtDREB1C transgenic plants. Further investigation of the effects of the AtDREB1C TF on the crosstalk between different plant hormones, including biosynthesis and signaling, as well as their roles in plant growth and responses to abiotic stresses, will help elucidate the mechanism underlying the effects of AtDREB1C-mediated plant hormone homeostasis on plant growth and stress responses.
Author Contributions
TW, KD, CC, and YZ designed the experiments, analyzed the data and wrote the manuscript. TW and KD performed the main experiments in this study. QZ, YG, and YL contributed reagents, materials, and helped in the stress experiments. QZ, YG, MY, and LZ helped with the detection of photosynthesis and enzyme activities. XZ, CW, and ZL contributed to data analyses and discussion. All authors read and approved the final 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.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (31371682, 31371249), the Fundamental Research Funds for the Central Universities (ZYGX2014J081) and the National Natural Science Foundation of Tianjin (13JCZDJC29000).
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.00052/full#supplementary-material
References
Bhardwaj, A. R., Joshi, G., Kukreja, B., Malik, V., Arora, P., Pandey, R., et al. (2015). Global insights into high temperature and drought stress regulated genes by RNA-Seq in economically important oilseed crop Brassica juncea. BMC Plant Biol. 15:9. doi: 10.1186/s12870-014-0405-1
Bihmidine, S., Lin, J., Stone, J. M., Awada, T., Specht, J. E., and Clemente, T. E. (2013). Activity of the Arabidopsis RD29A and RD29B promoter elements in soybean under water stress. Planta 237, 55–64. doi: 10.1007/s00425-012-1740-9
Byun, M. Y., Lee, J., Cui, L. H., Kang, Y., Oh, T. K., Park, H., et al. (2015). Constitutive expression of DaCBF7, an Antarctic vascular plant Deschampsia antarctica CBF homolog, resulted in improved cold tolerance in transgenic rice plants. Plant Sci. 236, 61–74. doi: 10.1016/j.plantsci.2015.03.020
Cabello, J. V., Lodeyro, A. F., and Zurbriggen, M. D. (2014). Novel perspectives for the engineering of abiotic stress tolerance in plants. Curr. Opin. Biotech. 26, 62–70. doi: 10.1016/j.copbio.2013.09.011
Chen, H. Y., Hsieh, E. J., Cheng, M. C., Chen, C. Y., Hwang, S. Y., and Lin, T. P. (2016). ORA47 (octadecanoid-responsive AP2/ERF-domain transcription factor 47) regulates jasmonic acid and abscisic acid biosynthesis and signaling through binding to a novel cis-element. New Phytol. 211, 599–613. doi: 10.1111/nph.13914
Cramer, G. R., Urano, K., Delrot, S., Pezzotti, M., and Shinozaki, K. (2011). Effects of abiotic stress on plants: a systems biology perspective. BMC Plant Biol. 11:163. doi: 10.1186/1471-2229-11-163
Cui, G., Duan, L., Jin, B., Qian, J., Xue, Z., Shen, G., et al. (2015). Functional divergence of diterpene 18 syntheses in the medicinal plant Salvia miltiorrhiza Bunge. Plant Physiol. 169, 1607–1618. doi: 10.1104/pp.15.00695
Datta, K., Baisakh, N., Ganguly, M., Krishnan, S., Yamaguchi Shinozaki, K., and Datta, S. K. (2012). Overexpression of Arabidopsis and Rice stress genes’ inducible transcription factor confers drought and salinity tolerance to rice. Plant Biotechnol. J. 10, 579–586. doi: 10.1111/j.1467-7652.2012.00688.x
de Paiva Rolla, A. A., de Fátima Corrêa Carvalho, J., Fuganti-Pagliarini, R., Engels, C., do Rio, A., Marin, S. R., et al. (2014). Phenotyping soybean plants transformed with rd29A: AtDREB1A for drought tolerance in the greenhouse and field. Transgenic Res. 23, 75–87. doi: 10.1007/s11248-013-9723-6
Depuydt, S., and Hardtke, C. S. (2011). Hormone Signalling Crosstalk in Plant Growth Regulation. Curr. Biol. 21, 365–373. doi: 10.1016/j.cub.2011.03.013
Ding, X., Cao, Y., Huang, L., Zhao, J., Xu, C., Li, X., et al. (2008). Activation of the Indole-3-Acetic Acid-Amido Synthetase GH3-8 Suppresses Expansin Expression and Promotes Salicylate- and Jasmonate-Independent Basal Immunity in Rice. Plant Cell 20, 228–240. doi: 10.1105/tpc.107.055657
Dou, H., Xv, K., Meng, Q., Li, G., and Yang, X. (2015). Potato plants ectopically expressing Arabidopsis thaliana CBF3 exhibit enhanced tolerance to high-temperature stress. Plant Cell Environ. 38, 61–72. doi: 10.1111/pce.12366
Du, H., Wu, N., Fu, J., Wang, S., Li, X., Xiao, J., et al. (2012). A GH3 family member, OsGH3-2, modulates auxin and abscisic acid levels and differentially affects drought and cold tolerance in rice. J. Exp. Bot. 63, 6467–6480. doi: 10.1093/jxb/ers300
Fang, Y., and Xiong, L. (2015). General mechanisms of drought response and their application in drought resistance improvement in plants. Cell. Mol. Life Sci. 72, 673–689. doi: 10.1007/s00018-014-1767-0
Ferdous, J., Hussain, S. S., and Shi, B. J. (2015). Role of microRNAs in plant drought tolerance. Plant Biotechnol. J. 13, 293–305. doi: 10.1111/pbi.12318
Gilmour, S. J., Fowler, S. G., and Thomashow, M. F. (2004). Arabidopsis transcriptional activators CBF1, CBF2, and CBF3 have matching functional activities. Plant Mol. Biol. 54, 767–781.
Guttikonda, S. K., Valliyodan, B., Neelakandan, A. K., Tran, L. S. P., Kumar, R., Quach, T. N., et al. (2014). Overexpression of AtDREB1D transcription factor improves drought tolerance in soybean. Mol. Biol. Rep. 41, 7995–8008. doi: 10.1007/s11033-014-3695-3
Hichri, I., Muhovski, Y., Clippe, A., Žižková, E., Dobrev, P. I., Motyka, V., et al. (2015). SlDREB2, a tomato dehydration-responsive element-binding 2 transcription factor, mediates salt stress tolerance in tomato and Arabidopsis. Plant Cell Environ. 39, 1–17. doi: 10.1111/pce.12591
Jaillais, Y., and Chory, J. (2010). Unraveling the paradoxes of plant hormone signaling integration. Nat. Struct. Mol. Biol. 17, 642–645. doi: 10.1038/nsmb0610-642
Jain, M., and Khurana, J. P. (2009). Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. FEBS J. 276, 3148–3162. doi: 10.1111/j.1742-4658.2009.07033.x
Jia, Y., Ding, Y., Shi, Y., Zhang, X., Gong, Z., and Yang, S. (2016). The cbfs triple mutants reveal the essential functions of CBFs in cold acclimation and allow the definition of CBF regulons in Arabidopsis. New Phytol. 212, 345–353. doi: 10.1111/nph.14088
Jones, A. M. (2016). A new look at stress: abscisic acid patterns and dynamics at high-resolution. New Phytol. 210, 38–44. doi: 10.1111/nph.13552
Kang, J., Zhang, H., Sun, T., Shi, Y., Wang, J., Zhang, B., et al. (2013). Natural variation of C-repeat-binding factor (CBFs) genes is a major cause of divergence in freezing tolerance among a group of Arabidopsis thaliana populations along the Yangtze River in China. New Phytol. 199, 1069–1080. doi: 10.1111/nph.12335
Kidokoro, S., Watanabe, K., Ohori, T., Moriwaki, T., Maruyama, K., Mizoi, J., et al. (2015). Soybean DREB1/CBF-type transcription factors function in heat and drought as well as cold stress-responsive gene expression. Plant J. 81, 505–518. doi: 10.1111/tpj.12746
Kong, L., Cheng, J., Zhu, Y., Ding, Y., Meng, J., Chen, Z., et al. (2015). Degradation of the ABA co-receptor ABI1 by PUB12/13 U-box E3 ligases. Nat. Commun. 6:8630. doi: 10.1038/ncomms9630
Kovalchuk, N., Jia, W., Eini, O., Morran, S., Pyvovarenko, T., Fletcher, S., et al. (2013). Optimization of TaDREB3 gene expression in transgenic barley using cold-inducible promoters. Plant Biotechnol. J. 11, 659–670. doi: 10.1111/pbi.12056
Langmead, B., Trapnell, C., Pop, M., and Salzberg, S. L. (2009). Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 10:R25. doi: 10.1186/gb-2009-10-3-r25
Lata, C., and Prasad, M. (2011). Role of DREBs in regulation of abiotic stress responses in plants. J. Exp. Bot. 62, 4731–4748. doi: 10.1093/jxb/err210
Lee, J. T., Prasad, V., Yang, P. T., Wu, J. F., David Ho, T. H., Charng, Y. Y., et al. (2003). Expression of Arabidopsis CBF1 regulated by an ABA/stress inducible promoter in transgenic tomato confers stress tolerance without affecting yield. Plant Cell Environ. 26, 1181–1190. doi: 10.1046/j.1365-3040.2003.01048.x
Lee, S. C., and Luan, S. (2012). ABA signal transduction at the crossroad of biotic and abiotic stress responses. Plant Cell Environ. 35, 53–60. doi: 10.1111/j.1365-3040.2011.02426.x
Li, B., and Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12:323. doi: 10.1186/1471-2105-12-323
Li, J., Sima, W., Ouyang, B., Wang, T., Ziaf, K., Luo, Z., et al. (2012). Tomato SlDREB gene restricts leaf expansion and internode elongation by downregulating key genes for gibberellin biosynthesis. J. Exp. Bot. 63, 6407–6420. doi: 10.1093/jxb/ers295
Li, X., Cheng, X., Liu, J., Zeng, H., Han, L., and Tang, W. (2011). Heterologous expression of the Arabidopsis DREB1A/CBF3 gene enhances drought and freezing tolerance in transgenic Lolium perenne plants. Plant Biotechnol. Rep. 5, 61–69. doi: 10.1007/s11816-010-0157-9
Licausi, F., Ohme-Takagi, M., and Perata, P. (2013). APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors: mediators of stress responses and developmental programs. New Phytol. 199, 639–649. doi: 10.1111/nph.12291
Luo, J., Song, W., Yang, G., Xu, H., and Chen, K. (2015). Compound Danshen (Salvia miltiorrhiza) dripping pill for coronary heart disease: an overview of systematic reviews. Am. J. Chin. Med. 43, 25–43. doi: 10.1142/S0192415X15500020
Ma, J. T., Yin, C. C., Guo, Q. Q., Zhou, M. L., Wang, Z. L., and Wu, Y. M. (2015). A novel DREB transcription factor from Halimodendron halodendron leads to enhance drought and salt tolerance in Arabidopsis. Biol. Plant. 59, 74–82. doi: 10.1007/s10535-014-0467-9
Mao, H., Yu, L., Han, R., Li, Z., and Liu, H. (2016). ZmNAC55, a maize stress-responsive NAC transcription factor, confers drought resistance in transgenic Arabidopsis. Plant Physiol. Biochem. 105, 55–66. doi: 10.1016/j.plaphy.2016.04.018
Mendiondo, G. M., Gibbs, D. J., Szurman-Zubrzycka, M., Korn, A., Marquez, J., Szarejko, I., et al. (2016). Enhanced waterlogging tolerance in barley by manipulation of expression of the N-end rule pathway E3 ligase PROTEOLYSIS6. Plant Biotechnol. J. 14, 40–50. doi: 10.1111/pbi.12334
Morran, S., Eini, O., Pyvovarenko, T., Parent, B., Singh, R., Ismagul, A., et al. (2011). Improvement of stress tolerance of wheat and barley by modulation of expression of DREB/CBF factors. Plant Biotechnol. J. 9, 230–249. doi: 10.1111/j.1467-7652.2010.00547.x
Oono, Y., Yazawa, T., Kawahara, Y., Kanamori, H., Kobayashi, F., Sasaki, H., et al. (2014). Genome-wide transcriptome analysis reveals that cadmium stress signaling controls the expression of genes in drought stress signal pathways in rice. PLoS ONE 9:e96946. doi: 10.1371/journal.pone.0096946
Osakabe, Y., Yamaguchi-Shinozaki, K., Shinozaki, K., and Tran, L. S. P. (2014). ABA control of plant macroelement membrane transport systems in response to water deficit and high salinity. New Phytol. 202, 35–49. doi: 10.1111/nph.12613
Park, J. E., Park, J. Y., Kim, Y. S., Staswick, P. E., Jeon, J., Yun, J., et al. (2007). GH3-mediated auxin homeostasis links growth regulation with stress adaptation response in Arabidopsis. J. Biol. Chem. 282, 10036–10046. doi: 10.1074/jbc.M610524200
Park, S., Lee, C.-M., Doherty, C. J., Gilmour, S. J., Kim, Y., and Thomashow, M. F. (2015). Regulation of the Arabidopsis CBF regulon by a complex low-temperature regulatory network. Plant J. 82, 193–207. doi: 10.1111/tpj.12796
Patel, M., Milla-Lewis, S., Zhang, W., Templeton, K., Reynolds, W. C., Richardson, K., et al. (2015). Overexpression of ubiquitin-like LpHUB1 gene confers drought tolerance in perennial ryegrass. Plant Biotechnol. J. 13, 689–699. doi: 10.1111/pbi.12291
Peleg, Z., and Blumwald, E. (2011). Hormone balance and abiotic stress tolerance in crop plants. Curr. Opin. Plant Biol. 14, 290–295. doi: 10.1016/j.pbi.2011.02.001
Rai, G. K., Rai, N. P., Rathaur, S., Kumar, S., and Singh, M. (2013). Expression of rd29A::AtDREB1A/CBF3 in tomato alleviates drought-induced oxidative stress by regulating key enzymatic and non-enzymatic antioxidants. Plant Physiol. Biochem. 69, 90–100. doi: 10.1016/j.plaphy.2013.05.002
Ravikumar, G., Manimaran, P., Voleti, S., Subrahmanyam, D., Sundaram, R., Bansal, K., et al. (2014). Stress-inducible expression of AtDREB1A transcription factor greatly improves drought stress tolerance in transgenic indica rice. Transgenic Res. 23, 421–439. doi: 10.1007/s11248-013-9776-6
Rehman, S., and Mahmood, T. (2015). Functional role of DREB and ERF transcription factors: regulating stress-responsive network in plants. Acta Physiol. Plant. 37, 1–14. doi: 10.1007/s11738-015-1929-1
Reis, R. R., da Cunha, B. A., Martins, P. K., Martins, M. T., Alekcevetch, J. C., Chalfun, A. Jr., et al. (2014). Induced over-expression of AtDREB2A CA improves drought tolerance in sugarcane. Plant Sci. 22, 59–68. doi: 10.1016/j.plantsci.2014.02.003
Seiler, C., Harshavardhan, V. T., Rajesh, K., Reddy, P. S., Strickert, M., Rolletschek, H., et al. (2011). ABA biosynthesis and degradation contributing to ABA homeostasis during barley seed development under control and terminal drought-stress conditions. J. Exp. Bot. 62, 2615–2632. doi: 10.1093/jxb/erq446
Shan, D. P., Huang, J. G., Yang, Y. T., Guo, Y. H., Wu, C. A., Yang, G. D., et al. (2007). Cotton GhDREB1 increases plant tolerance to low temperature and is negatively regulated by gibberellic acid. New Phytol. 176, 70–81. doi: 10.1111/j.1469-8137.2007.02160.x
Sharabi-Schwager, M., Lers, A., Samach, A., Guy, C. L., and Porat, R. (2010). Overexpression of the CBF2 transcriptional activator in Arabidopsis delays leaf senescence and extends plant longevity. J. Exp. Bot. 61, 261–273. doi: 10.1093/jxb/erp300
Shavrukov, Y., Baho, M., Lopato, S., and Langridge, P. (2015). The TaDREB3 transgene transferred by conventional crossings to different genetic backgrounds of bread wheat improves drought tolerance. Plant Biotechnol. J. 14, 313–322. doi: 10.1111/pbi.12385
Shi, Y., Ding, Y., and Yang, S. (2015). Cold signal transduction and its interplay with phytohormones during cold acclimation. Plant Cell Physiol. 56, 7–15. doi: 10.1093/pcp/pcu115
Shi, Y., Tian, S., Hou, L., Huang, X., Zhang, X., Guo, H., et al. (2012). Ethylene signaling negatively regulates freezing tolerance by repressing expression of CBF and type-A ARR genes in Arabidopsis. Plant Cell 24, 2578–2595. doi: 10.1105/tpc.112.098640
Song, S., Huang, H., Gao, H., Wang, J., Wu, D., Liu, X., et al. (2014). Interaction between MYC2 and ETHYLENE INSENSITIVE3 modulates antagonism between jasmonate and ethylene signaling in Arabidopsis. Plant Cell 26, 263–279. doi: 10.1105/tpc.113.120394
Sreenivasulu, N., Harshavardhan, V. T., Govind, G., Seiler, C., and Kohli, A. (2012). Contrapuntal role of ABA: does it mediate stress tolerance or plant growth retardation under long-term drought stress? Gene 506, 265–273. doi: 10.1016/j.gene.2012.06.076
Storani, L., Hernando, C. E., Staneloni, R. J., Ploschuk, E., Rugnone, M. L., Striker, G. G., et al. (2015). AtCBF1 overexpression confers tolerance to high light conditions at warm temperatures in potato plants. Am. J. Potato Res. 92, 619–635. doi: 10.1007/s12230-015-9476-2
Thomas, W. (2015). Drought-resistant cereals: impact on water sustainability and nutritional quality. Proc. Nutr. Soc. 74, 191–197. doi: 10.1017/S0029665115000026
Tillett, R. L., Wheatley, M. D., Tattersall, E. A., Schlauch, K. A., Cramer, G. R., and Cushman, J. C. (2012). The Vitis vinifera C-repeat binding protein 4 (VvCBF4) transcriptional factor enhances freezing tolerance in wine grape. Plant Biotechnol. J. 10, 105–124. doi: 10.1111/j.1467-7652.2011.00648.x
Todaka, D., Shinozaki, K., and Yamaguchi-Shinozaki, K. (2015). Recent advances in the dissection of drought-stress regulatory networks and strategies for development of drought-tolerant transgenic rice plants. Front. Plant Sci. 6:84. doi: 10.3389/fpls.2015.00084
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Baren, M. J. V., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Vahala, J., Felten, J., Love, J., Gorzsás, A., Gerber, L., Lamminmäki, A., et al. (2013). A genome-wide screen for ethylene-induced Ethylene Response Factors (ERFs) in hybrid aspen stem identifies ERF genes that modify stem growth and wood properties. New Phytol. 200, 511–522. doi: 10.1111/nph.12386
Wang, C., Liu, S., Dong, Y., Zhao, Y., Geng, A., Xia, X., et al. (2016). PdEPF1 regulates water-use efficiency and drought tolerance by modulating stomatal density in poplar. Plant Biotechnol. J. 14, 849–860. doi: 10.1111/pbi.12434
Wang, D., Yao, W., Song, Y., Liu, W., and Wang, Z. (2012). Molecular characterization and expression of three galactinol synthase genes that confer stress tolerance in Salvia miltiorrhiza. J. Plant Physiol. 169, 1838–1848. doi: 10.1016/j.jplph.2012.07.015
Wang, H., Tian, C., Duan, J., and Wu, K. (2008). Research progresses on GH3s, one family of primary auxin-responsive genes. Plant Growth Regul. 56, 225–232. doi: 10.1007/s10725-008-9313-4
Wei, T., Deng, K., Gao, Y., Liu, Y., Yang, M., Zhang, L., et al. (2016a). Arabidopsis DREB1B in transgenic Salvia miltiorrhiza increased tolerance to drought stress without stunting growth. Plant Physiol. Biochem. 104, 17–28. doi: 10.1016/j.plaphy.2016.03.003
Wei, T., Deng, K., Liu, D., Gao, Y., Liu, Y., Yang, M., et al. (2016b). Ectopic expression of DREB transcription factor, AtDREB1A, confers tolerance to drought in transgenic Salvia miltiorrhiza. Plant Cell Physiol. 57, 1593–1609. doi: 10.1093/pcp/pcw084
Wu, J., Folta, K. M., Xie, Y., Jiang, W., Lu, J., and Zhang, Y. (2017). Overexpression of Muscadinia rotundifolia CBF2 gene enhances biotic and abiotic stress tolerance in Arabidopsis. Protoplasma 254, 239–251. doi: 10.1007/s00709-015-0939-6
Wu, Y., Liu, C., Kuang, J., Ge, Q., Zhang, Y., and Wang, Z. (2014). Overexpression of SmLEA enhances salt and drought tolerance in Escherichia coli and Salvia miltiorrhiza. Protoplasma 251, 1191–1199. doi: 10.1007/s00709-014-0626-z
Xiu, Y., Iqbal, A., Zhu, C., Wu, G., Chang, Y., Li, N., et al. (2015). Improvement and transcriptome analysis of root architecture by overexpression of Fraxinus pennsylvanica DREB2A transcription factor in Robinia pseudoacacia L.‘Idaho.’ Plant Biotechnol. J. 14, 1456–1469. doi: 10.1111/pbi.12509
Xu, Z., Peters, R. J., Weirather, J., Luo, H., Liao, B., Zhang, X., et al. (2015). Full-length transcriptome sequences and splice variants obtained by a combination of sequencing platforms applied to different root tissues of Salvia miltiorrhiza and tanshinone biosynthesis. Plant J. 82, 951–961. doi: 10.1111/tpj.12865
Yamasaki, K., Kigawa, T., Seki, M., Shinozaki, K., and Yokoyama, S. (2013). DNA-binding domains of plant-specific transcription factors: structure, function, and evolution. Trends Plant Sci. 18, 267–276. doi: 10.1016/j.tplants.2012.09.001
Yamasaki, Y., and Randall, S. K. (2016). Functionality of soybean CBF/DREB1 transcription factors. Plant Sci. 246, 80–90. doi: 10.1016/j.plantsci.2016.02.007
Zeppel, M. J., Harrison, S. P., Adams, H. D., Kelley, D. I., Li, G., Tissue, D. T., et al. (2014). Drought and resprouting plants. New Phytol. 206, 583–589. doi: 10.1111/nph.13205
Zhang, N., Liu, B., Ma, C., Zhang, G., Chang, J., Si, H., et al. (2014). Transcriptome characterization and sequencing-based identification of drought-responsive genes in potato. Mol. Biol. Rep. 41, 505–517. doi: 10.1007/s11033-013-2886-7
Zhang, Y., Yan, Y. P., Wu, Y. C., Hua, W. P., Chen, C., Ge, Q., et al. (2014). Pathway engineering for phenolic acid accumulations in Salvia miltiorrhiza by combinational genetic manipulation. Metab. Eng. 21, 71–80. doi: 10.1016/j.ymben.2013.10.009
Zhang, X., Dong, J., Liu, H., Wang, J., Qi, Y., and Liang, Z. (2016). Transcriptome sequencing in response to salicylic acid in Salvia miltiorrhiza. PLoS ONE 11:e0147849. doi: 10.1371/journal.pone.0147849
Zhao, S., Zhang, J., Tan, R., Yang, L., and Zheng, X. (2015). Enhancing diterpenoid concentration in Salvia miltiorrhiza hairy roots through pathway engineering with maize C1 transcription factor. J. Exp. Bot. 66, 7211–7226. doi: 10.1093/jxb/erv418
Zhou, M., Xu, M., Wu, L., Shen, C., Ma, H., and Lin, J. (2014). CbCBF from Capsella bursa-pastoris enhances cold tolerance and restrains growth in Nicotiana tabacum by antagonizing with gibberellin and affecting cell cycle signaling. Plant Mol. Biol. 85, 259–275. doi: 10.1007/s11103-014-0181-1
Zhou, Y., Sun, W., Chen, J., Tan, H., Xiao, Y., Li, Q., et al. (2016). SmMYC2a and SmMYC2b played similar but irreplaceable roles in regulating the biosynthesis of tanshinones and phenolic acids in Salvia miltiorrhiza. Sci. Rep. 6:22852. doi: 10.1038/srep22852
Keywords: antioxidative metabolism, AtDREB1C, drought, photosynthetic capacity, Salvia miltiorrhiza, transcriptome
Citation: Wei T, Deng K, Zhang Q, Gao Y, Liu Y, Yang M, Zhang L, Zheng X, Wang C, Liu Z, Chen C and Zhang Y (2017) Modulating AtDREB1C Expression Improves Drought Tolerance in Salvia miltiorrhiza. Front. Plant Sci. 8:52. doi: 10.3389/fpls.2017.00052
Received: 30 September 2016; Accepted: 10 January 2017;
Published: 24 January 2017.
Edited by:
Agnieszka Ludwików, Adam Mickiewicz University in Poznañ, PolandReviewed by:
Prakash Lakshmanan, Sugar Research Australia, AustraliaDominique Job, Centre National de la Recherche Scientifique (CNRS), France
Copyright © 2017 Wei, Deng, Zhang, Gao, Liu, Yang, Zhang, Zheng, Wang, Liu, Chen and Zhang. 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: Chengbin Chen, Y2hlbmNiQG5hbmthaS5lZHUuY24= Yong Zhang, emhhbmd5b25nOTE2QHVlc3RjLmVkdS5jbg==
†These authors have contributed equally to this work.