Comprehensive analysis of differentially expressed genes and transcriptional regulation induced by salt stress in two contrasting cotton genotypes
© Peng et al.; licensee BioMed Central Ltd. 2014
Received: 30 April 2014
Accepted: 4 August 2014
Published: 5 September 2014
Cotton (Gossypium spp.) is one of the major fibre crops of the world. Although it is classified as salt tolerant crop, cotton growth and productivity are adversely affected by high salinity, especially at germination and seedling stages. Identification of genes and miRNAs responsible for salt tolerance in upland cotton (Gossypium hirsutum L.) would help reveal the molecular mechanisms of salt tolerance. We performed physiological experiments and transcriptome sequencing (mRNA-seq and small RNA-seq) of cotton leaves under salt stress using Illumina sequencing technology.
We investigated two distinct salt stress phases—dehydration (4 h) and ionic stress (osmotic restoration; 24 h)—that were identified by physiological changes of 14-day-old seedlings of two cotton genotypes, one salt tolerant and the other salt sensitive, during a 72-h NaCl exposure. A comparative transcriptomics was used to monitor gene and miRNA differential expression at two time points (4 and 24 h) in leaves of the two cotton genotypes under salinity conditions. The expression patterns of differentially co-expressed unigenes were divided into six groups using short time-servies expression miner software. During a 24-h salt exposure, 819 transcription factor unigenes were differentially expressed in both genotypes, with 129 unigenes specifically expressed in the salt-tolerant genotype. Under salt stress, 108 conserved miRNAs from known families were differentially expressed at two time points in the salt-tolerant genotype. We further analyzed the predicted target genes of these miRNAs along with the transcriptome for each time point. Important expressed genes encoding membrane receptors, transporters, and pathways involved in biosynthesis and signal transduction of calcium-dependent protein kinase, mitogen-activated protein kinase, and hormones (abscisic acid and ethylene) were up-regulated. We also analyzed the salt stress response of some key miRNAs and their target genes and found that the expressions of five of nine target genes exhibited significant inverse correlations with their corresponding miRNAs. On the basis of these results, we constructed molecular regulatory pathways and a potential regulatory network for these salt-responsive miRNAs.
Our comprehensive transcriptome analysis has provided new insights into salt-stress response of upland cotton. The results should contribute to the development of genetically modified cotton with salt tolerance.
KeywordsCotton Salt stress Leaf transcriptome Transcription factor MicroRNA
High soil salinity, a devastating environmental stress, typically causes major reductions in crop productivity and quality . Over 800 million hectares, equivalent to 6.5% of the world’s total land area, are currently estimated to be impacted by high salt concentration . In China, approximately 100 million hectares distributed over 16 provinces have been affected by increased salinity . To stabilize global crop production, the problem of salinity must be urgently addressed.
Salt stress generally induces a combination of dehydration/osmotic-related effects and damage as a consequence of excess sodium ions that greatly affect plant growth and crop production . Osmotic stress and Na+ stress are considered to be the two major components of the plant salt-stress response . Plants employ various mechanisms to deal with salt stress; these mechanisms include minimization of the amount of salt taken up by roots and its partitioning at tissue and cellular levels to avoid buildup of toxic concentrations in the cytosol of functional leaves [6, 7]. Much effort has been devoted to revealing the molecular mechanisms of plant salt tolerance, with the ultimate goal of improving salt tolerance of crop plants.
Cotton (Gossypium hirsutum L.) is not only the world’s leading textile fiber, but is also a major oil crop. Although cotton is the second most salt-tolerant herbaceous crop [8, 9], its growth and productivity are adversely affected by high salinity, especially at germination and the young seedling stage . Salinity suppresses primary root growth , inhibits the length and numbers of secondary roots , and limits photosynthesis and respiration, flowering, boll and fiber quality, and ion uptake in cotton, resulting in significant yield losses . Salt stress has also been found to regulate the expression level of many genes in different bio-processes and pathways, including morphological adaptation, maintenance of ion homeostasis, cell signal transduction, and oxidative stress mitigation [14–16]. Identifying salt-tolerance genes is an important component of the breeding of salt-tolerant crop plants through genetic engineering. Although many genes controlling response to high salinity have been identified in model plants, only a few salt stress-inducible genes, such as NHX1 (GhNHX1) , DREB (GhDREB1) , ERF (GhERF2- GhERF6) [19–21], NAC (GhNAC1-GhNAC6) , metallothionein (GhMT3a) , GhMPK2, GhMKK1, and CCCH-type zinc finger (GhZFP1)  have been documented in cotton. With recent advances in genomic sequencing and transcriptome mapping (microarray and high-throughput sequencing), some salt-related genes and regulatory factors have been identified in cotton on a large scale at the genome-wide level [27–33]. Nevertheless, the molecular basis of cotton tolerance to salt stress remains to be discovered.
Small RNAs are important post-transcriptional regulators of gene expression. miRNAs are a highly conserved class of endogenous, non-coding RNAs that range in length from 19 to 25 nucleotides (nt) . In plants, cleavage or translational repression of target mRNA appears to be the prevalent method of post-transcriptional regulation [35, 36]. Current findings that some plant miRNAs respond to stress conditions and that some miRNA targets are stress-related genes suggest that miRNAs play important roles in plant stress response. miRNA expression profiles in response to salt stress have been analyzed in Arabidopsis thaliana, Oryza sativa, Zea mays , Gossypium hirsutum, and Caragana intermedia. In addition to studying the effects of stress on miRNAs, the identification of miRNA targets is important. Among cotton cultivars, however, some miRNAs related to salt tolerance have been found to differ [42–45]. Little is known regarding their expression profiles in response to salt stress, and their roles in salt adaptation remain unclear.
To obtain insights into the initial transcriptional regulation and differential expression of functional genes in response to high salinity, we exposed two cotton genotypes (salt-tolerant Earlistaple 7 and salt-sensitive Nan Dan Ba Di Da Hua) to salt stress for 0–72 h. We then compared dynamic changes in ion accumulation and osmotic adjustment between Earlistaple 7 and Nan Dan Ba Di Da Hua and analyzed the expression profiles of responsive genes and miRNAs. According to their expression patterns, salt stress-inducible genes could be broadly discerned as early-, late-, and sustained-induced regulated genes, and their functional significance was characterized. We also identified and characterized conserved miRNAs in the two genotypes at two time points (4 and 24 h) under salt stress. Putative target genes of these miRNAs were predicted from the transcripts by computational methods. We established a model of the miRNA-mediated regulatory network by combining the mRNA-seq and small RNA-seq information. Our comprehensive physiological- and molecular-level analysis of salt stress responses should facilitate illumination of the mechanisms of salt stress tolerance in G. hirsutum.
Physiological differences between salt-sensitive and salt-tolerant genotypes under salinity stress
Genotypes Nan Dan Ba Di Da Hua and Earlistaple 7 were subjected to short-term salinity stress (200 mM NaCl) at the seedling stage for 4, 24, 48, and 72-h. Seedling response to salt was evaluated based on leaf relative water content (RWC), relative electrical conductivity (REC), and chlorophyll and ion (Na+ and K+) content of leaves and roots (Figure 1).
Highly significant (P < 0.001) differences were observed in leaf RWC and chlorophyll content between the two genotypes after 48- and 24-h salt treatments, respectively, although these values were sharply decreased in both genotypes after 4 h of exposure to NaCl compared with controls (Figure 1A, C). After 4 and 24 h of treatment, the REC of Nan Dan Ba Di Da Hua was highly significantly (P < 0.001) increased throughout the duration of stress, while that of Earlistaple 7 was maintained at a lower level. After 48 h, however, the RECs of the two genotypes were similar (Figure 1B). As shown in Figure 1D–F, leaf and root Na+ concentrations increased in both genotypes under salinity stress. Characteristically, Nan Dan Ba Di Da Hua roots accrued less Na+ as compared with Earlistaple 7 (Figure 1D), whereas Earlistaple 7 contained a significantly lower level of leaf Na+ (Figure 1E, F). This difference can probably be attributed to the higher Na+ absorption capacity of roots of salt-tolerant Earlistaple 7 compared with salt-sensitive Nan Dan Ba Di Da Hua, which prevents the transport of excessive Na+ to the seedling shoots.
Visual damage owing to salinity stress appeared on the leaves of Nan Dan Ba Di Da Hua beginning approximately 24 to 72 h after 200-mM NaCl treatment. After 0.5 h, distinct wilting and dehydration caused by osmotic stress  was observed on leaves of both genotypes. By 4 h, both genotypes showed more severe wilting, which was followed by gradual recovery to an upright position (24 h). The restoration of Earlistaple 7 was obviously better than that of Nan Dan Ba Di Da Hua, which experienced slight wilting of cotyledons and old leaves, and dry necrosis of leaf edges. After 48 h, these phenomena were markedly more noticeable in both genotypes, although plants of Earlistaple 7 were in much better condition than those of Nan Dan Ba Di Da Hua. According to the above results, significant physiological and morphological differences were observed at different salt-treatment time points in response to salt stress, which was accordingly divided into two typical phases: dehydration stress (4 h) and ionic stress (24 h). We therefore selected these two time points for further study.
De novomRNA-seq assembly and annotation of non-redundant unigenes
Six cDNA and six small RNA libraries were constructed from total RNA extracted from 14-day-old seedling leaves of two upland cottons (cultivars Nan Dan Ba Di Da Hua and Earlistaple 7) treated with 200 mM NaCl (or water as a control) for 4 and 24 h.
The data quality of mRNA-seq and the transcripts in two genotypes
Nan Dan Ba Di Da Hua
Length distribution of the transcripts and unigenes clustered from the de novo assembly
All (> = 200 bp)
> = 500 bp
> = 1000 bp
The numbers and distribution rate of unigenes in the databases of NR, SWISS-PROT, TrEMBL, CDD, PRAM, KOG and KEEG
Identification of salt-responsive, differentially expressed unigenes
The numbers of all differentially expressed unigenes (DEUs) at different salt-stressed time points of two genotypes
Nan Dan Ba Di Da Hua
0 versus 4 h
0 versus 24 h
0 versus 4 h
0 versus 24 h
Dynamic change analysis of co-expressed unigenes
To understand dynamic expression patterns of DEUs under salt stress at different time points, the various gene expression modes in Earlistaple 7 were identified. The most representative significantly DEU expression mode was then selected for further analysis. To uncover the mechanisms underlying salt sensitivity in Nan Dan Ba Di Da Hua, we used Short Time-series Expression Miner (STEM) software to analyze abundance changes of differentially expressed TF unigenes and other significant functional unigenes filtered by different genotypes at different time points. The analytical procedure used is outlined in Additional file 3.
Summary of the numbers of differentially co-expressed unigenes from the cluster 1 and cluster 2 by STEM analysis
ck, 4 h
ck, 24 h
ck, 4 h, 24 h
72,919(78.2% ) d
The classification and expression pattern of the identified differentially co-expressed unigenes
I ck,4 h
II ck,24 h
III ck,4 h,24 h
IV ck,4 h
V ck,24 h
VI ck,4 h,24 h
Early induced up-regulated
Late induced up-regulated
Continuous induced up-regulated
Early induced down-regulated
Late induced down-regulated
Continuous induced down-regulated
N-E common unigenesa no.
E specific unigenesb no.
Analysis of salt stress-responsive TFs and other salt tolerance-related functional genes
Although expression patterns of many unigenes were similar between Earlistaple 7 and Nan Dan Ba Di Da Hua, 78.2% of salt-responsive unigenes identified in Earlistaple 7 were more or less differentially expressed (Nominal false discovery rate (FDR) < 0.001) in Nan Dan Ba Di Da Hua at one or more time points (Table 5). TFs are important regulators of gene expression. The related the expression abundance and annotation information of TF unigenes in the six data sets by STEM analysis were obtained (see Additional file 2C-H).
The specific transcription factor unigenes response to salt stress in Earlistaple 7
Earlistaple 7 Log2 ratio
NAC domain-containing protein 7
Zinc finger protein 2
AP2/ERF and B3 domain transcription factor RAV
Probable WRKY transcription factor 27
Zinc finger protein 6
CCAAT-binding transcription factor. Putative
Zinc finger CCCH domain-containing protein 33
Zinc finger CCCH domain-containing protein 29
Homeodomain-leucine zipper protein HD2 GhHB2
Putative dehydration responsive GhDREB1A
Probable protein phosphatase 2C 1
Ethylene-responsive transcription factor ERF105
Zinc finger CCCH domain-containing protein 29
Nuclear transcription factor Y subunit C-2 NFYC2
Probable WRKY transcription factor 11
WRKY transcription factor. Putative
MADS-box transcription factor 50
WUSCHEL-related homeobox 4
Probable transcription factor KAN4
Auxin response factor 3 GhARF3
Zinc finger protein CONSTANS-LIKE 5
Homeobox protein GhHB1
NAC domain-containing protein 42
Zinc finger protein ZAT5
Scarecrow-like protein 9
Transcription factor. Putative
TT2 like MYB transcription factor
Homeobox-leucine zipper protein ATHB-14
Dof zinc finger protein DOF4.6
ZF-HD homeobox protein At4g24660
Zinc finger protein-like Ser/Thr protein kinase-like protein
Probable protein phosphatase 2C 38
Transcription factor. Putative
WRKY transcription factor PmWRKY110 (Fragment)
BZIP-like protein Gossypium hirsutum
Under NaCl stress, most of the salt-responsive TF unigenes exhibited genotype specificity. We identified 129 TF unigenes whose expressions were specific to Earlistaple 7, of which 85 were up-regulated and 44 were suppressed under salt stress (see Additional file 2 F-H). These 129 TF unigenes, which corresponded to 1.6% of the 8,127 salt-responsive genes of Earlistaple 7, were only expressed during either early or late stages of salt stress. We note that some TF unigene families expressed under both 4- and 24-h salt stress were identical to the above-mentioned TF families commonly expressed in both genotypes, but the expressed TF was different. The expression trend of the WRKY TF unigene (homologous to PtWRKY25) was different at the two time points—decreasing at the beginning and then increasing, while that of others, such as AP2-EREBP, MYB, HB, NAC, WRKY, and bHLH family members, experienced a continual rise.
Many other salt tolerance-related unigenes encoding transporters, Ca2+ binding protein, heat-shock proteins (HSPs), detoxificants, and dehydration-response proteins were differentially expressed under 4- and 24-h salt stress (see Additional file 4). Several differentially expressed genes encoded transporters contributing to the re-establishment of ionic and osmotic homeostasis under salt stress, such as ABC transporters (i.e., ABC transporter B, C, and G family members), sodium/hydrogen exchanger (NHX1), ATPases (Ca2+-ATPase, H+-ATPase, and V-type H+-transporting ATPase), potassium transporter, and nitrate transporter. Genes encoding some Ca2+ sensors, namely CML and CPK, and one CBL-interacting protein kinase (CIPK) gene were identified in the different libraries. Two cotton CBL genes (GhCBL1 and GhCBL8) were highly expressed under 4- and/or 24-h salt stress in Earlistaple 7. Unigenes encoding CDPK and CIPK members were up-regulated under 4- and/or 24-h salt stress in both genotypes. Genes from several HSP families (e.g., HSP20, HSP70, Hsp90, and DanJ) as well as genes related to reactive oxygen species (ROS) scavenging and detoxification (CAT, POD, GPX, ALDH, and RBOH) and dehydration response (LEAs and RING/U-box domain-containing protein/XERICO) were highly up-regulated at different salt-treatment time points either in Earlistaple 7 or Nan Dan Ba Di Da Hua (see Additional file 4). Finally, genes involved in salt-stress tolerance pathways, such as CDPK, starch and sucrose metabolism, carbohydrate metabolism, and MAPK signal pathways, were found to be induced at different time points. These identified unigenes are potentially salt-induced genes and may serve as a supplementary resource for salt-tolerance gene mining.
Deep sequencing of small RNAs and identification of conserved miRNAs under salt stress
To characterize differences between Nan Dan Ba Di Da Hua and Earlistaple 7 at the miRNA regulatory level and to identify conserved miRNAs in cotton under salt stress, we constructed and analyzed six small RNA libraries of the two genotypes. A total of 48,585,151 sequence reads ranging in length from 17 to 35 nt were generated. After removing adapter sequences, low-quality tags, and reads shorter than 18 nt, we obtained 22,681,645 clean reads, which corresponded to more than 2 million unique sequences for each sample (see Additional file 6A). As shown in Additional file 6B, most small RNAs were 20–24 nt in size; the 24-nt class was the most heavily represented (50.56%) in the Nan Dan Ba Di Da Hua small RNA control library, followed by 23-nt (7.37%), 21-nt (10.27%), and 22-nt classes (7.14%).
Distribution of sequence reads among different RNA categories in six libraries of cotton
Unique sRNA number
Total sRNA number
The number of identified miRNAs and the abundance of miRNA families are listed in Additional file 6C–D and Additional file 7. The reads for each miRNA family were normalized as “reads per million mapped reads” (RPM). Among miRNA families, miR166 exhibited the highest abundance (3,909,410 RPM). Other highly expressed families, including miR396, miR482, miR159, miR1697, miR394, and miR156, were represented by approximately 100,000 RPM, whereas miR319, miR2118, miR858, miR828, miR829 and miR1310 families had less than 100 RPM. Differences among miRNA families were also revealed by the number of members they contained (see Additional file 6D and Additional file 7). miR166 was the largest family, with 29 members, while miR159, miR396, miR171, and miR156 possessed 28, 26, 23, and 20 members, respectively. Different family members also displayed drastically different expression levels (see Additional file 7). In the miR156 family, for instance, RPM of ghr-miR156a was 308 in the control Nan Dan Ba Di Da Hua library, while the RPM of cca-miR156b was only 2.6.
Combined analysis of conserved miRNAs and their target genes under salt stress
Prerequisites for identification of candidate salt tolerance-related genes and miRNAs in cotton
Two main factors are responsible for salt stress-induced inhibition of plant growth: osmotic stress and ionic stress . Different mechanisms control development at different time points in plants exposed to salinity. In this study, we monitored the physiological responses of two cotton genotypes with differing salt tolerances at two distinct phases. We found that these phases were characterized by increases or decreases of clusters of differentially expressed genes and miRNAs.
The positive relationship between water content and Na+ concentration causes ions to accumulate in vacuoles, which lowers water potential, increases the driving force for water uptake, and increases turgor . Maintenance of K+ and Na+ homeostasis is crucial under salt stress . Our study results suggest that after being taken up by the roots, Na + is transported to shoots through the stems and the roots of salt-tolerant genotype Earlistaple 7 were able to retain more Na+ than Nan Dan Ba Di Da Hua (Figure 1D, E). To avoid excessive ion accumulation in roots owing to continued Na+ uptake, compensatory decreases in other cations occurred in the roots at 24 h. Shoot leaves accumulated Na+ in addition to the cations already present, with the tolerant genotype thereby experiencing the strongest ionic stress beginning at 24 h (Figure 1F).
By examining these dynamics, we gained insights into physiological responses of the NaCl-acclimation process, which consisted of an initial dehydration phase in leaves (0–4 h) and subsequent NaCl accumulation (4–24 h) followed by restoration of osmotic homeostasis at a new ionic stress level (24 h) and either final adjustment to a steady ion balance or ion-induced damage (24–72 h and beyond). The timing of these adaptive responses differed between roots and leaves. In previous studies, either roots or leaves were used to identify differences between early and late transcriptional responses to salt treatment in cotton [31–34]. Our study of two genotypes with contrasting salt tolerances revealed that Na+ accumulated mainly in shoots, and that significant physiological and morphological differences occurred at different time points under salt stress: 4 h (owing to dehydration stress) and 24 h (owing to ionic stress). These two important time points thus corresponded to the two typical phases of plant growth response to salt stress.
In some studies, genes differentially expressed in response to salt treatment have been identified using a single plant genotype [47–49]. Two or more genotypes with contrasting salt tolerance have only been used in a few reported investigations [46, 55–57]. To identify genes responsible for salt tolerance, analysis of genotypes with similar genetic backgrounds and contrasting salt tolerances is preferable. Identification of common genes that are differentially expressed under stress conditions between tolerant and sensitive genotypes with different genetic backgrounds has been suggested as an alternative approach . As pointed out in , it is impossible to isolate stress tolerance-related genes from stress-responsive ones unless transcriptional-level differences are compared between tolerant and sensitive genotypes under stress conditions . Genes related to salt tolerance can thus be identified in genotypes with different genetic background and contrasting tolerances by comparative analysis of differentially expressed genes and miRNAs at different salt-treatment time points. We therefore used a pair of salt-tolerant and salt-sensitive genotypes in this study to identify differentially expressed genes and miRNAs.
Abundance of differentially expressed TF unigenes in response to salt stress in tolerant and sensitive genotypes
To govern necessary transcriptional regulation in response to developmental and environmental conditions, plant genomes encode a large number of TFs. In the present study, we identified 3,202 TF unigenes (see Additional file 2A) by global TF classification of DEUs. Many TFs in the AP2/EREBP family, followed to a lesser extent by WRKY, MYB, NAC, DBP, bHLH, C2H2, C3H, C2C2, and CCAAT family members, were up-regulated during salt stress (Figure 2 and Additional file 2B). Some members of these families have been reported to participate in abiotic stress response, suggesting their roles in plant adaptation to stress [18, 59–62]. Seven unigene clusters were significantly assigned as comprising co-expressed genes according to STEM (Figure 3A). The mechanism responsible for these gene co-expression clusters is unclear, although the co-expression of clustered genes may be partially regulated by the interaction of common elements and TFs . A previous promoter analysis has suggested that some of these co-expressed genes may be regulated by common TFs . In light of the above views, we separated out the differentially co-expressed unigenes from all differentially expressed genes encoding TFs to identify the types of TFs regulating response to 4- and 24-h salt stress (Table 5 and Additional file 2C–H).
Salt-responsive TF unigenes in the two cotton genotypes
Dramatic changes were observed in TF unigene abundance and types (common and specific) at 4- and 24-h salt-stress time points in the two contrasting genotypes. Among these changes, 478 TF unigenes (groups III and IV in Table 6) exhibited similar patterns of altered expression between the two genotypes. This dramatic variation is consistent with the plant switching from normal growth and development to stress-specific responses. Our observation is the first reported up-regulation of GhbZIP, GhDREB1C, GhMYB, GhMYB2, and GhWAKY under 4- and 24-h salt-stress treatment. AtRAP2-12 and AtRAP2-3, homologs of cotton genes GhERT and GbERF2, are reported to increase rapidly in Arabidopsis under salt stress [65, 66]. Our results illuminate the roles of these genes in response to salt stress (in addition to abiotic stresses). In an earlier study, bHLH family member GhbHLH1 was transitorily induced in leaves by abscisic acid (ABA) and polyethylene glycol treatments . In our study, GhbHLH1 was significantly up-regulated (see Additional file 2E), indicating its additional involvement in cotton salt-stress response. Seven TF genes (three GhNAC, two GhHB and two GhMYB genes) were identified at 4- and 24-h salt-stress time points, suggesting that their protein products may play an important role in salt-stress response in addition to the involvement in ABA signaling uncovered in previous research [22, 68, 69]. Transcripts encoding C2H2-type zinc finger proteins (ZAT10) that have been implicated in salt stress tolerance were up-regulated under 4- and 24-h salt stress. Overexpression of C2H2-type zinc finger proteins (STZ/ZAT10 homologs) has been shown to induce the expression of several stress genes related to salt tolerance, dehydration, and cold stress [70, 71]. GhMADS7 exhibited higher expression levels in the two cotton genotypes upon 4 and 24 h of salt treatment. Although GhMADS7 is a well-known participant in cotton fiber elongation, salt-stress response is a novel function for this TF; alternative splicing may have altered its cellular role .
Although 115 “osmotic stress”-specific TF unigenes were observed after 4 h of stress treatment, most of the 226 “salt stress”-specific TF unigenes appeared to be ion-specific at 24 h (groups I, II, IV, and IV in Table 5). These two salt-stress treatment phases clearly triggered very different responses in the two contrasting genotypes. Another finding of our study was that 37 unigenes belonging to the AP2/EREBP family of TFs were also related to salt tolerance. DREB/CBF subfamily genes, such as DREB2A[73, 74] and DREB2B, are reportedly involved in osmotic (4 h) and ionic (24 h) stress. Interestingly, we also found that some type-B ARRs (e.g., AtPRR7) were induced or repressed at 4- and 24-h time points (Additional file 2C, D). These general salt-responsive genes identified at 4 h (osmotic stress) are thus similar to those identified at 24 h (ionic stress). Because those unigenes were expressed in both tolerant and sensitive genotypes, they may not be directly responsible for salt tolerance.
Salt-responsive TF unigenes in the salt-tolerant cotton genotype
We detected 129 unigenes—48 at 4 h, 42 at 24 h, and 39 at both time points (Table 5 and Additional file 2 F, G, H)—that were differentially expressed under NaCl treatment in the salt-tolerant genotype but not in the salt-sensitive one. In addition, unigene expression in the tolerant line Earlistaple 7 differed dramatically between 24- vs. 4-h salt stress (Table 5). One possible explanation for this temporal variation is that the expression of upstream response genes (including TFs) was gradually reduced during osmotic stress recovery, with the upstream cascade genes then activated during ionic stress. Other putative TFs related to salt tolerance that were identified in this study included AP2-EREBP, MYB, C2C2-Dof, bZIP, and bHLH (see Additional file 2). Previous research has also demonstrated that some differentially expressed TF genes are associated with plant salt tolerance, including GhHB2, GhDREB1A, AtC3H29, AtPP2C66, and OsC3H33. The discovery in this study of many TFs in various families indicates that a complicated transcriptional regulation network is involved in response to high salinity in cotton.
Differentially expressed regulators and functional unigenes in tolerant and sensitive genotypes under salt stress
CDPK, SOS, and MAPK pathways are important signaling pathways associated with salt stress response . CaMs/CMLs, CDPKs, and CIPKs represent the three major classes of plant Ca2+ sensors that are involved in many developmental and stress-induced signal transduction pathways . In the present study, CDPKs (homologous to AtCPK7and AtCPK30) were mainly up-regulated at 4 and 24 h in the salt-tolerant genotype, with their expression remaining at low levels over 24 h of salt stress in the salt-sensitive genotype. We also found that a CPK5 gene, homologous to PtCPK5, was up-regulated at 4 and 24 h. GhCDPK32, homologous to AtCPK32, was identified and exhibited peak expression at 24 h of treatment (see Additional file 4). AtCPK32 overexpression has been reported to affect the expression of several ABF4-regulated genes and to be involved in ABA/salt-stress response . Furthermore, a GhCPK5 gene was cloned and found to be responsive to salt tolerance . Overexpression of OsCPK7 and AtCPK30 enhances tolerance to cold and to salinity and drought [85, 86]. The observed up-regulation of MEKK1-MPK4/6 at 4 and 24 h of salt treatment indicates that the MAPK signaling pathway mediates salt-stress signaling in cotton (Figure 11). In addition, GhMPK6 was up-regulated under salinity at two time points in the salt-tolerant genotype, while its expression remained unchanged in the salt-sensitive genotype. GhCAT1 was also over-expressed both at 4 and 24 h in the salt-tolerant genotype (see Additional file 4). The genes in the MAPK signaling pathway likely exist in cotton as well. The role of the MEKK1-MKK2-MPK4/MPK6 cascade in cold- and salt-stress signaling has been demonstrated in Arabidopsis . Luo et al.  have suggested that GhMPK6 plays a major role in ABA-induced CAT1 expression and H2O2 production. Although GhMPK7 has an important function in salicylic acid (SA)-regulated broad-spectrum resistance to pathogen infection , it was a novel MAPK involved in salt-stress response. The expression of genes involved in MAPK and Ca2+ signaling pathways, along with other kinases, at various time points under salt stress suggests that the induction of signaling pathways can activate responsive targets for adaptation to salt stress in cotton.
Plant hormones such as ABA, ethylene (ET), SA, and jasmonate (JA) are essential for plant adaptation to abiotic stress conditions . In our study, cotton unigenes associated with ABA and ET signaling pathways showed differential expression at different time points under salt stress (Figure 11 and Additional file 4). Four genes homologous to AtNCED3, AtHVA22A, AtHVA22E, and AtEDL3, which are related to ABA synthesis or regulation, were up-regulated after 4 and 24 h of salt treatment in the salt-tolerant line. Up-regulation of PLC2, ABI1/2, and HAB1/2 at 4 and 24 h, PLC4 at 4 h, and PLC1 at 24 h in the salt-tolerant genotype implied that ABA signaling pathway induction occurred under salt stress (Figure 11). As one of the main plant hormones, ABA regulates numerous developmental processes and adaptive stress responses in plants. Recent studies have shown that 9-cisepoxycarotenoid dioxygenase (NCED3) is the key enzyme involved in ABA biosynthesis . AtHVA22 is differentially regulated by ABA, cold, dehydration, and salt stresses . In A. thaliana, strong and rapid induction of the EDL3 gene has been observed under osmotic stress and high salinity . Phospholipase C (PLC) and protein phosphatase 2C (PP2C) are involved in the ABA signaling pathway . During salt stress (between 4 and 24 h), genes related to ET biosynthesis (homologs of ACS1, ACS12, ACO1, and ACO3) were also up-regulated. ET receptor (ETR1, ETR2, and EIN4), ET signaling pathway (EIN3, ERF1, ERF2, and MEKK1-MKK2-MPK4/6 kinases), and feedback mechanism (CTR1) genes were also up-regulated at 4- and/or 24-h time points (Figure 11 and Additional file 4). ET regulates stress and defense responses along with many key events of plant growth and development. It binds to membrane protein receptors (ETR1, ETR2, ERS1, ERS2, and EIN4) and activates the MAPK cascade (MKK9-MPK3/MPK6) which, in turn, stabilizes the EIN3 protein. EIN3 activates the expression of ERF1 involved in secondary transcription in the signaling pathway, while CTR1 (belonging to MAPKKK) acts as a negative regulator and degrades EIN3 . In our study, only AtMPK6 was detected. We suspect that additional genes exist that were undetected. An earlier study also revealed that ET is involved in drought-stress tolerance as well as cotton fiber development . Our results suggest the occurrence of multiple-hormone cross-talk in response to osmotic and salt stresses at different time points.
Some salt-responsive downstream genes, especially those involved in antioxidant system and transport functions (ion translocations across plasma and vacuolar membranes), may also provide a well-characterized clue to understanding the mechanisms responsible for salt tolerance. Under salt stress, several transporters of toxic metal ions, lipids, water, and sugar molecules were differentially expressed, such as ABC transporters, H+-ATPase, sucrose transport, inositol transporter, gated channel, sulfate transporter, NHX1, H+ antiporter, chloride channel, and potassium channel proteins. Another aspect of salt tolerance is the regulation of ROS to protect membranes and macromolecules. In our study, several genes related to detoxification were up-regulated at various time points under salt stress. These expressed genes were homologous to AOMI, ALDH, NMT1, ANN1, ALA1, AAP3, VAKP722, POD, CAT, RBOH, GR/GPX, and GST[97–99].
Mediation of a potential regulatory network via conserved miRNAs in the G. hirsutumsalt-tolerant genotype
Over the past few years, miRNAs have been reported to be prominent gene regulatory factors in plant tolerance to environmental stresses such as drought, cold, heat, and high salinity . Much effort has been devoted to understanding their role in salt-stress responses of various plants including Zea mays, Populus euphratica, Gossypium hirsutum, and Thellungiella salsuginea.
In our study, salt-responsive miRNAs were detected in both cotton genotypes. A majority of the responsive miRNAs were either detected only in the salt-tolerant genotype or had distinct and inverse expression trends at 4- or 24-h time points between the two cotton genotypes (Figure 8B; Additional file 9A). This dichotomy suggests that common-regulated miRNAs may be part of the fundamental mechanism of salt-stress adaptation, while differentially expressed miRNAs are responsible for the distinct salt sensitivities between the two cotton genotypes (see Additional file 10). Our results are consistent with a previous hypothesis that highly expressed miRNAs are primarily responsible for control of basic cellular and developmental pathways common to most eukaryotes, whereas lower-expressed miRNAs are involved in regulation of lineage-specific pathways and functions .
Because miRNAs regulate the specific genes targeting mRNAs for degradation or translation inhibition, identification of the potential target unigenes of miRNAs is crucially important for understanding miRNA-mediated processes such as salt tolerance in plants. According to the PCC and linear regression analysis, the candidate negative correlation miRNA/target genes were further analyzed (Figure 7). It is noted that there is four clumps along the X-axis obviously: miRNA strongly down, slightly down, slightly up, and very strongly up (Figure 7B, E). Since some miRNAs were strongly repressed in salt stress than that in control condition,and some miRNAs were strongly inducted in salt stress compared with that in control condition, the log2(E4/E0) and log2 (E24/E0) value of some miRNAs in different treatment were more than 10 or less than -10. These miRNAs may be play an important role in response to salt stress. The identified differentially expressed miRNAs which their log2 (E4/E0) and log2 (E24/E0) value range of -2 ~ -1 and 1 ~ 4 can regard as slightly down and slightly up-regulated miRNAs. Therefore, the data points are separated into four clumps along the X-axis. There is no doubt that all these candidate negative correlation miRNA/target genes should be emphasized. Here, we found that the differentially expressed, conserved miRNA targets of salt-tolerant Earlistaple 7 generally encoded TFs, enzymes, and stress response, transport, and signal transduction proteins (Figure 9).
Many miRNA target genes encode mRNAs of TFs, indicating upstream regulation of miRNAs during developmental processes and environmental responses [104, 105]. Some miRNAs, including miR156, miR164, miR166, miR171, miR172, miR393, miR396, and miR397, exhibited altered expression profiles in at least one genotype in our study. Based on transcriptome data prediction results, they regulate the expression of TFs (Figure 9). In this study, three miR156s (bna-miR156a, osa-miR156f-3p, and gma-miR156q) and one miR172 (zma-miR172a) were up-regulated and down-regulated, respectively, under 24-h salt treatment in the tolerant line Earlistaple 7. In addition, three of the above miRNAs may function by changing normal levels of downstream gene transcripts through modulation of SPL5 and RAP2-7 expression (see Additional file 11). miR156 and miR172 are also involved in the regulation of flowering time and floral development in Arabidopsis, which they accomplish by negatively regulating SBP-LIKE proteins (SPLs) and AP2-like factors, respectively . These results suggest that the three conserved miRNAs may play significant regulatory roles in shoot maturation and the vegetative- to reproductive-phase transition. The NAC family is involved in embryo and shoot meristem development, lateral root formation, auxin signaling, and defense and abiotic stress responses . The conserved miRNA ptc-miR164f, which was verified by stem-loop RT-PCR to target the NAC domain-containing protein (NAC021/AtNAC1), was highly down-regulated at 4 and 24 h of salt stress in Earlistaple 7 (see Additional file 11 and Figure 10B). We observed that the CUP-SHAPED COTYLEDON 2 gene (comp81541_c0_seq1; homologous to AtNAC098), which is targeted and negatively regulated by at-miR164a, was especially down-regulated only at 24 h in Earlistaple 7 (see Additional file 11 and Figure 9). This result suggests that the different miRNA164s may negatively regulate different NACs during at the early stage of salt stress of salt stress. The No Apical Meristem (NAC) protein is involved in plant hormonal control and defense [108, 109]. In Arabidopsis, AtNAC1-overexpressing lines are bigger than the wild type, with larger leaves, thicker stems, and more abundant roots , suggesting that NAC1 may be an early auxin-responsive gene. CUP-SHAPED COTYLEDON 2 (a NAC family member) and at-miR164A are transcribed in overlapping domains at the margins of young leaf primordia, with transcription gradually restricted to the sinus where the Arabidopsis leaf margins become serrated . Most often, the regulatory pathway consisting of miR393 that targets TIR1, a negative regulator in auxin signaling, may participate in lateral root development in Arabidopsis . The other target was a bHLH transcription factor, which was involved in many biological processes, but there were no phenotypic changes and functional mechanisms owing to miR393/bHLH interaction in Arabidopsis . In our study, all miR393s (ath-miR393a, osa-miR393a, csi-miR393, and osa-miR393b-5p), which were highly expressed at 4 and 24 h of salt stress in salt-tolerant Earlistaple 7 but were unaffected or down-regulated in salt-sensitive Nan Dan Ba Di Da Hua, targeted and tested as negatively related to a bHLH78 homolog unigene (comp105131_c0_seq1) by stem-loop RT-PCR (see Additional file 11 and Figure 10B). Induction of miR393 at the early stage of salt stress (4 and 24 h in the salt-tolerant genotype) suppressed the expression of its target AtbHLH78.In our proposed miRNA-mediated regulatory network based on miRNA-TF interactions (Figure 9), one miRNA may regulated more than one TF gene. In addition, differentially expressed conserved miRNAs may also target non-TF genes. Comparison of target genes of conserved miRNAs illustrated in the regulatory network suggests that the same miRNA response to salt stress at different time points may regulate the same or different genes to exert different functions. Nevertheless, these putative interactions need to be confirmed by experimental data.
In this study, we first generated a global transcription map of genes and miRNAs expressed in leaves of salt-tolerant cotton under salt stress. We detected 29,136 and 33,492 unigenes exhibiting differential co-expression after 4 and 24 h of salt stress, respectively. GO annotations demonstrated that unigenes belonging to a wide range of functional categories are involved in salinity defense response at the two different time points. We also identified 320 miRNAs conserved across the two genotypes. Of these, 108 members of miR156, miR159, miR164, miR167, miR171, miR172, miR393, miR396, miR397, and miR482 families were differentially expressed at 4 and 24 h under salt stress in the salt-tolerant genotype. Some miRNA targets were predicted to encode abiotic stress-responsive TFs and enzymes. Expressions of some selected differentially regulated miRNAs and their targets were found to be negatively regulated. Comparative transcriptome analysis of salt-tolerant and salt-sensitive genotypes indicated that essential salt tolerance-related genes and miRNAs were differentially expressed and were only weakly or non-responsive to salt stress in the salt-sensitive cotton. Finally, we have provided a list of salt tolerance-related genes that will facilitate candidate gene discovery and molecular marker development for salt-tolerant breeding in cotton.
Plant materials and NaCl treatments
Plant materials comprised two upland cotton (G. hirsutum) genotypes: salt-tolerant Earlistaple 7, which was selected from Earlistaple 808 and introduced by the AR-SEA-USDA Pee Dee Experiment Station in the United States in 1982, and salt-sensitive Nan Dan Ba Di Da Hua, introduced from Guangxi Province, China, before 1978. Seeds of the two cotton genotypes were available in the National Mid-term Genebank of the Institute of Cotton Research, Chinese Academy of Agricultural Sciences (ICR-CAAS) after signing the Material Transfer Agreement (MTA) . Seedling cultivation and all experiments were carried out at the State Key Laboratory of Cotton Biology, Anyang, China.
Hand-selected seeds (120 of each genotype) of Nan Dan Ba Di Da Hua and Earlistaple 7 were surface-sterilized in 70% (v/v) ethanol for 15 s and 4% (w/v) sodium hypochlorite for 15 min, and then rinsed several times with sterile distilled water. The seeds were submerged in sterile water for 12 h at room temperature and then sown in sterile silica sand. After 3 days, germinated seedlings were transferred to hydroponic containers (600 × 350 × 120 mm) containing half-strength Hoagland’s solution (pH 6.0)  and grown in a phytotron (KR-III; Zhengzhou Henan Ke-rui Inc.). Growth conditions were 28/22 day/night temperature, 60–80% relative humidity, and a 14 h/10 h light/dark cycle under 450 μmol⋅m-2⋅s light intensity. The culture solutions were changed after about 14 days, when seedlings had produced three leaves. Seedlings showing normal growth were randomly divided into two groups; one group was placed into tanks filled with a 200 mM solution of NaCl, and the remaining seedlings were transferred to tanks filled with plain water to serve as control. After exposure to the two solutions for 4, 24, 48, and 72 h, seedlings of control and treated Nan Dan Ba Di Da Hua and Earlistaple 7 were harvested directly into liquid nitrogen and stored at -80°C until used for RNA extraction.
Leaf RWC, REC, chlorophyll content, and ion measurements
Plant fresh weight (FW) was measured immediately after harvest. The leaves were floated on deionized water for 8 h at 4°C. The turgid leaves were quickly weighed (TW) and their dry mass (DW) was measured after oven-drying at 105°C for 10 min followed by 80°C for 24 h. RWC was calculated as follows: RWC (%) = (FW - DW) / (TW - DW) × 100 . For REC, 0.5 g of fresh leaves were rinsed 3 times with ddH2O and were placed in an Erlenmeyer flask containing 40 mL of ddH2O incubated at room temperature for 4 h. The electrical conductivity of the solution (C1) was measured using an EM38 conductivity meter (ICT international, Armidale, NSW, Australia). The solution was boiled for 10 min and cooled to room temperature, and electrical conductivity (C2) was re-measured. REC was calculated as C1 / C2 × 100% . To measure leaf chlorophyll, 0.2 g of seedling leaves were incubated in 80% acetone in darkness at 4°C overnight. After centrifugation at 5000 × g for 5 min at 4°C, the absorbance of the supernatant was measured at 663 and 645 nm with a DU800 spectrophotometer (Beckman Coulter, Brea, CA, USA). Chlorophyll content was calculated as described by Porra (1989) . Concentrations of Na+ and K+ in roots and leaves were determined by inductively coupled plasma-optical emission spectrometry (ICP-OES) (Optima 2100 DV; Perkin-Elmer, Wellesley USA) according to the manufacturer’s instructions. Fifteen shoots per genotype were analyzed (3 from control and 12 from salt-treated plants). To ensure more accurate measurements, five independent plants were pooled as a sample (about 0.2–0.6 g) to reduce the effects of biological variation. For the Na+ and K+ determinations, the samples were completely digested in 12 mL 65% HNO3 and then analyzed on the ICP-OES spectrometer against National Institute of Standards and Technology traceable standards and controls. For the physiological experiments, three or more independent biological replicates of each control and time-treated sample were performed.
RNA extraction and qRT-PCR
Total RNA was isolated from cotton leaf samples treated with either 200 mM NaCl or water (as a control) for 4 or 24 h using a modified CTAB method and purified using Qiagen RNeasy columns (Qiagen, Hilden, Germany). First-strand cDNA was synthesized using a PrimeScript RT reagent kit with gDNA eraser (Takara, Dalian, China) according to the manufacturer’s instructions. Primers are listed in Additional file 5A. qRT-PCR was performed on a 7500 Fast Real-Time PCR system (Applied Biosystems, Inc., California USA). The 20-μL reaction solutions contained 10 μL SYBR Green Premix Ex Taq (Takara), 1 ng cDNA sample, and 0.2 μM gene-specific primers. Three replicate PCRs were run per cDNA sample, and three cDNAs from different samples were used as biological replicates. Relative expression levels were calculated by the 2-ΔΔT method  using GhActin (GI: AY305733) as an internal standard to normalize cDNA content .
The expression profiles of 10 salt-responsive mature miRNAs were assayed by stem-loop RT-PCR. Stem-loop reverse-transcription primers were designed following the method described by Chen et al.  and Varkonyi-Gasic et al. . Their 3′ ends, which were complementary to the 6-nt 3′ end of the miRNA, were combined with the 44-nt sequence 5′-GTCGTATCCAGTGCAGGGTCCGAGGTATTCGCACTGGATACGAC-3′ (Additional file 5B) to constitute the primer. Subsequent PCRs used a 5′ primer matching the ~18 nt at the 5′ end of the target miRNA. The 3′ primer was a universal reverse-transcription primer. The stem-loop reverse transcription reactions were performed using a PrimeScript RT reagent kit with gDNA eraser (Takara) according to the supplier’s manual. U6 snRNA was used as a reference for miRNA expression validation. Stem-loop RT primers and miRNA-specific PCR primers are listed in Additional file 5B, C. As in the qRT-PCR analysis, the 2-ΔΔT method was used to calculate relative expression levels of salt-responsive mature miRNAs.
Transcriptome sequencing and de novoassembly
Poly(A) mRNA was isolated from 10 μg total RNA of each leaf sample using oligo (dT) beads. The mRNA was sheared into short fragments with fragmentation buffer and selected for construction of Illumina RNA-Seq libraries according to the described protocol. Each library had an insert size of approximately 200 bp and was sequenced on a HiSeq 2000 system (Illumina, California USA). Leaf-sample read lengths were 90 bp. The raw image data was transformed by base-calling into sequence data and stored in FASTQ format. To control for the influence of error rate in the Solexa data results, quality pretreatment of the raw data was performed: the sliding window method was first used to remove low-quality fragments (quality threshold = 20, error rate = 1%, window size = 5 bp, and length threshold = 35 bp), followed by removal of partial (<35 bp) sequences containing “N” from the reads. After filtering of low-quality and dirty raw reads, transcriptome de novo assembly was carried out using the paired-end assembly method as implemented in Trinity (http://sourceforge.net/projects/trinityrnaseq/files/; trinityrnaseq_r2012-10-05) . After removal of repeats from the spliced sequences, transcripts were obtained of the expected length and size. The longest transcript in each locus (“comp * _ * _c”; Chrysalis Clusters module) was taken as the unigene. To assign possible annotations, all unigene sequences were searched using Blastx (e-value < 10-5) against the following protein databases: NR (in NCBI), SWISS-PROT (UniProt), TrEMBL (http://www.bioinfo.pte.hu/more/TrEMBL.htm), CDD (in NCBI http://www.ncbi.nlm.nih.gov/cdd/), PFAM (http://pfam.janelia.org/), and KOG (http://www.ncbi.nlm.nih.gov/COG/). GO functional annotations were obtained using the Blast2GO program , and all unigenes were assigned into GO functional categories of molecular function, biological process, and cellular components using WEGO software . To identify possible TFs in the transcripts, we used the BLAST tool, sequence information in the PlantTFDB database (http://planttfdb.cbi.pku.edu.cn), and annotations from the SWISS-PROT and NR databases.
Identification and functional analysis of DEUs
To identify genes displaying significant expression changes during NaCl treatment, differentially expressed transcripts were analyzed by comparing 4- and 24-h libraries with the control library, with the number of reads in all six libraries normalized to reads per kilobase of exon model per million mapped reads (RPKM) . For better accuracy, only unigenes longer than 200 bp were included in the differential unigene expression analysis. Using Bowtie v0.12.8 software and the single-end mapping method, each read per unigene was counted as 1/n by permitting alignment of one read on multiple unigenes. RPKM were calculated according to the formula RPKM = 106 × C / 10-3 × N × L, where C is the number of reads uniquely aligned to the unigene, N is the total number of reads uniquely aligned to all unigenes in the specific sample, and L is the number of nt in the unigene. The P-value corresponding to a differentially expressed unigene in two samples was determined from Audic’s algorithm , with the FDR method applied to determine the P-value threshold in multiple tests. To determine whether a unigene was differentially expressed between two samples, we therefore used threshold criteria based on the FDR statistical method and the ratio of RPKMs, i.e., FDR < 0.001 and |log2 Ratio| ≥ 1.
Gene expression patterns uncovered during the investigated two time points were analyzed and interpreted using STEM v1.3.8 . GO enrichment analyses for differentially co-expressed unigenes were performed with the differential and full unigene sets treated as prospect and background, respectively, using the hypergeometric distribution algorithm “phyper” in STEM. P-values (<0.05 or 0.01) were calculated between the prospective unigene and selected GO subcategories, followed by FDR adjustment.
Small RNA sequencing and identification of conserved miRNAs
For small RNA library construction, total RNA was separately extracted from the leaves of the same six samples used for generation of mRNA-seq libraries. After confirming RNA quality by agarose gel electrophoresis (28S:18S > 1.5) and on an Agilent 2100 bioanalyzer (RNA Integrity Number ≥ 8.0), the extracted RNA was separated by 15% denaturing polyacrylamide gel electrophoresis (PAGE) to recover the population of small RNAs (size range 17–35 nt) present. The small RNAs were ligated at their 5′ and 3′ ends to a pair of Solexa adapters, and the resulting ligation products were gel-purified by 10% denaturing PAGE and reverse transcribed. The cDNAs obtained in this fashion were sequenced on an Illumina HiSeq 2000 system (Illumina, California USA). After removal of low-quality reads and adapter sequences from the Solexa sequencing reads, unique sequences of 17–35 nt were used for further analysis. First, the unique sequences were queried against ribosomal and transfer RNAs from tRNAdb (http://trna.bioinf.uni-leipzig.de/DataOutput/Search), SILVArRNA (http://www.arb-silva.de/), and NONCODE v3.0 (http://www.noncode.org/NONCODERv3/; plant database only, single-nt mismatches allowed) databases to obtain matching rRNA, tRNA, snRNA, and snoRNA sequences. Any small RNA sequences having exact matches to these sequences were removed from the analysis. The remaining unique small RNA sequences were then BLASTn-searched against the conserved plant miRNAs in the miRBase database (miRBase 19.0; http://www.mirbase.org/) to identify conserved miRNAs. Only perfectly matched sequences were considered to be conserved miRNAs. Because we were only conducting a combined analysis of transcriptome and conserved miRNA data, we did not perform a novel miRNA prediction on the content.
Analyses of conserved miRNA differential expression, target unigene prediction, and miRNA-target pair functional regulatory relationships
A high-throughput sequencing abundance profile analysis was performed based on the sequence reads in each library. miRNA expression abundance in each library was calculated as RPM (reads per million) according to the formula RPM = mapped reads × 106 / total reads. Calculation of P-values for comparing the miRNA expression between salt-treated samples and the control sample was based on previously established methods [126, 128]. Specifically, we used the following log2 Ratio formula: log2 Ratio = log2 (miRNA reads in the salt treatment / miRNA reads in the control). Differential expression of miRNAs in the comparisons N4/N0 (where N is Nan Dan Ba Di Da Hua and E is Earlistaple 7), N24/N0, E4/E0, and E24/E0 were screened based on the chi-square test with filters of P < 0.05 and |log2 Ratio| ≥ 1.0.
To search target genes, the conserved miRNA sequences were matched against the cotton transcriptome unigene database using psRNATarget (http://plantgrn.noble.org/psRNATarget/). Sequences with less than 3 nt mismatches compared with the query miRNA sequences were selected.
Each data point represented the mean of three or more replicated treatments, with each treatment consisting of at least five plants. Physiological data statistical analyses were performed by Tukey’s method of two-way ANOVA using IBM SPSS Statistics v19.0 (SPSS Inc., Chicago, IL, USA). P-values less than 0.05 were considered as statistically significant. R software (http://www.r-project.org) was used for the calculation of Pearson correlation coefficients between miRNAs and their targets and for the construction of a heatmap of co-expressed TF unigenes and miRNA expression profiles.
Nan Dan Ba Di Da Hua
Relative water content
Relative electric conductivity
Differentially expressed unigene
False discovery rate
Reads per kilobase of gene model per million mapped reads
Reads per million
Real-time quantitative RT-PCR
Stem-loop reverse transcription-PCR
Short time-servies expression miner
Pearson correlation coefficient.
This work was supported by the grant of the “Twelfth Five-Year Plan” National Science and Technology Support Program of China [No. 2013BAD01B03].
- Wang X, Fan P, Song H, Chen X, Li X, Li Y: Comparative proteomic analysis of differentially expressed proteins in shoots of Salicornia europaea under different salinity. J Proteome Res. 2009, 8: 3331-3345. 10.1021/pr801083a.PubMedGoogle Scholar
- Rozema J, Flowers T: Crops for a salinized world. Science. 2008, 322: 1478-1480. 10.1126/science.1168572.PubMedGoogle Scholar
- Zhao KF, Fan H, Song J: Species, types, vegetation of halophytes in China and its economic potential. Utilization of halophytes and sustainable development of local agriculture. 2002, Beijing: China Meteorological Press, 1-9.Google Scholar
- Hasegawa PM, Bressan RA, Zhu J-K, Bohnert HJ: Plant cellular and molecular responses to high salinity. Annu Rev Plant Physiol Plant Mol Biol. 2000, 51: 463-499. 10.1146/annurev.arplant.51.1.463.PubMedGoogle Scholar
- Munns R: Comparative physiology of salt and water stress. Plant Cell Environ. 2002, 25: 239-250. 10.1046/j.0016-8025.2001.00808.x.PubMedGoogle Scholar
- Munns R: Genes and salt tolerance: bringing them together. New Phytol. 2005, 167: 645-663. 10.1111/j.1469-8137.2005.01487.x.PubMedGoogle Scholar
- Loescher W, Chan Z, Grumet R: Options for Developing Salt-tolerant Crops. HortScience. 2011, 46: 1085-1092.Google Scholar
- Maas E: Salt tolerance of plants. Handb Plant Sci Agric. 1984, 2: 57-75.Google Scholar
- Leidi EO, Saiz J: Is salinity tolerance related to Na accumulation in upland cotton (Gossypium hirsutum) seedlings?. Plant and Soil. 1997, 190: 67-75. 10.1023/A:1004214825946.Google Scholar
- Ahmad S, Khan N, Iqbal M, Hussain A, Hassan M: Salt tolerance of cotton (Gossypium hirsutum L.). Asian J Plant Sci. 2002, 1: 715-719.Google Scholar
- Reinhardt D, Rost T: Primary and lateral root development of dark-and light-grown cotton seedlings under salinity stress. Botanica Acta (Germany). 1995, 108: 457-465. 10.1111/j.1438-8677.1995.tb00521.x.Google Scholar
- Leidi E: Genotypic variation of cotton in response to stress by NaCl or PEG. REUR Technical Series. 1994Google Scholar
- Brugnoli E, Lauteri M: Effects of salinity on stomatal conductance, photosynthetic capacity, and carbon isotope discrimination of salt-tolerant (Gossypium hirsutum L.) and salt-sensitive (Phaseolus vulgaris L.) C3 non-halophytes. Plant Physiol. 1991, 95: 628-635. 10.1104/pp.95.2.628.PubMed CentralPubMedGoogle Scholar
- Flowers T: Improving crop salt tolerance. J Exp Biol. 2004, 55: 307-319.Google Scholar
- Chinnusamy V, Zhu J, Zhu J-K: Salt stress signaling and mechanisms of plant salt tolerance. Genetic engineering. 2006, Springer, 141-177.Google Scholar
- Roy SJ, Tucker EJ, Tester M: Genetic analysis of abiotic stress tolerance in crops. Curr Opin Plant Biol. 2011, 14: 232-239. 10.1016/j.pbi.2011.03.002.PubMedGoogle Scholar
- Wu C-A, Yang G-D, Meng Q-W, Zheng C-C: The cotton GhNHX1 gene encoding a novel putative tonoplast Na+/H+ antiporter plays an important role in salt stress. Plant Cell Physiol. 2004, 45: 600-607. 10.1093/pcp/pch071.PubMedGoogle Scholar
- Gao S-Q, Chen M, Xia L-Q, Xiu H-J, Xu Z-S, Li L-C, Zhao C-P, Cheng X-G, Ma Y-Z: A cotton (Gossypium hirsutum) DRE-binding transcription factor gene, GhDREB, confers enhanced tolerance to drought, high salt, and freezing stresses in transgenic wheat. Plant Cell Rep. 2009, 28: 301-311. 10.1007/s00299-008-0623-9.PubMedGoogle Scholar
- Jin LG, Li H, Liu JY: Molecular characterization of three ethylene responsive element binding factor genes from cotton. J Integr Plant Biol. 2010, 52: 485-495.PubMedGoogle Scholar
- Jin L-G, Liu J-Y: Molecular cloning, expression profile and promoter analysis of a novel ethylene responsive transcription factor gene GhERF4 from cotton (Gossypium hirsutum). Plant Physiol Biochem. 2008, 46: 46-53. 10.1016/j.plaphy.2007.10.004.PubMedGoogle Scholar
- Champion A, Hébrard E, Parra B, Bournaud C, Marmey P, Tranchant C, Nicole M: Molecular diversity and gene expression of cotton ERF transcription factors reveal that group IXa members are responsive to jasmonate, ethylene and Xanthomonas. Mol Plant Pathol. 2009, 10: 471-485. 10.1111/j.1364-3703.2009.00549.x.PubMedGoogle Scholar
- Meng C, Cai C, Zhang T, Guo W: Characterization of six novel NAC genes and their responses to abiotic stresses in Gossypium hirsutum L. Plant Sci. 2009, 176: 352-359. 10.1016/j.plantsci.2008.12.003.Google Scholar
- Xue T, Li X, Zhu W, Wu C, Yang G, Zheng C: Cotton metallothionein GhMT3a, a reactive oxygen species scavenger, increased tolerance against abiotic stress in transgenic tobacco and yeast. J Exp Biol. 2009, 60: 339-349.Google Scholar
- Zhang L, Xi D, Li S, Gao Z, Zhao S, Shi J, Wu C, Guo X: A cotton group C MAP kinase gene, GhMPK2, positively regulates salt and drought tolerance in tobacco. Plant Mol Biol. 2011, 77: 17-31. 10.1007/s11103-011-9788-7.PubMedGoogle Scholar
- Lu W, Chu X, Li Y, Wang C, Guo X: Cotton GhMKK1 Induces the Tolerance of Salt and Drought Stress, and Mediates Defence Responses to Pathogen Infection in Transgenic Nicotiana benthamiana. PLoS One. 2013, 8: e68503-10.1371/journal.pone.0068503.PubMed CentralPubMedGoogle Scholar
- Guo YH, Yu YP, Wang D, Wu CA, Yang GD, Huang JG, Zheng CC: GhZFP1, a novel CCCH‒type zinc finger protein from cotton, enhances salt stress tolerance and fungal disease resistance in transgenic tobacco by interacting with GZIRD21A and GZIPR5. New Phytol. 2009, 183: 62-75. 10.1111/j.1469-8137.2009.02838.x.PubMedGoogle Scholar
- Kawasaki S, Borchert C, Deyholos M, Wang H, Brazille S, Kawai K, Galbraith D, Bohnert HJ: Gene expression profiles during the initial phase of salt stress in rice. Plant Cell. 2001, 13: 889-905. 10.1105/tpc.13.4.889.PubMed CentralPubMedGoogle Scholar
- Dai Yin C, LUO YH, Min S, Da L, LIN HX: Salt-responsive genes in rice revealed by cDNA microarray analysis. Cell Res. 2005, 15: 796-810. 10.1038/sj.cr.7290349.Google Scholar
- Ueda A, Kathiresan A, Bennett J, Takabe T: Comparative transcriptome analyses of barley and rice under salt stress. Theor Appl Genet. 2006, 112: 1286-1294. 10.1007/s00122-006-0231-4.PubMedGoogle Scholar
- Rodriguez-Uribe L, Higbie SM, Stewart JM, Wilkins T, Lindemann W, Sengupta-Gopalan C, Zhang J: Identification of salt responsive genes using comparative microarray analysis in Upland cotton (Gossypium hirsutum L.). Plant Sci. 2011, 180: 461-469. 10.1016/j.plantsci.2010.10.009.PubMedGoogle Scholar
- Yao D, Zhang X, Zhao X, Liu C, Wang C, Zhang Z, Zhang C, Wei Q, Wang Q, Yan H: Transcriptome analysis reveals salt-stress-regulated biological processes and key pathways in roots of cotton (Gossypium hirsutum L). Genomics. 2011, 98: 47-55.PubMedGoogle Scholar
- Wang G, Zhu Q, Meng Q, Wu C: Transcript profiling during salt stress of young cotton (Gossypium hirsutum) seedlings via Solexa sequencing. Acta Physiol Plant. 2012, 34: 107-115. 10.1007/s11738-011-0809-6.Google Scholar
- Xu P, Liu Z, Fan X, Gao J, Zhang X, Zhang X, Shen X: De novo transcriptome sequencing and comparative analysis of differentially expressed genes in Gossypium aridum under salt stress. Gene. 2013, 525: 26-34. 10.1016/j.gene.2013.04.066.PubMedGoogle Scholar
- Carrington JC, Ambros V: Role of microRNAs in plant and animal development. Science. 2003, 301: 336-338. 10.1126/science.1085242.PubMedGoogle Scholar
- Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D: Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005, 8: 517-527. 10.1016/j.devcel.2005.01.018.PubMedGoogle Scholar
- Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O: Widespread translational inhibition by plant miRNAs and siRNAs. Science. 2008, 320: 1185-1190. 10.1126/science.1159151.PubMedGoogle Scholar
- Liu H-H, Tian X, Li Y-J, Wu C-A, Zheng C-C: Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA. 2008, 14: 836-843. 10.1261/rna.895308.PubMed CentralPubMedGoogle Scholar
- Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu J-K: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Bio. 2008, 8 (1): 25-10.1186/1471-2229-8-25.Google Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- Yin Z, Li Y, Yu J, Liu Y, Li C, Han X, Shen F: Difference in miRNA expression profiles between two cotton cultivars with distinct salt sensitivity. Mol Biol Rep. 2012, 39: 4961-4970. 10.1007/s11033-011-1292-2.PubMedGoogle Scholar
- Zhu J, Li W, Yang W, Qi L, Han S: Identification of microRNAs in Caragana intermedia by high-throughput sequencing and expression analysis of 12 microRNAs and their targets under salt stress. Plant Cell Rep. 2013, 32: 1339-1349. 10.1007/s00299-013-1446-x.PubMedGoogle Scholar
- Zhang B, Wang Q, Wang K, Pan X, Liu F, Guo T, Cobb GP, Anderson TA: Identification of cotton microRNAs and their targets. Gene. 2007, 397: 26-37. 10.1016/j.gene.2007.03.020.PubMedGoogle Scholar
- Ruan M-B, Zhao Y-T, Meng Z-H, Wang X-J, Yang W-C: Conserved miRNA analysis in Gossypium hirsutum through small RNA sequencing. Genomics. 2009, 94: 263-268. 10.1016/j.ygeno.2009.07.002.PubMedGoogle Scholar
- Pang M, Woodward AW, Agarwal V, Guan X, Ha M, Ramachandran V, Chen X, Triplett BA, Stelly DM, Chen ZJ: Genome-wide analysis reveals rapid and dynamic changes in miRNA and siRNA sequence and expression during ovule and fiber development in allotetraploid cotton (Gossypium hirsutum L). Genome Biol. 2009, 10: R122-10.1186/gb-2009-10-11-r122.PubMed CentralPubMedGoogle Scholar
- Xue W, Wang Z, Du M, Liu Y, Liu J-Y: Genome-wide analysis of small RNAs reveals eight fiber elongation-related and 257 novel microRNAs in elongating cotton fiber cells. BMC Genomics. 2013, 14: 629-10.1186/1471-2164-14-629.PubMed CentralPubMedGoogle Scholar
- Zhang X, Yao D, Wang Q, Xu W, Wei Q, Wang C, Liu C, Zhang C, Yan H, Ling Y: mRNA-seq Analysis of the Gossypium arboreum transcriptome reveals tissue selective signaling in response to water stress during seedling stage. PLoS One. 2013, 8: e54762-10.1371/journal.pone.0054762.PubMed CentralPubMedGoogle Scholar
- Brinker M, Brosché M, Vinocur B, Abo-Ogiala A, Fayyaz P, Janz D, Ottow EA, Cullmann AD, Saborowski J, Kangasjärvi J: Linking the salt transcriptome with physiological responses of a salt-resistant Populus species as a strategy to identify genes important for stress acclimation. Plant Physiol. 2010, 154: 1697-1709. 10.1104/pp.110.164152.PubMed CentralPubMedGoogle Scholar
- Zhang X, Zhen J, Li Z, Kang D, Yang Y, Kong J, Hua J: Expression profile of early responsive genes under salt stress in upland cotton (Gossypium hirsutum L.). Plant Mol Biol Rep. 2011, 29: 626-637. 10.1007/s11105-010-0269-y.Google Scholar
- Z-h D, Zheng L-l, Wang J, Gao Z, Wu S-b, Qi Z, Wang Y-c: Transcriptomic profiling of the salt-stress response in the wild recretohalophyte Reaumuria trigyna. BMC Genomics. 2013, 14: 29-10.1186/1471-2164-14-29.Google Scholar
- Munns R, Tester M: Mechanisms of salinity tolerance. Annu Rev Plant Biol. 2008, 59: 651-681. 10.1146/annurev.arplant.59.032607.092911.PubMedGoogle Scholar
- Hey SJ, Byrne E, Halford NG: The interface between metabolic and stress signalling. Ann Bot. 2010, 105: 197-203. 10.1093/aob/mcp285.PubMed CentralPubMedGoogle Scholar
- Lopez-Gomollon S, Mohorianu I, Szittya G, Moulton V, Dalmay T: Diverse correlation patterns between microRNAs and their targets during tomato fruit development indicates different modes of microRNA actions. Planta. 2012, 236: 1875-1887. 10.1007/s00425-012-1734-7.PubMedGoogle Scholar
- Flowers T, Hajibagherp M, Yeo A: Ion accumulation in the cell walls of rice plants growing under saline conditions: evidence for the Oertli hypothesis. Plant Cell Environ. 1991, 14: 319-325. 10.1111/j.1365-3040.1991.tb01507.x.Google Scholar
- Zhu J-K: Regulation of ion homeostasis under salt stress. Curr Opin Plant Biol. 2003, 6: 441-445. 10.1016/S1369-5266(03)00085-2.PubMedGoogle Scholar
- Walia H, Wilson C, Condamine P, Liu X, Ismail AM, Zeng L, Wanamaker SI, Mandal J, Xu J, Cui X: Comparative transcriptional profiling of two contrasting rice genotypes under salinity stress during the vegetative growth stage. Plant Physiol. 2005, 139: 822-835. 10.1104/pp.105.065961.PubMed CentralPubMedGoogle Scholar
- Ouyang B, Yang T, Li H, Zhang L, Zhang Y, Zhang J, Fei Z, Ye Z: Identification of early salt stress response genes in tomato root by suppression subtractive hybridization and microarray analysis. J Exp Biol. 2007, 58: 507-520.Google Scholar
- Wang Y, Yang L, Zheng Z, Grumet R, Loescher W, Zhu J-K, Yang P, Hu Y, Chan Z: Transcriptomic and Physiological Variations of Three Arabidopsis Ecotypes in Response to Salt Stress. PLoS One. 2013, 8: e69036-10.1371/journal.pone.0069036.PubMed CentralPubMedGoogle Scholar
- Guo P, Baum M, Grando S, Ceccarelli S, Bai G, Li R, Von Korff M, Varshney RK, Graner A, Valkoun J: Differentially expressed genes between drought-tolerant and drought-sensitive barley genotypes in response to drought stress during the reproductive stage. J Exp Biol. 2009, 60: 3531-3544.Google Scholar
- Sun J, Jiang H, Xu Y, Li H, Wu X, Xie Q, Li C: The CCCH-type zinc finger proteins AtSZF1 and AtSZF2 regulate salt stress responses in Arabidopsis. Plant Cell Physiol. 2007, 48: 1148-1158. 10.1093/pcp/pcm088.PubMedGoogle Scholar
- Zhou QY, Tian AG, Zou HF, Xie ZM, Lei G, Huang J, Wang CM, Wang HW, Zhang JS, Chen SY: Soybean WRKY-type transcription factor genes, GmWRKY13, GmWRKY21, and GmWRKY54, confer differential tolerance to abiotic stresses in transgenic Arabidopsis plants. Plant Biotechnol J. 2008, 6: 486-503. 10.1111/j.1467-7652.2008.00336.x.PubMedGoogle Scholar
- Zou M, Guan Y, Ren H, Zhang F, Chen F: A bZIP transcription factor, OsABI5, is involved in rice fertility and stress tolerance. Plant Mol Biol. 2008, 66: 675-683. 10.1007/s11103-008-9298-4.PubMedGoogle Scholar
- Meng X, Li F, Liu C, Zhang C, Wu Z, Chen Y: Isolation and characterization of an ERF transcription factor gene from cotton (Gossypium barbadense L.). Plant Mol Biol Rep. 2010, 28: 176-183. 10.1007/s11105-009-0136-x.Google Scholar
- Michalak P: Coexpression, coregulation, and cofunctionality of neighboring genes in eukaryotic genomes. Genomics. 2008, 91: 243-248. 10.1016/j.ygeno.2007.11.002.PubMedGoogle Scholar
- Ge Y, Li Y, Lv D-K, Bai X, Ji W, Cai H, Wang A-X, Zhu Y-M: Alkaline-stress response in Glycine soja leaf identifies specific transcription factors and ABA-mediated signaling factors. Funct Integr Genomic. 2011, 11: 369-379. 10.1007/s10142-010-0191-2.Google Scholar
- Büttner M, Singh KB: Arabidopsis thaliana ethylene-responsive element binding protein (AtEBP), an ethylene-inducible, GCC box DNA-binding protein interacts with an ocs element binding protein. Proc Natl Acad Sci U S A. 1997, 94: 5961-5966. 10.1073/pnas.94.11.5961.PubMed CentralPubMedGoogle Scholar
- Licausi F, Kosmacz M, Weits DA, Giuntoli B, Giorgi FM, Voesenek LA, Perata P, van Dongen JT: Oxygen sensing in plants is mediated by an N-end rule pathway for protein destabilization. Nature. 2011, 479: 419-422. 10.1038/nature10536.PubMedGoogle Scholar
- Meng C-M, Zhang T-Z, Guo W-Z: Molecular cloning and characterization of a novel Gossypium hirsutum L. bHLH gene in response to ABA and drought stresses. Plant Mol Biol Rep. 2009, 27: 381-387. 10.1007/s11105-009-0112-5.Google Scholar
- Qin Y-F, Li D-D, Wu Y-J, Liu Z-H, Zhang J, Zheng Y, Li X-B: Three cotton homeobox genes are preferentially expressed during early seedling development and in response to phytohormone signaling. Plant Cell Rep. 2010, 29: 1147-1156. 10.1007/s00299-010-0901-1.PubMedGoogle Scholar
- Ni Y, Wang X, Li D, Wu Y, Xu W, Li X: Novel cotton homeobox gene and its expression profiling in root development and in response to stresses and phytohormones. Acta Bioch Bioph Sin. 2008, 40: 78-84. 10.1111/j.1745-7270.2008.00371.x.Google Scholar
- Sakamoto H, Maruyama K, Sakuma Y, Meshi T, Iwabuchi M, Shinozaki K, Yamaguchi-Shinozaki K: Arabidopsis Cys2/His2-type zinc-finger proteins function as transcription repressors under drought, cold, and high-salinity stress conditions. Plant Physiol. 2004, 136: 2734-2746. 10.1104/pp.104.046599.PubMed CentralPubMedGoogle Scholar
- Huang J, Sun S-J, Xu D-Q, Yang X, Bao Y-M, Wang Z-F, Tang H-J, Zhang H: Increased tolerance of rice to cold, drought and oxidative stresses mediated by the overexpression of a gene that encodes the zinc finger protein ZFP245. Biochem Bioph Res Co. 2009, 389: 556-561. 10.1016/j.bbrc.2009.09.032.Google Scholar
- Lightfoot DJ, Malone KM, Timmis JN, Orford SJ: Evidence for alternative splicing of MADS-box transcripts in developing cotton fibre cells. Molecul Genet Genomics. 2008, 279: 75-85. 10.1007/s00438-007-0297-y.Google Scholar
- Sakuma Y, Maruyama K, Qin F, Osakabe Y, Shinozaki K, Yamaguchi-Shinozaki K: Dual function of an Arabidopsis transcription factor DREB2A in water-stress-responsive and heat-stress-responsive gene expression. Proc Natl Acad Sci U S A. 2006, 103: 18822-18827. 10.1073/pnas.0605639103.PubMed CentralPubMedGoogle Scholar
- Sakuma Y, Maruyama K, Osakabe Y, Qin F, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Functional analysis of an Arabidopsis transcription factor, DREB2A, involved in drought-responsive gene expression. Plant Cell. 2006, 18: 1292-1309. 10.1105/tpc.105.035881.PubMed CentralPubMedGoogle Scholar
- Nakashima K, Shinwari ZK, Sakuma Y, Seki M, Miura S, Shinozaki K, Yamaguchi-Shinozaki K: Organization and expression of two Arabidopsis DREB2 genes encoding DRE-binding proteins involved in dehydration-and high-salinity-responsive gene expression. Plant Mol Biol. 2000, 42: 657-665. 10.1023/A:1006321900483.PubMedGoogle Scholar
- Shan DP, Huang JG, Yang YT, Guo YH, Wu CA, Yang GD, Gao Z, Zheng CC: Cotton GhDREB1 increases plant tolerance to low temperature and is negatively regulated by gibberellic acid. New Phytol. 2007, 176: 70-81. 10.1111/j.1469-8137.2007.02160.x.PubMedGoogle Scholar
- Servet C, Benhamed M, Latrasse D, Kim W, Delarue M, Zhou D-X: Characterization of a phosphatase 2C protein as an interacting partner of the histone acetyltransferase GCN5 in Arabidopsis. Biochim Biophys Acta. 2008, 1779: 376-382. 10.1016/j.bbagrm.2008.04.007.PubMedGoogle Scholar
- Jamil M, Iqbal W, Bangash A, Rehman SU, Imran QM, Rha ES: Constitutive expression of OSC3H33, OSC3H50 and OSC3H37 genes in rice under salt stress. Pak J Bot. 2010, 42: 4003-4009.Google Scholar
- Gobert A, Park G, Amtmann A, Sanders D, Maathuis FJ: Arabidopsis thaliana cyclic nucleotide gated channel 3 forms a non-selective ion transporter involved in germination and cation transport. J Exp Biol. 2006, 57: 791-800.Google Scholar
- Mahajan S, Pandey GK, Tuteja N: Calcium-and salt-stress signaling in plants: shedding light on SOS pathway. Arch Biochem Biophys. 2008, 471: 146-158. 10.1016/j.abb.2008.01.010.PubMedGoogle Scholar
- Xiong L, Schumaker KS, Zhu J-K: Cell signaling during cold, drought, and salt stress. Plant Cell. 2002, 14 (suppl 1): S165-S183.PubMed CentralPubMedGoogle Scholar
- Weinl S, Kudla J: The CBL-CIPK Ca2+-decoding signaling network: function and perspectives. New Phytol. 2009, 184: 517-528. 10.1111/j.1469-8137.2009.02938.x.PubMedGoogle Scholar
- H-i C, Park H-J, Park JH, Kim S, Im M-Y, Seo H-H, Kim Y-W, Hwang I, Kim SY: Arabidopsis calcium-dependent protein kinase AtCPK32 interacts with ABF4, a transcriptional regulator of abscisic acid-responsive gene expression, and modulates its activity. Plant Physiol. 2005, 139: 1750-1761. 10.1104/pp.105.069757.Google Scholar
- Zhen J, Zhang X, Wang Y, Zhai Z, Li Z, Hua J: Isolation and sequence analysis of salt responsive gene GhCPK5 from Gossypium hirsutum L. J Chin Agric Univ. 2011, 16: 8-14.Google Scholar
- Saijo Y, Hata S, Kyozuka J, Shimamoto K, Izui K: Over expression of a single Ca2+-dependent protein kinase confers both cold and salt/drought tolerance on rice plants. Plant J. 2000, 23: 319-327. 10.1046/j.1365-313x.2000.00787.x.PubMedGoogle Scholar
- Sheen J: Ca2+-dependent protein kinases and stress signal transduction in plants. Science. 1996, 274: 1900-1902. 10.1126/science.274.5294.1900.PubMedGoogle Scholar
- Teige M, Scheikl E, Eulgem T, Dóczi R, Ichimura K, Shinozaki K, Dangl JL, Hirt H: The MKK2 Pathway Mediates Cold and Salt Stress Signaling in Arabidopsis. Mol Cell. 2004, 15: 141-152. 10.1016/j.molcel.2004.06.023.PubMedGoogle Scholar
- Luo J, Zhao L-L, Gong S-Y, Sun X, Li P, Qin L-X, Zhou Y, Xu W-L, Li X-B: A cotton mitogen-activated protein kinase (GhMPK6) is involved in ABA-induced CAT1 expression and H2O2 production. J Genet Genomics. 2011, 38: 557-565. 10.1016/j.jgg.2011.10.003.PubMedGoogle Scholar
- Shi J, An H-L, Zhang L, Gao Z, Guo X-Q: GhMPK7, a novel multiple stress-responsive cotton group C MAPK gene, has a role in broad spectrum disease resistance and plant development. Plant Mol Biol. 2010, 74: 1-17. 10.1007/s11103-010-9661-0.PubMedGoogle Scholar
- Peleg Z, Blumwald E: Hormone balance and abiotic stress tolerance in crop plants. Curr Opin Plant Biol. 2011, 14: 290-295. 10.1016/j.pbi.2011.02.001.PubMedGoogle Scholar
- Seo M, Koshiba T: Complex regulation of ABA biosynthesis in plants. Trends Plant Sci. 2002, 7: 41-48.PubMedGoogle Scholar
- Chen C-N, Chu C-C, Zentella R, Pan S-M, Ho T-HD: AtHVA22 gene family in Arabidopsis: phylogenetic relationship, ABA and stress regulation, and tissue-specific expression. Plant Mol Biol. 2002, 49: 631-642. 10.1023/A:1015593715144.Google Scholar
- Koops P, Pelser S, Ignatz M, Klose C, Marrocco-Selden K, Kretsch T: EDL3 is an F-box protein involved in the regulation of abscisic acid signalling in Arabidopsis thaliana. J Exp Biol. 2011, 62: 5547-5560.Google Scholar
- Hirayama T, Shinozaki K: Perception and transduction of abscisic acid signals: keys to the function of the versatile plant hormone ABA. Trends Plant Sci. 2007, 12: 343-351. 10.1016/j.tplants.2007.06.013.PubMedGoogle Scholar
- Yoo S-D, Cho Y-H, Tena G, Xiong Y, Sheen J: Dual control of nuclear EIN3 by bifurcate MAPK cascades in C2H4 signalling. Nature. 2008, 451: 789-795. 10.1038/nature06543.PubMed CentralPubMedGoogle Scholar
- Padmalatha KV, Dhandapani G, Kanakachari M, Kumar S, Dass A, Patil DP, Rajamani V, Kumar K, Pathak R, Rawat B: Genome-wide transcriptomic analysis of cotton under drought stress reveal significant down-regulation of genes and pathways involved in fibre elongation and up-regulation of defense responsive genes. Plant Mol Biol. 2012, 78: 223-246. 10.1007/s11103-011-9857-y.PubMedGoogle Scholar
- Vanlerberghe GC, Cvetkovska M, Wang J: Is the maintenance of homeostatic mitochondrial signaling during stress a physiological role for alternative oxidase?. Physiol Plantarum. 2009, 137: 392-406. 10.1111/j.1399-3054.2009.01254.x.Google Scholar
- Kirch H-H, Schlingensiepen S, Kotchoni S, Sunkar R, Bartels D: Detailed expression analysis of selected genes of the aldehyde dehydrogenase (ALDH) gene superfamily in Arabidopsis thaliana. Plant Mol Biol. 2005, 57: 315-332. 10.1007/s11103-004-7796-6.PubMedGoogle Scholar
- Kotchoni SO, Kuhns C, Ditzer A, KIRCH HH, Bartels D: Over expression of different aldehyde dehydrogenase genes in Arabidopsis thaliana confers tolerance to abiotic stress and protects plants against lipid peroxidation and oxidative stress. Plant Cell Environ. 2006, 29: 1033-1048. 10.1111/j.1365-3040.2005.01458.x.PubMedGoogle Scholar
- Sunkar R: MicroRNAs with macro-effects on plant stress responses. Seminars in cell & developmental biology. 2010, Academic Press, 21: 805-811Google Scholar
- Li B, Qin Y, Duan H, Yin W, Xia X: Genome-wide characterization of new and drought stress responsive microRNAs in Populus euphratica. J Exp Biol. 2011, 62: 3765-3779.Google Scholar
- Zhang Q, Zhao C, Li M, Sun W, Liu Y, Xia H, Sun M, Li A, Li C, Zhao S: Genome-wide identification of Thellungiella salsuginea microRNAs with putative roles in the salt stress response. BMC Plant Biol. 2013, 13: 180-10.1186/1471-2229-13-180.PubMed CentralPubMedGoogle Scholar
- Glazov EA, Cottee PA, Barris WC, Moore RJ, Dalrymple BP, Tizard ML: A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach. Genome Res. 2008, 18: 957-964. 10.1101/gr.074740.107.PubMed CentralPubMedGoogle Scholar
- Sunkar R, Chinnusamy V, Zhu J, Zhu J-K: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12: 301-309. 10.1016/j.tplants.2007.05.001.PubMedGoogle Scholar
- Lelandais-Brière C, Sorin C, Declerck M, Benslimane A, Crespi M, Hartmann C: Small RNA diversity in plants and its impact in development. Curr Genomics. 2010, 11: 14-23. 10.2174/138920210790217918.PubMed CentralPubMedGoogle Scholar
- Wu G, Park MY, Conway SR, Wang J-W, Weigel D, Poethig RS: The Sequential Action of miR156 and miR172 Regulates Developmental Timing in Arabidopsis. Cell. 2009, 138: 750-759. 10.1016/j.cell.2009.06.031.PubMed CentralPubMedGoogle Scholar
- Olsen AN, Ernst HA, Leggio LL, Skriver K: NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 2005, 10: 79-87.PubMedGoogle Scholar
- Xie Q, Frugis G, Colgan D, Chua N-H: Arabidopsis NAC1 transduces auxin signal downstream of TIR1 to promote lateral root development. Gene Dev. 2000, 14: 3024-3036. 10.1101/gad.852200.PubMed CentralPubMedGoogle Scholar
- Duval M, Hsieh T-F, Kim SY, Thomas TL: Molecular characterization of AtNAM: A member of the Arabidopsis NAC domain superfamily. Plant Mol Biol. 2002, 50: 237-248. 10.1023/A:1016028530943.PubMedGoogle Scholar
- Nikovics K, Blein T, Peaucelle A, Ishida T, Morin H, Aida M, Laufs P: The balance between the MIR164A and CUC2 genes controls leaf margin serration in Arabidopsis. Plant Cell. 2006, 18: 2929-2945. 10.1105/tpc.106.045617.PubMed CentralPubMedGoogle Scholar
- Meng Y, Ma X, Chen D, Wu P, Chen M: MicroRNA-mediated signaling involved in plant root development. Biochem Bioph Res Co. 2010, 393: 345-349. 10.1016/j.bbrc.2010.01.129.Google Scholar
- Toledo-Ortiz G, Huq E, Quail PH: The Arabidopsis basic/helix-loop-helix transcription factor family. Plant Cell. 2003, 15: 1749-1770. 10.1105/tpc.013839.PubMed CentralPubMedGoogle Scholar
- X-M DU, J-L SUN, Z-L ZHOU, Y-H GU, Z-E PAN, S-P HE, B-Y PANG, L-R WANG: Current situation and the future in collection, preservation, evaluation and utilization of cotton germplasm in China. J Plant Genet Resour. 2012, 13 (2): 163-168.Google Scholar
- Hoagland DR, Arnon DI: The water-culture method for growing plants without soil. Calif Agric Exp Stn Circ. 1950, 347 (2nd edit):
- Smart RE, Bingham GE: Rapid estimates of relative water content. Plant Physiol. 1974, 53: 258-260. 10.1104/pp.53.2.258.PubMed CentralPubMedGoogle Scholar
- Cao W-H, Liu J, He X-J, Mu R-L, Zhou H-L, Chen S-Y, Zhang J-S: Modulation of ethylene responses affects plant salt-stress responses. Plant Physiol. 2007, 143: 707-719.PubMed CentralPubMedGoogle Scholar
- Porra R, Thompson W, Kriedemann P: Determination of accurate extinction coefficients and simultaneous equations for assaying chlorophylls a and b extracted with four different solvents: verification of the concentration of chlorophyll standards by atomic absorption spectroscopy. BBA-Bioenergetics. 1989, 975: 384-394. 10.1016/S0005-2728(89)80347-0.Google Scholar
- Livak KJ, Schmittgen TD: Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2- ΔΔCT Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedGoogle Scholar
- Kosmas SA, Argyrokastritis A, Loukas MG, Eliopoulos E, Tsakas S, Kaltsikes PJ: Isolation and characterization of drought-related trehalose 6-phosphate-synthase gene from cultivated cotton (Gossypium hirsutum L). Planta. 2006, 223: 329-339. 10.1007/s00425-005-0071-5.PubMedGoogle Scholar
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-e179. 10.1093/nar/gni178.PubMed CentralPubMedGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP: Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007, 3: 12-10.1186/1746-4811-3-12.PubMed CentralPubMedGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.PubMed CentralPubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (suppl 2): W293-W297.PubMed CentralPubMedGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedGoogle Scholar
- Audic S, Claverie J-M: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- Ernst J, Bar-Joseph Z: STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006, 7: 191-10.1186/1471-2105-7-191.PubMed CentralPubMedGoogle Scholar
- Man MZ, Wang X, Wang Y: POWER_SAGE: comparing statistical tests for SAGE experiments. Bioinformatics. 2000, 16: 953-959. 10.1093/bioinformatics/16.11.953.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.