Next Article in Journal
The Importance of G Protein-Coupled Receptor Kinase 4 (GRK4) in Pathogenesis of Salt Sensitivity, Salt Sensitive Hypertension and Response to Antihypertensive Treatment
Next Article in Special Issue
De Novo Transcriptome Sequencing of Low Temperature-Treated Phlox subulata and Analysis of the Genes Involved in Cold Stress
Previous Article in Journal
Whole Exome Sequencing for a Patient with Rubinstein-Taybi Syndrome Reveals de Novo Variants besides an Overt CREBBP Mutation
Previous Article in Special Issue
A Ribosomal Protein AgRPS3aE from Halophilic Aspergillus glaucus Confers Salt Tolerance in Heterologous Organisms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 16
Issue 3
10.3390/ijms16035714
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

Genome-Wide Identification and Analysis of Drought-Responsive Genes and MicroRNAs in Tobacco

by
Fuqiang Yin
1,2,†,
Cheng Qin
2,3,†,
Jian Gao
2,†,
Ming Liu
1,
Xirong Luo
3,
Wenyou Zhang
1,
Hongjun Liu
2,
Xinhui Liao
4,
Yaou Shen
2,
Likai Mao
4,
Zhiming Zhang
2,
Haijian Lin
2,
Thomas Lübberstedt
5 and
Guangtang Pan
2,*
1
School of Agricultural Sciences, Xichang College, Xichang 615000, China
2
Maize Research Institute of Sichuan Agricultural University/Key Laboratory of Biology and Genetic Improvement of Maize in Southwest Region, Ministry of Agriculture, Chengdu 611130, China
3
Zunyi Academy of Agricultural Sciences, Zunyi 563102, China
4
Beijing Genomics Institute-Shenzhen, Shenzhen 518083, China
5
Department of Agronomy, Iowa State University, Ames, IA 50011, USA
*
Author to whom correspondence should be addressed.
†
These authors are equally contributed to this work.
Int. J. Mol. Sci. 2015, 16(3), 5714-5740; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms16035714
Submission received: 7 December 2014 / Revised: 19 January 2015 / Accepted: 29 January 2015 / Published: 12 March 2015
(This article belongs to the Special Issue Abiotic Stress and Gene Networks in Plants)
Browse Figures
Versions Notes

Abstract

:
Drought stress response is a complex trait regulated at transcriptional and post-transcriptional levels in tobacco. Since the 1990s, many studies have shown that miRNAs act in many ways to regulate target expression in plant growth, development and stress response. The recent draft genome sequence of Nicotiana benthamiana has provided a framework for Digital Gene Expression (DGE) and small RNA sequencing to understand patterns of transcription in the context of plant response to environmental stress. We sequenced and analyzed three Digital Gene Expression (DGE) libraries from roots of normal and drought-stressed tobacco plants, and four small RNA populations from roots, stems and leaves of control or drought-treated tobacco plants, respectively. We identified 276 candidate drought responsive genes (DRGs) with sequence similarities to 64 known DRGs from other model plant crops, 82 were transcription factors (TFs) including WRKY, NAC, ERF and bZIP families. Of these tobacco DRGs, 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. An integrated gene regulatory network is proposed for the molecular mechanisms of tobacco root response to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and their target transcripts. This network analysis serves as a reference for future studies on tobacco response stresses such as drought, cold and heavy metals.

Graphical Abstract

1. Introduction

Tobacco, Nicotiana benthamiana, is an important agricultural and economic crop in China [1]. It is one of the most commonly used species to study molecular plant-microbe interactions [2,3]. Numerous efforts are underway to improve plant productivity, including reduction of yield loss to environmental stress, such as drought, water logging, salt loading, and freezing. These stresses have adverse effects on plant growth, and plant stress responses are regulated by multiple signaling pathways [4,5]. In general, plants respond and adapt to stress through various biochemical and physiological processes, thereby acquiring stress tolerance [6].
Decades of research into the effects of drought on plant physiology and development have generated a plethora of information. It has been shown that regulatory networks of gene expression under drought and cold stress are active [6], and that particular transcription factors are employed to enhance drought tolerance in plants [7]. Drought stress increases endogenous abscisic acid (ABA) levels and induces ABA-dependent and ABA-independent transcriptional regulatory networks [8]. Many of these responses can be mimicked by external application of ABA or polyethylene glycol (PEG) [6,9]. Drought-responsive transcriptional networks have been primarily developed from related studies in Arabidopsis, and some important promoter elements were confirmed experimentally, such as ABA-responsive elements (ABREs) and coupling elements (CE) [10,11,12]. However, the networks that underlie these responses in tobacco have not been fully characterized.
In addition to drought responsive genes (DRGs) and those transcription factors mentioned, microRNAs (miRNAs) are short (20–22 nt), endogenously expressed, non-translated RNAs that function in posttranscriptional gene regulation. The mature miRNA form is processed from a larger primary transcript (pri-miRNA) [13] and incorporated into the RNA-induced silencing complex (RISC). RISC directs binding to complementary mRNA transcripts, followed by target cleavage or reversible inhibition of translation [14]. MiRNAs act in a variety of ways and they can spatially restrict, temporally regulate, dampen or mutually exclude the expression of their target gene [14]. MiRNAs play critical roles in different developmental pathways, including root [15,16] and shoot development [17]. Additionally, they respond to plant hormones [18], and a variety of stresses like drought, salt, cold and oxidative stress [19,20,21,22]. For instance, ABA treatment and drought stress induces miR159 to accumulate and target MYB TFs that positively regulate ABA responses during seed germination in Arabidopsis, with miR159 as a negative feedback regulator of ABA responses [23]. Zhao et al. [21] reported that miRNA169g accumulates during drought stress and miR169g is regulated by dehydration responsive element (DREs) in rice, suggesting this miRNA plays a role in drought stress. There is a complex interplay between transcriptional and posttranscriptional regulation of drought response in Arabidopsis, but this has not been extensively characterized in tobacco.
The recent draft genome sequence of Nicotiana benthamiana [3] has provided a framework for the identification and functional characterization of genes and genetic networks for tobacco crop improvement and basic research. The availability of next-generation sequencing technologies, such as Digital Gene Expression (DGE) [24] and small RNA sequencing [25,26], together with the genome sequence, offer the opportunity to understand patterns of transcription in the context of plant responses to environmental stresses. Although several studies have documented individual genes that responded under drought-stress treatment of tobacco [27,28,29,30], a global study of transcription patterns of single or diverse organs using a single platform is lacking.
The root is generally the first organ in terrestrial plants that comes into contact with the drought stress. In this study, we used DGE [24,31], small RNA sequencing [32], and computational approaches to develop and map the transcriptional and posttranscriptional gene regulatory network that operates in response to drought stress in tobacco root tissue. This network analysis serves as a reference for future studies on tobacco response to various stresses, such as to drought, cold and heavy metals.

2. Results

2.1. Physiological Analysis of PRO, SOD, and MDA under Drought Conditions

To determine the optimal time to analysis gene expression for drought stress, uniform seedlings of tobacco with six leaves were challenged with 20% PEG6000 to simulate drought and then sampled at six durations of stress (3, 6, 12, 24, 48, and 96 h) (see “Experimental Section” for details). As shown in Figure 1A, the plants grew regularly under normal condition but leaves gradually turned yellow under longer duration of stress. After 48 h, plants were seriously damaged, and after 96 h they were wilted and dying. Superoxide dismutase (SOD) activities increased significantly (p < 0.01), and reached their peak after 48 h at about four times the control level (Figure 1B). This increase was paralleled by increased content of proline (PRO) and malondialdehyde (MDA) (Figure 1C), indicating that these three components in tobacco roots may act synergistically under drought stress. The general trend of the activity of SOD was different from that of PRO and MDA contents in that the former reached its peak at 48 h and then leveled off. PRO and MDA continued to ascend slowly and reached their peak after 96 h. These results suggest that SOD activity, PRO and MDA content all significantly increased at 6 and 48 h relative to 3 and 24 h, and the optimal times for drought stress assays were 6 and 48 h. Accordingly, our experiments of tobacco roots at 0 h (control) and two drought-treatment times were designated as NCK, N6H, and N48H.
Ijms 16 05714 g001 1024
Figure 1. The growth status of tobacco under drought stress and corresponding physiological analyses. (A) The growth status of control plant (0 h) and drought-treatment plants at six time-points (3, 6, 12, 24, 48, and 96 h) with 20% PEG6000. The graphs (B) and (C) depicted the average change of superoxide dismutase (SOD) activities and contents of proline (PRO) and malondialdehyde (MDA), respectively, in the actual value ± standard error of the mean from control (time 0) across time for each treatment. Different lower or upper cases above standard error denoted the significantly difference of every two time –points (one-way ANOVA and Fisher’s LSD, p < 0.01, n = 3), whereas the same letter indicated no highly significant difference (p < 0.01). The gray bars in graph C represent the contents of PRO (µmol/g) and the black bars, the contents of MDA (µmol/g) at time 0 and six drought-treatment time-points. Time 0 represented the level obtained from the normal seedlings with six leaves.
Figure 1. The growth status of tobacco under drought stress and corresponding physiological analyses. (A) The growth status of control plant (0 h) and drought-treatment plants at six time-points (3, 6, 12, 24, 48, and 96 h) with 20% PEG6000. The graphs (B) and (C) depicted the average change of superoxide dismutase (SOD) activities and contents of proline (PRO) and malondialdehyde (MDA), respectively, in the actual value ± standard error of the mean from control (time 0) across time for each treatment. Different lower or upper cases above standard error denoted the significantly difference of every two time –points (one-way ANOVA and Fisher’s LSD, p < 0.01, n = 3), whereas the same letter indicated no highly significant difference (p < 0.01). The gray bars in graph C represent the contents of PRO (µmol/g) and the black bars, the contents of MDA (µmol/g) at time 0 and six drought-treatment time-points. Time 0 represented the level obtained from the normal seedlings with six leaves.
Ijms 16 05714 g001

2.2. Generation and Quality Assessment of the DGE Dataset

To document gene expression profiles of tobacco roots responding to drought stress, three DGE libraries (NCK, N6H, and N48H) were sequenced and high-quality sequences were obtained. The total sequence reads of all libraries ranged from 3.33 M (million) to 3.39 M (Supplementary Figure S1A and Table S1) with an average of 3.37 M. The N48H library had the highest number of distinct tags (242,065), followed by the NCK (229,344) and N6H (221,248) libraries (Figure S1B). Copy number of most of the distinct tags (over 67%) ranged from 2 to 5 (2 ≤ frequency ≤ 5). However, a small number of distinct tags (less than 2.2%) with frequency higher than 100 covered over 43% of all the clean tags in all three libraries (Supplementary Figure S1B).
When total tag number in NCK reached 1 million, the number of identified genes started to level off, and stabilized when the number of tags reached 3 million (Supplementary Figure S2A). The N6H and N48H data show a similar trend (Supplementary Figure S2A). This suggests additional distinct genes would not be identified after total clean tag number reached a certain value. Moreover, after sequencing depths reached 1 million tags, the number of distinct tags discovered dropped dramatically in all three libraries (Supplementary Figure S2B). In other words, increasing sequence depth resulted in a slow and stable accumulation of new distinct tags up to a point of saturation.

2.3. Mapping DGE Tags to the Tobacco Reference Genome Database

We used SOAP2 software [33] to map all distinct tags to the draft sequence of the Nicotiana benthamiana genome (v.0.4.4) [3], and the clean tags that matched perfectly or with only one mismatch were analyzed further. Mapping results showed that 73.40%, 72.81% and 70.80% of total clean tags mapped to the reference database (sense or anti-sense) for NCK, N6H and N48H, respectively, and 63.20%, 60.76% and 62.82% of reference genes were identified in NCK, N6H and N48H, respectively (Supplementary Table S1). In addition, around 30% of reads could not be mapped to the whole genome in three samples. In total, we detected 13,101 genes expressed in all three samples and 21,128 genes collectively (73.6% of the annotated transcriptome for tobacco) (Supplementary Figure S3).

2.4. Expression of Tobacco Root Genes under Drought Stress

We examined the dynamics of gene expression under drought stress using DGE data, and only 1887 out of 21,128 genes that were differentially expressed among NCK, N6H, and N48H, represented 8.9% of the root transcriptome (Supplementary Table S2). A sample of 22 genes with significant differences in gene expression was randomly selected for validation via qRT-PCR. For 10 of these genes, the expression dynamics of genes detected in the qRT-PCR experiments were indistinguishable from that based on three DGE data, and other 10 genes at least consistent with that obtained from these DGE data (Figure 2).
Ijms 16 05714 g002 1024
Figure 2. Expression of 22 genes in control and drought-treated roots of tobacco. (A) The levels of 22 genes detected by Solexa sequencing in three Digital Gene Expression (DGE) libraries (NCK, N6H and N48H); (B) The levels of these 22 genes detected by quantitative Real-Time PCR (qRT-PCR) in NCK, N6H and N48H (normalized to actin; n = 3). The green, purple and blue bars in graph B depicted the qRT-PCR relative expression level ± standard error of three replicates for each gene in NCK, N6H and N48H, respectively. NCK, N6H and N48H represented three samples of tobacco roots at 0 h (control, without drought-treatment) and two drought-treatment time-points (6 and 48 h), respectively.
Figure 2. Expression of 22 genes in control and drought-treated roots of tobacco. (A) The levels of 22 genes detected by Solexa sequencing in three Digital Gene Expression (DGE) libraries (NCK, N6H and N48H); (B) The levels of these 22 genes detected by quantitative Real-Time PCR (qRT-PCR) in NCK, N6H and N48H (normalized to actin; n = 3). The green, purple and blue bars in graph B depicted the qRT-PCR relative expression level ± standard error of three replicates for each gene in NCK, N6H and N48H, respectively. NCK, N6H and N48H represented three samples of tobacco roots at 0 h (control, without drought-treatment) and two drought-treatment time-points (6 and 48 h), respectively.
Ijms 16 05714 g002
We used Gene Ontology annotation to assign genes to functional categories and grouped genes by expression dynamics using the K-Means clustering algorithm [34]. We identified six clusters (Figure 3A,B and Supplementary Table S2). Genes that encode enzymes for fatty acid metabolism, transferase (transferring acyl groups), oxidoreductase, ethanol metabolism, primary alcohol metabolism and transferase were greatly enriched in cluster K1, including 426 genes, representing genes that were expressed at the highest levels in N48H. Genes with peak expression in N6H (clusters K3, 206 genes) were found to be required for transferase activity, carbohydrate catabolic process, cellular lipid metabolism, hydrolase activity, isoprenoid metabolic process, lipid metabolism, organelle membrane biosynthesis, aminoglycan metabolic process, and polysaccharide catabolic process, which suggests the root transcriptome undergoes drought stress first. Root gene expression was lowest in N48H (clusters K4, 539 genes) and was mostly restricted to genes involved in lipid metabolic processes. Genes encoding enzymes for polyamine metabolism, flavin containing compounds, polysaccharide, riboflavin, vitamin and water soluble vitamins showed a decline in expression after 6 h but reversed direction and their expressions were greatly enriched after 48 h (clusters K5, 275 genes).
Ijms 16 05714 g003 1024
Figure 3. The dynamic profiles of all differentially expressed genes in tobacco roots on response to drought stress. (A) K-means clustering showing the expression profile of 1887 differentially expressed genes. Six clusters were identified in NCK, N6H and N48H from these genes above. These clusters are presented in A. Error bars show standard deviation; (B) Analysis of functional category enrichment among the six major clusters; (C) Analysis of pathways enrichment among the six major clusters.
Figure 3. The dynamic profiles of all differentially expressed genes in tobacco roots on response to drought stress. (A) K-means clustering showing the expression profile of 1887 differentially expressed genes. Six clusters were identified in NCK, N6H and N48H from these genes above. These clusters are presented in A. Error bars show standard deviation; (B) Analysis of functional category enrichment among the six major clusters; (C) Analysis of pathways enrichment among the six major clusters.
Ijms 16 05714 g003
In addition, biological pathways influenced by drought were evaluated by enrichment analysis of all differentially expressed genes. Significantly enriched metabolic pathways and signal transduction pathways were also identified. A total of 17 pathways were affected based on K1 thru K6 clusters (p < 0.05), and differentially expressed genes with pathway annotation were listed according to enrichment priority (Figures 3A,C and Supplementary Table S2). The expression level of genes in cluster K1 was lowest in the control but dramatically increased with time. These genes were associated with glutathione metabolism, fatty acid elongation, stilbenoid, diarylheptanoid and gingerol synthesis, synthesis of secondary metabolites, phenylpropanoid synthesis and arginine and proline metabolism (Figure 3C). Genes in the K2 cluster were mainly annotated as synthesis of secondary metabolites, flavonoid synthesis, stilbenoid, diarylheptanoid and gingerol synthesis, metabolic pathways, phenylpropanoid synthesis, glutathione metabolism, arginine and proline metabolism and sesquiterpenoid and triterpenoid synthesis (Figure 3C). The expression of genes in K3 cluster began to increase after 6 h and leveled off after 48 h. The K3 cluster with most of the enriched biological pathways affected by drought, included synthesis of secondary metabolites, terpenoid backbone synthesis, synthesis and degradation of ketone bodies, metabolic pathways, amino sugar and nucleotide sugar metabolism, phenylalanine, tyrosine and tryptophan synthesis, cysteine and methionine metabolism, diterpenoid synthesis, and isoquinoline alkaloid synthesis (Figure 3C). Although the genes in K4 and K5 showed a very high expression in the control, no significantly enriched pathways were found in K4, K5, and K6 clusters.

2.5. Resolving Transcription Factors (TFs) among Differentially Expressed Genes

A primary objective was to identify genes that encode TFs and elucidate the dynamics of accumulation of TFs under drought stress in our DGE data. To do this, we retrieved putative orthologs of tobacco genes based on information from the EnsemblCompara gene trees [35] at Sol Genomics network (solgenomics.net), PlantGDB (plantgdb.org), and NIBC (www.ncbi.nlm.nih.gov/). We then queried known plant TFs in the Plant Transcription Factor Database (v2.0, http://planttfdb.cbi.edu.cn/) and matched 609 tobacco TFs with sequence similarities to TFs of known plants. All TFs were detected in roots of tobacco seedlings in this study (Figure 4A and Supplementary Table S3A). Of these, 82 TFs, belonging to 24 families, were differentially expressed during stress (Figure 4B and Supplementary Table S3B). These TFs associated with functions in drought tolerance (MYB, NAC, and ERF) [7], development and meristem maintenance or identity (HD-ZIP, NF-YA, NAC, GRAS, and TCP) [36,37], defense/stress signaling pathways (HSP, WRKY, and bZIP) [38], hormone-mediated or stress-mediated signaling by auxin (AUX/IAA), brassinosteroids (BES), or ethylene and stress (AP2/ERF) [38]. Functional investigation of the roles of TFs was analyzed in detail in following sections.
Ijms 16 05714 g004 1024
Figure 4. Dynamics of transcription factors in response to drought stress. (A) Dendrogram of the transcription factors. 609 expressed transcription factors from NCK, N6H and N48H were identified and clustered into three lineages (G1, G2 and G3) by K-means clustering (see “Materials and methods” for details); (B) Distribution of 82 differentially expressed transcription factor families among G1, G2 and G3.
Figure 4. Dynamics of transcription factors in response to drought stress. (A) Dendrogram of the transcription factors. 609 expressed transcription factors from NCK, N6H and N48H were identified and clustered into three lineages (G1, G2 and G3) by K-means clustering (see “Materials and methods” for details); (B) Distribution of 82 differentially expressed transcription factor families among G1, G2 and G3.
Ijms 16 05714 g004
The abundance of most of these TFs (53%) was highest in the control (G1), whereas 10% was highest after 6 h (G2). The reminder (37%) had peak expressions after 48 h (G3). We also identified family-specific expression trends (Figure 4B and Table S3B). Members of the C2H2 (three genes), MYB (3), bHLH (4), WRKY (3), ERF (13) and Dof (3) families of transcriptional regulators were highly expressed in the control. Several GATA (1), MYB (3) and ERF (3) TFs accumulated their highest levels after 3 h stress. Transcriptional regulators including MYB (6), NAC (6), MYB-related (2), NF-YA (2), HD-ZIP (2) and ERF (2) were preferentially expressed after 48 h of drought stress.

2.6. Identifying Candidate Drought Responsive Genes (DRGs) in Tobacco Roots

The 64 DRGs amalgamated from various drought studies on Arabidopsis, rice, tobacco, and maize (Table S4), were selected for this analysis. Most of these were signaling component genes acting upstream of TFs and the transcription factor encoding genes in the regulatory network of gene expression in drought and cold stress responses [6,7]. We identified 276 candidate DRGs in tobacco with sequence similarity to known genes and Nicotiana benthamiana annotation [3] (Table S4). About 40% (110 out of 276 genes) were TFs including WRKY, NAC, ERF, and bZIP families.
In our study, we also investigated the roles of these candidate DRGs, and found 46 differentially expressed DRGs under drought stress (Table 1 and Supplementary Figure S4). Of these, 21 (46%) were TFs which belonged to NAC (6), MYB (4), ERF (10), and bZIP (1) families. Other DRGs, such as GRF6, ABF1, APX2, SIPK, and ZPT2, have different expression patterns in response to drought stress.
Table
Table 1. Possible drought responsive genes (DRGs) with different expression in tobacco roots.
Table 1. Possible drought responsive genes (DRGs) with different expression in tobacco roots.
Gene Name 1FunctionHomolog_Niben_IDPlantGDB_or_unigenes_IDTranfactor_IDFamilyTPM-NCKTPM-N6HTPM-N48H
GRF6 (14-3-3 protein GF14 lambda)Drought toleranceNbS00060175g0001.1---161.7480.67123.31
NbS00044217g0013.1---5.898.2718.9
NbS00020600g0007.1---22.981333.3
NbS00007737g0010.1---158.5277.48121.51
ABF1Drought toleranceNbS00021408g0025.1---0.8818.030.6
ANAC055Drought and salt toleranceNbS00023955g0006.1gnl|UG|Nta#S45443622Nta005482NAC108.7158.5139
NbS00008328g0013.1gnl|UG|Nta#S45488905Nta010847NAC35.0660.2873.51
ANAC072Drought and salt toleranceNbS00007567g0016.1gnl|UG|Nta#S33559198Nta013275NAC5.33.2513.8
NbC26152828g0003.1gnl|UG|Nta#S45463754Nta012057NAC23.5748.7699.01
NbS00025931g0004.1gnl|UG|Nta#S33566495Nta005129NAC7.957.6823.4
NbS00028594g0003.1gnl|UG|Nta#S45443622Nta005482NAC19.1521.8756.7
APX2Drought resistanceNbS00018810g0008.1---0.882.966.6
bZIP1Drought, salt and disease toleranceNbS00007503g0111.1gnl|UG|Nta#S33578426Nta007642bZIP84.5540.4828.5
CBF4Drought and freezing toleranceNbS00018604g0003.1gnl|UG|Nta#S45471037Nta012107ERF10.90.010.01
CpMYB10Desiccation and salinity toleranceNbS00016700g0002.1gnl|UG|Nta#S45484072Nta012264MYB2.3610.051.8
NbS00007512g0018.1---11.4910.343.3
DREB1B/CBF1Increased tolerance to drought, cold and salinityNbS00009125g0101.1gnl|UG|Nta#S33537030Nta005632ERF114.3111.5216.2
DREB2ADrought resistanceNbS00032542g0006.1gnl|UG|Nta#S33577665Nta005987ERF160.5614.7831.5
NbS00025865g0001.1gnl|UG|Nta#S33537030Nta005632ERF13.260.897.5
NbS00002425g0103.1gnl|UG|Nta#S33577665Nta005987ERF5.30.010.6
DREB3Drought toleranceNbS00016353g0001.1---20.033.252.7
FAD3Drought resistanceNbS00037943g0006.1---17.9727.4842
NbS00033197g0010.1---16.29.4629.7
GmERF3Drought toleranceNbS00020925g0011.1gnl|UG|Nta#S33555957Nta004926ERF27.127.7812.6
HSFA2Resistance to environmental stressesNbS00009669g0213.1---0.012.0717.1
MYB60Drought toleranceNbS00046172g0003.1gnl|UG|Nta#S40644996Nta005483MYB1.771.7712.3
NbS00009475g0003.1---5.893.250.6
NbS00009284g0003.1gnl|UG|Nta#S33578151Nta004281MYB65.4108.7528.8
NbS00007772g0004.1gnl|UG|Nta#S33578094Nta008536MYB1.770.8910.2
NCED1Drought and salinity resistanceNbS00014845g0012.1---51.2687.1732.7
NtERD10BImproved drought- and low temperature stress toleranceNbS00041294g0010.1---0.010.017.5
NbS00029572g0007.1---0.010.0113.2
NbS00019308g0009.1---2.365.6186.41
NtERD10CImproved drought- and low-temperature stress toleranceNbS00018964g0003.1---75.1212.4121.9
SRK2E/OST1/SnRK2.6Response to dehydration stressNbS00019609g0014.1---3.837.9818
NbS00018358g0015.1---2.362.969.9
PP2Ac-1Drought resistance; maintain RWC and membrane stabilityNbS00043074g0007.1---12.678.5721.9
NbS00020903g0004.1---81.0252.01110.41
RD22Response to drought stressNbS00017616g0005.1---3.549.1622.8
RD29A and RD29Blow-temperature-responsive and desiccation-responsiveNbS00001669g0005.1---9.433.841.8
SIPKosmotic stress/pathogen resistanceNbS00060107g0004.1---36.5342.8578.91
WIPKosmotic stress/pathogen resistanceNbS00041241g0001.1gnl|UG|Nta#S37452582Nta004233ERF7.660.591.5
NbS00028162g0009.1gnl|UG|Nta#S37452581Nta012156ERF40.076.89.9
NbS00017618g0003.1gnl|UG|Nta#S37452582Nta004233ERF58.335.912.1
ZPT2-1 (renamed from EPF1)Drought toleranceNbS00017486g0002.1---15.324.434.5
NbS00002494g0018.1---44.7833.119.8
1 For a version of the table with references for each gene, see Table S4. Abbreviations: Nicotiana benthamiana, Niben; TPM, transcripts per million clean tags.

2.7. Small RNA Sequencing

To examine miRNAs drought stress-response, two small RNA libraries were constructed based on the result of the physiological index measurement, and another two libraries from leaves and stems of the control plant. Four libraries were then sequenced using the Illumina high-through put sequencing technology. The resulting raw sequence reads (more than fourteen million for each library, 15–35 nt) were processed computationally to remove the 3'adapter and this yielded a total of 29,088,506 genome-matching reads (>18 nt) from the four libraries (6,939,963, 7,514,006, 7,331,406 and 7,303,131 reads from leaf, stem, and root controls, and treated root libraries, respectively) (Supplementary Table S5). The outstanding sequences that could not be mapped onto the tobacco genome were probably derived from unsequenced Nicotiana benthamiana genomic regions, sequence differences between Nicotiana species, sequence errors or contamination, and thus were excluded from subsequent analyses. More than 80% of these genome-mapped small RNAs ranged 20–24 nt reads in length with 24 and 21 nt as the modes which were consistent with the size of Dicer-like (DCL) cleavage products and similar to previous studies in other plant species (Supplementary Figure S5).
The genome-matched small RNA sequences were clustered into several RNA classes such as known miRNAs, repeats, rRNA, tRNA, snRNA, snoRNA and others (Supplementary Table S5). Known tobacco miRNAs account for 21.0% of all sequence reads for the leaf library (highest percentage) and 14.1% for the stem library (lowest percentage), suggesting that mature miRNAs were highly enriched in our small RNA libraries. However, after analyzing the number of unique sequences, the proportion of small RNA sequences derived from known miRNAs represented only a very small fraction (0.12%–0.16%) of the total number (Supplementary Table S5). The highest fraction of unique sequences (>93%) was unclassified small RNA sequences, which probably include novel miRNA candidates and other classes of regulatory RNAs.

2.8. Known Tobacco miRNAs Expressed in Different Tissues and Drought Responsive miRNAs in Tobacco Roots

Currently, miRBase (release 21 [39] contains 165 mature miRNAs sequences in tobacco. Of these, 122 were detected in our sequencing datasets (Supplementary Table S6) and all of the 18 miRNA families conserved in Arabidopsis and tobacco were included. Conserved miRNAs were far more abundant than non-conserved miRNAs in our libraries as reported previously [32,39,40]. MiR166 and miR168 were the most abundant miRNA families which accounted for about 57% and 16% respectively of the total sequence reads from the known miRNAs datasets. However, 44 experimentally identified tobacco miRNA families (33 miRNAs) were not detected in our dataset [25]. A sample of 18 expressed miRNAs was randomly selected for validation by stem-loop qRT-PCR. The expression trend of these miRNAs in control library relative to drought-treated library detected by Solexa small RNA sequencing were basically consistent with that detected by stem-loop qRT-PCR (Figure 5).
Ijms 16 05714 g005 1024
Figure 5. Expression of 18 miRNAs randomly selected in control and drought-treated roots of tobacco. The black bars represent the fold change (log2) in control library relative to drought-treated library detected by Solexa small RNA sequencing, while the gray bars represent the fold change (log2) in control roots relative to drought-treated roots detected by stem-loop qRT-PCR (normalized to 5S rRNA; n = 3).
Figure 5. Expression of 18 miRNAs randomly selected in control and drought-treated roots of tobacco. The black bars represent the fold change (log2) in control library relative to drought-treated library detected by Solexa small RNA sequencing, while the gray bars represent the fold change (log2) in control roots relative to drought-treated roots detected by stem-loop qRT-PCR (normalized to 5S rRNA; n = 3).
Ijms 16 05714 g005
As reported [41], deep sequencing with good reproducibility is a powerful tool for analyzing genome-wide patterns of miRNA expression. The change in frequency of miRNAs between drought-treated and control libraries suggest that their expression was regulated in response to drought stress. As shown in Figure 6, comparison of the normalized sequence reads of the miRNAs between the two libraries indicated that 5 known tobacco miRNA families had relative changes (log2root-ck/root-treat) greater than other 18 miRNA families and thus might be differentially or extremely differentially expressed. However, miR159, miR169, miR402, and miR408 sequence reads displayed no meaningful changes between the two libraries (Figure 6) even though their expression had been reported to be affected by drought stress treatments in other plants [21,23,42,43].
Ijms 16 05714 g006 1024
Figure 6. The expression profiles of the normalized sequence reads of the known miRNA families in the drought-treated library relative to the control library. * p-value < 0.01 is considered as significant; ** p-value < 0.001 is considered as extreme significant.
Figure 6. The expression profiles of the normalized sequence reads of the known miRNA families in the drought-treated library relative to the control library. * p-value < 0.01 is considered as significant; ** p-value < 0.001 is considered as extreme significant.
Ijms 16 05714 g006

2.9. Potential Targets of Drought Responsive miRNAs

miRNAs regulate gene expression at the posttranscriptional level by base paring with complementary sequences of mRNA thus inducing gene silencing or target degradation. Some miRNA families have been reported to respond to salt or drought stress in plants [42]. To understand the function of miRNAs’ response to drought in tobacco, we compiled a complete list of miRNA targets in tobacco identified by degradome sequencing [25]. Only 27 target GSS sequences, corresponding to 87 transcripts, could be mapped to the tobacco reference genome (Supplementary Table S7). Only targets of two drought-responsive miRNAs were obtained and listed in Table 2. Of these targets, miR160 and miR395 are highly conserved miRNA families and have been validated to target those mRNAs encoding ARF transcription factors and AST68 (Sulfate transporter 2.1) proteins, in Arabidopsis [44] and/or tomato [45], respectively. Previous studies showed that freeing ARF17 or ARF16 from miR160 regulation by transgenic Arabidopsis expressing miRNA-resistant targets resulted in dramatic morphological changes and alters basal levels of auxin-induced transcripts [46,47]. The expression of the validated targets of drought-responsive miRNAs was investigated by DGE data (Table 2). We found that there was a negative correlation for miR395 in the expression patterns of the miRNAs and their targets. This is consistent with miRNA function in guiding the cleavage of target mRNAs.
Table
Table 2. Drought responsive miRNA families and their targets in tobacco.
Table 2. Drought responsive miRNA families and their targets in tobacco.
miRNA_familyRegulatedTarget_GSS_IDMapped_Niben_IDNCK_TPMN6H_TPMN48H_TPMDescription
miR160downET792241NbS00018737g0003.10.591.181.20ARF Transcription factor
miR160downET792303NbS00005871g0012.10.010.010.60ARF Transcription factor
miR160downFH226342NbS00059497g0003.14.713.844.20ARF Transcription factor
miR160downFH226342NbS00026905g0003.16.783.844.50ARF Transcription factor
miR160downFH226421NbS00028965g0003.13.542.661.20ARF Transcription factor
miR160downFH226421NbS00036438g0004.10.010.590.01ARF Transcription factor
miR395upET937275NbS00014405g0001.12.361.480.90unknown protein
miR395upET937275NbS00009197g0109.112.373.252.40unknown
miR395upET937275NbS00006001g0005.111.7817.432.10unknown protein

3. Discussion

We first performed DGE analysis to determine genome-wide patterns of gene expression of tobacco roots under drought stress. Given the nature of the DGE system, we pooled biological replicates for each group to make representative samples for deep sequencing analysis based on physiological analysis of PRO, SOD, and MDA. We detected expression of 21,128 genes, which accounted for 73.6% of predicted genes in the tobacco genome. Approximately 30% of the clean tags could not be mapped to any genes in all three samples. These non-mapped tags most likely represent regions where the tobacco reference sequence is incomplete [3], or there are allelic sequence differences between the reference genome and the cultivar used in this study. Alternatively, differential mRNA processing events exist for most tobacco genes, such as alternative splicing [34]. Another reason may be that RNA-seq data for tobacco genome annotation should represent all major tissue types, developmental stages and responses to abiotic and biotic stresses [48]. If the annotated transcripts contained SNP, RNA editing and InDel which are located in CATG sites by chance, and, therefore, the reads could not be mapped to the reference genome sequence.
Through the use of DGE, this global analysis of gene expression provided a comprehensive dataset responding to drought stress in roots of tobacco seedlings. We identified six clusters for all differentially expressed genes and coarsely assigned them to 21 functional categories (p < 0.05) (Figures 3A,B and Supplementary Table S2). Interestingly, there were no other overlapping GO functional enrichments between clusters except transferase, transferring acyl groups, indicating that these genes of different clusters were predicted to be involved in many plant biological processes, including defense [6].
Pathway enrichment analysis revealed 17 pathways were significantly affected by drought stress (Figure 3A,C and Supplementary Table S2). These pathways were enriched in K1/K2/K3 clusters. Glutathione metabolism and biosynthesis of secondary metabolites pathways were the most affected in K1 and K2/K3 clusters, respectively. The former finding implied that Glutathione (GSH) was an abundant and ubiquitous antioxidant with proposed roles in the maintenance of tissue antioxidant defenses and in the regulation of redox-sensitive signal transduction. The size of glutathione pool and its redox status were tightly correlated with the tolerance of plants to drought stress [49,50]. Genes that function in fatty acid (FA) biosynthesis, elongation and lipid beta-oxidation were also enriched in FA elongation of K1, suggesting an elaboration of plastid machinery under drought stress. A number of genes annotated as biosynthesis-related (e.g., NAD(P)-binding Rossmann-fold superfamily protein (NbS00002389g0014.1), Glucose/ribitol dehydrogenase (NbS00007017g0016.1) and ADH-like UDP-glucose dehydrogenase (NbS00019080g0006.1)) were enriched in biosynthesis of secondary metabolite pathway in cluster K2/K3, which was previously described as plant defense-related pathway [51]. We also noted two differentially expressed genes (Aquaporin (NbS00007456g0005.1) and Glutathione transferase, (NbS00049661g0007.1)), whose homologs in other plants play an important role in drought tolerance [52,53,54] and oxidative stress [55], respectively.
Previous studies have uncovered gene expression, transcriptional regulation, and signal transduction in plant responses to drought [6]. In the signal transduction network that leads from the perception of stress signals to the expression of stress-responsive genes, transcription factors (TFs) play an essential role [7]. In our study, TFs detected in DGE analysis included 82 members of 24 TF families with differential expression, of which 21 were DRGs, indicating that these TFs probably played important roles in response to drought.
In addition, some miRNAs have been shown to be stress regulated at the posttranscriptional level by repressing mRNA expression and could be involved in cell responses to various abiotic stresses such as salinity and drought [42]. Some miRNA target genes are stress-responsive TFs or functional genes [42], indicating that miRNA-dependent posttranscriptional regulation may play a role in plant stress response. In this study, we used high-throughput sequencing to confirm expression of the known miRNAs and their changes upon drought treatment. We detected the expression of 122 known miRNAs (39 families) in tobacco roots, stems and leaves. Upon drought treatment, five conserved miRNA families showed differential regulation (Figure 6). One of these drought-responsive conserved miRNA families, miR395, has been shown to be stress responsive in other plant species [42], suggesting that these miRNAs might also play a positive regulatory role in drought-tolerance in tobacco roots. Additionally, most targets of the new miRNAs identified by degradome sequencing [25] were mapped to the reference tobacco genome and 87 potential target genes were found. Ten targets responsive to drought stress suggest the importance of these miRNAs in stress responses in tobacco (Supplementary Table S7). However, expression of some miRNA families (miR159, miR169, miR402 and miR408) had no meaningful changes between libraries thus indicating that different stress conditions might evoke diverse plant responses and physiological adaption.
With the availability of regulatory networks of gene expression in drought and cold stress responses [6], an integrated gene regulatory network has been proposed for the molecular mechanisms of the response of tobacco roots to drought stress using differentially expressed DRGs, the changed expression profiles of miRNAs and subsequent target transcripts as a basis [56] (Figure 7). Two pathways (ABA-dependent and ABA-independent) can shed light on cell mechanisms involved in stress signaling and/or adaptation at transcriptional regulation. In the ABA-dependent pathway, NCED1 was involved in rapid emergency response to drought and Figure 7 shows transcription cascades involved in slow and adaptive processes in stress responses, such as those involving AREB/ABF, MYB, bZIP, NAC and SnRK2.6 protein kinases which are involved in ABA signaling [57,58]. In the ABA-independent pathway, unknown proteins were thought to function as an osmo-sensor and function upstream of the ERF system [59]. In addition, the responsive miR395 also showed a transitory expression model. Right side of Figure 7 shows the proposed regulation cascades in tobacco roots after initiation of drought stress. Our study provides valuable information for future studies of the molecular mechanisms underlying drought tolerance in tobacco roots and other plants.
Ijms 16 05714 g007 1024
Figure 7. A probable model for regulatory network being involved in drought responses in tobacco roots. The integrated gene regulatory network has been proposed on basis of previous regulatory network model of Arabidopsis in the drought and cold stress responses. Abbreviations: ABA, abscisic acid; NCED, 9-cis-epoxycarotenoid dioxygenase gene; SnRK2.6, SNF1-related protein kinase2.6/OST1 kinase; NAC, NAM, ATAF, and CUC transcription factor; bZIP, basic-domain leucine-zipper; ABF, ABRE-binding factor; AREB, ABA-responsive-element-binding protein; ERF, ethylene-responsive-element-binding factor; DRE, dehydration-responsive element; DREB, DRE-binding factor; DREB1B, DREB2A, DREB3 and GmERF3 are the members of ERF family.
Figure 7. A probable model for regulatory network being involved in drought responses in tobacco roots. The integrated gene regulatory network has been proposed on basis of previous regulatory network model of Arabidopsis in the drought and cold stress responses. Abbreviations: ABA, abscisic acid; NCED, 9-cis-epoxycarotenoid dioxygenase gene; SnRK2.6, SNF1-related protein kinase2.6/OST1 kinase; NAC, NAM, ATAF, and CUC transcription factor; bZIP, basic-domain leucine-zipper; ABF, ABRE-binding factor; AREB, ABA-responsive-element-binding protein; ERF, ethylene-responsive-element-binding factor; DRE, dehydration-responsive element; DREB, DRE-binding factor; DREB1B, DREB2A, DREB3 and GmERF3 are the members of ERF family.
Ijms 16 05714 g007

4. Experimental Section

4.1. Plant Materials, Growth Conditions and Stress Treatments

A flue-cured tobacco (Nicotiana tabacum L.) cultivar named Honghua Dajinyuan, known for drought-tolerance and widely grown in Southwest China, was used in our study. The seed was sterilized and germinated in an incubator as reported earlier [60]. The germinated seed was then sown in spots with soil matrix and grown until the seedlings developed six leaves at 28/21 °C day/night temperatures (16 h day/ 8 h night) and relative humidity of 70%. To mimic drought, seedlings were treated by adding 20% PEG6000 to the matrix. Samples were taken 0, 3, 6, 12, 24, 48, and 96 h after treatment [60]. The control was not treated and is the time zero treatment. Roots were collected from all treatments. In addition, leaf and stem samples were collected from the control. Tissues were snap-frozen in liquid nitrogen, and then stored at −80 °C for physiological index measurements [61,62,63], RNA extraction, and quantitative real-time PCR (qRT-PCR) validation.

4.2. Physiological Analyses

Ultraviolet-visible transmittance spectrophotometers (UV-VIS) (Shimadzu, UV-2450, Kyoto, Japan) were used to measure three physiological indexes: superoxide dismutase (SOD) activity, and proline (PRO) and malondialdehyde (MDA) content of the root samples [61,62,63]. A physiological assay kit was used and tests were performed according to the manufacturer’s protocol (Nanjing Jiancheng Bioengineering Institute, Nanjing, China). All tests were replicated three times. Single factor analysis of variance was performed using SPSS16.0 (SPSS Inc., Chicago, IL, USA). Fisher’s Least Significance Difference (LSD) with multiple comparisons was used to separate treatments at (p < 0.01, n = 3) level of significance.

4.3. DGE, Small RNA Library Construction, and Solexa Sequencing

Total RNA was isolated from the frozen roots samples by using Trizol Reagent (Invitrogen, Carlsbad, CA, USA) according to manufacturer’s instructions. DGE library preparation was performed in parallel by using the Illumina gene expression sample preparation kit as described by [31,64]. Each tunnel will generate millions of raw reads with sequencing length of 35 bp. Raw data (tag sequences) were deposited in the GEO database (No. GES43058).
The samples of our small RNA libraries were used based on the result of physiological index measurement as follows: equal quantities (10 µg) of total RNA isolated from tobacco roots treated with two time points (6 and 48 h) were mixed together to construct the drought-treated small RNA library (Root-treat), and total RNA prepared from the control roots sample was used to construct the control small RNA library (Root-ck). In addition, we constructed two libraries from leaves and stems of the control plants. These libraries were constructed using the Small RNA Sample Prep Kit (Illumina, San Diego, CA, USA) as described by Tang et al. [25] and Yin et al. [32].

4.4. Analysis and Mapping of Digital Gene Expression Tags

Raw sequencing image data were transformed by base calling into sequences. Raw data reads were stored in FASTA format, and analysis was conducted as described by Shen et al. [31,64]. Prior to mapping to the reference database, all sequences were filtered to trim the 3' adaptor sequence, filter empty tags (reads with only 3' adaptor sequences but no tags) and low-quality tags containing Ns, and remove tags which are too long or too short. Because of the lack of EST sequence annotation of Nicotiana tabacum L., we chose the genome of its close relative Nicotiana benthamiana (v.0.4.4) [3] as reference, which is a virtual library containing all possible CATG + 17 base-length sequences and present high degree of sequence similarity (>90%) with Nicotiana tabacum L. (Supplementary Figure S6). All clean tags were mapped to the reference sequences and a mismatch of only 1 bp was considered. Clean tags that were mapped to the tobacco genome reference sequences from multiple genes were filtered. The remaining clean tags were designated unambiguous clean tags. The expression level of each gene was estimated by the frequency of clean tags and then normalized to TPM (transcripts per million clean tags) [65], which is a standard method and extensively used in DGE analysis [66]. The expression level of each gene was measured by the normalized number of matched clean tags. KOG functional classification, Gene Ontology (GO), pathway annotation and enrichment analyses were based on the Clusters of Orthologous Groups of proteins (COGs) (http://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/COG) [67], Gene Ontology Database (http://www.geneontology.org/) [68] and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway (http://www.genome.jp/kegg/) [69], respectively. When we investigated pathways in which different genes were involved and enriched, q-value was used to aid identification according to the previous description [70].

4.5. Identification of Differentially Expressed Genes and Cluster Analysis

For differential expression analysis between control and treatment samples, the probability that one gene G is equally expressed in two samples was expressed by p values [71]. Thresholds of p values for multiple tests were determined by False Discovery Rates (FDR) [72]. In addition, we applied the R package [73] to identify differentially expressed genes with the random sampling model based on the read count for each gene under drought stress. We used “FDR ≤ 0.001 and the absolute value of log2-fold change ≥ 1” [74] as the threshold for judging the significance of different gene expressions. More stringent criteria with a smaller FDR and higher fold-change value were used to identify differentially expressed genes. GO functional enrichment analysis of differentially expressed genes was carried out using Blast2GO (version 2.3.5, https://www.blast2go.com/). KEGG pathway analyses of differentially expressed genes were performed using Cytoscape software (version 3.1.1, http://www.softpedia.com/get/Science-CAD/Cytoscape.shtml) with the ClueGO plug (http://www.ici.upmc.fr/cluego/cluegoDownload.shtml) [75]. GO annotations were performed using AgriGO (version 1.2; http://bioinfo.cau.edu.cn/agriGO/). K-means clustering was previously described by Li et al. [34].

4.6. Analysis of Small RNA Sequencing Data and Identification of Drought Responsive miRNAs

The overall procedure for analyzing Solexa small RNA libraries was performed as previously described [25,26]. The Illumina Gerald pipeline was used to process and extract the first 36 bases of each read. Adaptor sequences were identified and trimmed from each read using a customized Perl script. Reads in which the adaptor could not be identified were discarded. We used the Nicotiana benthamiana genome (v.0.4.4) database (ftp://ftp.solgenomics.net/genomes/Nicotiana_benthamiana) [3], and Sanger Rfam data (ftp://ftp.sanger.ac.uk/pub/databases/Rfam/) [76] to identify sequences originating from protein-coding genes, repeats, rRNA, tRNA, snRNA, and snoRNA. SOAP2 software [33] was used to align the trimmed reads to the high-confidence set of 165 pre-miRNAs from miRBase (version 21; http://www.mirbase.org/) [25]. For each library, we counted the number of trimmed reads within the 18–25nt ranges that were mapping to each pre-miRNA and normalized by the total number of 19–24 nt trimmed reads in the library. Trimmed reads that were < 18 or ≥ 25 nt were not considered in this analysis. The relative change of individual miRNA families of detected reads for treatments and p-values of the t-test were calculated; differentially expressed miRNA families were those with p-values <0.01.

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

In order to verify the DGE sequencing results, 22 differentially expressed genes were randomly selected and qRT-PCR performed (Supplementary Table S8A). Specific primers for an internal gene, actin (GenBank no. X63603, Forward: 5'-CGCGAAAAGATGACTCAAATC-3', Reverse: 5'-AGATCCTTTCTGATATCCACG-3'), were used for the normalization of reactions [77]. qRT-PCRs were performed using SYBR Premix Ex Taq protocol (TaKaRa, Shiga, Japan) on an Applied Biosystems 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). RNA samples were the same as those for the DGE experiments from biological replicates, and triplicate experiments were performed using independent plant materials, from which results were calculated using the 2−ΔΔCt method [78].
Expression trends of 18 drought responsive mature miRNAs were assayed by stem-loop reverse transcription-PCR (RT–PCR). 200 ng of total RNA was used for the initiation of the reverse transcription reaction. The stem–loop reverse transcription primers were designed following the method described by Chen et al. [79] and Varkonyi-Gasic et al. [80]. The reverse transcription product was amplified using a miRNA-specific forward primer and a universal reverse primer. The stem–loop reverse transcription reactions were performed using One Step PrimeScript® miRNA cDNA Synthesis Kit (TaKaRa, Japan), and PCR primers were then added to perform the PCR on the Applied Biosystems 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). One of the uniformly expressed 5S rRNAs was used as the internal control for stem-loop qRT-PCR [81]. Supplementary Table S8B shows the sequences of stem-loop RT primers and miRNA-specific PCR primers.

5. Conclusions

In summary, this study provided a comprehensive analysis of drought-responsive genes and microRNAs expression profiles of tobacco roots by combined DGE, small RNA sequencing, and computational approaches, and revealed 54 differentially expressed DRGs included 21 TFs, which belonged to 4 TF families such as NAC (6), MYB (4), ERF (10), and bZIP (1). Additionally, we confirmed expression of 39 known miRNA families (122 members) and five conserved miRNA families, which showed differential regulation under drought stress. Targets of miRNAs were further surveyed based on a recently published study, of which ten targets were DRGs. Finally, we developed and mapped the transcriptional and posttranscriptional gene regulatory network that operated in response to drought stress in tobacco root tissue. This network analysis provided new clues for future studies on tobacco response to various stresses, such as to drought, cold and heavy metals.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/16/03/5714/s1.

Acknowledgments

We appreciated the helpful comments on the manuscript from Dan Meckenstock. This work was supported by grants from the National Natural Science Foundation of China (No. 30900901), the Sichuan Province Special Sci-Tech Projects (No. 2006J13-039), the Ministry of Agriculture in China (No. 2011ZX08003-003) and the Open Project Program of Key Laboratory of Tobacco Genetic and Breeding, State Tobacco Monopoly Administration (Molecular Mechanism of Drought Resistant in Tobacco).

Author Contributions

Guangtang Pan, Cheng Qin and Fuqiang Yin designed research; Fuqiang Yin, Cheng Qin and Jian Gao performed research, and they contributed equally to this study. Cheng Qin, Xirong Luo, Zhiming Zhang and Likai Mao analyzed data; Ming Liu and Xirong Luo conducted physiological analyses and qRT-PCR experiments; Fuqiang Yin, Cheng Qin and Thomas Lübberstedt wrote the paper. Wenyou Zhang, Xinhui Liao, Yaou Shen, Zhiming Zhang and Hongjun Liu contributed reagents/materials/analysis tools. All authors have read and approved the manuscript for publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hu, T.-W.; Mao, Z. Economic Analysis of Tobacco and Options for Tobacco Control: China Case Study; Center for Tobacco Control Research and Education: UC San Francisco, CA, USA, 2002. [Google Scholar]
  2. Goodin, M.M.; Zaitlin, D.; Naidu, R.A.; Lommel, S.A. Nicotiana benthamiana: Its history and future as a model for plant-pathogen interactions. Mol. Plant-Microbe Interact. 2008, 21, 1015–1026. [Google Scholar] [CrossRef] [PubMed]
  3. Bombarely, A.; Rosli, H.G.; Vrebalov, J.; Moffett, P.; Mueller, L.A.; Martin, G.B. A draft genome sequence of nicotiana benthamiana to enhance molecular plant-microbe biology research. Mol. Plant-Microbe Interact. 2012, 25, 1523–1530. [Google Scholar] [CrossRef] [PubMed]
  4. Glazebrook, J. Genes controlling expression of defense responses in Arabidopsis—2001 status. Curr. Opin. Plant Biol. 2001, 4, 301–308. [Google Scholar] [CrossRef] [PubMed]
  5. Knight, H.; Knight, M.R. Abiotic stress signalling pathways: Specificity and cross-talk. Trends Plant Sci. 2001, 6, 262–267. [Google Scholar] [CrossRef] [PubMed]
  6. Shinozaki, K.; Yamaguchi-Shinozaki, K.; Seki, M. Regulatory network of gene expression in the drought and cold stress responses. Curr. Opin. Plant Biol. 2003, 6, 410–417. [Google Scholar] [PubMed]
  7. Hussain, S.S.; Kayani, M.A.; Amjad, M. Transcription factors as tools to engineer enhanced drought stress tolerance in plants. Biotechnol. Progr. 2011, 27, 297–306. [Google Scholar] [CrossRef]
  8. Shinozaki, K.; Yamaguchi-Shinozaki, K. Gene networks involved in drought stress response and tolerance. J. Exp. Bot. 2007, 58, 221–227. [Google Scholar] [CrossRef]
  9. Zhu, J.K. Salt and drought stress signal transduction in plants. Annu. Rev. Plant Biol. 2002, 53, 247–273. [Google Scholar] [CrossRef] [PubMed]
  10. Choi, H.; Hong, J.; Ha, J.; Kang, J.; Kim, S.Y. ABFs, a family of ABA-responsive element binding factors. J. Biol. Chem. 2000, 275, 1723–1730. [Google Scholar] [CrossRef] [PubMed]
  11. Uno, Y.; Furihata, T.; Abe, H.; Yoshida, R.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Arabidopsis basic leucine zipper transcription factors involved in an abscisic acid-dependent signal transduction pathway under drought and high-salinity conditions. Proc. Natl. Acad. Sci. USA 2000, 97, 11632–11637. [Google Scholar] [CrossRef] [PubMed]
  12. Narusaka, Y.; Nakashima, K.; Shinwari, Z.K.; Sakuma, Y.; Furihata, T.; Abe, H.; Narusaka, M.; Shinozaki, K.; Yamaguchi-Shinozaki, K. Interaction between two cis-acting elements, ABRE and DRE, in ABA-dependent expression of Arabidopsis rd29A gene in response to dehydration and high-salinity stresses. Plant J. 2003, 34, 137–148. [Google Scholar] [CrossRef] [PubMed]
  13. Lee, Y.; Kim, M.; Han, J.; Yeom, K.H.; Lee, S.; Baek, S.H.; Kim, V.N. MicroRNA genes are transcribed by RNA polymerase II. EMBO J. 2004, 23, 4051–4060. [Google Scholar] [CrossRef] [PubMed]
  14. Voinnet, O. Origin, biogenesis, and activity of plant microRNAs. Cell 2009, 136, 669–687. [Google Scholar] [CrossRef] [PubMed]
  15. Boualem, A.; Laporte, P.; Jovanovic, M.; Laffont, C.; Plet, J.; Combier, J.P.; Niebel, A.; Crespi, M.; Frugier, F. MicroRNA166 controls root and nodule development in Medicago truncatula. Plant J. 2008, 54, 876–887. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, H.; Nussbaum-Wagler, T.; Li, B.; Zhao, Q.; Vigouroux, Y.; Faller, M.; Bomblies, K.; Lukens, L.; Doebley, J.F. The origin of the naked grains of maize. Nature 2005, 436, 714–719. [Google Scholar] [CrossRef] [PubMed]
  17. Golz, J.F. Signalling between the shoot apical meristem and developing lateral organs. Plant Mol. Biol. 2006, 60, 889–903. [Google Scholar] [CrossRef] [PubMed]
  18. Liu, Q.; Zhang, Y.C.; Wang, C.Y.; Luo, Y.C.; Huang, Q.J.; Chen, S.Y.; Zhou, H.; Qu, L.H.; Chen, Y.Q. Expression analysis of phytohormone-regulated microRNAs in rice, implying their regulation roles in plant hormone signaling. FEBS Lett. 2009, 583, 723–728. [Google Scholar] [CrossRef] [PubMed]
  19. Ding, D.; Zhang, L.; Wang, H.; Liu, Z.; Zhang, Z.; Zheng, Y. Differential expression of miRNAs in response to salt stress in maize roots. Ann. Bot. 2009, 103, 29–38. [Google Scholar] [CrossRef] [PubMed]
  20. Lu, S.; Sun, Y.H.; Chiang, V.L. Stress-responsive microRNAs in Populus. Plant J. 2008, 55, 131–151. [Google Scholar] [CrossRef] [PubMed]
  21. Zhao, B.; Liang, R.; Ge, L.; Li, W.; Xiao, H.; Lin, H.; Ruan, K.; Jin, Y. Identification of drought-induced microRNAs in rice. Biochem. Biophys. Res. Commun. 2007, 354, 585–590. [Google Scholar] [CrossRef] [PubMed]
  22. Zhou, X.; Wang, G.; Sutoh, K.; Zhu, J.K.; Zhang, W. Identification of cold-inducible microRNAs in plants by transcriptome analysis. Biochim. Biophys. Acta 2008, 1779, 780–788. [Google Scholar] [CrossRef] [PubMed]
  23. Reyes, J.L.; Chua, N.H. ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J. 2007, 49, 592–606. [Google Scholar] [CrossRef] [PubMed]
  24. Eveland, A.L.; Satoh-Nagasawa, N.; Goldshmidt, A.; Meyer, S.; Beatty, M.; Sakai, H.; Ware, D.; Jackson, D. Digital gene expression signatures for maize development. Plant Physiol. 2010, 154, 1024–1039. [Google Scholar] [CrossRef] [PubMed]
  25. Tang, S.; Wang, Y.; Li, Z.; Gui, Y.; Xiao, B.; Xie, J.; Zhu, Q.-H.; Fan, L. Identification of wounding and topping responsive small RNAs in tobacco (Nicotiana tabacum). BMC Plant Biol. 2012, 12, 28. [Google Scholar] [CrossRef] [PubMed]
  26. Li, T.; Li, H.; Zhang, Y.-X.; Liu, J.-Y. Identification and analysis of seven H2O2-responsive miRNAs and 32 new miRNAs in the seedlings of rice (Oryza sativa L. ssp. indica). Nucleic Acids Res. 2011, 39, 2821–2833. [Google Scholar] [CrossRef] [PubMed]
  27. Rizhsky, L.; Liang, H.; Mittler, R. The combined effect of drought stress and heat shock on gene expression in tobacco. Plant Physiol. 2002, 130, 1143–1151. [Google Scholar] [CrossRef] [PubMed]
  28. Ma, X.; Wang, Y.; Xie, S.; Wang, C.; Wang, W. Glycinebetaine application ameliorates negative effects of drought stress in tobacco. Russ. J. Plant Physiol. 2007, 54, 472–479. [Google Scholar] [CrossRef]
  29. Karakas, B.; Ozias-Akins, P.; Stushnoff, C.; Suefferheld, M.; Rieger, M. Salinity and drought tolerance of mannitol-accumulating transgenic tobacco. Plant Cell Environ. 1997, 20, 609–616. [Google Scholar] [CrossRef]
  30. Cho, E.; Hong, C. Over-expression of tobacco NtHSP70–1 contributes to drought-stress tolerance in plants. Plant Cell Rep. 2006, 25, 349–358. [Google Scholar] [CrossRef] [PubMed]
  31. Ding, H.; Qin, C.; Luo, X.; Li, L.; Chen, Z.; Liu, H.; Gao, J.; Lin, H.; Shen, Y.; Zhao, M.; et al. Heterosis in early maize ear inflorescence development: A genome-wide transcription analysis for two maize inbred lines and their hybrid. Int. J. Mol. Sci. 2014, 15, 13892–13915. [Google Scholar] [CrossRef] [PubMed]
  32. Yin, F.; Gao, J.; Liu, M.; Qin, C.; Zhang, W.; Yang, A.; Xia, M.; Zhang, Z.; Shen, Y.; Lin, H.; et al. Genome-wide analysis of water-stress-responsive microRNA expression profile in tobacco roots. Funct. Integr. Genomics 2014, 14, 319–332. [Google Scholar] [CrossRef] [PubMed]
  33. Li, R.; Yu, C.; Li, Y.; Lam, T.W.; Yiu, S.M.; Kristiansen, K.; Wang, J. SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics 2009, 25, 1966–1967. [Google Scholar] [CrossRef] [PubMed]
  34. Li, P.; Ponnala, L.; Gandotra, N.; Wang, L.; Si, Y.; Tausta, S.L.; Kebrom, T.H.; Provart, N.; Patel, R.; Myers, C.R.; et al. The developmental dynamics of the maize leaf transcriptome. Nat. Genet. 2010, 42, 1060–1067. [Google Scholar] [CrossRef] [PubMed]
  35. Vilella, A.J.; Severin, J.; Ureta-Vidal, A.; Heng, L.; Durbin, R.; Birney, E. EnsemblCompara GeneTrees: Complete, duplication-aware phylogenetic trees in vertebrates. Genome Res. 2009, 19, 327–335. [Google Scholar] [CrossRef] [PubMed]
  36. Riechmann, J.L.; Heard, J.; Martin, G.; Reuber, L.; Jiang, C.Z.; Keddie, J.; Adam, L.; Pineda, O.; Ratcliffe, O.J.; Samaha, R.R.; et al. Arabidopsis transcription factors: Genome-wide comparative analysis among eukaryotes. Science 2000, 290, 2105–2110. [Google Scholar] [CrossRef] [PubMed]
  37. Cartolano, M.; Castillo, R.; Efremova, N.; Kuckenberg, M.; Zethof, J.; Schwarz-Sommer, Z.; Gerats, T.; Vandenbussche, M. A conserved microRNA module exerts homeotic control over Petunia hybrida and Antirrhinum majus floral organ identity. Nat. Genet. 2007, 39, 901–905. [Google Scholar] [CrossRef] [PubMed]
  38. Singh, K.B.; Foley, R.C.; Oñate-Sánchez, L. Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol. 2002, 5, 430–436. [Google Scholar] [CrossRef] [PubMed]
  39. Jones-Rhoades, M.W.; Bartel, D.P.; Bartel, B. MicroRNAs and their regulatory roles in plants. Annu. Rev. Plant Biol. 2006, 57, 19–53. [Google Scholar] [CrossRef] [PubMed]
  40. Gao, J.; Yin, F.; Liu, M.; Luo, M.; Qin, C.; Yang, A.; Yang, S.; Zhang, Z.; Shen, Y.; Lin, H.; Pan, G. Identification and characterisation of tobacco microRNA transcriptome using high-throughput sequencing. Plant Biol. 2015. [Google Scholar] [CrossRef]
  41. Carthew, R.W.; Sontheimer, E.J. Origins and mechanisms of miRNAs and siRNAs. Cell 2009, 136, 642–655. [Google Scholar] [CrossRef] [PubMed]
  42. Covarrubias, A.A.; Reyes, J.L. Post-transcriptional gene regulation of salinity and drought responses by plant microRNAs. Plant Cell Environ. 2010, 33, 481–489. [Google Scholar] [CrossRef] [PubMed]
  43. Lu, S.; Sun, Y.H.; Shi, R.; Clark, C.; Li, L.; Chiang, V.L. Novel and mechanical stress-responsive MicroRNAs in Populus trichocarpa that are absent from Arabidopsis. Plant Cell 2005, 17, 2186–2203. [Google Scholar] [CrossRef] [PubMed]
  44. Addo-Quaye, C.; Eshoo, T.W.; Bartel, D.P.; Axtell, M.J. Endogenous siRNA and miRNA Targets Identified by Sequencing of the Arabidopsis Degradome. Curr. Biol. CB 2008, 18, 758–762. [Google Scholar] [CrossRef]
  45. Itaya, A.; Bundschuh, R.; Archual, A.J.; Joung, J.-G.; Fei, Z.; Dai, X.; Zhao, P.X.; Tang, Y.; Nelson, R.S.; Ding, B. Small RNAs in tomato fruit and leaf development. Biochim. Biophys. Acta (BBA) 2008, 1779, 99–107. [Google Scholar] [CrossRef]
  46. Wang, J.-W.; Wang, L.-J.; Mao, Y.-B.; Cai, W.-J.; Xue, H.-W.; Chen, X.-Y. Control of Root Cap Formation by MicroRNA-Targeted Auxin Response Factors in Arabidopsis. Plant Cell Online 2005, 17, 2204–2216. [Google Scholar] [CrossRef]
  47. Mallory, A.C.; Bartel, D.P.; Bartel, B. MicroRNA-directed regulation of arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell Online 2005, 17, 1360–1375. [Google Scholar] [CrossRef]
  48. Potato Genome Sequencing Consortium; Xu, X.; Pan, S.; Cheng, S.; Zhang, B.; Mu, D.; Ni, P.; Zhang, G.; Yang, S.; Li, R.; et al. Genome sequence and analysis of the tuber crop potato. Nature 2011, 475, 189–195. [Google Scholar]
  49. Dhindsa, R.S. Drought stress, enzymes of glutathione metabolism, oxidation injury, and protein synthesis in tortula ruralis. Plant Physiol. 1991, 95, 648–651. [Google Scholar] [CrossRef] [PubMed]
  50. May, M.J.; Vernoux, T.; Leaver, C.; Montagu, M.V.; Inzé, D. Glutathione homeostasis in plants: Implications for environmental sensing and plant development. J. Exp. Bot. 1998, 49, 649–667. [Google Scholar]
  51. Ma, S.; Gong, Q.; Bohnert, H.J. Dissecting salt stress pathways. J. Exp. Bot. 2006, 57, 1097–1107. [Google Scholar] [CrossRef] [PubMed]
  52. Knipfer, T.; Besse, M.; Verdeil, J.L.; Fricke, W. Aquaporin-facilitated water uptake in barley (Hordeum vulgare L.) roots. J. Exp. Bot. 2011, 62, 4115–4126. [Google Scholar]
  53. Heinen, R.B.; Ye, Q.; Chaumont, F. Role of aquaporins in leaf physiology. J. Exp. Bot. 2009, 60, 2971–2985. [Google Scholar] [CrossRef] [PubMed]
  54. Chaumont, F.; Moshelion, M.; Daniels, M.J. Regulation of plant aquaporin activity. Biol. Cell 2005, 97, 749–764. [Google Scholar] [CrossRef] [PubMed]
  55. Kilili, K.G.; Atanassova, N.; Vardanyan, A.; Clatot, N.; Al-Sabarna, K.; Kanellopoulos, P.N.; Makris, A.M.; Kampranis, S.C. Differential roles of tau class glutathione S-transferases in oxidative stress. J. Biol. Chem. 2004, 279, 24540–24551. [Google Scholar] [CrossRef] [PubMed]
  56. Lo, K.; Raftery, A.; Dombek, K.; Zhu, J.; Schadt, E.; Bumgarner, R.; Yeung, K. Integrating external biological knowledge in the construction of regulatory networks from time-series expression data. BMC Syst. Biol. 2012, 6, 101. [Google Scholar] [CrossRef] [PubMed]
  57. Yoshida, R.; Hobo, T.; Ichimura, K.; Mizoguchi, T.; Takahashi, F.; Aronso, J.; Ecker, J.R.; Shinozaki, K. ABA-activated SnRK2 protein kinase is required for dehydration stress signaling in arabidopsis. Plant Cell Physiol. 2002, 43, 1473–1483. [Google Scholar] [CrossRef] [PubMed]
  58. Mustilli, A.-C.; Merlot, S.; Vavasseur, A.; Fenzi, F.; Giraudat, J. Arabidopsis OST1 protein kinase mediates the regulation of stomatal aperture by abscisic acid and acts upstream of reactive oxygen species production. Plant Cell Online 2002, 14, 3089–3099. [Google Scholar] [CrossRef]
  59. Shinozaki, K.; Yamaguchi-Shinozaki, K. Molecular responses to dehydration and low temperature: differences and cross-talk between two stress signaling pathways. Curr. Opin. Plant Biol. 2000, 3, 217–223. [Google Scholar] [CrossRef] [PubMed]
  60. Frazier, T.; Sun, G.; Burklew, C.; Zhang, B. Salt and drought stresses induce the aberrant expression of microRNA genes in tobacco. Mol. Biotechnol. 2011, 49, 159–165. [Google Scholar] [CrossRef] [PubMed]
  61. Cakmak, I.; Horst, W.J. Effect of aluminium on lipid peroxidation, superoxide dismutase, catalase, and peroxidase activities in root tips of soybean (Glycine max). Physiol. Plantarum 1991, 83, 463–468. [Google Scholar] [CrossRef]
  62. Zhang, J.; Kirkham, M.B. Drought-stress-induced changes in activities of superoxide dismutase, catalase, and peroxidase in wheat species. Plant Cell Physiol. 1994, 35, 785–791. [Google Scholar]
  63. Delauney, A.J.; Verma, D.P.S. Proline biosynthesis and osmoregulation in plants. Plant J. 1993, 4, 215–223. [Google Scholar] [CrossRef]
  64. Shen, Y.; Jiang, Z.; Yao, X.; Zhang, Z.; Lin, H.; Zhao, M.; Liu, H.; Peng, H.; Li, S.; Pan, G. Genome Expression profile analysis of the immature maize embryo during dedifferentiation. PLoS One 2012, 7, e32237. [Google Scholar] [CrossRef] [PubMed]
  65. AC’t Hoen, P.; Ariyurek, Y.; Thygesen, H.H.; Vreugdenhil, E.; Vossen, R.H.A.M.; de Menezes, R.X.; Boer, J.M.; van Ommen, G.J.B.; den Dunnen, J.T. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36, e141. [Google Scholar] [CrossRef] [PubMed]
  66. Morrissy, A.S.; Morin, R.D.; Delaney, A.; Zeng, T.; McDonald, H.; Jones, S.; Zhao, Y.; Hirst, M.; Marra, M.A. Next-generation tag sequencing for cancer gene expression profiling. Genome Res. 2009, 19, 1825–1835. [Google Scholar] [CrossRef] [PubMed]
  67. Tatusov, R.L.; Fedorova, N.D.; Jackson, J.D.; Jacobs, A.R.; Kiryutin, B.; Koonin, E.V.; Krylov, D.M.; Mazumder, R.; Mekhedov, S.L.; Nikolskaya, A.N.; et al. The COG database: An updated version includes eukaryotes. BMC Bioinf. 2003, 4, 41. [Google Scholar] [CrossRef]
  68. Consortium, T.G.O. The gene ontology project in 2008. Nucleic Acids Res. 2008, 36, D440–D444. [Google Scholar] [CrossRef] [PubMed]
  69. Kanehisa, M. The KEGG Database. In “In Silico” Simulation of Biological Processes; John Wiley & Sons, Ltd: Hoboken, NJ, USA, 2002; pp. 91–103. [Google Scholar]
  70. Storey, J.D.; Tibshirani, R. Statistical significance for genomewide studies. Proc. Natl. Acad. Sci. USA 2003, 100, 9440–9445. [Google Scholar] [CrossRef] [PubMed]
  71. Audic, S.; Claverie, J.M. The significance of digital gene expression profiles. Genome Res. 1997, 7, 986–995. [Google Scholar] [PubMed]
  72. Reiner, A.; Yekutieli, D.; Benjamini, Y. Identifying differentially expressed genes using false discovery rate controlling procedures. Bioinformatics 2003, 19, 368–375. [Google Scholar] [CrossRef] [PubMed]
  73. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [PubMed]
  74. Benjamini, Y.; Yekutieli, D. The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 2001, 29, 1165–1188. [Google Scholar] [CrossRef]
  75. Bindea, G.; Mlecnik, B.; Hackl, H.; Charoentong, P.; Tosolini, M.; Kirilovsky, A.; Fridman, W.-H.; Pagès, F.; Trajanoski, Z.; Galon, J. ClueGO: A Cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009, 25, 1091–1093. [Google Scholar] [CrossRef]
  76. Gardner, P.P.; Daub, J.; Tate, J.G.; Nawrocki, E.P.; Kolbe, D.L.; Lindgreen, S.; Wilkinson, A.C.; Finn, R.D.; Griffiths-Jones, S.; Eddy, S.R.; et al. Rfam: Updates to the RNA families database. Nucleic Acids Res. 2009, 37, D136–D140. [Google Scholar] [CrossRef] [PubMed]
  77. Pignocchi, C.; Fletcher, J.M.; Wilkinson, J.E.; Barnes, J.D.; Foyer, C.H. The function of ascorbate oxidase in tobacco. Plant Physiol. 2003, 132, 1631–1641. [Google Scholar] [CrossRef] [PubMed]
  78. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time Quantitative PCR and the 2−ΔΔCt method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  79. Chen, C.; Ridzon, D.A.; Broomer, A.J.; Zhou, Z.; Lee, D.H.; Nguyen, J.T.; Barbisin, M.; Xu, N.L.; Mahuvakar, V.R.; Andersen, M.R.; et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33, e179. [Google Scholar] [CrossRef] [PubMed]
  80. Varkonyi-Gasic, E.; Wu, R.; Wood, M.; Walton, E.; Hellens, R. Protocol: A highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 2007, 3, 1–12. [Google Scholar] [CrossRef] [PubMed]
  81. Lang, Q.; Jin, C.; Lai, L.; Feng, J.; Chen, S.; Chen, J. Tobacco microRNAs prediction and their expression infected with Cucumber mosaic virus and Potato virus X. Mol. Biol. Rep. 2011, 38, 1523–1531. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Yin, F.; Qin, C.; Gao, J.; Liu, M.; Luo, X.; Zhang, W.; Liu, H.; Liao, X.; Shen, Y.; Mao, L.; et al. Genome-Wide Identification and Analysis of Drought-Responsive Genes and MicroRNAs in Tobacco. Int. J. Mol. Sci. 2015, 16, 5714-5740. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms16035714

AMA Style

Yin F, Qin C, Gao J, Liu M, Luo X, Zhang W, Liu H, Liao X, Shen Y, Mao L, et al. Genome-Wide Identification and Analysis of Drought-Responsive Genes and MicroRNAs in Tobacco. International Journal of Molecular Sciences. 2015; 16(3):5714-5740. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms16035714

Chicago/Turabian Style

Yin, Fuqiang, Cheng Qin, Jian Gao, Ming Liu, Xirong Luo, Wenyou Zhang, Hongjun Liu, Xinhui Liao, Yaou Shen, Likai Mao, and et al. 2015. "Genome-Wide Identification and Analysis of Drought-Responsive Genes and MicroRNAs in Tobacco" International Journal of Molecular Sciences 16, no. 3: 5714-5740. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms16035714

Article Metrics

Back to TopTop