Genome-wide identification of Thellungiella salsuginea microRNAs with putative roles in the salt stress response
- Quan Zhang†1,
- Chuanzhi Zhao†2,
- Ming Li2,
- Wei Sun1,
- Yan Liu1,
- Han Xia2,
- Mingnan Sun2,
- Aiqin Li2,
- Changsheng Li2,
- Shuzhen Zhao2,
- Lei Hou2,
- Jean-François Picimbon2,
- Xingjun Wang1, 2Email author and
- Yanxiu Zhao1Email author
© Zhang et al.; licensee BioMed Central Ltd. 2013
Received: 6 August 2013
Accepted: 4 November 2013
Published: 15 November 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.
KeywordsThellungiella salsuginea Salt stress miRNA identification Solexa sequencing Expression analysis
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
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
Known miRNA families identified in Thellungiella
Counts in CL
Counts in TL
Conserved miRNA families
Non-conserved miRNA families
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.
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
Novel miRNAs identified in Thellungiella
Novel miRNA sequence
Genomic location of pre-miRNA
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
Target gene prediction of known and novel miRNAs
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
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.
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).
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.View Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- Bartel DP: MicroRNAs: Genomics, biogenesis, mechanism, and function. Cell. 2004, 116: 281-297. 10.1016/S0092-8674(04)00045-5.PubMedView Article
- 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.View Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- Chen XM: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303: 2022-2025. 10.1126/science.1088060.PubMedView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.View Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedPubMed CentralView Article
- 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.View Article
- 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.View Article
- 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.PubMedPubMed CentralView Article
- 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-37View Article
- Ambros V, Bartel B, Bartel DP: A uniform system for microRNA annotation. RNA. 2003, 9: 277-279. 10.1261/rna.2183803.PubMedPubMed CentralView Article
- 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.View Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.092111View Article
- 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.PubMedView Article
- 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.View Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.View Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- Amtmann A, Sanders D: Mechanisms of Na+ uptake by plant cells. Adv Bot Res. 1999, 29: 75-112.View Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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.PubMedView Article
- 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/btn025PubMedView Article
- 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.PubMedView Article
- 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.PubMedPubMed CentralView Article
- 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.PubMedView Article
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.