- Research article
- Open Access
Genome-wide identification of Thellungiella salsuginea microRNAs with putative roles in the salt stress response
BMC Plant Biology volume 13, Article number: 180 (2013)
MicroRNAs are key regulators of plant growth and development with important roles in environmental adaptation. The microRNAs from the halophyte species Thellungiella salsuginea (salt cress), which exhibits extreme salt stress tolerance, remain to be investigated. The sequenced genome of T. salsuginea and the availability of high-throughput sequencing technology enabled us to discover the conserved and novel miRNAs in this plant species. It is interesting to identify the microRNAs from T. salsuginea genome wide and study their roles in salt stress response.
In this study, two T. salsuginea small RNA libraries were constructed and sequenced using Solexa technology. We identified 109 miRNAs that had previously been reported in other plant species. A total of 137 novel miRNA candidates were identified, among which the miR* sequence of 26 miRNAs was detected. In addition, 143 and 425 target mRNAs were predicted for the previously identified and Thellungiella-specific miRNAs, respectively. A quarter of these putative targets encode transcription factors. Furthermore, numerous signaling factor encoding genes, defense-related genes, and transporter encoding genes were amongst the identified targets, some of which were shown to be important for salt tolerance. Cleavage sites of seven target genes were validated by 5’ RACE, and some of the miRNAs were confirmed by qRT-PCR analysis. The expression levels of 26 known miRNAs in the roots and leaves of plants subjected to NaCl treatment were determined by Affymetrix microarray analysis. The expression of most tested miRNA families was up- or down-regulated upon NaCl treatment. Differential response patterns between the leaves and roots were observed for these miRNAs.
Our results indicated that diverse set of miRNAs of T. salsuginea were responsive to salt stress and could play an important role in the salt stress response.
Thellungiella salsuginea, previously named Thellungiella halophila, separated from Arabidopsis around 15–30 million years ago and occupies an intermediate position between the Brassica and Arabidopsis clades [1–3]. Whereas this crucifer exhibits tolerance to high salt, drought, and cold stress, its close relative Arabidopsis thaliana is sensitive to these stresses. Comparative physiology and molecular biology studies demonstrated that different strategies are employed by these two species to adapt to stress conditions. T. salsuginea accumulates less sodium under salt stress than does Arabidopsis by using different ion uptake and translocation strategies. In contrast to Arabidopsis, the concentrations of sodium and potassium are not negatively correlated in T. salsuginea. Under salt stress, T. salsuginea plants undergo relatively high rates of transpiration and thus maintain water uptake and ion transport, which are crucial for plant survival . The accumulation of high concentrations of compatible osmolytes in the cytosol, including sugars, polyols, and amino acids, is considered to be an effective strategy for acquiring tolerance to salt stress . T. salsuginea contains higher levels of proline than does Arabidopsis under unstressed conditions, and accumulates 19-fold more proline in the shoot than does Arabidopsis under salt stress . Microarray, qRT-PCR, and metabolite profiling analyses defined distinct stress-responsive genes and pathways that were active during the adaptation to stress conditions in T. salsuginea and Arabidopsis [3, 4, 7]. The expression patterns of many transcription factors diverged between these two species . Given that the target genes of microRNAs (miRNAs) are predominantly transcription factors , we thought that it would be interesting to investigate the role of miRNAs in the regulation of salt tolerance in T. salsuginea.
Gene regulation by small RNAs represents a recently discovered layer of gene expression regulation in eukaryotes, and miRNAs are the best characterized of the small RNAs in plants . miRNAs are ~21 nt non-coding RNAs that negatively regulate gene expression by mediating mRNA cleavage or translational inhibition of its target mRNA. Mature miRNA is preferentially incorporated into the RNA-Induced Silencing Complex (RISC) to direct the post-transcriptional repression of target mRNAs [9–11]. Previous studies demonstrated that miRNA regulates various aspects of plant growth and development. For example, miR165/miR166 is a critical regulator of adaxial/abaxial patterning in developing leaves , miR172 governs floral organ development [13, 14], and miR159/miR319 regulates the level of TEOSINTE BRANCHED/CYCLOIDEA/PCF (TCP)-family transcripts, which control leaf development . Increasing data from Arabidopsis, Oryza sativa (rice), and Populus indicate that miRNAs are stress-regulated and may play critical roles in the plant’s adaptation to environmental stress [11, 16–19]. However, the roles of miRNAs in the extreme stress tolerance exhibited by the halophyte T. salsuginea have not been investigated. The release of the entire genome sequence of T. salsuginea[20, 21] has greatly stimulated functional and comparative genomics studies of this species, and has enabled us to perform whole genome identification of T. salsuginea miRNAs, to predict their target genes, and to investigate how the salt-stress response differs from that in Arabidopsis at the level of miRNA regulation.
Although T. salsuginea can withstand exposure to high concentrations of NaCl [7, 22] or sudden osmotic shock (i.e., transfer from 500 mM to 0 mM NaCl), the germination of T. salsuginea seeds is hypersensitive to NaCl. Shoot succulence of T. salsuginea was affected and growth was inhibited by 200 mM NaCl [4, 23]. A previous study showed that 150 and 250 mM NaCl represent a minor irritation and moderate stress, respectively, for T. salsuginea. To identify T. salsuginea miRNAs and investigate their possible roles in salt tolerance, we constructed two small RNA libraries. One library, which was used as the control, was derived from the tissues of 40-d-old plants grown under non-stress conditions, and the other library was constructed using 40-d-old plants treated with 200 mM NaCl for 24 h. We identified 109 known miRNAs, 26 novel miRNAs, and 111 plausible novel miRNA candidates from T. salsuginea. Target gene analysis and an investigation of the patterns of miRNA accumulation revealed important information about miRNAs in T. salsuginea that may elucidate salt tolerance mechanisms in this halophyte plant and others.
Comparison of sRNA populations in CL and TL
To identify miRNAs in T. salsuginea, we constructed two small RNA libraries. The control library (CL) was constructed using plants grown in Hoagland solution without NaCl, and the treatment library (TL) used plants treated with 200 mM NaCl for 24 h. Previous studies indicated that the shoot succulence of T. salsuginea was affected and the growth was inhibited by 200 mM NaCl [4, 23]. Total RNA was isolated from the roots and leaves 40 d after germination. Solexa technology (Illumina, Shenzheng, China) was used to sequence the small RNA enriched libraries, and 12,932,436 and 13,608,695 raw reads were generated from the CL and TL, respectively (Additional file 1: Table S1). The length of the small RNA sequences ranged from 10–30 nt (Figure 1). In both libraries, 24-nt sequences represented the largest group, followed by a group of 21 nt. Sequences of 20 nt and 23 nt represented the third largest group. The abundance of most small RNAs from 10–30 nt in the TL was higher than that in the CL. Only 20, 21, 28, 29, and 30-nt small RNAs were less abundant in the TL than in the CL (Figure 1).
After removal of the adaptor sequences, polyA sequences, and sequences of less than 18 nt, 12,010,658 and 12,330,771 clean reads were obtained from the CL and TL, respectively, and used for further analysis and miRNA discovery (Additional file 1: Table S1). The two libraries shared 831,074 unique small RNAs, which accounted for 14.16% of the total unique sequences. We detected 2,593,871 (44.18%) and 2,446,118 (41.66%) unique sequences in the CL and TL, respectively. Together, the two libraries shared 76.04% (18,509,850 out of 24,341,429) of the total small RNA sequences. We identified 3,060,900 (12.57%) and 2770679 (11.38%) small RNAs that were specifically present in the CL and TL, respectively (Additional file 2: Table S2). About 60% of these sequences could be successfully mapped to a particular region of the T. salsuginea genome (Additional file 3: Table S3).
Based on a comparison with the reference genome and relevant bioinformatics analysis, these small RNAs fell into the following categories: rRNA, tRNA, snRNA, exon_antisense, exon_sense, intron_antisense, intron_sense, unannotated, and miRNA candidate sequences. The miRNA candidate sequences represented only about 1.4% of the unique small RNAs (Additional file 4: Table S4). Most sequences (85% of the total unique sequences in both libraries) were unannotated (Additional file 4: Table S4). Although the total number of candidate miRNA sequences from the TL (1502495) was less than that from the CL (1904244), more clean reads were sequenced from the TL. Similarly, fewer unique candidate miRNA sequences were identified in the TL (47015) than in the CL (47793) (Additional file 3: Table S3 and Additional file 4: Table S4). These results indicate that miRNA biogenesis is down regulated in T. salsuginea when the plants are treated with 200 mM NaCl.
Identification of known miRNAs and their variants
To identify miRNAs in T. salsuginea that have previously been reported in other plant species, small RNAs were compared with published miRNAs by conducting a Blastn search against miRbase 19.0 (April, 2013), in which 5141 plant miRNAs were registered. We identified 109 known miRNAs belonging to 23 conserved miRNA families and 9 known miRNA families (Table 1, Additional file 5: Table S5). The pre-miRNAs of 70 of these known miRNAs were discovered, and the secondary structures of these miRNAs are listed in Additional file 6: Figure S1. Previous studies indicated that 24 conserved miRNA families are present in both moncot and dicot model species [11, 24]. miR397 is the only conserved miRNA family that was not detected in these two T. salsuginea libraries. We only considered sequences that perfectly matched previously identified mature miRNA sequences as known miRNAs. Sequences with one or two mismatches or gaps relative to the known miRNAs and that were represented in a minimum of 10 reads in at least one library were considered as miRNA variants (Additional file 7: Table S6). We identified a total of 1022 miRNA variants, which we named tsa-miRv1 to tsa-miRv1022. Many variants of some miRNA families were identified, highlighting the complexity of the miRNA variants. For example, we identified 395 variants of miR156.
We identified the Thellungiella homolog of most miRNAs that were conserved in Arabidopsis, Populus trichocarpa, Vitis vinifera, Oryza sativa, and Brassica napus, including miR156a, miR157a, d, miR164a, miR166a, and miR167a (Additional file 5: Table S5). A few miRNAs, such as miR824, miR827, and miR845, were only present in one to three of the five species. miR5139 and miR4995 were not present in these five species, while miR894 was only reported in Populus trichocarpa (Additional file 5: Table S5). In agreement with previous bioinformatics studies, these results indicate the varied extent to which different miRNA families are conserved . In most cases, multiple members from one conserved miRNA family were detected; for example, we discovered 20 different members of miR156/157 and 13 different members of miR165/miR166. Interestingly, some miRNAs, including miR158, miR319, and miR398, were scarce, while their variants were abundant. For example, miR158 was sequenced about 20 times, while its variants (108 members) were sequenced more than 130,000 times in these two libraries. tsa-miRv396 and tsa-miRv397, the two variants of miR158, were the most abundant variants in T. salsuginea, being detected 54,501 and 24,028 times in CL and 31,356 and 16,098 times in the TL, respectively (Additional file 7: Table S6). These results demonstrate that tsa-miRv396 and tsa-miRv397 are the dominant mature sequences of miR158 in T. salsuginea.
The abundance of identified conserved miRNAs was significantly different, ranging from only a few to more than half a million counts. The most abundant expressed miRNA in T. salsuginea was miR156a, which was detected 665,846 and 394,884 times in the CL and TL, respectively (Figure 2, Additional file 5: Table S5). miR156q, miR157a, miR157d, miR166a, miR167a, miR168a, and miR172a were highly expressed, with more than 10,000 reads sequenced in each library. About 10 other miRNA members showed a moderate expression level, with 1,000 to 9,000 reads in each library. The abundance of most other miRNAs (> 80%) was low, ranging from a few to a hundred counts in each library (Figure 2, Additional file 5: Table S5.
A comparison of the 109 known miRNAs in the two libraries revealed that the abundance of 26 miRNAs decreased under salt treatment (ratio ≤ 0.67), including three miRNAs detected only in the CL. Only 17 miRNAs exhibited more than a 1.5-fold upregulation upon 200 mM NaCl treatment. However, the abundance of some miRNAs was extremely low in both the CL and TL, suggesting that the differential expression of these miRNAs may not be linked to the salt response. Overall, the abundance of miRNAs that exhibited increased expression in the TL was low. However, miR161, miR165a, miR166a, and miR168a were strongly expressed in both libraries and were up-regulated in the TL (Figure 2, Additional file 5: Table S5). The abundance of 380 variants was down-regulated (ratio ≤ 0.67), including that of 7 variants detected only in the CL. The NaCl treatment up-regulated the expression of 220 variants by more than 1.5-fold, including 4 variants detected only in the TL. The expression patterns of most variants were in accordance with those of known miRNAs (Additional file 7: Table S6). T. salsuginea miR894 and its variants tsa-miRv977 and tsa-miRv978 were highly expressed in both libraries, particularly after the NaCl treatment. T. salsuginea miR408 was strongly induced by the NaCl treatment, while Arabidopsis miR408 was induced by cold and mannitol, but not by salt . The differential response of these miRNAs to salt stress suggests that they may have a key role in salt tolerance.
Identification of novel miRNAs
We used previously described criteria to predict specific miRNAs of T. salsuginea[26, 27]. The precursor of each novel miRNA was identified and could form a proper secondary hairpin structure with a maximal free energy of -18 kcal/mol. We identified 26 novel miRNA families and 111 plausible candidate miRNAs (PCs) from the two libraries (Table 2, Additional file 8: Figure S2 and Additional file 7: Table S7). The free energy of the identified novel miRNAs and PCs ranged from -18.3 kcal mol-1 to -173.5.kcal mol-1. Blastn searches against all nucleotide sequences in NCBI databases were performed using these miRNA sequences as queries. No homolog was found in any other plant species, suggesting that these newly identified candidate miRNAs are all Thellungiella-specific. The abundance of most miR* sequences identified in this study was less than 10 counts, except for tsa-miRn3, tsa-miRn8, tsa-miRn11, and tsa-miRn20 (Table 2). Previous studies indicated that both the 5p and 3p sequences of some miRNAs could be recruited into the silencing complex and function as miRNAs . The abundance of the majority of novel miRNAs and plausible candidates was relatively low compared with that of conserved miRNAs. However, some novel miRNAs were abundant; for example, miRn13 (3777 counts in the CL, 2846 counts in the TL), tsa-miRn20-3p (639 counts in the CL, 523 counts in the TL), miRn7 (1927 counts in the CL), PC003 (397 counts in the CL, 652 counts in the TL), and PC070 (646 counts in the CL, 3201 counts in the TL) (Table 2, Additional file 7: Table S7).
Many miRNAs were detected in only one library; for instance, tsa-miRn6, tsa-miRn21, and tsa-miRn22 were only present in the TL, and tsa-miRn1, tsa-miRn7, tsa-miRn10, and tsa-miRn14 only in the CL (Table 2). Thirty-one novel candidate miRNAs were only detected in the CL and thirty were only found in the TL (Additional file 7: Table S7). Due to the low abundance of most novel miRNAs, the specificity of these miRNAs in one library may not truly reflect a specific expression pattern under a particular growth condition. However, a few putative miRNAs were expressed only in one library at relatively high abundance (Table 2, Additional file 7: Table S7). These miRNAs may function as key regulators of the salt stress response in T. salsuginea. Previous functional studies of the predicted target genes of tsa-miRn7 and PC069 support this hypothesis [29, 30] (Additional file 7: Table S8). Considering the environmental, developmental, and tissue-specific regulation of miRNA expression, the miRNAs identified in our study probably represent only part of the T. salsuginea miRNA population. It remains to be determined how much overlap exists between the miRNAs identified in this study and those that are responsive to cold or drought stress.
Expression analysis of miRNAs by microarray and qRT-PCR analysis
The expression of some conserved miRNAs in leaves and roots was examined in the presence or absence of NaCl treatment using an Affymetrix miRNA chip. About 400 probes on this chip corresponded with 45 plant miRNAs. These probes (miRNAs) shared identical sequences with the Thellungella miRNAs. Twenty-six of these miRNAs were represented by more than three probes on the chip and were used in the expression analysis. The remaining 19 miRNAs, which had less than three probes on the chip, were not used in the expression analysis. In this study, leaf and root tissues were sampled separately after 3, 9, and 24 h of NaCl treatment. The expression of some conserved miRNAs was up- or down-regulated (1.5-fold) by the NaCl treatment (Figure 3). Differential expression patterns in the roots and leaves were observed for some miRNAs (Figure 3). For example, the expression of tsa-miR164a and tsa-miR171c decreased in roots, but increased in leaves after salt treatment.
To verify the expression of the newly identified Thellungella miRNAs, a stem-loop qRT-PCR assay was carried out. The expression of nine conserved miRNAs and two candidate novel miRNAs (PC065 and PC073) was analyzed. The results were consistent with the Solexa sequencing data (Figure 4).
Target gene prediction of known and novel miRNAs
We predicted more than 143 potential target genes for all of the 32 known T. salsuginea miRNA families. These target genes included 42 (29%) transcription factors (Figure 5A, Additional file 7: Table S10). Interestingly, three stress-responsive transcription factors, one ABA-responsive factor, and two ethylene-responsive factors were the predicted targets of Thellungiella miRNAs. Genes involved in a broad range of developmental and physiological processes were also predicted to be miRNA targets in this study. Several genes encoding ion transporters or ion channels were identified as target genes of Thellungiella miRNAs, including an inward rectifying potassium channel, high-affinity nickel-transport family protein, Ca2+ transporting ATPase, plasma membrane sulphate transporter, boron transporter, calcium-binding protein annexin 7, Fe(II) transport protein IRT1, putative cadmium/zinc-transporting ATPase 3, and sulfite exporter TauE. In addition, many target genes identified in this study were directly or indirectly involved in the stress response or stress tolerance, including two salt-inducible proteins and two E3 ubiquitin ligases .
Following the published criteria for target gene prediction , 425 target genes were predicted from 26 novel and 111 PC Thellungilla miRNAs. These predicted target genes had a wide range of biological functions, and many encode transcription factors. One-hundred-and-nine targets encoded transcription factors or splicing factors thought to affect gene expression directly. Many of these transcription factors (e.g., auxin response factor, no apical meristem (NAM) protein, NAC transcription factor, and WRKY transcription factor) are key regulators in the salt, drought, or cold stress response [32–35]. Interestingly, ion and macromolecular transporters, which may play key roles in ion homeostasis or the accumulation of compatible molecules under stress conditions, were significantly enriched in targets of these novel miRNAs. Many miRNAs were predicted to target a variety of transporters and ion channels, such as an anion transporter, ion channel CLC-3, ABC transporter, and efflux pump (MatE). Furthermore, a number of other proteins known to be key players in the stress response were predicted to be targets of the Thellungiella novel miRNAs or miRNA candidates (Figure 5, Additional file 7: Table S8). The discovery that genes directly or indirectly involved in stress response pathways are potential targets of the newly identified Thellungiella miRNAs suggests that miRNA regulation is an important layer of stress tolerance in this plant.
Validating the cleavage of target mRNA by 5’ RACE
To test whether the identified miRNAs could mediate the cleavage of the predicted target mRNAs, a 5’ RACE experiment was performed. The targets of the conserved miRNAs, miR156/157, miR160, miR160-3p, miR164, miR171a, and miR171c, were cleaved at a site corresponding to the 10th miRNA nucleotide from the 5’ end (Figure 6). Mismatches between the miRNA molecules and their binding sites on the target were rare. miR171a was 100% complementary to the binding site of the target. One or two mismatches were indentified in the other miRNA/target pairs. No mismatches were found at the 10th and 11th nucleotides in any of these miRNA/target pairs.
Due to high levels of sequence similarity, the different members within an miRNA family may target the same gene. For example, several miR156/157 members could target the SBP transcription factor at the same site on the target mRNA. We analyzed an SBP gene with high sequence similarity to Arabidopsis SPL10 by 5’ RACE. The target was cleaved at the position between the 9th and 10th nucleotides by Thellungiella miR156a, miR156e, and miR157b. However, this position was corresponding to the 10th and 11th nucleotide of Thellungiella miR156n, miR156o, and miR157a (Figure 6, Table 1, Additional file 5: Table S5). In Phaseolus vulgaris, the cleavage site of the SBP-like related gene was located between the 10th and 11th nucleotide of miR156a . An SBP transcription factor, a PPR repeat protein, and a MADS box transcription factor were predicted to be targets of tsa-miR156/157, tsa-miR158, and tsa-miR824, respectively. However, the 5’ RACE experiment only detected cleavage of these genes upstream or downstream of the miRNA/target-binding region, but not at the predicted sites (data not shown). Thus, the target cleavage positions were predominantly at around the 10th nucleotide from the 5’ end of the miRNA, but some might mediate cleavage at other locations.
Accumulating evidence suggests that miRNAs play essential roles in both abiotic and biotic stress responses. However, a systematic identification of miRNAs from the stress-tolerant model plant T. salsuginea had hitherto not been accomplished. Taking advantage of the sequenced genome of T. salsuginea and high-throughput sequencing technology, the present study reported the first large-scale identification of T. salsuginea miRNAs. A total of 109 known miRNAs belonging to 23 conserved miRNA families and 9 non-conserved known miRNA families, as well as 26 novel miRNAs and 111 plausible candidates were discovered. Some of the miRNAs were differentially expressed in the presence or absence of 200 mM NaCl. Interestingly, three novel and 30 plausible candidates of novel miRNAs were only detected in NaCl-treated samples.
Most of the conserved miRNAs exhibited decreased abundance upon treatment with 200 mM NaCl. Only 17 miRNAs showed more than a 1.5-fold increase in abundance after NaCl treatment. The overall abundance of these differentially expressed miRNAs was low in both libraries, with the exception of miR166f, miR168a, tsa-miR408, and tsa-miR408-5p. The abundance of these miRNAs was high in both libraries. In agreement with our results, previous studies consistently demonstrated that miR166, miR168a, and miR408 were all responsive to stress conditions, including salt and drought stress [19, 37–40].
The expression of miR389 was inhibited and the expression of miR397 and miR402 was induced by salt, cold, dehydration, and ABA in Arabidopsis , while these three miRNAs and their variants were not detected in T. salsuginea. Phaseolus miR393 was induced by NaCl, dehydration, and ABA, but Thellungiella miR393a was only detected once in each library, and miR393b-3p showed decreased expression upon salt treatment. In agreement with results from Phaseolus, Thellungiella miR159b exhibited increased expression levels under salt stress (Additional file 5: Table S5). MiR167c was strongly up-regulated by 24 h of NaCl treatment in Arabidopsis, while Thellungiella miR167a-h showed similar or decreased expression levels under NaCl stress conditions  (Additional file 5: Table S5). miR398 expression was induced by salt treatment in Populus, while it was repressed in Arabidopsis . Very low levels of miR398 expression were detected in both control and NaCl-treated Thellungiella libraries (1 vs. 5); however, its variant, tsa-miRv878, was sequenced with relatively high abundance (98 vs. 237), indicating that NaCl induced miR398 in Thellungiella. The conservation and divergence of miRNA expression under different stress conditions illustrates the complexity of miRNA regulation.
The target genes of all identified known miRNA families, 26 novel miRNAs, and 109 plausible candidates were predicted. The identified target genes are involved in multiple biological processes. Many transport proteins or ion channels were predicted as targets of conserved and Thellungiella-specific miRNAs. Some of these proteins have been well studied and their roles in salt tolerance or the stress response have been established. Maintenance of Na+, K+ homeostasis is required for plant survival in a high salinity environment. HKT1, a K+ transporter, contributed significantly to Na+ uptake in Arabidopsis and rice and plays key roles in salt tolerance [42–44]. Substantial transcriptional regulation of the high affinity K+ transporter was observed under salt stress . The K+ transporter gene was identified as a novel miRNA target in Thellungiella. An inward rectifying potassium channel, a member of the AKT protein family, exhibits non-specific Na+ uptake activity under high concentrations of Na+. In the halophytic plant M. crystallinum, the expression of the AKT1 homolog decreased in response to salt stress . We found that the inward rectifying potassium channel was the target of Thellungiella miRNA.
Besides tolerance to salinity, Thellungiella exhibits high tolerance to dehydration and cold stresses [23, 48–51]. The novel miRNA candidates identified in our study may represent small RNAs that are expressed in a relatively high level under non-stress or salt stress conditions. A previous study revealed that gene expression in response to cold, drought, and salinity showed little overlap . Therefore, we do not expect much overlap in Thellungiella miRNAs involved in the plant’s response to different stresses. It would be interesting to determine which miRNAs from Thellungiella are activated upon exposure to cold or drought stresses.
About thirteen million high-quality small RNAs were sequenced from each library. Totally, we identified 109 known miRNAs and 137 novel miRNA candidates from T. salsuginea. In addition, 568 target mRNAs were predicted. Differentially expressed miRNAs between NaCl treated and the control samples were analyzed using Solexa sequence data and part of them were validated by microarray and together with qRT-PCR. This study lay the basis for functional studies of miRNA in T. salsuginea and enables us to access the salt tolerance mechanism of this plant species from the layer of small RNA regulation.
Thellungiella salsuginea (Shandong ecotype) seeds were surface-sterilized, sown on 1/2 MS agar plates, stratified at 4°C for one week, and then moved to a growth chamber for one week at 26°C/22°C (day/night) and with 16 h of illumination per day. Seedlings were transferred to a planting box containing Hoagland solution, which was placed at a relative humidity of 75% and 26/20°C day/night temperature with a light intensity of 3000 lx. Forty-day-old plants treated with 200 mM NaCl for 24 h and non-treated plants were used to prepare small RNA libraries.
Small RNA library construction and sequencing
Total RNA was extracted from T. salsuginea using Trizol Reagent (Invitrogen, USA), following the manufacturer’s protocol. RNA quality was evaluated by electrophoresis on a 1% agarose gel. RNA concentration was quantified using BioSpectrometer fluorescence (Eppendorf, GER). Two small RNA libraries of T. salsuginea were generated and were sequenced by Solexa technology. The sequencing procedure was conducted as described in our previous work . Briefly, after PAGE purification of small RNA molecules and ligation of a pair of Solexa adaptors to their 5' and 3' ends, the small RNAs were amplified for 17 cycles using the adaptor primers, for cluster generation and sequencing analysis. The image files generated by the sequencer were then processed to produce digital-quality data by the BGI small RNA pipeline. After masking the adaptor sequences and removing contaminated reads, clean reads were processed for computational analysis. A computational pipeline was developed to process the sequencing data and to distinguish miRNA candidates from other types of small RNAs.
The raw reads were processed using PHRED and CROSS MATCH. Vector sequences, low quality sequences, and sequences that are shorter than 10 nt were removed. The clean reads were used for further analysis. Small RNA tags were mapped to the T. salsuginea genomes using the SOAP program . rRNA, scRNA, snoRNA, snRNA, and tRNA were removed by conducting a BlastN search against GenBank and Rfam. The remaining small RNAs were aligned with exons and introns of genes to identify degraded mRNA fragments.
The T. salsuginea genome sequence was downloaded from the T. salsuginea Sequencing Resource website (http://www.phytozome.net/thellungiella.php; T. salsuginea genome reference sequence ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Thalophila/assembly/). Small RNAs were mapped to the genome using PatMaN. The hairpin structure was predicted using online software mfold (http://www.bioinfo.rpi.edu/applications/mfold/). MiRNA candidates were predicted using miRCat (http://srna-tools.cmp.uea.ac.uk/plant)  with default parameters. T. salsuginea transcript reference sequences (ftp://ftp.jgi-psf.org/pub/compgen/phytozome/v9.0/Thalophila/annotation/ annotation/Thalophila_173_transcript.fa.gz) and NCBI dbEST (http://www.ncbi.nlm.nih.gov/nucest/?term=Thellungiella) were used for the target prediction. The number of small RNAs was normalized by TP10M (tags per ten million), which represents the abundance of small RNAs. The TL/CL ratio was used to calculate the fold changes in expression level of known and novel miRNAs.
Identification of known miRNAs and their variants
The mature miRNAs of miRBase 19.0 (http://www.mirbase.org/) were downloaded and aligned with small RNA tags by BlastN: blastall-p blastn-F F -e 0.01. Perfectly matched sequences were considered as known miRNAs in T. salsuginea, and the other sequences, with no more than two mismatches or gaps, were considered as variants. In addition, the abundance of the variant was no less than 10 counts in at least one miRNA library. The known miRNAs* were predicted in the hairpin structures using Mireap software.
Up to April 2013, 25141 mature miRNAs from 193 species are available in miRBase 19.0. An miRNA gene always has different variants due to imprecise processing or various biochemical modifications, and the miRNA annotation is defined by the dominant sequence if it is not the same as sequences in miRBase. Deep sequencing could detect miRNA* sequences that are too rare to be cloned by traditional cloning and sequencing methods. MiRNA genes show a characteristic “mapping pattern”, which is the result of sequential processing of Drosha and Dicer. According to this feature, false miRNAs can be removed from the known small RNA dataset by manual check. The number of reads sequenced was used to measure the abundance of a specific miRNA species in the library. Due to inaccurate processing and base modification, miRNA always has more than one sequence form, and the sum of all variants of an miRNA gene could be used as a digital measure of its expression.
Identification of novel miRNAs and plausible candidates of novel miRNAs
Novel miRNAs and plausible candidates of novel miRNAs were predicted based on the hairpin structure of miRNA precursors using Mireap software with the following parameters: dominant mature sequences reside at the stem region with a size of 18–25 nt; the abundance is at least ten counts; the precursor of each novel miRNA was identified and could form a proper secondary hairpin structure; and the distance between miRNA and miRNA* is 16–300 nt. The flank sequence length of the miRNA precursor is 20 nt and the maximal free energy allowed to form a hairpin structure is -18 kcal/mol. miRNAs were considered novel if the miRNA* sequence was detected in at least one library. miRNAs for which the miRNA* was not detected in the sequenced libraries were considered as plausible novel miRNA candidates.
Identification and functional annotation of the target genes
Target genes of conserved and novel miRNAs were predicted using perl software pick_plant_target.pl with the following parameters: (1) no more than four mismatches were permitted between the sRNA and target (G-U bases count as 0.5 mismatches); (2) no more than two adjacent mismatches in the miRNA/target duplex; (3) no adjacent mismatches in positions 2–12 of the miRNA/target duplex (5' of miRNA); (4) no mismatches in positions 10–11 of the miRNA/target duplex; (5) no more than 2.5 mismatches in positions 1–12 of the miRNA/target duplex (5' of miRNA); and (6) the minimum free energy (MFE) of the miRNA/target duplex should be ≥ 75% of the MFE of miRNA bound to its perfect complement. The coding sequences of predicted target genes were downloaded (http://www.phytozome.net/search.php?show=text&method=Org_Tsalsuginea, http://www.ncbi.nlm.nih.gov/nucest/?term=Thellungiella) and submitted to NCBI for BlastX analysis with an e-value of 1e-10 to identify putative gene homologs and the functions of potential targets.
MiRNA expression profiling using microarray analysis
In Affymetrix genechip Sanger miRNA V11.0, about 400 probes corresponding to 45 unique plant miRNAs with sequences identical to T. salsuginea miRNA sequences. Total RNAs were extracted from leaves and roots using Trizol Reagent (Invitrogen, USA). One microgram of total RNA was used to isolate low molecular weight RNA using polyethylene glycol solution precipitation. Low molecular weight RNAs were labeled with Cy3 and hybridized with miRNA microarrays (Shanghai BioCorp., China). Arrays were scanned using a GeneChip®Scanner 3000. Hybridization signals were extracted from the TIFF images using GCOS1.4 software. miRNA QC Tool software was used for data summarization, normalization, and quality control. Low signal probes and probes with fewer than three replicates on the chip were removed.
Detection of miRNA expression using stem-loop qRT-PCR
The stem-loop qRT-PCR primers (Additional file 9: Table S10) were designed and the reverse transcription was carried out according to previous studies [55, 56]. Stem-loop qRT-PCR was performed in an ABI StepOne Real-Time PCR System (Applied Biosystems) using SYBR Green I (TOYOBO) with the following program: 94°C for 10 min, followed by 40 cycles of 94°C for 15 s, 60°C for 10 s, and 72°C for 25 s. All reactions were conducted in triplicate, U6 was used as an internal control, and the sample without template was used as a negative control. The relative expression of the miRNA was calculated using ABI StepOne Real-Time PCR System software and the 2-ΔΔCt method.
Verification of the miRNA cleavage site by 5’ RACE
Total RNA was isolated from NaCl-treated and control plants using the CTAB method. Equal amounts of total RNA were pooled together and used for cDNA synthesis following the protocols of the SMARTer™ RACE cDNA Amplification Kit (Clontech). Gene-specific reverse nest primers were designed from the predicted targets and used in combination with the 5’ adapter primers (UPM, NUP) to amplify the cleaved transcripts. PCR reactions were separated by agarose gel electrophoresis, and distinct bands of the expected size were cloned into pMD18 simple-T vectors (Takara, Dalian, China) and sequenced using the general primer M13.
Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res. 2003, 13: 137-144. 10.1101/gr.751803.
Koch MA, Haubold B, Mitchell-Olds T: Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). Mol Biol Evol. 2000, 17: 1483-1498. 10.1093/oxfordjournals.molbev.a026248.
Gong Q, Li P, Ma S, Rupassara SI, Bohnert HJ: Salinity stress adaptation competence in the extremophile Thellungiella halophila in comparison with its relative Arabidopsis thaliana. Plant J. 2005, 44: 826-839. 10.1111/j.1365-313X.2005.02587.x.
Volkov V, Wang B, Dominy PJ, Fricke W, Amtmann A: Thellungiella halophila, a salt-tolerant relative of Arabidopsis thaliana, possesses effective mechanisms to discriminate between potassium and sodium. Plant Cell Environ. 2003, 27: 1-14.
Xiong L, Zhu JK: Molecular and genetic aspects of plant responses to osmotic stress. Plant Cell Environ. 2002, 25: 131-139. 10.1046/j.1365-3040.2002.00782.x.
Kant S, Kant P, Raveh E, Barak S: Evidence that differential gene expression between the halophyte, Thellungiella halophila, and Arabidopsis thaliana is responsible for higher levels of the compatible osmolyte proline and tight control of Na + uptake in T. halophila. Plant Cell Environ. 2006, 29: 1220-1234. 10.1111/j.1365-3040.2006.01502.x.
Bressan RA, Zhang C, Zhang H, Hasegawa PM, Bohnert HJ, Zhu JK: Learning from the Arabidopsis experience. The next gene search paradigm. Plant physiology. 2001, 127: 1354-1360. 10.1104/pp.010752.
Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress induced miRNA. Mol Cell. 2004, 14: 787-799. 10.1016/j.molcel.2004.05.027.
Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.
Subramanian S, Fu Y, Sunkar R, Barbazuk WB, Zhu JK, Yu O: Novel and nodulation-regulated microRNAs in soybean roots. BMC Genomics. 2008, doi:10.1186/1471-2164-9-160
Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Bio. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.
Kidner CA, Martienssen RA: The developmental role of microRNA in plants. Curr Opin Plant Biol. 2005, 8: 38-44. 10.1016/j.pbi.2004.11.008.
Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003, 15: 2730-2741. 10.1105/tpc.016238.
Chen XM: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303: 2022-2025. 10.1126/science.1088060.
Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425: 257-263. 10.1038/nature01958.
Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16: 2001-2019. 10.1105/tpc.104.022830.
Lu W, Li J, Liu F, Gu J, Guo C, Xu L, Zhang H, Xiao K: Expression pattern of wheat miRNAs under salinity stress and prediction of salt-inducible miRNAs targets. Front Agric China. 2011, 5: 413-422. 10.1007/s11703-011-1133-z.
Navarro L, Dunoyer P, Jay F, Arnold B, Dharmasiri N, Estelle M, Voinnet O, Jones JD: A plant miRNA contributes to antibacterial resistance by repressing auxin signaling. Science. 2006, 312: 436-439. 10.1126/science.1126088.
Liu HH, Tian X, Li YJ, Wu CA, Zheng CC: Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA. 2008, 14: 836-843. 10.1261/rna.895308.
Dassanayake M, Oh DH, Haas JS, et al: The genome of the extremophile crucifer Thellungiella parvula. Nat Genet. 2011, 43: 913-918. 10.1038/ng.889.
Wu HJ, Zhang Z, Wang JY, et al: Insights into salt tolerance from the genome of Thellungiella salsuginea. USA: Proc. Natl Acad. Sci; 2012. doi:10.1073/pnas.1209954109
Zhu JK: Plant salt tolerance. Trends Plant Sci. 2011, 6: 66-71.
Inan G, Zhang Q, Li P, et al: Salt cress. A halophyte and cryophyte Arabidopsis relative model system and its applicability to molecular genetic analyses of growth and development of extremophiles. Plant physiology. 2004, 135: 1-20. 10.1104/pp.900109.
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20: 3407-3425. 10.1101/gad.1476406.
Sunkar R, Jagadeeswaran G: In silica identification of conserved microRNAs in large number of diverse plant species. BMC Plant Biol. 2008, 8: 13-10.1186/1471-2229-8-13. doi:10.1186/1471-2229-8-37
Ambros V, Bartel B, Bartel DP: A uniform system for microRNA annotation. RNA. 2003, 9: 277-279. 10.1261/rna.2183803.
Meyer BC, Axtell MJ, Bartel B, et al: Criteria for annotation of plant microRNAs. Plant Cell. 2008, 20: 3186-3190. 10.1105/tpc.108.064311.
Tong L, Lin L, Wu S, et al: MiR-10a* up-regulates coxsackievirus B3 biosynthesis by targeting the 3D-coding sequence. Nucleic Acids Res. 2013, 41: 3760-3771. 10.1093/nar/gkt058.
Huang MD, Wu WL: Overexpression of TMAC2., a novel negative regulator of abscisic acid and salinity response., has pleiotropic effects in Arabidopsis thaliana. Plant Mol Biol. 2007, 63: 557-569. 10.1007/s11103-006-9109-8.
Perera IY, Hung CY, Moore CD, Stevenson-Paulik J, Boss WF: Transgenic Arabidopsis plants expressing the type 1 inositol 5-phoshatase exibit increased drought tolerance and altered abscisic acid signaling. Plant Cell. 2008, 20: 2876-2893. 10.1105/tpc.108.061374.
Lee JJ, Kim WT: Regulation of abiotic stress signal transduction by E3 ubiquitin ligases in Arabidopsis. Mol Cells. 2011, 31: 201-208. 10.1007/s10059-011-0031-9.
Singh K, Foley RC, Oñate-Sánchez L: Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002, 5: 430-436. 10.1016/S1369-5266(02)00289-3.
Hu H, Dai M, Yao J, Xiao B, Li X, Zhang Q, Xiong L: Overexpressing a NAM, ATAF, and CUC (NAC) transcription factor enhances drought resistance and salt tolerance in rice. Proc Natl Acad Sci USA. 2006, 103: 12987-12992. 10.1073/pnas.0604882103.
Ruiz-Ferrer V, Voinnet O: Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol. 2006, 60: 485-510. doi:10.1146/annurev.arplant.043008.092111
Nakashima K, Takasaki H, Mizoi J, Shinozaki K, Yamaguchi-Shinozaki K: NAC transcription factors in plant abiotic stress responses. Biochim Biophys Acta. 1819, 2012: 97-103. doi:10.1016/j.bbagrm.2011.10.005
Arenas-Huertero C, Pérez B, Rabanal F, Blanco-Melo D, De la Rosa C, Estrada-Navarrete G, Sanchez F, Covarrubias AA, Reyes JL: Conserved and novel miRNAs in the legume Phaseolus vulgaris in response to stress. Plant Mol Biol. 2009, 70: 385-401. 10.1007/s11103-009-9480-3.
Kong YM, Elling AA, Chen B, Deng XW: Differential expression of microRNAs in maize inbred and hybrid lines during salt and drought stress. Am J Plant Sci. 2010, 1: 69-76. 10.4236/ajps.2010.12009.
Zhou L, Liu Y, Liu Z, Kong D, Duan M, Luo L: Genome wide identification and analysis of drought-responsive microRNAs in Oryza sativa. J Exp Bot. 2010, 61: 4157-4168. 10.1093/jxb/erq237.
Li W, Cui X, Meng Z, Huang X, Xie Q, Wu Heng, Jin H, Zhang D, Liang W: Transcritional regulation of Arabidopsis miR168a and ARGONAUTE1 homeostasis in abscisic acid and abiotic stress responses. Plant Physiol. 2012, 158: 1279-1292. 10.1104/pp.111.188789.
Trindada I, Capitão C, Dalmay T, Fevereiro MP, Santos DM: miR398 and miR408 are up-regulated in response to water deficit in Medicago truncatula. Planta. 2010, 231: 705-716. 10.1007/s00425-009-1078-0.
Jia X, Wang WX, Ren L, Chen QJ, Mendu V, Willcut B, Dinkins R, Tang X, Tang G: Differential and dynamic regulation of miR398 in response to ABA and salt stress in Populus tremula and Arabidopsis thaliana. Plant Mol Biol. 2009, 71: 51-59. 10.1007/s11103-009-9508-8.
Uozumi N, Kim EJ, Rubio F, Yamaguchi T, Muto S, Tsuboi A, Bakker EP, Nakamura T, Schroeder JI: The Arabidopsis HKT1 gene homolog mediates inward Na + currents in Xenopus laevis oocytes and Na + uptake in Saccharomyces cerevisiae. Plant Physiol. 2000, 122: 1249-1259. 10.1104/pp.122.4.1249.
Rus A, Yokoi S, Sharkuu A, Reddy M, Lee BH, Matsumoto TK, Koiwa H, Zhu JK, Bressan RA, Hasegawa PM: AtHKT1 is a salt tolerance determinant that controls Na+ entry into plant roots. Proc Natl Acad Sci USA. 2001, 98: 14150-14155. 10.1073/pnas.241501798.
Horie T, Yoshida K, Nakayama H, Yamada K, Oiki S, Shinmyo A: Two types of HKT transporters with different properties of Na + and K + transport in Oryza sativa. Plant J. 2001, 27: 129-138. 10.1046/j.1365-313x.2001.01077.x.
Walia H, Wilson C, Condamine P, Liu X, Ismail AM, Close TJ: Comparative transcriptional profiling of two contrasting rice genotype under salinity stress during the vegetative growth stage. Plant Physiol. 2005, 139: 822-835. 10.1104/pp.105.065961.
Amtmann A, Sanders D: Mechanisms of Na+ uptake by plant cells. Adv Bot Res. 1999, 29: 75-112.
Su H, Golldack D, Katsuhara M, Zhao C, Bohnert HJ: Expression and stress dependent induction of potassium channel transcripts in the common ice plant. Plant Physiol. 2001, 125: 604-614. 10.1104/pp.125.2.604.
Griffith M, Timonin M, Wong AC, Gray GR, Akhter SR, Saldanha M, Rogers MA, Weretilnyk EA, Moffatt B: Thellungiella: an Arabidopsis-related model plant adapted to cold temperatures. Plant Cell Environ. 2007, 30: 529-538. 10.1111/j.1365-3040.2007.01653.x.
Wong CE, Li Y, Whitty BR, Diaz-Camino C, Akhter SR, Brandle JE, Golding GB, Weretilnyk EA, Moffatt BA, Griffith M: Expressed sequence tags from the Yukon ecotype of Thellungiella reveal that gene expression in response to cold, drought and salinity showed little overlap. Plant Mol Biol. 2005, 58: 561-574. 10.1007/s11103-005-6163-6.
Arbona V, Argamasilla R, Gómez-Cadenas A: Common and divergent physiological, hormonal and metabolic responses of Arabidopsis thaliana and Thellungiella halophila to water and salt stress. J Plant Physiol. 2010, 167: 1342-1350. 10.1016/j.jplph.2010.05.012.
Lugan R, Niogret MF, Leport L, Guegan JP, Larher FR, Savoure A, Kopka J, Bouchereau A: Metabolome and water homeostasis analysis of Thellungiella salsuginea suggests that dehydration tolerance is a key response to osmotic stress in this halophyte. Plant Journal. 2010, 64: 215-229. 10.1111/j.1365-313X.2010.04323.x.
Zhao CZ, Xia H, Frazier TP, Yao YY, Bi YP, Li AQ, Li MJ, Li CS, Zhang BH, Wang XJ: Deep sequencing identifies novel and conserved microRNAs in peanuts (Arachis hypogaea L.). BMC Plant Biol. 2010, 5: 10:3-doi:10.1186/1471-2229-10-3
Li R, Li Y, Kristiansen K, Wang J: SOAP: short oligonucleotide alignment program. Bioinformatics. 2008, 24: 713-714. 10.1093/bioinformatics/btn025. doi:10.1093/bioinformatics/btn025
Moxon S, Schwach F, Maclean D, Dalmay T, Studholme DJ, Moulton V: A toolkit for analysing large-scale plant small RNA datasets. Bioinformatics. 2008, 24: 2252-2253. 10.1093/bioinformatics/btn428.
Chen C, Ridzon DA, Broomer AJ, et al: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33 (20): e179-10.1093/nar/gni178.
Shen J, Xie K, Xiong L: Global expression profiling of rice microRNAs by one-tube stem-loop reverse transcription quantitative PCR revealed important roles of microRNAs in abiotic stress responses. Mol Genet Genomics. 2010, 284: 477-488. 10.1007/s00438-010-0581-0.
This work is supported by the National Basic Research Program of China (2012CB114204), National Natural Science Foundation of China (30870196, 31101427), the Initial Special Research for 973 Program of China (2012CB126313), Shandong Province Taishan Scholar Foundation (tshw20100416), grants from Shandong Province (2008GG10002002, ZR2010CQ008, 201103023).
The authors declare that they have no competing interests.
QZ, XW and YZ designed the study. QZ and CZ carried out most of the experiments, data analysis and wrote the material and method part of the manuscript. XW, JP and YZ wrote the manuscript, made the figures and finalized the tables. WS, HX, YL, AL, CL, ML, MS, SZ, LH, performed experiments and took care of the plants. All authors read and approved the final manuscript.
Quan Zhang, Chuanzhi Zhao contributed equally to this work.