Next Article in Journal
The Antimicrobial Activity of Annona emarginata (Schltdl.) H. Rainer and Most Active Isolated Compounds against Clinically Important Bacteria
Previous Article in Journal
Study on the Polar Extracts of Dendrobium nobile, D. officinale, D. loddigesii, and Flickingeria fimbriata: Metabolite Identification, Content Evaluation, and Bioactivity Assay
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 23
Issue 5
10.3390/molecules23051186
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analysis in Chinese Cabbage (Brassica rapa ssp. pekinensis) Provides the Role of Glucosinolate Metabolism in Response to Drought Stress

by
Seung Hee Eom
1,
Seung-A Baek
2,
Jae Kwang Kim
2 and
Tae Kyung Hyun
1,*
1
Department of Industrial Plant Science and Technology, College of Agricultural, Life and Environmental Sciences, Chungbuk National University, Cheongju 28644, Korea
2
Division of Life Sciences and Convergence Research Center for Insect Vectors, Incheon National University, Incheon 22012, Korea
*
Author to whom correspondence should be addressed.
Molecules 2018, 23(5), 1186; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules23051186
Submission received: 15 April 2018 / Revised: 12 May 2018 / Accepted: 12 May 2018 / Published: 15 May 2018
(This article belongs to the Section Chemical Biology)
Browse Figures
Versions Notes

Abstract

:
Although drought stress is one of the most limiting factors in growth and production of Chinese cabbage (Brassica rapa L. ssp. pekinensis), the underlying biochemical and molecular causes are poorly understood. In the present study, to address the mechanisms underlying the drought responses, we analyzed the transcriptome profile of Chinese cabbage grown under drought conditions. Drought stress transcriptionally activated several transcription factor genes, including AP2/ERFs, bHLHs, NACs and bZIPs, and was found to possibly result in transcriptional variation in genes involved in organic substance metabolic processes. In addition, comparative expression analysis of selected BrbZIPs under different stress conditions suggested that drought-induced BrbZIPs are important for improving drought tolerance. Further, drought stress in Chinese cabbage caused differential acclimation responses in glucosinolate metabolism in leaves and roots. Analysis of stomatal aperture indicated that drought-induced accumulation of glucosinolates in leaves directly or indirectly controlled stomatal closure to prevent water loss, suggesting that organ-specific responses are essential for plant survival under drought stress condition. Taken together, our results provide information important for further studies on molecular mechanisms of drought tolerance in Chinese cabbage.

1. Introduction

Because the increasing world population and worldwide climate change affects agriculture in several ways, an intergovernmental panel on climate change has concluded that increased concentrations of greenhouse gases will lead to dry conditions in the subtropics, creating widespread drought stress in agricultural regions [1]. Drought is considered as the dominant factor reducing crop growth and yield, and is expected to cause serious problems in plant growth and crop production over more than 50% of all agricultural lands by 2050 [2]. In addition, it has been shown that drought sensitivity in crops has increased over the past twenty years [3]. Therefore, improving drought tolerance is a major goal for crop breeders. Under drought conditions, plants display various physiological and biochemical responses at the cellular and whole-plant levels, leading to a range of specific and nonspecific phenotypes [4]. These responses include stomatal movement, repression of cell growth, development, and photosynthesis, and alteration in biosynthetic pathways, antioxidant pathways, and the respiration pathway [5], suggesting that drought tolerance is an outcome of a series of molecular, cellular, physiological, and biochemical processes mediated via induction and/or repression of genes and their regulation through complex transcriptional networks [6]. This indicates that particular attention to drought-stress response genes and drought-stress-induced transcriptional networks will be required for successful yield protection against drought.
Chinese cabbage (Brassica rapa L. ssp. pekinensis) is a subspecies of B. rapa and is considered an economically important cruciferous vegetable in Asia, particularly in China, Korea, and Japan. Given its significant economic value and its genetically close relationship with Arabidopsis thaliana, Chinese cabbage has been used as a model crop for studies on functional genomics of Brassica species and the evolution of polyploid genomes [7]. In addition, high-throughput techniques have been widely used for understanding molecular mechanisms in Chinese cabbage in response to environmental stresses. For example, transcriptome analysis using digital gene expression profiling suggested that genes encoding transcription factors (TFs) including NAC, MYB, HSF (heat shock factor), WRKY, bHLH (basic helix-loop-helix), and ERF (ethylene-responsive factors), antioxidant proteins (superoxide dismutase, peroxidase, catalase, and glutathione S-transferase), and proteins involved in osmolyte synthesis contribute to salt tolerance [8]. Further, comparative transcriptome analysis of different varieties of Chinese cabbage has revealed common and variety-specific responsive transcripts, which can serve as a helpful resource to explore novel candidate genes for improving stress tolerance [9]. Although further functional characterization of these genes will be essential to address how they modulate stress-tolerance, analysis of transcriptome changes help in understanding the molecular basis of plant adaptation to environmental stresses. Morphological, physiological, and biochemical changes in Chinese cabbage under drought conditions have been well-defined and extensively studied [10,11,12], yet little is known about the genome-wide responses of transcripts to drought stress.
Recently, with the development of next-generation sequencing (NGS) technology, RNA-sequencing (RNA-seq) for transcriptome analysis has been widely used to identify differentially expressed genes (DEGs) among different treatment, tissues or growth periods, suggesting that RNA-seq is powerful tool to obtain an overall view of gene expression profiles [13]. In the present study, to investigate the molecular basis of response to drought stress in Chinese cabbage, the transcriptomes of leaves and roots under drought conditions were analyzed using digital gene expression profiling. Drought stress-inducible or -repressible genes were identified, and further classified as common or specifically regulated. In addition, comparative analysis of DEGs in leaves and roots revealed the importance of glucosinolate metabolism in controlling the response to drought conditions. Taken together, our results provide an overview of molecular mechanisms triggered by drought stress in plants, and will be helpful in unraveling the basic mechanisms of environmental stress tolerance.

2. Results and Discussion

2.1. Physiological Response to Drought Stress in Chinese Cabbage

Drought stress is a major abiotic stress, inducing accelerated production of several reactive oxygen species including superoxide, singlet oxygen, hydroxyl radicals, and H2O2 causing oxidative damage to proteins, DNA, and lipids in different cellular compartments [14]. Therefore, oxidative stress markers including lipid peroxidation in terms of MDA, ROS accumulation, and protein carbonylation have been analyzed to identify variations in physiological response to drought stress [15,16,17]. We analyzed the physiological responses of Chinese cabbage subjected to drought stress, to evaluate efficacy of the treatment. As shown in Figure 1B, drought stress induced significant changes in relative water content. The leaves of the plant showed accumulation of H2O2 (Figure 1C), MDA (Figure 1D), and protein carbonylation (Figure 1E), with their levels increasing 1.9-, 2.7-, and 2.7-fold respectively after four days of drought stress treatment (Stage 3), compared with those in the control plants (Stage 1) (Figure 1A). These physiological changes suggested that the drought stress treatment was effective, and the leaves and roots of Chinese cabbage plants grown under control conditions (Stage 1) and drought stress conditions (Stage 2 and Stage 3) were harvested to unravel the molecular basis of the response.

2.2. RNA-Seq and Identification of Drought-Responsive Genes

During the last decade, several exciting studies have reported the development of new disease-resistant crops using genetic engineering tools, but efforts to generate drought-tolerant crops has been less successful [18]. This might be due to the complex responses of crops to drought stress. In fact, drought stress induces various morphological, physiological, and biochemical changes in crops controlled by numerous small-effect loci and hundreds of genes [19], indicating that knowledge of the biochemical and molecular responses to drought should contribute to improving plant tolerance under water-limited conditions. Therefore, to obtain a global overview of molecular mechanisms involved in drought stress response, total RNA was extracted from roots and leaves of control and drought-treated plants and used for generation of transcriptome libraries on an Illumina HiSeq™ 2500 sequencing platform. After removing low-quality reads, 46 to 54 million clear reads (4.64 to 5.49 Gb) from each sample were acquired for further analysis (Table 1).
All reads were deposited in the National Agricultural Biotechnology Information Center (NABIC, http://nabic.rda.go.kr). To identify the DEGs in the control and drought-treated groups, we constructed six comparison groups, and DEGs were defined as genes with value of |log2 (fold change)| ≥ 1. As shown in Figure 2A, a total of 7535 DEGs (3557 up-regulated and 3978 down-regulated) were detected comparing the L1 and L2 libraries. Comparative analysis of the L2 and L3 libraries showed that the expression of 2666 genes was increased and that of 2608 genes was down-regulated in L3. Upon drought stress exposure, more DEGs were down-regulated in roots than in leaves (Figure 2A), indicating that the response mechanism against drought stress differs between roots and leaves. A total of 1280 genes were exclusively up-regulated in L2, and the expression levels of 1282 genes were increased in L3 but not in L2. In roots, transcripts of 1135 genes were down-regulated solely in R2, whereas 1752 transcripts were down-regulated exclusively in R3 but not in R2 (Figure 2B). Further, 467 DEGs from leaves and 240 DEGs from roots were found to be up-regulated in both tissues, whereas 424 DEGs from leaves and 1181 DEGs from roots were found to be down-regulated in both tissues (Figure 2B). Non-specific lipid-transfer protein 3 (XP_009132002.1), dehydrin Rab18-like (XP_013609020.1), and sucrose synthase 3 (XP_009111303.1) were among those identified to be up-regulated in both roots and leaves, whereas aquaporin TIP2-1 (XP_013653350.1) was down-regulated in both tissues.
To understand the universal response of Chinese cabbage to drought stress, all screened DEGs in the six comparison groups were subjected to GO enrichment analysis, and 60 functional groups classified into three categories comprising “cellular component,” “molecular function,” and “biological process” were identified (Table S1). The category “biological process,” consisting of 27 functional groups, exhibited the highest number of annotations alongside “cellular process”, followed by “metabolic process”, and “single-organism process”. In addition, the drought-induced DEGs revealed various affected metabolic processes, mainly the organic substance metabolic processes, cellular metabolic processes and primary metabolic processes (Figure 2C). Photosynthesis is one of the important processes affected by drought, due to changes including decrease in turgor pressure, reduction in CO2 assimilation and diffusion, and impairment in the photosynthetic apparatus [20,21,22], suggesting that alterations in primary metabolic processes under drought conditions (Figure 2C) are mediated via drought-induced photosynthesis inhibition.
Further, transcript abundances of genes involved in the abscisic acid (ABA) biosynthesis and degradation pathway were dynamically changed in response to drought stress, whereas a number of genes involved in biosynthesis of auxin, gibberellin, or ethylene were down-regulated (Figure S1), suggesting that the hormonal regulation of drought stress response in Chinese cabbage is primarily mediated via ABA.

2.3. Transcription Factors and MapMan Analysis of Chinese Cabbage Genes Associated with Drought-Stress Response

Transcription factors (TFs) are known as major players in various transcriptional regulatory mechanisms in stress-induced signal transduction pathways. Therefore, understanding the role of stress-inducible TFs is important for engineering stress-tolerant plants by modulating a large set of genes [23]. In Chinese cabbage transcriptome libraries generated in the above mentioned RNA-seq, a total of 3663 TFs mainly distributed into 45 TF families was found (Figure S2). Among these TFs, the most abundant differentially expressed TF family was the AP2/ERF family (Figure S3). Based on the difference in copy numbers and similarity of AP2/ERF domains, this family is further divided into five subfamilies, the Apetala2 (AP2), the ethylene responsive element binding factor (ERF), the Related to ABI3/VP1 (RAV), the dehydration-responsive element-binding (DREB) protein, and the Soloist [24], However, our transcriptome libraries did not contain members of the Soloist subfamily (Figure S4), probably due to its tissue-specific expression. In fact, expression of the Soloist gene was detected only in the bud of Chinese cabbage plants [25]. The ERF subfamily, which plays significant roles in regulating gene expression in response to biotic and abiotic stresses [6], was the most abundant class of the AP2/ERF family differentially expressed in response to drought (Figure S4).
The bZIP TF family, which contain a basic region for DNA-binding and a leucine zipper motif for dimerization, is one of the largest and most diverse TF families in plants [26], and plays diverse functions in plant growth and development and environmental stress responses, especially in ABA signaling related to stress responses [27]. Out of 136 bZIP TFs identified in B. rapa (BrbZIP) [28], 132 bZIP TFs were found in Chinese cabbage transcriptome libraries, and exhibited differential expression patterns (Figure 3A). Among them, 64 bZIP TFs were mainly differentially expressed in drought-treated plants. A total of 17 bZIP TFs were increased in both leaves and roots of drought-treated plants, and five bZIP TFs were only up-regulated in leaves. In addition, seven bZIP TFs were only down-regulated in leaves, six bZIP TFs only in roots, and three bZIP TFs were down-regulated in both roots and leaves. Drought stress-responsive gene expression is regulated by TFs operating in ABA-dependent and ABA-independent signaling pathways. Among these TFs, members of the bZIP TF family regulate the transcriptional induction of genes directly or indirectly linked with stress tolerance in plants [29]. In Arabidopsis, accumulated evidence has suggested that group-A bZIP including AtbZIP35, AtbZIP36, AtbZIP37, AtbZIP38, AtbZIP39, AtbZIP40, and AtbZIP66 play crucial roles in activating plant ABA signaling or in abiotic stress response [30,31]. We also found that BrbZIP23 (XP_009112062.1), BrbZIP15 (XP_009101738.1), BrbZIP17 (XP_009123528.1), BrbZIP21 (XP_009135624.1), BrbZIP22 (XP_009145807.1), and BrbZIP10 (XP_009104023.1), belonging to group-A BrbZIPs, were involved in drought response (Figure 3B). In addition, BrbZIP101 (XP_009111876.1), BrbZIP114 (XP_009138305.1), and BrbZIP115 (XP_009111076.1), belonging to group-S bZIPs, were also found to be induced under drought conditions. Although the physiological function of group-S bZIPs remains poorly understood, the increased transcription levels of group-S bZIPs due to abiotic stresses might suggest their regulatory roles in stress response in plants [32,33,34,35]. Further, poplar PtabZIP1-like (PtabZIP1L), homologous to the Arabidopsis bZIP1 belonging to group-S bZIPs, has been known to be a positive modulator of drought resistance acting through modulation of multiple metabolic pathways including flavonoid biosynthesis [36]. These findings indicate that drought-induced group-S BrbZIPs (Figure 3C) might be interesting candidate genes for improving drought tolerance.
To validate the gene expression results from RNA-seq, the expression of nine bZIP TFs up-regulated in both roots and leaves was analyzed using qRT-PCR. As shown in Figure 3B, the expression pattern of these genes was found to vary in response to drought. In addition, exposure of Chinese cabbage to NaCl or wounding also significantly induced or repressed the expression of these genes, although the changes were relatively lower than those due to drought treatment (Figure S5). Further, selected bZIP TFs were found to show complex expression patterns in response to drought, NaCl, and wounding, as evidenced by cluster analysis results represented as a heatmap (Figure 3C). In cluster I, BrbZIP (XP_009111876.1), also known as BrbZIP101 [26] exhibited relative higher transcript accumulation in drought-treated leaves and was transiently down-regulated in response to NaCl and wounding stresses, whereas BrbZIPs in cluster II only showed transient down-regulation due to wounding stress. These results indicate that BrbZIPs have divergent functions in response to environmental stresses, and the increased expression of BrbZIP101 might suggest its role as a marker of drought-stress response.
Functional enrichment analysis with respect to known gene ontologies using MapMan can be performed to obtain a global overview of high-throughput data in the context of metabolic pathways and biological processes [37]. To investigate the metabolic pathways affected by drought stress, MapMan analysis was performed on L2/L1, L3/L1, R2/R1, and R3/R1 comparisons based on DEGs. In leaves and roots of drought-treated plants, several metabolic pathways including those involving cell wall synthesis, light reactions, and the Calvin cycle were found to be down-regulated, whereas those involving fermentation, starch synthesis, and amino acid synthesis were up-regulated (Figure 4). However, the impact of drought stress on the transcript levels of genes involved in these metabolic pathways in roots was higher than that in leaves. Noticeably, the genes involved in the sulfur-containing secondary metabolites (S-misc) category pathways were up-regulated in the leaves but down-regulated in the roots under drought-stress conditions, indicating differences between organs of Chinese cabbage in response to drought stress.

2.4. Alteration in Glucosinolate Metabolism in Rresponse to Drought Stress

In higher plants, a variety of sulfur-containing secondary metabolites are synthesized, which often play a crucial role in the survival of plants under biotic and abiotic stresses [38]. The alteration in sulfur metabolism under drought conditions is consistent with the promotion of root growth [39]. The down-regulation of glucosinolate metabolism in roots upon drought suggests a switch in sulfur metabolism from glucosinolate to other metabolic processes, which are required for enhancing root growth. Interestingly, we found that the genes involved in sulfur-containing secondary metabolites have opposing expression patterns in leaves and roots during response to drought stress (Figure 4). With respect to the S-misc BIN category, the genes mainly represented were those responsible for glucosinolate biosynthesis including branched-chain amino acid aminotransferase 4 (BCAT4) and methylthioalkylmalate synthase 1 (MAMI1) (Figure S6, Table S2). To test whether the increased expression of these genes correlated with the accumulation of glucosinolates, we analyzed drought-induced variation in glucosinolate composition and content. As shown in Figure 5, the total glucosinolate content was increased in response to drought stress in leaves of Chinese cabbage, suggesting that the accumulation of glucosinolates in L3 was due to increase transcription of glucosinolate biosynthesis genes.
In addition, variation in the levels of identified glucosinolate compounds between non-stressed plants (L1) and drought-treated plants (L3) was observed. The compound 4-methoxyglucobrassicin (0.73 ± 0.01 nmol/mg of D.W) was found to be the major constituent in Chinese cabbage cultivars (Chunkwang), followed by progoitrin (0.55 ± 0.01 nmol/mg of D.W) and glucobrassicin (0.42 ± 0.001 nmol/mg of D.W), whereas glucobrassicanapin (1.50 ± 0.02 nmol/mg of D.W) and 4-methoxyglucobrassicin (1.02 ± 0.002 nmol/mg of D.W) were the abundant glucosinolates in L3. Further, drought-stress treatment led to an approximately six-fold increase in levels of gluconapin and 4-hydroxyglucobrassicin, indicating that drought stress increased the transcription level of genes involved in the glucosinolate biosynthetic pathway (Figure S6) and glucosinolate accumulation in the leaves of Chinese cabbage (Figure 5), consistent with the protein competition model which predicts the effects of environmental factors on secondary metabolism [40,41]. Similarly, the accumulation of glucosinolates has been observed in Brassica species as a part of the plant’s response to drought, although this response was highly variable depending on the duration and intensity of drought [40]. These findings indicate that drought stress significantly influenced the glucosinolate composition and content through modulating the expression of a range of glucosinolate biosynthesis genes.
Glucosinolate degradation products generated by the enzymatic action of myrosinases have attracted attention for several years due to their important role in plant defense against pests and pathogens [42]. In addition, the multiple roles of glucosinolates and their degradation products in heat tolerance, water transport, and transcriptional reprogramming suggest the physiological importance of these compounds in plant response to environmental stresses [43,44,45]. Exogenous application of glucosinolate degradation products, isothiocyanates, nitriles, and thiocyanates has been found to induce stomatal closure, allowing suppression of water loss [46,47]. Here, we hypothesized that drought-induced accumulation of glucosinolates may regulate stomatal closure leading to suppression of water loss. To test this hypothesis, glucosinolate-derived products were extracted from L1 and L3 samples and were applied to leaf disks floated on the opening medium. Application of glucosinolate-derived products significantly reduced stomatal apertures compared with application of mock-control (solvent control) (Figure 6).
The movement of the stomatal guard cells was found to be more responsive to glucosinolate-derived products obtained from L3 than to those obtained from L1 (Figure 6), indicating that the accumulation of glucosinolates induced drought tolerance through stimulating stomatal closure. Taken together, the increased concentration of sulfate transported from the roots to the leaves might be required for enhancing the anti-transpiration effect of ABA in the stomata in leaves [48,49], and glucosinolate degradation products function analogously to ABA in stomatal aperture regulation [50]. These findings suggest that sulfur metabolism plays significant roles in drought stress signaling and responses [51] and indicate that variation in glucosinolate metabolism between roots and leaves (Figure 5) is required for successful defense against drought conditions.

3. Materials and Methods

3.1. Plant Material and Stress Treatment

Seeds of Chinese cabbage (B. rapa L. subsp. pekinensis) cultivar Chunkwang were obtained from the Sakata Korea Seed Co. Ltd. (Seoul, South Korea), and were germinated and grown in a growth chamber with a long photoperiod (16 h light/8 h dark) at 24 °C. Soil water content was measured using HH2 Moisture Meter (Delta-T Devices, Cambridge, UK). The water content in the soil was controlled to be 60 to 65% of saturated soil water content. For drought stress treatment, 4-week-old plants (S1) were subjected to water withholding for 2 d or 4 d, resulting in a soil water content of 20% (S2) or 5% (S3), respectively. In addition, 4-week-old plants were watered with 250 mM NaCl or effectively wounded using sterile syringe needles to analyze the expression pattern of bZIP TFs. Each treatment was carried out three times, resulting in three biological replicates.

3.2. Determination of Relative Water Content, Malondialdehyde (MDA) Content, H2O2 Accumulation, and Protein Carbonyl Content

Relative water content was determined according to Soltys-Kalina et al. [52], and MDA content was determined using the thiobarbituric acid (TBA) reaction as described by Spanic et al. [53]. About 200 mg of ground material obtained from each sample (leaves from S1, S2, and S3) was used. MDA concentration was estimated by subtracting the non-specific turbidity at 600 nm from the absorption at 532 nm, and was expressed as MDA in nmol per mg of fresh weight using an absorbance coefficient of extinction (155 mM−1 cm−1).
As described by Spanic et al. [53], the supernatant obtained via 0.1% ice-cooled TCA extraction was used to analyze H2O2 concentration. After reaction with 1 M KI, the absorbance at 390 nm was determined via comparison to a standard curve prepared using H2O2. H2O2 concentration in each sample was expressed as µM per mg of fresh weight.
For protein carbonyl determination, 300 mg of frozen dried samples were incubated with protein extraction buffer containing 50 mM Tris-HCl, pH 7.5, 1 mM EDTA, 5 mM MgCl2, and 10 mM β-mercaptoethanol. After centrifugation at 12,000 g for 15 min at 4 °C, the protein content in the supernatant was determined following the Bradford method with bovine serum albumin as standard. Protein carbonyl content was determined using a fluorometric protein carbonyl content assay kit (BioVision, Milpitas, CA, USA) following the manufacturer’s instructions.

3.3. Identification and Functional Annotation of DEGs

Total RNA was extracted from each sample using the FavorPrep Plant Total RNA Mini Kit (Favorgen Biotech Corp, Pingtung, Taiwan), and quantified using an DeNovix DS-11 (DeNovix Inc., Wilmington, DE, USA). Equal amounts of total RNA from different stages of each sample were pooled to generate a cDNA library. Equal amounts of total RNA from different stages of each sample were pooled to generate a cDNA library. After DNaseI treatment, poly(A) mRNA was isolated with oligo(dT) magnetic beads and fragmented to short fragments. Then, the cDNA library was synthesized as described by Choi et al. [54] and sequenced using an Illumina HiSeq™ 2500 sequencing platform. The quality of raw data was controlled by using the FastQC tool, and adapter sequences, empty reads, low-quality reads (with ambiguous sequence, N), and reads with more than 10% Q < 20 bases (i.e., with a base quality of less than 20) were discarded using Trimmomatic v.0.33. The clean reads were mapped to the B. rapa reference sequence (BrapaV2.1, http://brassicadb.org/brad/datasets/pub/Genomes/Brassica_rapa/V2.0/V2.1/Chr/) using the HISAT2 aligner [55]. Transcript levels were calculated using SAMtools (http://samtools.sourceforge.net/) and relative transcript abundances were analyzed using DEseq [56]. DEGs were determined by combining a p value cutoff of 0.01 and adjusting to |log2 (fold change)| ≥ 1.For functional annotation, gene ontology (GO) enrichment analysis of DEGs was performed using Blast2GO (v.4.1.0) In addition, MapMan metabolism overview maps were drawn as described by Yu et al. [57].

3.4. Quantitative Real-Time PCR (qRT-PCR) Analysis

To further verify the expression profiles observed in the RNA-Seq data, nine bZIP TFs were selected for qRT-PCR analysis performed using the SYBR® Green Real-time PCR Master Mix (TOYOBO, Co., Ltd., Osaka, Japan) in the CFX96TM Real-time system (Bio-Rad, Hercules, CA, USA) with default parameters. The conditions for all reactions were as follows: 95 °C for 60 s, 40 cycles of 95 °C for 5 s, followed by 60 °C for 15 s, and 72 °C for 20 s. Melting curve analysis was performed to confirm the PCR specificity. The relative expression levels of selected bZIP TFs were normalized to expression of the internal reference gene actin [58]. Specific primer pairs used in qRT-PCR are listed in Table S3.

3.5. HPLC Analysis of Glucosinolates

Extraction and desulfation of glucosinolates were performed following the method described by Kim et al. [59]. Desulfated glucosinolates were analyzed using an HPLC system (1200 series, Agilent Technologies, Waldbronn, Germany) with an Inertsil ODS-3 (C18) column (4.6 mm × 250 mm, particle size 5 μm, GL Sciences, Tokyo, Japan). The mobile phase consisted of solvents A (water) and B (acetonitrile) with a gradient elution comprising 0 min (99% A), 0–18 min (1–30% B, linear gradient), and 18–30 min (30% B, isocratic). Glucosinolate content was calculated as sinigrin equivalents [60].

3.6. Measurement of Stomatal Aperture

Stomatal apertures were measured as described by Hossain et al. [47]. Briefly, leaf disks were floated on the opening medium (5 mM KCl, 50 μM CaCl2 and 10 mM MES-Tris, pH 6.15) for 2 h under light conditions. To analyze the effect of glucosinolate-derived products on stomatal movement, glucosinolate-derived products was extracted from 50 mg of frozen dried samples (leaves from S1 and S3) as described by Lola-Luz et al. [61]. Each extract was re-suspended in 5 mL HPLC-grade methanol and filtered using a membrane filter (0.2 μm pore size). Then, leaf disks floated on the opening medium were treated with 5% of extract (50 μL of each extract suspended in 5 mL MeOH was mixed with 1 mL of opening medium). Two hours after treatment, stomatal apertures were captured using an CKX53 microscope (Olympus, Tokyo, Japan) and the width and length of 100 stomatal apertures were measured using the image processing software ImageJ (National Institute of Mental Health, Bethesda, MD, USA).

4. Conclusions

Because of the complexity of drought as a stress signal, understanding the biochemical and molecular mechanisms of plant response has remained a major challenge in plant biology and breeding. The present study provides an overview of the changes in transcript abundance in leaves and roots under drought conditions. DEGs found to be responsive to drought suggest differential regulation patterns in the leaf and root, and drought-induced group-S BrbZIPs are interesting candidates for improving plant tolerance. Further, we suggest that organ-specific metabolic responses to drought are part of the survival strategies of plants. Taken together, our dataset provides a crucial starting point for future efforts in understanding the biochemical and molecular mechanisms underlying response to drought stress in Chinese cabbage.

Supplementary Materials

The following are available online.

Author Contributions

S.H.E. and T.K.H. conceived and designed the experiments, analyzed the data, and wrote the manuscript. S.-AB. and J.K.K. performed the experiments. The authors read and approve the final manuscript.

Acknowledgments

This work was carried out with the support of “Cooperative Research Program for Agriculture Science and Technology Development (Project No. PJ01261303)” Rural Development Administration, Republic of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, S.; Vanderbeld, B.; Wan, J.; Huang, Y. Narrowing down the targets: Towards successful genetic engineering of drought-tolerant crops. Mol. Plant 2010, 3, 469–490. [Google Scholar] [CrossRef] [PubMed]
  2. Vurukonda, S.S.; Vardharajula, S.; Shrivastava, M.; SkZ, A. Enhancement of drought stress tolerance in crops by plant growth promoting rhizobacteria. Microbiol. Res. 2016, 184, 13–24. [Google Scholar] [CrossRef] [PubMed]
  3. Lobell, D.B.; Roberts, M.J.; Schlenker, W.; Braun, N.; Little, B.B.; Rejesus, R.M.; Hammer, G.L. Greater sensitivity to drought accompanies maize yield increase in the U.S. Midwest. Science 2014, 344, 516–519. [Google Scholar] [CrossRef] [PubMed]
  4. Beck, E.H.; Fettig, S.; Knake, C.; Hartig, K.; Bhattarai, T. Specific and unspecific responses of plants to cold and drought stress. J. Biosci. 2007, 32, 501–510. [Google Scholar] [CrossRef] [PubMed]
  5. Shinozaki, K.; Yamaguchi-Shinozaki, K. Gene networks involved in drought stress response and tolerance. J. Exp. Bot. 2007, 58, 221–227. [Google Scholar] [CrossRef] [PubMed]
  6. Joshi, R.; Wani, S.H.; Singh, B.; Bohra, A.; Dar, Z.A.; Lone, A.A.; Pareek, A.; Singla-Pareek, S.L. Transcription factors and plants response to drought stress: Current understanding and future directions. Front. Plant Sci. 2016, 7, 1029. [Google Scholar] [CrossRef] [PubMed]
  7. Brassica rapa Genome Sequencing Project Consortium. The genome of the mesopolyploid crop species Brassica rapa. Nat. Genet. 2011, 43, 1035–1039. [Google Scholar]
  8. Qiu, N.; Liu, Q.; Li, J.; Zhang, Y.; Wang, F.; Gao, J. Physiological and transcriptomic responses of Chinese cabbage (Brassica rapa L. ssp. pekinensis) to salt stress. Int. J. Mol. Sci. 2017, 18, 1953. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, A.; Hu, J.; Huang, X.; Li, X.; Zhou, G.; Yan, Z. Comparative transcriptome analysis reveals heat-responsive genes in Chinese cabbage (Brassica rapa ssp. chinensis). Front. Plant Sci. 2016, 28, 939. [Google Scholar] [CrossRef] [PubMed]
  10. Sun, C.; Johnson, J.M.; Cai, D.; Sherameti, I.; Oelmüller, R.; Lou, B. Piriformospora indica confers drought tolerance in Chinese cabbage leaves by stimulating antioxidant enzymes, the expression of drought-related genes and the plastid-localized CAS protein. J. Plant Physiol. 2010, 167, 1009–1017. [Google Scholar] [CrossRef] [PubMed]
  11. Lee, S.G.; Lee, H.J.; Kim, S.K.; Choi, C.S.; Park, S.T.; Jang, Y.A.; Do, K.R. Effects of vernalization, temperature, and soil drying periods on the growth and yield of Chinese cabbage. Korean J. Hortic. Sci. 2015, 33, 820–828. [Google Scholar] [CrossRef]
  12. Yan, M. Seed priming stimulate germination and early seedling growth of Chinese cabbage under drought stress. S. Afr. J. Bot. 2015, 99, 88–92. [Google Scholar] [CrossRef]
  13. Hrdlickova, R.; Toloue, M.; Tian, B. RNA-Seq methods for transcriptome analysis. Wiley Interdiscip. Rev. RNA 2017, 8, e1364. [Google Scholar] [CrossRef] [PubMed]
  14. Miller, G.; Suzuki, N.; Ciftci-Yilmaz, S.; Mittler, R. Reactive oxygen species homeostasis and signalling during drought and salinity stresses. Plant Cell Environ. 2010, 33, 453–467. [Google Scholar] [CrossRef] [PubMed]
  15. Furlan, A.; Llanes, A.; Luna, V.; Castro, S. Physiological and biochemical responses to drought stress and subsequent rehydration in the symbiotic association peanut-Bradyrhizobium sp. ISRN Agron. 2012, 2012, 318083. [Google Scholar] [CrossRef]
  16. Toscano, S.; Farieri, E.; Ferrante, A.; Romano, D. Physiological and biochemical responses in two ornamental shrubs to drought stress. Front. Plant Sci. 2016, 7, 645. [Google Scholar] [CrossRef] [PubMed]
  17. Cao, Y.; Luo, Q.; Tian, Y.; Meng, F. Physiological and proteomic analyses of the drought stress response in Amygdalus Mira (Koehne) Yü et Lu roots. BMC Plant Biol. 2017, 17, 53. [Google Scholar] [CrossRef] [PubMed]
  18. Des Marais, D.L. Into the fourth dimension. eLife 2017, 6, e31328. [Google Scholar] [CrossRef] [PubMed]
  19. Hu, H.; Xiong, L. Genetic engineering and breeding of drought-resistant crops. Annu. Rev. Plant Biol. 2014, 65, 715–741. [Google Scholar] [CrossRef] [PubMed]
  20. Ohashi, Y.; Nakayama, N.; Saneokai, H.; Fujita, K. Effects of drought stress on photosynthetic gas exchange, chlorophyll fluorescence and stem diameter of soybean plants. Biol. Plant 2006, 50, 138–141. [Google Scholar] [CrossRef]
  21. Hu, L.; Wang, Z.; Huang, B. Diffusion limitations and metabolic factors associated with inhibition and recovery of photosynthesis from drought stress in a C perennial grass species. Physiol. Plant. 2010, 139, 93–106. [Google Scholar] [CrossRef] [PubMed]
  22. Pinheiro, C.; Chaves, M.M. Photosynthesis and drought: Can we make metabolic connections from available data? J. Exp. Bot. 2011, 62, 869–882. [Google Scholar] [CrossRef] [PubMed]
  23. Singh, D.; Laxmi, A. Transcriptional regulation of drought response: A tortuous network of transcriptional factors. Front. Plant Sci. 2015, 6, 895. [Google Scholar] [CrossRef] [PubMed]
  24. Lakhwani, D.; Pandey, A.; Dhar, Y.V.; Bag, S.K.; Trivedi, P.K.; Asif, M.H. Genome-wide analysis of the AP2/ERF family in Musa species reveals divergence and neofunctionalisation during evolution. Sci. Rep. 2016, 6, 18878. [Google Scholar] [CrossRef] [PubMed]
  25. Song, X.; Li, Y.; Hou, X. Genome-wide analysis of the AP2/ERF transcription factor superfamily in Chinese cabbage (Brassica rapa ssp. pekinensis). BMC Genom. 2013, 14, 573. [Google Scholar] [CrossRef] [PubMed]
  26. Bai, Y.; Zhu, W.; Hu, X.; Sun, C.; Li, Y.; Wang, D.; Wang, Q.; Pei, G.; Zhang, Y.; Guo, A.; et al. Genome-wide analysis of the bZIP gene family identifies two ABI5-like bZIP transcription factors, BrABI5a and BrABI5b, as positive modulators of ABA signalling in Chinese cabbage. PLoS ONE 2016, 11, e0158966. [Google Scholar] [CrossRef] [PubMed]
  27. Ji, L.; Wang, J.; Ye, M.; Li, Y.; Guo, B.; Chen, Z.; Li, H.; An, X. Identification and characterization of the Populus AREB/ABF subfamily. J. Integr. Plant Biol. 2013, 55, 177–186. [Google Scholar] [CrossRef] [PubMed]
  28. Hwang, I.; Jung, H.J.; Park, J.I.; Yang, T.J.; Nou, I.S. Transcriptome analysis of newly classified bZIP transcription factors of Brassica rapa in cold stress response. Genomics 2014, 104, 194–202. [Google Scholar] [CrossRef] [PubMed]
  29. Agarwal, P.K.; Jha, B. Transcription factors in plants and ABA dependent and independent abiotic stress signalling. Biol. Plant. 2010, 54, 201–212. [Google Scholar] [CrossRef]
  30. Nakashima, K.; Yamaguchi-Shinozaki, K.; Shinozaki, K. The transcriptional regulatory network in the drought response and its crosstalk in abiotic stress responses including drought, cold, and heat. Front. Plant Sci. 2014, 5, 170. [Google Scholar] [CrossRef] [PubMed]
  31. Julkowska, M.M.; Testerink, C. Tuning plant signaling and growth to survive salt. Trends Plant Sci. 2015, 20, 586–594. [Google Scholar] [CrossRef] [PubMed]
  32. Lee, S.S.; Yang, S.H.; Bererich, T.; Miyazaki, A.; Kusano, T. Characterization of AtbZIP2, AtbZIP11 and AtbZIP53 from the group S basic region-leucine zipper family in Arabidopsis thaliana. Plant Biotechnol. 2006, 23, 249–258. [Google Scholar] [CrossRef]
  33. Obertello, M.; Krouk, G.; Katari, M.S.; Runko, S.J.; Coruzzi, G.M. Modeling the global effect of the basic leucine zipper transcription factor 1 (bZIP1) on nitrogen and light regulation in Arabidopsis. BMC Syst. Biol. 2010, 4, 111. [Google Scholar] [CrossRef] [PubMed]
  34. Sun, X.; Li, Y.; Cai, H.; Bai, X.; Ji, W.; Ding, X.; Zhu, Y. The Arabidopsis AtbZIP1 transcription factor is a positive regulator of plant tolerance to salt, osmotic and drought stresses. J. Plant Res. 2012, 125, 429–438. [Google Scholar] [CrossRef] [PubMed]
  35. Moon, S.J.; Han, S.Y.; Kim, D.Y.; Yoon, I.S.; Shin, D.; Byun, M.O.; Kwon, H.B.; Kim, B.G. Ectopic expression of a hot pepper bZIP-like transcription factor in potato enhances drought tolerance without decreasing tuber yield. Plant Mol. Biol. 2015, 89, 421–431. [Google Scholar] [CrossRef] [PubMed]
  36. Dash, M.; Yordanov, Y.S.; Georgieva, T.; Tschaplinski, T.J.; Yordanova, E.; Busov, V. Poplar PtabZIP1-like enhances lateral root formation and biomass growth under drought stress. Plant J. 2017, 89, 692–705. [Google Scholar] [CrossRef] [PubMed]
  37. Bisaga, M.; Lowe, M.; Hegarty, M.; Abberton, M.; Ravagnani, A. Deep sequencing of suppression subtractive hybridisation drought and recovery libraries of the non-model crop Trifolium repens L. Front. Plant Sci. 2017, 8, 213. [Google Scholar] [CrossRef] [PubMed]
  38. Rausch, T.; Wachter, A. Sulfur metabolism: A versatile platform for launching defence operations. Trends Plant Sci. 2005, 10, 503–509. [Google Scholar] [CrossRef] [PubMed]
  39. Ahmad, N.; Malagoli, M.; Wirtz, M.; Hell, R. Drought stress in maize causes differential acclimation responses of glutathione and sulfur metabolism in leaves and roots. BMC Plant Biol. 2016, 16, 247. [Google Scholar] [CrossRef] [PubMed]
  40. Del Carmen Martínez-Ballesta, M.; Moreno, D.A.; Carvajal, M. The physiological importance of glucosinolates on plant response to abiotic stress in Brassica. Int. J. Mol. Sci. 2013, 14, 11607–11625. [Google Scholar] [CrossRef] [PubMed]
  41. Jones, C.G.; Hartley, S.E. A protein competition model of phenolic allocation. Oikos 1999, 86, 27–44. [Google Scholar] [CrossRef]
  42. Brader, G.; Mikkelsen, M.D.; Halkier, B.A.; Tapio Palva, E. Altering glucosinolate profiles modulates disease resistance in plants. Plant J. 2006, 46, 758–767. [Google Scholar] [CrossRef] [PubMed]
  43. Hara, M.; Harazaki, A.; Tabata, K. Administration of isothiocyanates enhances heat tolerance in Arabidopsis thaliana. Plant Growth Regul. 2013, 69, 71–77. [Google Scholar] [CrossRef]
  44. Martínez-Ballesta Mdel, C.; Muries, B.; Moreno, D.Á.; Dominguez-Perles, R.; García-Viguera, C.; Carvajal, M. Involvement of a glucosinolate (sinigrin) in the regulation of water transport in Brassica oleracea grown under salt stress. Physiol. Plant. 2014, 150, 145–160. [Google Scholar] [CrossRef] [PubMed]
  45. Kissen, R.; Øverby, A.; Winge, P.; Bones, A.M. Allyl-isothiocyanate treatment induces a complex transcriptional reprogramming including heat stress, oxidative stress and plant defence responses in Arabidopsis thaliana. BMC Genom. 2016, 17, 740. [Google Scholar] [CrossRef] [PubMed]
  46. Khokon, M.A.; Jahan, M.S.; Rahman, T.; Hossain, M.A.; Muroyama, D.; Minami, I.; Munemasa, S.; Mori, I.C.; Nakamura, Y.; Murata, Y. Allyl isothiocyanate (AITC) induces stomatal closure in Arabidopsis. Plant Cell Environ. 2011, 34, 1900–1906. [Google Scholar] [CrossRef] [PubMed]
  47. Hossain, M.S.; Ye, W.; Hossain, M.A.; Okuma, E.; Uraji, M.; Nakamura, Y.; Mori, I.C.; Murata, Y. Glucosinolate degradation products, isothiocyanates, nitriles, and thiocyanates, induce stomatal closure accompanied by peroxidase-mediated reactive oxygen species production in Arabidopsis thaliana. Biosci. Biotechnol. Biochem. 2013, 77, 977–983. [Google Scholar] [CrossRef] [PubMed]
  48. Gallardo, K.; Courty, P.E.; Le Signor, C.; Wipf, D.; Vernoud, V. Sulfate transporters in the plant’s response to drought and salinity: Regulation and possible functions. Front. Plant Sci. 2014, 5, 580. [Google Scholar] [CrossRef] [PubMed]
  49. Ernst, L.; Goodger, J.Q.; Alvarez, S.; Marsh, E.L.; Berla, B.; Lockhart, E.; Jung, J.; Li, P.; Bohnert, H.J.; Schachtman, D.P. Sulphate as a xylem-borne chemical signal precedes the expression of ABA biosynthetic genes in maize roots. J. Exp. Bot. 2010, 61, 3395–3405. [Google Scholar] [CrossRef] [PubMed]
  50. Zhu, M.; Assmann, S.M. Metabolic signatures in response to abscisic acid (ABA) treatment in Brassica napus guard cells revealed by metabolomics. Sci. Rep. 2017, 7, 12875. [Google Scholar] [CrossRef] [PubMed]
  51. Chan, K.X.; Wirtz, M.; Phua, S.Y.; Estavillo, G.M.; Pogson, B.J. Balancing metabolites in drought: The sulfur assimilation conundrum. Trends Plant Sci. 2013, 18, 18–29. [Google Scholar] [CrossRef] [PubMed]
  52. Soltys-Kalina, D.; Plich, J.; Strzelczyk-Żyta, D.; Śliwka, J.; Marczewski, W. The effect of drought stress on the leaf relative water content and tuber yield of a half-sib family of ‘Katahdin’-derived potato cultivars. Breed. Sci. 2016, 66, 328–331. [Google Scholar] [CrossRef] [PubMed]
  53. Spanic, V.; Viljevac Vuletic, M.; Abicic, I.; Marcek, T. Early response of wheat antioxidant system with special reference to Fusarium head blight stress. Plant Physiol. Biochem. 2017, 115, 34–43. [Google Scholar] [CrossRef] [PubMed]
  54. Choi, J.H.; Kim, H.; Hyun, T.K. Transcriptome analysis of Abeliophyllum distichum NAKAI reveals potential molecular markers and candidate genes involved in anthocyanin biosynthesis pathway. S. Afr. J. Bot. 2018, 116, 34–41. [Google Scholar] [CrossRef]
  55. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed]
  56. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  57. Yu, Y.; Zhen, S.; Wang, S.; Wang, Y.; Cao, H.; Zhang, Y.; Li, Y.; Yan, Y. Comparative transcriptome analysis of wheat embryo and endosperm responses to ABA and H2O2 stresses during seed germination. BMC Genom. 2016, 17, 1–18. [Google Scholar] [CrossRef] [PubMed]
  58. Kim, A.-R.; Lim, H.; Park, S.H.; Ji, H.; Yoon, U.-H.; Kim, T.-H.; Cho, J.-I.; Kim, C.-K.; Lee, G.-S. Molecular characterization of a putatively stress-related gene BrTSR53 isolated from Chinese cabbage (Brassica rapa). Korean J. Breed. Sci. 2014, 46, 143–151. [Google Scholar] [CrossRef]
  59. Kim, J.K.; Chu, S.M.; Kim, S.J.; Lee, D.J.; Lee, S.Y.; Lim, S.Y.; Ha, S.H.; Kweon, S.J.; Cho, H.S. Variation of glucosinolates in vegetable crops of Brassica rapa L. ssp. pekinensis. Food Chem. 2010, 119, 423–428. [Google Scholar] [CrossRef]
  60. Baek, S.A.; Jung, Y.H.; Lim, S.H.; Park, S.U.; Kim, J.K. Metabolic profiling in Chinese cabbage (Brassica rapa L. subsp. pekinensis) cultivars reveals that glucosinolate content is correlated with carotenoid content. J. Agric. Food Chem. 2016, 64, 4426–4434. [Google Scholar] [CrossRef] [PubMed]
  61. Lola-Luz, T.; Hennequart, F.; Gaffney, M. Effect on yield, total phenolic, total flavonoid and total isothiocyanate content of two broccoli cultivars (Brassica oleraceae var italica) following the application of a commercial brown seaweed extract (Ascophyllum nodosum). Agric. Food Sci. 2014, 23, 28–37. [Google Scholar]
Sample Availability: Samples of the compounds are available from the authors.
Molecules 23 01186 g001 550
Figure 1. Physiological response to drought stress in Chinese cabbage. (A) Phenotypes of Chinese cabbage plants after exposure to drought stress for four days. S1 (Stage 1), normal water supply; S2 (Stage 2), water withholding for 2 d (soil water content at 20%); S3 (Stage 3), water withholding for 4 d (soil water content at 5%). Relative water content (B) and changes in levels of H2O2 (C), MDA (D), and protein carbonylation (E) after drought treatment were determined. Values are averages from three biological replicates consisting three plants per sample were used. The data are presented as mean ± standard error (SE). * p < 0.05, ** p < 0.01, and *** p < 0.001 represent the significant differences in comparison with Stage 1 (S1).
Figure 1. Physiological response to drought stress in Chinese cabbage. (A) Phenotypes of Chinese cabbage plants after exposure to drought stress for four days. S1 (Stage 1), normal water supply; S2 (Stage 2), water withholding for 2 d (soil water content at 20%); S3 (Stage 3), water withholding for 4 d (soil water content at 5%). Relative water content (B) and changes in levels of H2O2 (C), MDA (D), and protein carbonylation (E) after drought treatment were determined. Values are averages from three biological replicates consisting three plants per sample were used. The data are presented as mean ± standard error (SE). * p < 0.05, ** p < 0.01, and *** p < 0.001 represent the significant differences in comparison with Stage 1 (S1).
Molecules 23 01186 g001
Molecules 23 01186 g002 550
Figure 2. Analysis of differentially expressed genes (DEGs). (A) Distribution of up- and down-regulated DEGs in each comparison; (B) Venn diagram analysis of DEGs in different comparisons among groups; (C) Classification of DEGs based on metabolism categories. L1, L2, and L3, leaf samples obtained from Stage 1, Stage 2, and Stage 3 plants, respectively; R1, R2, and R3, root samples obtained from Stage 1, Stage 2, and Stage 3 plants, respectively.
Figure 2. Analysis of differentially expressed genes (DEGs). (A) Distribution of up- and down-regulated DEGs in each comparison; (B) Venn diagram analysis of DEGs in different comparisons among groups; (C) Classification of DEGs based on metabolism categories. L1, L2, and L3, leaf samples obtained from Stage 1, Stage 2, and Stage 3 plants, respectively; R1, R2, and R3, root samples obtained from Stage 1, Stage 2, and Stage 3 plants, respectively.
Molecules 23 01186 g002
Molecules 23 01186 g003 550
Figure 3. Change in bZIP transcription factor (TF) expression during the response to abiotic stresses. (A) The heatmap shows the expression profiles of bZIP TFs in different comparisons among groups. The grey bar indicates bZIP TFs selected for qRT-PCR analysis; (B) The expression patterns of the selected bZIP TFs were analyzed using qRT-PCR. Transcript levels of the selected bZIP TFs were normalized to those of Chinese cabbage actin, and gene expression is relative to Stage 1 plants (L1 or R1) set to a value of 1 for each biological replicate. Values are averages from three independent biological experiments. Data are means ± SE. * p < 0.05, ** p < 0.01, and *** p < 0.001 represent the significant differences in comparison with L1 or R1; (C) Differential expression of the selected bZIP TFs under different stress treatments. Expression is indicated as a log2 ratio of experimental treatments relative to control samples and visualized in heatmaps showing hierarchical clustering of the selected bZIP TFs.
Figure 3. Change in bZIP transcription factor (TF) expression during the response to abiotic stresses. (A) The heatmap shows the expression profiles of bZIP TFs in different comparisons among groups. The grey bar indicates bZIP TFs selected for qRT-PCR analysis; (B) The expression patterns of the selected bZIP TFs were analyzed using qRT-PCR. Transcript levels of the selected bZIP TFs were normalized to those of Chinese cabbage actin, and gene expression is relative to Stage 1 plants (L1 or R1) set to a value of 1 for each biological replicate. Values are averages from three independent biological experiments. Data are means ± SE. * p < 0.05, ** p < 0.01, and *** p < 0.001 represent the significant differences in comparison with L1 or R1; (C) Differential expression of the selected bZIP TFs under different stress treatments. Expression is indicated as a log2 ratio of experimental treatments relative to control samples and visualized in heatmaps showing hierarchical clustering of the selected bZIP TFs.
Molecules 23 01186 g003
Molecules 23 01186 g004 550
Figure 4. Mapman metabolism overview showing the differentially expressed genes (DEGs) in each comparison. The different colors represent the log2 values of the gene expression levels in response to drought stress.
Figure 4. Mapman metabolism overview showing the differentially expressed genes (DEGs) in each comparison. The different colors represent the log2 values of the gene expression levels in response to drought stress.
Molecules 23 01186 g004
Molecules 23 01186 g005 550
Figure 5. Glucosinolate content in control (L1) and drought-treated (L3) Chinese cabbage leaves. Values are averages from three biological replicates. Error bars indicate S.E. * p < 0.05 and ** p < 0.01 indicate significantly different from L1.
Figure 5. Glucosinolate content in control (L1) and drought-treated (L3) Chinese cabbage leaves. Values are averages from three biological replicates. Error bars indicate S.E. * p < 0.05 and ** p < 0.01 indicate significantly different from L1.
Molecules 23 01186 g005
Molecules 23 01186 g006 550
Figure 6. Stomatal movement in Chinese cabbage in response to glucosinolate-derived products obtained from control (L1) and drought-treated (L3) Chinese cabbage leaves. Apertures were analyzed as width/length after 2 h of treatment. Values are means ± S.E., *** p < 0.001, n > 100.
Figure 6. Stomatal movement in Chinese cabbage in response to glucosinolate-derived products obtained from control (L1) and drought-treated (L3) Chinese cabbage leaves. Apertures were analyzed as width/length after 2 h of treatment. Values are means ± S.E., *** p < 0.001, n > 100.
Molecules 23 01186 g006
Table
Table 1. Summary of RNA sequencing data from six RNA libraries of control and drought-stressed leaves and roots.
Table 1. Summary of RNA sequencing data from six RNA libraries of control and drought-stressed leaves and roots.
Sample No.TissueDrought StressClean ReadsClean Bases (Gb)Accession Number (NABIC)
L1LeafStage 153,744,0645.38NN-4956
L2LeafStage 246,299,3124.64NN-4958
L3LeafStage 354,348,5425.45NN-4960
R1RootStage 154,764,1005.49NN-4957
R2RootStage 248,558,4684.87NN-4959
R3RootStage 347,117,4104.72NN-4961

Share and Cite

MDPI and ACS Style

Eom, S.H.; Baek, S.-A.; Kim, J.K.; Hyun, T.K. Transcriptome Analysis in Chinese Cabbage (Brassica rapa ssp. pekinensis) Provides the Role of Glucosinolate Metabolism in Response to Drought Stress. Molecules 2018, 23, 1186. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules23051186

AMA Style

Eom SH, Baek S-A, Kim JK, Hyun TK. Transcriptome Analysis in Chinese Cabbage (Brassica rapa ssp. pekinensis) Provides the Role of Glucosinolate Metabolism in Response to Drought Stress. Molecules. 2018; 23(5):1186. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules23051186

Chicago/Turabian Style

Eom, Seung Hee, Seung-A Baek, Jae Kwang Kim, and Tae Kyung Hyun. 2018. "Transcriptome Analysis in Chinese Cabbage (Brassica rapa ssp. pekinensis) Provides the Role of Glucosinolate Metabolism in Response to Drought Stress" Molecules 23, no. 5: 1186. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules23051186

Article Metrics

Back to TopTop