Next Article in Journal
Melatonin, a Full Service Anti-Cancer Agent: Inhibition of Initiation, Progression and Metastasis
Next Article in Special Issue
Overexpression of S-Adenosyl-l-Methionine Synthetase 2 from Sugar Beet M14 Increased Arabidopsis Tolerance to Salt and Oxidative Stress
Previous Article in Journal
Subtype Diagnosis of Primary Aldosteronism: Is Adrenal Vein Sampling Always Necessary?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 18
Issue 4
10.3390/ijms18040849
ijms-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 Sequencing of Dianthus spiculifolius and Analysis of the Genes Involved in Responses to Combined Cold and Drought Stress

by
Aimin Zhou
1,
Hongping Ma
1,
Enhui Liu
1,
Tongtong Jiang
1,
Shuang Feng
2,
Shufang Gong
1,* and
Jingang Wang
1,*
1
College of Horticulture and Landscape Architecture, Northeast Agricultural University, Harbin 150030, China
2
Key Laboratory of Saline-Alkali Vegetation Ecology Restoration in Oil Field (SAVER), Ministry of Education, Alkali Soil Natural Environmental Science Center (ASNESC), Northeast Forestry University, Harbin 150040, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(4), 849; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms18040849
Submission received: 18 March 2017 / Revised: 12 April 2017 / Accepted: 14 April 2017 / Published: 17 April 2017
(This article belongs to the Special Issue Selected Papers from the 6th National Plant Protein Research Congress)
Browse Figures
Versions Notes

Abstract

:
Dianthus spiculifolius, a perennial herbaceous flower and a member of the Caryophyllaceae family, has strong resistance to cold and drought stresses. To explore the transcriptional responses of D. spiculifolius to individual and combined stresses, we performed transcriptome sequencing of seedlings under normal conditions or subjected to cold treatment (CT), simulated drought treatment (DT), or their combination (CTDT). After de novo assembly of the obtained reads, 112,015 unigenes were generated. Analysis of differentially expressed genes (DEGs) showed that 2026, 940, and 2346 genes were up-regulated and 1468, 707, and 1759 were down-regulated in CT, DT, and CTDT samples, respectively. Among all the DEGs, 182 up-regulated and 116 down-regulated genes were identified in all the treatment groups. Analysis of metabolic pathways and regulatory networks associated with the DEGs revealed overlaps and cross-talk between cold and drought stress response pathways. The expression profiles of the selected DEGs in CT, DT, and CTDT samples were characterized and confirmed by quantitative RT-PCR. These DEGs and metabolic pathways may play important roles in the response of D. spiculifolius to the combined stress. Functional characterization of these genes and pathways will provide new targets for enhancement of plant stress tolerance through genetic manipulation.

1. Introduction

Adverse environmental conditions seriously affect the growth and development of plants. In particular, abiotic stresses such as cold, heat, drought and salinity have large detrimental effects on agricultural production. Investigations into the effects of stress have generally examined the responses to individual abiotic stresses [1,2]; however, abiotic stresses rarely occur individually and plants are often simultaneously exposed to a combination of stresses such as cold and drought. In the field, the combination of cold and drought is recognized as an important threat to plants during early spring and winter in temperate regions [3]. To date, the molecular mechanisms associated with plant responses to combined cold and drought stresses, such as changes in gene expression, signal transduction and regulatory networks, are largely unknown.
Transcriptome profiling has been carried out in various plant species to characterize their responses to individual stresses [4,5,6,7,8,9,10]. Only a few studies have investigated the molecular mechanisms underlying the tolerance to combined cold and drought stress. In Arabidopsis, plants with mutation in ein2 (ethylene insensitive) and aba1.6 (abscisic acid) show differential regulation of 2313 and 4131 transcripts, respectively, under cold and drought stress. Validation of some of the important differentially expressed genes (DEGs) showed that HSPs (heat shock proteins), COR47 (cold regulated proteins), and SPX1 (SPX domain containing proteins) were altered by stress in ein2 and aba1.6 mutants [11]. Furthermore, overexpression of an exogenous HSP and COR increases multiple abiotic stress tolerance in transgenic plants [12,13]. In tea (Camellia sinensis), 319 DEGs that enhance the resistance to combined cold and drought stress have been identified. Additionally, drought-induced leaf senescence is affected by the expression of senescence-associated genes and cold responsive genes and by the regulation of multiple metabolic pathways in response to a combined stress [14].
In contrast to model and economically important species, many wild plant species are able to grow in stressful environments because of their specialized genetic resources and adaptive mechanisms. Dianthus spiculifolius Schur, a perennial herbaceous flower and a member of Caryophyllaceae [15], shows strong resistance to cold and drought stresses and can survive winters in Russian Siberia and northeast China, where temperatures can fall to −40 °C. This strong adaptability suggests that D. spiculifolius may have distinct molecular mechanisms that allow it to tolerate extreme conditions. However, to date, no genomic information has been reported for D. spiculifolius.
In the present study, we sought to elucidate the adaptation mechanisms of D. spiculifolius to complex environmental conditions by examining the transcriptome profiles of seedlings grown under either cold or drought stress conditions, as well as a combination of the two stresses. We identified DEGs in cold-treated (CT), simulated drought-treated (DT), and combined cold and simulated drought-treated (CTDT) plants, and also identified numerous specific stress-related genes. We also analyzed the signal transduction and metabolic pathways associated with these DEGs. Our results provide a basis for understanding the mechanisms of the responses of D. spiculifolius to the combined stress and also provided a potentially valuable resource for the future development of stress-tolerant plants through genetic manipulation.

2. Results and Discussion

2.1. Transcriptome Sequencing of D. spiculifolius Subjected to Cold, Drought, and Combined Cold and Drought Stress

D. spiculifolius plants cultivated in an open field showed strong resistance to cold and drought stresses (Figure 1A,B). The plants grew to 8–10 cm in height and flowered in the autumn (Figure 1C), producing needle-like leaves and flowers, 2 cm in diameter (Figure 1D,E). In addition to asexual reproduction, D. spiculifolius can be propagated by seeds (Figure 1F–I). To investigate the changes in the patterns of transcription in response to combined cold and drought stress, we exposed 10-day-old seedlings to cold, simulated drought, and combined cold and simulated drought stress for 24 h. The Illumina HiSeq 2000 platform was used to analyze transcriptome sequences from the treated seedlings. In total, we obtained approximately 29–30 million raw reads for all the samples. After data cleaning and the removal of invalid reads, we obtained 23,624,696, 22,495,628, 23,317,141, and 22,965,345 clean reads for the control (CK), CT, DT, and CTDT samples, respectively. The reads had 46.45%, 46.32%, 46.62%, and 42.10% guanine and cytosine (GC) content in the CK, CT, DT, and CTDT samples, respectively. After clustering the high-quality reads, 184,782 transcripts and 112,015 unigenes were obtained for all the samples. The unigenes had an average length of 819 bp (Table 1).
To validate and annotate the assembled unigenes, sequence similarity searches were conducted using sequence- and domain-based alignments. In total, 36,142 (32.27% of all the unigenes), 29,149 (26.55%), 29,739 (26.54%), 25,073 (22.38%), and 10,845 (9.68%) unigenes were found in SwissProt, PFAM (Protein family), GO (Gene Ontology), KOG (euKaryotic Ortholog Groups), and KO (KEGG Ortholog) databases, respectively (Table 2).

2.2. Identification of Differentially Expressed Genes (DEGs) Responding to Combined Cold and Drought Stress

To investigate the gene expression during cold and drought stress, we used the reads per kilobase per million reads (RPKM) method to calculate the expression levels of the unigenes from CT, DT, and CTDT samples compared to the CK sample. As shown in Figure 2, 940 up-regulated and 707 down-regulated unigenes were identified in the DT sample. In the CT sample, 2026 unigenes were up-regulated and 1468 were down-regulated. In the CTDT sample, 2346 unigenes were up-regulated and 1759 were down-regulated. Among all the DEGs, 182 up-regulated and 116 down-regulated overlapping DEGs were identified in all the three treatment groups, respectively (Figure 2, Tables S1 and S2). This suggested that these overlapping DEGs may play an important role in the response of plants to combined cold and drought stress.

2.3. Functional Classification and Expression Analysis of DEGs

We used gene ontology (GO) annotation [16] to classify the possible functions of the overlapping DEGs. This analysis assigned the 182 up-regulated and 116 down-regulated DEGs to 39 and 18 functional groups, respectively; the three main categories were biological processes, cellular components, and molecular functions. In the molecular function category, the most frequent assignments of the up-regulated DEGs were transcription factor activity, sequence-specific DNA binding, and calcium binding. For cellular components, the DEGs were mostly assigned to nucleus, internal components of membrane, and cytoplasm. For biological processes, the majority of the DEGs were assigned to transcription, DNA-template, response to water dehydration, and cold (Figure 3A). The down-regulated DEGs were mainly assigned to the structural constituents of ribosome, chloroplast, and translation (Figure 3B).
We analyzed the protein–protein interactions (PPI) of the DEGs involved in five important GO terms using STRING (http://string.embl.de/). A PPI subnetwork including 47 nodes and 36 interactions was constructed (Figure 4). The degree of connectivity of each node was calculated; the top six nodes were ACC1 (acetyl-CoA carboxylase 1), PP2C16 (type-2C protein phosphatase 16), ATHB7 and ATHB12 (homeobox-leucine zipper proteins 7 and 12), NCED3 (9-cis-epoxycarotenoid dioxygenase 3), and MD37E (mediator of RNA polymerase II transcription subunit 37e). The transcription factors ATHB7 and ATHB12 were considered as key nodes that interacted with many proteins, including the transcription factors ERF79 (ethylene response factor), PP2C16, and NCED3. It has been reported that ATHB7 and ATHB12 modulate abscisic acid (ABA) signaling by regulating PP2C [17]. The ABA signaling pathway is one of the important signaling pathways in plants for response to cold and drought stresses [18,19]. Furthermore, PP2C16 interacted with the protein kinases MPK19 (mitogen-activated protein kinase 19), CIPKP (calcineurin B-like protein (CBL)-interacting serine/threonine-protein kinase 25), and CDL1 (serine/threonine-protein kinase CDL1), which can catalyze the phosphorylation of proteins and are an important part of cell signal recognition and transduction in the responses of plants to abiotic stress [18,20,21].
We produced a heat map based on the expression profiles to classify the key DEGs. The analysis showed that protein kinases and transcription factors constituted the largest number of DEGs (Figure 5). Among the 10 protein kinase genes, seven were up-regulated and three were down-regulated. The ten protein kinases included three CIPKs, two of which were up-regulated and the third of which was down-regulated. CIPKs are involved in the Ca2+-related signaling pathway that is central to the response of plants to both abiotic and biotic environmental stimuli [21,22,23]. Calcium is an important second messenger and it is thought that specific calcium ion signal patterns can result in the expression of targeted sets of genes [24]. CIPKs also interact with PP2C, which is involved in the ABA signal transduction pathway [25,26,27]. In the present study, two PP2C genes (PP2C16 and PP2C65) were identified: PP2C16 was up-regulated and PP2C65 was down-regulated. Among the 12 identified transcription factors, eight were up-regulated and the other four were down-regulated. Most transcription factors are associated with cold and drought stresses, such as DREB3 (dehydration-responsive element-binding protein 3), MYB39 (myeloblastosis 39), ERF79, and BBX32 (B-box zinc finger protein 32) [28,29,30,31,32]. Many stress response genes were identified, such as COR47 (cold-responsive gene 47), LEAs (late embryogenesis abundant) and HSPs (heat shock factor), and their expression was up-regulated. These genes are known to aid abiotic stress resistance in plants and help them cope with multiple stresses via ABA-dependent and ABA-independent pathways [33,34]. DREBs, also known as C-repeat binding factor (CBF), can regulate the expression of COR genes that confer freezing tolerance in plants [35]. The overexpression of ZmDREB2A results in the up-regulation of 44 genes, including LEAs, HSPs, and detoxification proteins, and increased drought tolerance in plants [36]. LEA proteins may have protective roles during cellular dehydration through physiological mechanisms such as hydration buffering, antioxidant protection and stabilization of sensitive processes [37,38]. As expected, the overexpression of LEAs from Oryza sativa enhances drought stress tolerance in transgenic rice [39,40]. In addition, ten transporter proteins and six auxin responsive proteins were identified. Most of the transporter proteins were up-regulated, whereas most of the auxin responsive proteins were down-regulated (Figure 5).

2.4. Validation of RNA-Seq Gene Expression

To confirm the reliability of the RNA-Seq data, quantitative real-time PCR (qRT-PCR) was performed on 16 randomly selected up-regulated genes. These DEGs included transcription factors (DsATHB7, DsDREB, DsAI5L3, DsBBX32, and DsMYB39), ABA signal-related proteins (DsASR1 and DsPP2C16), stress response proteins (DsHSP17, DsSTEP2, DsCOR47, and DsTIL), transporter proteins (DsPTR33, DsSWET1, and DsSWET16), and unknown proteins (DsUP1 and DsUP2). The 16 selected genes showed differential expression between control and stress-treated plants by qRT-PCR (Figure 6). The results of the qRT-PCR assay were correlated, with the RNA-Seq data (evaluated by RPKM), confirming the reproducibility of the RNA-Seq data.
Overall, we identified a large number of DEGs by transcriptome sequencing under cold and simulated drought stress, and annotated the regulatory pathways involving these DEGs. Our results provide further insights into the molecular mechanisms underlying the adaptation of D. spiculifolius to cold and drought conditions. However, there are limitations of the present study; for example, the simulation of drought by osmotic stress caused by mannitol could not completely replace the drought under natural conditions; the qPCR results of a few DEGs were not completely consistent with the RPKM data. Therefore, more reliable evidence from more experiments is desired for the validation of our results; for example, the increase in the type of stress and the experimental replications, qPCR validation at more time points, promoter analysis, examination of the vital physiological parameters, transgenic experiments, etc., are required and are being carried out for future reports.

3. Materials and Methods

3.1. Plant Growth Conditions and Stress Treatments

D. spiculifolius plants were grown in an open field at the Northeast Agricultural University (Harbin, China; 128.4° E, 45.0° N). D. spiculifolius seeds were surface sterilized and stratified at 4 °C for 2 days in the dark. After germination, the seedlings were grown on half-strength Murashige and Skoog (MS) medium (3% sucrose, 1% agar, pH 5.8) under a 12-h light/12-h dark photoperiod (100 μmol·m−2·s−1 light intensity) at 22 °C. Ten-day-old seedlings were assigned to control (no treatment, CK), cold (4 °C), simulated drought (300 mM mannitol), or combined cold and simulated drought (4 °C + 300 mM mannitol) treatment groups (10 seedlings per treatment) for 24 h. Three independent experiments were performed for each treatment condition. Seedlings from all the four groups were sampled simultaneously and were immediately stored at −80 °C until required for RNA extraction.

3.2. RNA-Seq

Total RNA from each sample was isolated using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. The quality of the RNA samples was confirmed using a NanoDrop ND-8000 spectrophotometer (NanoDrop, Wilmington, DE, USA); the A260/A280 ratios of all of the samples ranged from 1.8 to 2.0. The integrity of the RNA samples was assessed with an Agilent 2100 Bioanalyzer; no sign of degradation was found. The RNA samples (mixed RNA of seedlings from three independent experiments according to 1:1 proportion) was sent to the Bionova Biotech (Beijing, China) for RNA sequencing. cDNA library construction and Illumina sequencing were performed according to the manufacturer’s instructions (Illumina, San Diego, CA, USA). Three technical replicates were performed per sample. The RNA-seq read data were deposited in the NCBI Sequence Read Archive (NCBI SRA) under the accession number SRX2637828.

3.3. Transcriptome Data Processing

The clean reads were obtained from the raw data by removing the adaptor sequences, reads with ambiguous bases “N”, low quality reads (Q value < 10), and fragments smaller than 20 bp in length. The de novo assembly of the D. spiculifolius transcriptome in the absence of a reference genome was accomplished using Trinity [41]. Trinity combines the reads with a certain length of overlap to form longer fragments, known as contigs. These contigs were subjected to sequence clustering to form longer sequences. Such sequences were defined as unigenes.
All of the assembled unigenes were compared with the public protein databases, including Swiss-Prot (http://www.expasy.ch/sprot), PFAM (Protein family, http://www.sanger.ac.uk/Software/ Pfam/), GO (Gene Ontology, http://www.geneontology.org), KOG (euKaryotic Ortholog Groups, http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/COG/), and KO (KEGG Ortholog, http://www.genome.jp/kegg/), using BLAST2GO analysis with a cut-off E-value of 10−5. After obtaining GO annotations for the assembled unigenes, WEGO software [42] was used to determine GO functional classifications for understanding the distribution of unigene functions. The unigenes were aligned to the COG database to predict and classify their possible functions. In the final step, KEGG pathway [43] annotations were performed to analyze the metabolic pathways and functions of the unigenes.

3.4. Screening of Differentially Expressed Genes (DEGs)

The expression levels of the unigenes were calculated using the RPKM (Reads per kilobase per Million reads) method [44]. The RPKM method eliminated the influence of the length and sequencing level of the gene. The formula used for determining the RPKM value was the following:
RPKM = 10 6 C N L / 10 3
where C is the number of reads that uniquely aligned to one unigene; N is the total number of reads that uniquely aligned to all the unigenes; and L is the number of bases in the CDS of one unigene. All of the DEGs were mapped to each term of the GO and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases, and significant pathways were defined based on a corrected p-value ≤ 0.05. The DEGs were screened with an false discovery rate (FDR) threshold of 0.05 or less and an absolute log2 ratio of 1 or more.

3.5. Heatmap and Protein–Protein Interaction (PPI) Analysis of DEGs

A heat map demonstrating the gene expression data was created using Java TreeView [45]. For the PPI analysis, the DEGs were used to retrieve the interacting genes with STRING 9.1 [46]. All the interactions in STRING were provided with a probabilistic confidence score (combined score). The PPI network was constructed using STRING and was visualized using Cytoscape. The proteins in the network served as nodes and the degree of a node was assessed as the number of interactions with other proteins. The proteins with high degrees were considered as the hub nodes.

3.6. Quantitative Real-Time PCR of DEGs

Sixteen DEGs were randomly selected from among the total of DEGs and were used to validate the reliability of the libraries. The primers for these assays (Table S3) were designed using Primer 5.0 software and quantitative real-time PCR (qRT-PCR) was carried out as described by Ren et al. [47]. The reference sequence was DsActin7 (Ds66780). Three biological and three technical replicates were performed per sample.

4. Conclusions

In this study, we provide the first report of the transcriptome data of D. spiculifolius under cold, simulated drought, or combined cold and simulated drought stresses using the Illumina HiSeq 2000 platform. A large number of DEGs involved in cold and drought stresses were identified. These DEGs and their functional annotations provide a valuable resource for genetic and genomic studies in D. spiculifolius plants. Functional characterization of the DEGs and their pathways needs to be performed, as they could potentially be used as targets for marker-assisted selection or transgenics to enhance stress tolerance.

Supplementary Materials

Supplementary materials can be found at www.mdpi.com/1422-0067/18/4/849/s1.

Acknowledgments

This work was supported by grants from the National Key Research and Development Program of China (No. 2016YFC0500300), the Natural Science Foundation of Heilongjiang Province of China (No. C2016024), and the ‘Young Talent’ Project of Northeast Agricultural University of China (No. 16QC06).

Author Contributions

Aimin Zhou conceived the study, coordinated and performed experiments, and prepared the manuscript. Hongping Ma, Enhui Liu, Tongtong Jiang, and Shuang Feng performed data analysis. Shufang Gong and Jingang Wang analyzed the data with suggestions from Aimin Zhou and wrote the manuscript. All the authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chinnusamy, V.; Zhu, J.; Zhu, J.K. Cold stress regulation of gene expression in plants. Trends Plant Sci. 2007, 12, 444–451. [Google Scholar] [CrossRef] [PubMed]
  2. Huang, G.T.; Ma, S.L.; Bai, L.P.; Zhang, L.; Ma, H.; Jia, P.; Liu, J.; Zhong, M.; Guo, Z.F. Signal transduction during cold, salt, and drought stresses in plants. Mol. Biol. Rep. 2012, 39, 969–987. [Google Scholar] [CrossRef] [PubMed]
  3. Su, L.; Dai, Z.; Li, S.; Xin, H. A novel system for evaluating drought-cold tolerance of grapevines using chlorophyll fluorescence. BMC Plant Biol. 2015, 15, 82. [Google Scholar] [CrossRef] [PubMed]
  4. Qu, Y.; Zhou, A.; Zhang, X.; Tang, H.; Liang, M.; Han, H.; Zuo, Y. De novo transcriptome sequencing of low temperature-treated Phlox subulata and analysis of the genes involved in cold stress. Int. J. Mol. Sci. 2015, 16, 9732–9748. [Google Scholar] [CrossRef] [PubMed]
  5. Lu, X.; Zhou, X.; Cao, Y.; Zhou, M.; McNeil, D.; Liang, S.; Yang, C. RNA-seq analysis of cold and drought responsive transcriptomes of Zea mays ssp. mexicana L. Front. Plant Sci. 2017, 8, 136. [Google Scholar] [CrossRef] [PubMed]
  6. Calzadilla, P.I.; Maiale, S.J.; Ruiz, O.A.; Escaray, F.J. Transcriptome response mediated by cold stress in Lotus japonicus. Front. Plant Sci. 2016, 7, 374. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, L.; Zhang, F.; Zhang, F.; Wang, W.; Zhou, Y.; Fu, B.; Li, Z. Comparative transcriptome sequencing of tolerant rice introgression line and its parents in response to drought stress. BMC Genom. 2014, 15, 1026. [Google Scholar] [CrossRef] [PubMed]
  8. Tian, D.Q.; Pan, X.Y.; Yu, Y.M.; Wang, W.Y.; Zhang, F.; Ge, Y.Y.; Shen, X.L.; Shen, F.Q.; Liu, X.J. De novo characterization of the Anthurium transcriptome and analysis of its digital gene expression under cold stress. BMC Genom. 2013, 14, 827. [Google Scholar] [CrossRef] [PubMed]
  9. Pang, T.; Ye, C.Y.; Xia, X.; Yin, W. De novo sequencing and transcriptome analysis of the desert shrub, Ammopiptanthus mongolicus, during cold acclimation using Illumina/Solexa. BMC Genom. 2013, 14, 488. [Google Scholar] [CrossRef] [PubMed]
  10. An, D.; Yang, J.; Zhang, P. Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress. BMC Genom. 2012, 13, 64. [Google Scholar] [CrossRef] [PubMed]
  11. Kumar, D.; Hazra, S.; Datta, R.; Chattopadhyay, S. Transcriptome analysis of Arabidopsis mutants suggests a crosstalk between ABA, ethylene and GSH against combined cold and osmotic stress. Sci. Rep. 2016, 6, 36867. [Google Scholar] [CrossRef] [PubMed]
  12. Li, X.W.; Feng, Z.G.; Yang, H.M.; Zhu, X.P.; Liu, J.; Yuan, H.Y. A novel cold-regulated gene from Camellia sinensis, CsCOR1, enhances salt- and dehydration-tolerance in tobacco. Biochem. Biophys. Res. Commun. 2010, 394, 354–359. [Google Scholar] [CrossRef] [PubMed]
  13. Masand, S.; Yadav, S.K. Overexpression of MuHSP70 gene from Macrotyloma uniflorum confers multiple abiotic stress tolerance in transgenic Arabidopsis thaliana. Mol. Biol. Rep. 2016, 43, 53–64. [Google Scholar] [CrossRef] [PubMed]
  14. Zheng, C.; Wang, Y.; Ding, Z.; Zhao, L. Global transcriptional analysis reveals the complex relationship between tea quality, leaf senescence and the responses to cold-drought combined stress in Camellia sinensis. Front. Plant Sci. 2016, 7, 1858. [Google Scholar] [CrossRef] [PubMed]
  15. Holobiuc, M.; Blindu, R.; Mitoi, M.; Helepciuc, F.; Cristea, V. The establishment of an in vitro gene bank in Dianthus spiculifolius Schur. and D. glacialis ssp. gelidus (Schott Nym. et Kotschy) Tutin: I. The initiation of a tissue collection and the characterization of the cultures in minimal growth conditions. Ann. For. Res. 2009, 52, 117–128. [Google Scholar]
  16. Harris, M.A.; Clark, J.; Ireland, A.; Lomax, J.; Ashburner, M.; Foulger, R.; Eilbeck, K.; Lewis, S.; Marshall, B.; Mungall, C.; et al. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32, D258–D261. [Google Scholar] [PubMed]
  17. Valdes, A.E.; Overnas, E.; Johansson, H.; Rada-Iglesias, A.; Engstrom, P. The homeodomain-leucine zipper (HD-Zip) class I transcription factors ATHB7 and ATHB12 modulate abscisic acid signalling by regulating protein phosphatase 2C and abscisic acid receptor gene activities. Plant Mol. Biol. 2012, 80, 405–418. [Google Scholar] [CrossRef] [PubMed]
  18. Danquah, A.; de Zelicourt, A.; Colcombet, J.; Hirt, H. The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol. Adv. 2014, 32, 40–52. [Google Scholar] [CrossRef] [PubMed]
  19. Eremina, M.; Rozhon, W.; Poppenberger, B. Hormonal control of cold stress responses in plants. Cell. Mol. Life Sci. CMLS 2016, 73, 797–810. [Google Scholar] [CrossRef] [PubMed]
  20. Luan, S. The CBL-CIPK network in plant calcium signaling. Trends Plant Sci. 2009, 14, 37–42. [Google Scholar] [CrossRef] [PubMed]
  21. Yu, Q.; An, L.; Li, W. The CBL-CIPK network mediates different signaling pathways in plants. Plant Cell Rep. 2014, 33, 203–214. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, T.; Chaudhuri, S.; Yang, L.; Du, L.; Poovaiah, B.W. A calcium/calmodulin-regulated member of the receptor-like kinase family confers cold tolerance in plants. J. Biol. Chem. 2010, 285, 7119–7126. [Google Scholar] [CrossRef] [PubMed]
  23. Winfield, M.O.; Lu, C.; Wilson, I.D.; Coghill, J.A.; Edwards, K.J. Plant responses to cold: Transcriptome analysis of wheat. Plant Biotechnol. J. 2010, 8, 749–771. [Google Scholar] [CrossRef] [PubMed]
  24. Kim, M.C.; Chung, W.S.; Yun, D.J.; Cho, M.J. Calcium and calmodulin-mediated regulation of gene expression in plants. Mol. Plant 2009, 2, 13–21. [Google Scholar] [CrossRef] [PubMed]
  25. Sanchez-Barrena, M.J.; Martinez-Ripoll, M.; Albert, A. Structural Biology of a Major Signaling Network that Regulates Plant Abiotic Stress: The CBL-CIPK Mediated Pathway. Int. J. Mol. Sci. 2013, 14, 5734–5749. [Google Scholar] [CrossRef] [PubMed]
  26. Lan, W.Z.; Lee, S.C.; Che, Y.F.; Jiang, Y.Q.; Luan, S. Mechanistic analysis of AKT1 regulation by the CBL-CIPK-PP2CA interactions. Mol. Plant 2011, 4, 527–536. [Google Scholar] [CrossRef] [PubMed]
  27. Soon, F.F.; Ng, L.M.; Zhou, X.E.; West, G.M.; Kovach, A.; Tan, M.H.; Suino-Powell, K.M.; He, Y.; Xu, Y.; Chalmers, M.J.; et al. Molecular mimicry regulates ABA signaling by SnRK2 kinases and PP2C phosphatases. Science 2012, 335, 85–88. [Google Scholar] [CrossRef] [PubMed]
  28. Novillo, F.; Medina, J.; Salinas, J. Arabidopsis CBF1 and CBF3 have a different function than CBF2 in cold acclimation and define different gene classes in the CBF regulon. Proc. Natl. Acad. Sci. USA 2007, 104, 21002–21007. [Google Scholar] [CrossRef] [PubMed]
  29. Shavrukov, Y.; Baho, M.; Lopato, S.; Langridge, P. The TaDREB3 transgene transferred by conventional crossings to different genetic backgrounds of bread wheat improves drought tolerance. Plant Biotechnol. J. 2016, 14, 313–322. [Google Scholar] [CrossRef] [PubMed]
  30. Agarwal, M.; Hao, Y.; Kapoor, A.; Dong, C.H.; Fujii, H.; Zheng, X.; Zhu, J.K. A R2R3 type MYB transcription factor is involved in the cold regulation of CBF genes and in acquired freezing tolerance. J. Biol. Chem. 2006, 281, 37636–37645. [Google Scholar] [CrossRef] [PubMed]
  31. Licausi, F.; Ohme-Takagi, M.; Perata, P. APETALA2/ethylene responsive factor (AP2/ERF) transcription factors: Mediators of stress responses and developmental programs. New Phytol. 2013, 199, 639–649. [Google Scholar] [CrossRef] [PubMed]
  32. Gangappa, S.N.; Botto, J.F. The BBX family of plant transcription factors. Trends Plant Sci. 2014, 19, 460–470. [Google Scholar] [CrossRef] [PubMed]
  33. Yamaguchi-Shinozaki, K.; Shinozaki, K. Transcriptional regulatory networks in cellular responses and tolerance to dehydration and cold stresses. Ann. Rev. Plant Biol. 2006, 57, 781–803. [Google Scholar] [CrossRef] [PubMed]
  34. Yoshida, T.; Mogami, J.; Yamaguchi-Shinozaki, K. ABA-dependent and ABA-independent signaling in response to osmotic stress in plants. Curr. Opin. Plant Biol. 2014, 21, 133–139. [Google Scholar] [CrossRef] [PubMed]
  35. Zhou, M.Q.; Shen, C.; Wu, L.H.; Tang, K.X.; Lin, J. CBF-dependent signaling pathway: A key responder to low temperature stress in plants. Crit. Rev. Biotechnol. 2011, 31, 186–192. [Google Scholar] [CrossRef] [PubMed]
  36. Qin, F.; Kakimoto, M.; Sakuma, Y.; Maruyama, K.; Osakabe, Y.; Tran, L.S.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Regulation and functional analysis of ZmDREB2A in response to drought and heat stresses in Zea mays L. Plant J. 2007, 50, 54–69. [Google Scholar] [CrossRef] [PubMed]
  37. Tunnacliffe, A.; Wise, M.J. The continuing conundrum of the LEA proteins. Naturwissenschaften 2007, 94, 791–812. [Google Scholar] [CrossRef] [PubMed]
  38. Kosova, K.; Vitamvas, P.; Prasil, I.T. Wheat and barley dehydrins under cold, drought, and salinity—What can LEA-II proteins tell us about plant stress response? Front. Plant Sci. 2014, 5, 343. [Google Scholar] [CrossRef] [PubMed]
  39. Hu, T.Z.; Zhu, S.S.; Tan, L.L.; Qi, W.H.; He, S.; Wang, G.X. Overexpression of OsLEA4 enhances drought, high salt and heavy metal stress tolerance in transgenic rice (Oryza sativa L.). Environ. Exp. Bot. 2016, 123, 68–77. [Google Scholar] [CrossRef]
  40. Yu, J.; Lai, Y.; Wu, X.; Wu, G.; Guo, C. Overexpression of OsEm1 encoding a group I LEA protein confers enhanced drought tolerance in rice. Biochem. Biophys. Res. Commun. 2016, 478, 703–709. [Google Scholar] [CrossRef] [PubMed]
  41. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  42. Ye, J.; Fang, L.; Zheng, H.; Zhang, Y.; Chen, J.; Zhang, Z.; Wang, J.; Li, S.; Li, R.; Bolund, L.; et al. WEGO: A web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34, W293–W297. [Google Scholar] [CrossRef] [PubMed]
  43. Kanehisa, M.; Goto, S.; Sato, Y.; Furumichi, M.; Tanabe, M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40, D109–D114. [Google Scholar] [CrossRef] [PubMed]
  44. Mortazavi, A.; Williams, B.A.; McCue, K.; Schaeffer, L.; Wold, B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods 2008, 5, 621–628. [Google Scholar] [CrossRef] [PubMed]
  45. Saldanha, A.J. Java Treeview—Extensible visualization of microarray data. Bioinformatics 2004, 20, 3246–3248. [Google Scholar] [CrossRef] [PubMed]
  46. Franceschini, A.; Szklarczyk, D.; Frankild, S.; Kuhn, M.; Simonovic, M.; Roth, A.; Lin, J.; Minguez, P.; Bork, P.; von Mering, C.; et al. STRING v9.1: Protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013, 41, D808–D815. [Google Scholar] [CrossRef] [PubMed]
  47. Ren, D.; Madsen, J.S.; de la Cruz-Perera, C.I.; Bergmark, L.; Sorensen, S.J.; Burmolle, M. High-throughput screening of multispecies biofilm formation and quantitative PCR-based assessment of individual species proportions, useful for exploring interspecific bacterial interactions. Microb. Ecol. 2014, 68, 146–154. [Google Scholar] [CrossRef] [PubMed]
Ijms 18 00849 g001 550
Figure 1. Open field cultivation and phenotypic characteristics of D. spiculifolius. (A–C) Open field cultivation of D. spiculifolius in a dry summer (A), cold winter (B) and autumn (C); (D,E) Phenotype of D. spiculifolius (D) and its flower type (E). (F) Seeds of D. spiculifolius; (G–I) Seed germination (G) and seedling growth (H,I) of D. spiculifolius.
Figure 1. Open field cultivation and phenotypic characteristics of D. spiculifolius. (A–C) Open field cultivation of D. spiculifolius in a dry summer (A), cold winter (B) and autumn (C); (D,E) Phenotype of D. spiculifolius (D) and its flower type (E). (F) Seeds of D. spiculifolius; (G–I) Seed germination (G) and seedling growth (H,I) of D. spiculifolius.
Ijms 18 00849 g001
Ijms 18 00849 g002 550
Figure 2. Venn diagrams showing the number of differentially expressed genes (DEGs). The numbers of up-regulated (A) and down-regulated (B) genes in either cold-, drought-, or combined cold and drought-treated samples of D. spiculifolius seedling (compared to the control non-stressed plants) are shown. We used “false discovery rate ≤0.05 and the absolute value of log2Ratio ≥1” as the threshold to judge the significance of the differences in gene expression.
Figure 2. Venn diagrams showing the number of differentially expressed genes (DEGs). The numbers of up-regulated (A) and down-regulated (B) genes in either cold-, drought-, or combined cold and drought-treated samples of D. spiculifolius seedling (compared to the control non-stressed plants) are shown. We used “false discovery rate ≤0.05 and the absolute value of log2Ratio ≥1” as the threshold to judge the significance of the differences in gene expression.
Ijms 18 00849 g002
Ijms 18 00849 g003 550
Figure 3. Histogram of gene ontology (GO) classification of DEGs. (A) GO classification of up-regulated genes. (B) GO classification of down-regulated genes. The results are summarized in three main categories: biological processes (green), cellular component (yellow) and molecular function (blue). The y-axis on the left side indicates the number of DEGs.
Figure 3. Histogram of gene ontology (GO) classification of DEGs. (A) GO classification of up-regulated genes. (B) GO classification of down-regulated genes. The results are summarized in three main categories: biological processes (green), cellular component (yellow) and molecular function (blue). The y-axis on the left side indicates the number of DEGs.
Ijms 18 00849 g003
Ijms 18 00849 g004 550
Figure 4. PPI sub-network involving the DEGs. The square nodes represent GO terms; circular nodes represent genes; star-like nodes represent transcription factors. Gray lines indicate interactions between two proteins. The intensity of the color of the nodes is based on the p-value of the GO term or fold change in the DEGs. All information for each gene can be found in Tables S1 and S2.
Figure 4. PPI sub-network involving the DEGs. The square nodes represent GO terms; circular nodes represent genes; star-like nodes represent transcription factors. Gray lines indicate interactions between two proteins. The intensity of the color of the nodes is based on the p-value of the GO term or fold change in the DEGs. All information for each gene can be found in Tables S1 and S2.
Ijms 18 00849 g004
Ijms 18 00849 g005 550
Figure 5. Expression profiles of 58 DEGs. The bar represents the scale of the expression levels of each gene (log2 RPKM) in the heat map. The red rectangles represent the up-regulation of genes, and the green rectangles represent down-regulation. All information for each gene can be found in Tables S1 and S2.
Figure 5. Expression profiles of 58 DEGs. The bar represents the scale of the expression levels of each gene (log2 RPKM) in the heat map. The red rectangles represent the up-regulation of genes, and the green rectangles represent down-regulation. All information for each gene can be found in Tables S1 and S2.
Ijms 18 00849 g005
Ijms 18 00849 g006 550
Figure 6. The relative expression levels of representative DEGs from CK, CT, DT, or CTDT samples. The DsActin gene was used as an internal control, and the transcript level in non-stressed seedlings was set as 1.0. Error bars represent SE (n = 3).
Figure 6. The relative expression levels of representative DEGs from CK, CT, DT, or CTDT samples. The DsActin gene was used as an internal control, and the transcript level in non-stressed seedlings was set as 1.0. Error bars represent SE (n = 3).
Ijms 18 00849 g006
Table
Table 1. Summary of Illumina transcriptome sequencing of D. spiculifolius.
Table 1. Summary of Illumina transcriptome sequencing of D. spiculifolius.
ContentControl (CK)Cold-Treated (CT)Drought-Treated (DT)Combined-Treated (CTDT)
Total reads30,100,66629,709,67730,182,57229,577,457
Total clean reads23,624,69622,495,62823,317,14122,965,345
Guanine and cytosine (GC) percentage46.4546.3246.6246.10
Total number of transcripts184,782
Total number of unigenes112,015
Mean length of unigenes (bp)819
Table
Table 2. Number of functional annotations for all of the unigenes in public databases.
Table 2. Number of functional annotations for all of the unigenes in public databases.
Annotated DatabaseNumber of UnigenesPercentage (%)
SwissProt36,14232.27%
Protein family (PFAM)29,14926.55%
Gene Ontology (GO)29,73926.54%
euKaryotic Ortholog Groups (KOG)25,07322.38%
KEGG Ortholog (KO)10,8459.68%

Share and Cite

MDPI and ACS Style

Zhou, A.; Ma, H.; Liu, E.; Jiang, T.; Feng, S.; Gong, S.; Wang, J. Transcriptome Sequencing of Dianthus spiculifolius and Analysis of the Genes Involved in Responses to Combined Cold and Drought Stress. Int. J. Mol. Sci. 2017, 18, 849. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms18040849

AMA Style

Zhou A, Ma H, Liu E, Jiang T, Feng S, Gong S, Wang J. Transcriptome Sequencing of Dianthus spiculifolius and Analysis of the Genes Involved in Responses to Combined Cold and Drought Stress. International Journal of Molecular Sciences. 2017; 18(4):849. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms18040849

Chicago/Turabian Style

Zhou, Aimin, Hongping Ma, Enhui Liu, Tongtong Jiang, Shuang Feng, Shufang Gong, and Jingang Wang. 2017. "Transcriptome Sequencing of Dianthus spiculifolius and Analysis of the Genes Involved in Responses to Combined Cold and Drought Stress" International Journal of Molecular Sciences 18, no. 4: 849. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms18040849

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop