Identification and characterization of cold-responsive microRNAs in tea plant (Camellia sinensis) and their targets using high-throughput sequencing and degradome analysis
BMC Plant Biology volume 14, Article number: 271 (2014)
MicroRNAs (miRNAs) are approximately 19 ~ 21 nucleotide noncoding RNAs produced by Dicer-catalyzed excision from stem-loop precursors. Many plant miRNAs have critical functions in development, nutrient homeostasis, abiotic stress responses, and pathogen responses via interaction with specific target mRNAs. Camellia sinensis is one of the most important commercial beverage crops in the world. However, miRNAs associated with cold stress tolerance in C. sinensis remains unexplored. The use of high-throughput sequencing can provide a much deeper understanding of miRNAs. To obtain more insight into the function of miRNAs in cold stress tolerance, Illumina sequencing of C. sinensis sRNA was conducted.
Solexa sequencing technology was used for high-throughput sequencing of the small RNA library from the cold treatment of tea leaves. To align the sequencing data with known plant miRNAs, we characterized 106 conserved C. sinensis miRNAs. In addition, 215 potential candidate miRNAs were found, among, which 98 candidates with star sequences were chosen as novel miRNAs. Both congruously and differentially regulated miRNAs were obtained, and cultivar-specific miRNAs were identified by microarray-based hybridization in response to cold stress. The results were also confirmed by quantitative real-time polymerase chain reaction. To confirm the targets of miRNAs, two degradome libraries from two treatments were constructed. According to degradome sequencing, 455 and 591 genes were identified as cleavage targets of miRNAs from cold treatments and control libraries, respectively, and 283 targets were present in both libraries. Functional analysis of these miRNA targets indicated their involvement in important activities, such as development, regulation of transcription, and stress response.
We discovered 31 up-regulated miRNAs and 43 down-regulated miRNAs in ‘Yingshuang’, and 46 up-regulated miRNA and 45 down-regulated miRNAs in `Baiye 1’ in response to cold stress, respectively. A total of 763 related target genes were detected by degradome sequencing. The RLM-5′RACE procedure was successfully used to map the cleavage sites in six target genes of C. sinensis. These findings reveal important information about the regulatory mechanism of miRNAs in C. sinensis, and promote the understanding of miRNA functions during the cold response. The miRNA genotype-specific expression model might explain the distinct cold sensitivities between tea lines.
MicroRNAs (miRNAs) are a class of non-coding RNAs, approximately 19 ~ 21 nucleotides (nt) long, that function as post-transcriptional regulators in eukaryotes . The miRNA gene is processed by Dicer-like proteins into a stem-loop miRNA::miRNA* duplex, after transcription by Pol II or Pol III enzyme into primary miRNA . The miRNA::miRNA* duplex is then cleaved and transported from the nucleus into the cytoplasm. Single stranded miRNA then joins with Argonaute (AGO) to form an RNA-induced silencing complex (RISC) . Finally, the RISC down-regulates targets by either cleaving the target mRNAs or repressing translation . Plants miRNAs have important functions in response to biotic and abiotic stresses . In recent years, numerous miRNAs have been identified in plant genomes , suggesting that the identification of their target RNAs is essential for the functional analysis of miRNA.
Cold stress negatively affects plant growth and development by causing tissue injury and delayed growth, which significantly restrict the spatial distribution of plants and productivity of economic crops . Besides transcriptional regulation, miRNAs are also involved in cold-responsive gene regulatory networks. Sunkar and Zhu  showed that the expression levels of miR393 and miR319c are up-regulated by cold treatment. Microarray analysis revealed that approximately 17% of Arabidopsis miRNAs are up-regulated in response to low temperature at early stages of cold treatment . The abundance of miR169 and miR172 in Arabidopsis challenged with cold stress was determined via a computational, transcriptome-based approach and microarray analysis ,. Solexa sequencing analysis showed that the expression levels of three conserved miRNAs (miR169e, miR172b, and miR397) and 25 predicted miRNAs exhibit significant changes in response to cold stress in Brachypodium. Cold resistance of the plant depends on different regulatory gene expression types related to physiology, metabolism, and growth . In rice, 18 cold responsive rice miRNAs were identified using microarrays, and the members of the miR171 family showed diverse expression patterns . Deep sequencing led to the identification of 30 cold responsive miRNAs in Populous tomentosa. Although miRNAs have been extensively studied in Arabidopsis and other plant species, no systematic examination of miRNA has been performed on C. sinensis. Prabu  and Das  identified numerous conserved miRNAs and their targets in C. sinensis through in silico analysis. Six novel small RNA candidates were isolated and cloned; the small RNAs were validated through expression analysis in young and old leaves, during non-dormant and dormant growth phases of C. sinensis. However, further study is needed to elucidate the functions of miRNAs at a genome-wide level in response to cold stress in C. sinensis.
Tea plant (C. sinensis) is one of the most important commercial beverage crops in the world. Cold stress may negatively affect the growth, development, and spatial distribution of tea plant, decreasing its yield and quality. Generally, cultivar specific expression exhibits strong relevance to the physiological functions of the corresponding cultivars . Understanding cultivar-specific expression patterns of miRNA is necessary to gain insight into the functions of miRNA. Thus, ‘Yingshuang’ (YS, a cold-tolerant tea plant cultivar) and `Baiye 1’ (BY, a cold-sensitive tea plant cultivar) were chosen as two cultivars. In our study, high-throughput Solexa sequencing (Illumina Genome Analyzer) was employed to identify the C. sinensis miRNAs, which were responsive to cold stress, and 106 conserved miRNAs were obtained in the small RNA library. A selected number of cold-responsive and new miRNAs were then validated by Quantitative real-time polymerase chain reaction (qRT-PCR) combined with computational analysis. The identified miRNAs and their potential miRNA targets were predicted and confirmed by degradome sequencing. Abundantly conserved sequenced signatures were identified as the targets cleaved by conserved miRNAs, and novel miRNAs targeted different genes with various biological functions.
High-throughput sequencing of small RNAs in tea plant
Tea plants were stored at 4°C and 28°C for 1, 4, 8, 12, 24 and 48 h, respectively. A small RNA library of tea leaves, which was generated from a mixture of total RNAs from each cold-treatment stage, was subjected to high-throughput sequencing by the Illumina platform. Raw sequences were first subjected to an Illumina Pipeline filter provided by the supplier (Solexa 0.3). A total of 9,700,042 raw reads, representing 3,145,122 distinct sequences, were obtained. Reads without small RNA sequences, ranging from 15 nt to 30 nt in length, were filtered (Figure 1). The majority of the RNA sequences ranged from 19 nt to 25 nt in size. The most abundant small RNAs in the library were 24 nt long. The distribution of 24 nt small RNAs was approximately 45.96% and 69.67% in the total and unique sequences, respectively, whereas the distribution of 21 nt small RNAs in the total and unique sequences was approximately 13.68% and 5.67%, respectively. A total of 1,319,524 clean reads were obtained from the tea plant, including Rfam, rRNA, tRNA, snoRNA, snRNA, miRNA, other ncRNA, and repeats (Table 1). These clean reads were obtained by removing adaptor/acceptor sequences, filtering low quality tags, and cleaning the contaminants formed by adaptor-adaptor ligations and shorts RNAs less than 15 nt.
Conserved miRNAs in tea plant
To identify the conserved miRNAs in tea plant, we compared our dataset with known plant miRNAs, such as miRNA precursors and mature miRNAs, in miRBase 19.0 by miRAlign. Following the BLASTN searches and further sequence analysis, 106 unique sequences belonging to 25 families in the small RNA library were found to be orthologs of known miRNAs from other plant species, which were previously deposited in the miRBase database (Additional file 1: Table S1). Moreover, 57 miRNA*s were identified, which are considered to be strong evidence of bona fide miRNAs . Previous studies have predicted some miRNAs in tea plant. Six conserved miRNAs (csn-miR156a, csn-miR164, csn-miR169, csn-miR171a, csn-miR399, and csn-miR408) were identified and verified to those reported previously ,,,.
The number of different conserved miRNA family members was analyzed. The majority of the 25 miRNA families contained several members, and two families (miR166 and miR171) possessed multiple members, with 17 and 10 members, respectively, whereas seven miRNA families had only one member (Figure 2). The frequency of diverse members sequenced from the same or different miRNA families also varied drastically, ranging from one to 159,305 times. Among the 25 identified conserved csn-miRNAs, six miRNAs (csn-miR166a-1, csn-miR166a-2, csn-miR166a-3, csn-miR166a-4, csn-miR166a-5, and csn-miR166a-6) had the most reads, reaching up to 159,305. Moreover, four miRNA families were represented in tens of thousands, whereas some miRNAs (e.g., csn-miR156d, csn-miR395b, and csn-miR396f) had only one read sequenced. This large discrepancy in the expression levels of csn-miRNAs, deduced from the number of reads sequenced, could reflect the divergence of potential functions during the different stages of cold stress.
Novel miRNAs in tea plant
In our study, the stem-loop structure of miRNA precursors was used to predict novel miRNAs, and the secondary structures of novel miRNA precursors were obtained by Mfold . The secondary hairpin structures of the representative miRNAs are listed in Additional file 2: Figure S1. A total of 215 sequences were predicted to be potentially non-conserved miRNAs from the remaining unannotated sRNAs. Using updated plant miRNA annotation criteria , 98 sequences were recognized as novel miRNAs with high confidence and designated as novel C. sinensis miRNAs (Additional file 3: Table S2). The length of the mature miRNAs varied from 19 nt to 25 nt, with the majority being 24 nt long. Furthermore, the length of the novel miRNA precursors ranged from 72 nt to 264 nt, with an average length at 150 nt. The minimum free energy (MFE) of these novel miRNA precursors varied from -114.3 kcal mol−1 to −17.6 kcal mol−1, with an average of −61.79 kcal mol−1. The MFE index (MFEI) is a unique criterion to designate miRNAs. The MFEI is calculated by the equation: MFEI = (100 × MFE/L) / (G + C)% (L: the length of pre-miRNA). The sequence is most likely to be miRNA when the MFEI is more than 0.85 . These novel miRNAs showed lower abundance levels, compared with the conserved miRNAs, which was consistent with previous studies -. The novel miRNAs showed different expression levels, and their normalized reads were from one to 1644. Thus, csn-smR30 (1644), csn-smR65 (794), and csn-smR80 (668) were the most abundant miRNAs. Most novel miRNAs were observed in less than 100 times, whereas 43 (43.89%) csn-smRNAs were sequenced in less than 10 times. The nucleotide bias at each position of the 98 novel identified miRNA (Additional file 3: Table S2) shoved that the first nucleotide of the new miRNA genes tended to be (U) in general. As expected, miRNAs are loaded to the RISC assisted by AGO1. Research has shown that AGO1 proteins have more affinity with uracil in the 5′ terminus of miRNA, thus resulting in cloned miRNA sequences with uracil nucleotide bias in the first position .
miRNA microarray chip content and hybridization of arrays
Microarray-based hybridization was performed to analyze the expression of the newly identified miRNAs in tea plant. The sequences from miRBase (http://microrna.sanger.ac.uk/sequences/), and newly identified sequences in this study were used as probes for chip hybridization. The miRNAs showed different expression profiles under cold stress (4, 12, and 24 h) and non-treated conditions. Among the 3511 miRNA probes by microarray, a total of 303 and 349 conserved miRNAs were observed in ‘Yingshuang’ (YS, a cold-tolerant tea plant cultivar) and `Baiye 1’ (BY, a cold-sensitive tea plant cultivar) respectively (Additional file 4: Figure S2 and Additional file 5: Figure S3). The detected miRNAs were defined as a value of hybridization signal greater than 500, expression of miRNAs was significant difference when the signal ratio greater than 2 (|log|>1) and p was less than 0.01. Based on this principle, 158 tea plant miRNAs were differentially expressed compared with expression patterns under different cold stress stages in YS, and 159 miRNAs were differentially expressed in BY (Additional file 6: Table S3 and Additional file 7: Table S4), including 87 conserved miRNAs (p <0.01 and Signal >500) in both cultivars.
Both congruously and differentially regulated miRNAs were observed in our study, as well as cultivar-specific miRNAs. The majority of the differentially expressed miRNAs showed different expression patterns either among three cold stress stages, or between two tea cultivars. In YS, all of 31 miRNAs showed up-regulated trends, for example miR164, miR167, miR168, miR171, and so on, whereas 43 miRNAs (miR156, miR319, miR474, miR529, and the rest) showed down-regulated trends. By contrast, 46 miRNAs presented up-regulated trends, for instance miR168, miR474, miR1160, and so forth, while 45 miRNAs (miR159, miR166, miR171, miR529, etc.) presented down-regulated trends in BY. Three miRNA families (miR168, miR152, and miR2936) were uniformly regulated at four cold stress stages in the two plant cultivars. However, the expression level of miR171 and miR474 gradually increased in YS, but gradually declined in BY. These results strongly indicate that the regulatory patterns may be in accordance with delayed expression patterns in the cold-sensitive tea cultivar, which partly explains the distinct cold sensitivities between the two cultivars.
To confirm the microarray results, the abundance of several miRNAs was further analyzed by qRT-PCR. The result of qRT-PCR and abundance profiles of the microarray shared similar trends (Figure 3). Discrepancies were also found in the magnitude of response at different cold stress stages, which could be due to cross hybridization between the probe and other highly homologous miRNA family members. The discrepancies could also be due to data normalization between the two methods. The qRT-PCR data were normalized to the abundance of 5.8S rRNA, whereas the microarray data were normalized to the global abundance of all miRNAs detected by microarray. The use of 5.8S rRNA was technically unfeasible as the normalization standard for microarray data.
Targets identification for tea plant miRNA
We performed genome-wide analysis of miRNA-cleaved mRNAs to identify miRNA targets using high-throughput degradome sequencing technology ,. We sequenced 9,224,714 and 6,736820 signatures for each library (−C and + C). After removing duplications, 7,439,589 and 5,376,267 distinct reads were obtained for -C and + C libraries, respectively. Alignment of the distinct sequences to tea plant expressed sequence tag (EST) sequences yielded 37,088 and 37,011 unique signatures for -C and + C libraries, respectively. We identified sliced targets for known miRNA and novel miRNA candidates based on the method of CleaveLand pipeline . The abundance of sequences was plotted on each transcript (Additional file 8: Figure S4 and Additional file 9: Figure S5), and the sliced target transcripts were grouped into five classes according to the relative abundance of tags at the target sites . Based on this approach, Category 0 and Category 1 have more than one raw read at the position. The abundance at the position is equal to the maximum on the transcript, and only one maximum is on the transcript. Category 2 has more than one raw read at the position. The abundance at the position is less than the maximum but higher than the median for the transcript. Category 3 has more than one raw read at the position. The abundance at the position is equal or less than the median for the transcript. Category 4 has only one raw read at the position.
A total of 514 target transcripts were identified for 13 known miRNA families (Additional file 10: Table S5) based on our dataset, which shows most of the targets cleaved by the conserved miRNAs. A total of 332 targets in the + C library were identified for known conserved miRNA families, from which 9 (2.71%), 110 (33.13%), 98 (29.52%), 15 (4.52%) and 100 (30.12%) were grouped into categories 0, 1, 2, 3, and 4, respectively. For the -C library, 371 targets were identified, from which 19 (5.12%), 2 (0.54%), 96 (25.88%), 48 (12.94%), and 206 (55.52%) were grouped into category 0, 1, 2, 3, and 4, respectively (Figure 4A). Of these targets, 35.41% (182) were identified in the -C library, 27.82% (143) were identified in the + C library, and 36.77% (189) were present in both conditions (Figure 5A). Among the 13 conserved miRNA families, four (miR167, miR390, miR393, and miR398) were identified to have less than 10 targets, whereas the others target multiple transcripts. miR319 and miR160 had the highest number of targets, with 126 and 85 transcripts, respectively (Additional file 10: Table S5).
Forty novel miRNAs from 249 new candidates targets were identified (Additional file 11: Table S6). Profiling of the targets from the + C library showed that 3 (2.44%), 19 (15.45%), 46 (37.40%), 5 (4.06%), and 50 (40.65%) targets could be classified into categories 0, 1, 2, 3, and 4, respectively, and 2 (0.91%), 4 (1.82%), 93 (42.27%), 37 (16.82%), and 84 (38.18%) targets in the -C library fell into category 0, 1, 2, 3, and 4, respectively (Figure 4B). Of these targets, 50.60% (126) were identified in the -C library, 11.65% (29) were identified in the + C library, and 37.75% (94) were present in both conditions (Figure 5B). The distribution patterns of the two libraries differed, which suggests that the cleavage of targets by miRNAs was affected by cold stress.
Based on BLASTX analysis, 39.58% of the identified miRNA targets were generally homologous to conserved target genes that have already been found in Arabidopsis thaliana. Most of these conserved target genes were protein-coding genes, including zinc finger family protein (C2H2 and C3HC4 type), late embryogenesis abundant family protein (LEA), dormancy/auxin-associated family protein, and drought-responsive family protein, which are involved in plant growth, differentiation, development, and abiotic stress, respectively -. Among the identified miRNA targets, VQ motif-containing protein was a miR160 and miR408 target. The VQ motif represents the core of a protein-protein interaction domain, which is consistent with the interaction between another VQ motif protein with an RNA polymerase σ-factor ,. Thus, the identified tea plant miRNAs could regulate a wide range of genes in development and other physiological processes.
Identification of miRNA-guided cleavage of target mRNA using RLM-RACE
miRNAs, like small interfering RNA (siRNA), can direct the cleavage of their mRNA targets when these messages have extensive complementarity to the miRNAs -. This miRNA-directed cleavage can be detected by using a modified form of 5′ RNA ligase-mediated RACE (RLM-5′ RACE) because the 3′ product of the cleavage has two diagnostic properties: (1) a 5′ terminal phosphate, making it a suitable substrate for ligation to an RNA adaptor using T4 RNA ligase, and (2) a 5′ terminus that maps precisely to the nucleotide that pairs with the tenth nucleotide of the miRNA ,. To verify the nature the csn-miRNA target genes and study how the csn-miRNA regulate their target gene, RLM-5′ RACE experiment was employed, which was carried out in this study for further characterization of csn-miRNAs functions. All six of the csn-miRNAs 5′ end of the mRNA fragment mapped to the nucleotide that pairs to the tenth nucleotide of one of the miRNAs validated by PCR (Figure 6). CV014890.1, JK476458.1, FS943373.1, FS954022.1, GD254786.1, and FS955921.1 were confirmed as the real targets of csn-miR319b-1, csn-miR396b-2, csn-miR396c, csn-miR398, and csn-miR408 respectively, since all the 5′ ends of the mRNA fragments were mapped to the nucleotide the pairs to the eleventh nucleotide of miRNA with higher frequencies than depicted for each pairing oligo. From the precise sequences of the csn-miRNAs results, we know that the miRNA-guided cleavage in C. sinensis obeyed the principle that base-paring to the 5′ ‘seed’ region of the miRNA was the dominant factor for the miRNA target recognition, and that the cleavage site was mostly located at the eleventh nucleotide, just 3′ of the ‘seed’ sequence . All the six targets were found to have specific cleavage sites corresponding to the miRNA complementary sequences and might be regulated by the miRNAs in the style of siRNAs  directing the cleavage of mRNA targets with extensive complementarity to the miRNAs . FS943373.1 is similar to Arabidopsis proteins coded by plant calmodulin-binding protein-related, FS954022.1 coded for a protein highly homologous to rubredoxin-like superfamily protein, GD254786.1 coded for a transposable element gene, while FS955921.1 code for a protein highly homologous to VQ motif-containing protein (Table 2).
Gene ontology (GO) function analysis of targets
GO categories were assigned to all targets, including 514 known targets and 249 new candidates, according to three ontologies in GO: cellular component, molecular function, and biological process (Figure 7). Comparing the target gene functions of two libraries, more than 50% of the genes were classified into cellular component, of which 11 genes function belong to cell wall in + C library, however, there was no such target genes in -C library (Figure 7A). Based on the molecular function, genes were finally classified into eight classes, the three mainly represented GO terms were receptor activity (31%), other binding (31%), and kinase activity (10%) in + C library, while other binding accounted for 28%, followed by enzyme activity for 23%, receptor activity only 8% in -C library (Figure 7B). In the biological process, the target gene functions focused on the metabolic process (31%) and regulation of transcription (29%) in + C library, while the two class processes were only 10% and 11% in °C library, respectively (Figure 7C). This difference in the function of the target genes showed tea plant cell structure was severe damaged under cold stress. Moreover, stress-response genes were also identified as miRNA targets, including salt stress response, heat shock protein binding, and water deprivation response. The results imply the possible function of miRNAs in the regulation of biological processes involved in cold-stress.
Identification of miRNAs in tea plant
Many highly conserved miRNAs that exhibit particular expression patterns with specific timing and tissue specificity, have critical functions in growth, development, differentiation, apoptosis, metabolism and biotic and abiotic stress responses, regulating specific target mRNAs. Some tea plant miRNAs and their target genes have been identified using bioinformatical approaches in previous reports. Fourteen new C. sinensis miRNAs were recently identified from 47,452 available C. sinensis ESTs, and these miRNAs potentially target 51 mRNAs, which can act as transcription factors, and participate in transcription and signal transduction . Recent advances in high-throughput sequencing methods have revolutionized the identification of low-abundance, novel miRNAs in various species ,-. However, no comprehensive study on a novel miRNA discovery has been reported for tea plant. This study aimed to identify the evolutionary known and potentially novel tea plant-specific miRNAs recovered from cold stress tea plant libraries. The differential expression of miRNAs associated with cold stress response was also analyzed. Thus, approximately nine million sRNA raw reads were obtained from the sRNA library, in which 25 conserved miRNA families and 98 potentially novel miRNAs were successfully identified. The read number varied from one (miR156, miR395, and miR396) to 159,305 (miR166) (Additional file 1: Table S1), suggesting dramatically varied expression patterns among each miRNA family. However, only a small proportion of the conserved and novel tea miRNAs was detected, because of the unavailability of full genome sequences of tea plant. The number of miRNAs identified from tea plant appears to be far from saturation, and numerous unknown miRNAs remain to be discovered.
For a broader perspective of high-throughput sequencing of small RNAs from tea plant, we observed that small RNAs of 24 nt dominated the library of unique species, which was reported for other plant species, such as A. thaliana, Citrus trifoliata, Medicago truncatula, and Citrus sinensis. Length distribution analysis is effective in assessing the composition of small RNA samples. The overall distribution pattern of small RNAs (21 nt sRNAs =5.67%, 24 nt sRNAs =69.67%) in tea plant was significantly different from that in Populous trichocarpa, a model forest species, in which 21 nt RNAs are more abundant (37.16%) and 24 nt RNAs are less frequent (<5%) . A difference was also present between tea plant small RNAs and monocot species of maize . These results indicate that the small RNA transcriptome was complex across plant species, and could be significantly different between phylogenetically distant plant families .
Based on deep sequencing and hairpin structure prediction, we successfully identified 98 novel miRNAs, and 53 novel miRNAs with complementary miRNA* strands. Precursors of these miRNAs formed secondary hairpin structures with free energies ranging from −114.3 kcal mol−1 to −17.6 kcal mol−1 (average 61.79 kcal mol−1) (Additional file 3: Table S2). The secondary structures are not a unique feature of miRNAs because random inverted repeats may also form hairpin structures . Identification of anti-sense miRNAs (miRNAs*) of the candidate miRNAs provides credible evidence that these miRNAs are authentic miRNAs . One difference between conserved and novel miRNAs is that novel miRNAs are expressed at substantially lower levels or in a tissue-specific or environmentally inducible pattern ,.
Different expression profiles of cold-responsive miRNAs between two tea plant cultivars
The cultivar-specific and cold-responsive miRNAs were further determined using a customized microarray. A total of 18 and 14 conserved cold-responsive miRNA families were identified from YS and BY, respectively. Moreover, six and nine new cold-responsive miRNA families were identified from YS and BY, respectively. A quarter of miRNAs were down-regulated (27.2% in YS and 28.3% in BY), indicating that the expression of target genes controlled by these miRNAs was turned on to adapt to cold stress. We found that many members of families (miR156, miR159, and miR396) showed no consistent regulatory pattern for cold stress response, thereby suggesting different functions of miRNAs from the same family (Figures 8 and 9). Moreover, there was a certain degree of overlap of cold-stress pathway of miR168, miR529, and miR2936 by comparing miRNAs expression patterns between different cultivars. By contrast, the expression of miR164, miR408, miR1511, miR5368, miR172, miR482, miR529, and miR1160 was significantly different in the two cultivars (Figures 8 and 9).
miRNAs are known to silence genes post-transcriptionally by guiding target mRNAs for degradation or repressing translation . NADP-ME is a key enzyme that catalyzes the oxidative decarboxylation of L-malate to yield pyruvate, CO2, and NADPH in the presence of a divalent cation . In plants, malate can be decarboxylated in the mitochondrial matrix through the action of NAD-ME to produce pyruvate, which is oxidized by the tricarboxylic acid cycle. The genes for these two enzymes are predicted to be the target of miR474. The miR474 family was down-regulated in YS and up-regulated in BY. The result indicates that a cold induced, non-specific responding pathway was possibly involved in the maintenance of energy supply.
The miR171 family is a largely conserved miRNA family, and microarray assay revealed that miR171a in Arabidopsis, is induced at 6 h in response to cold stress . In our study, the miR171 family was significantly up-regulated in YS and down-regulated in BY, which suggests that miR171 members may perform different functions. AGO1 is regulated by miR168, AGO1-catalyzed mRNA cleavage, and is the only member of the Arabidopsis AGO family involved in the miRNA-directed mRNA cleavage processes. The co-adjustment of miR168 and AGO1 indicates negative feedback regulation in the miRNA pathway . In our study, miR168 was found to be a cold-responsive miRNA and the members of the csn-miR168 family were induced in the two tea cultivars after 12 h of cold treatment.
Degradome analysis of cold responsive targets for tea plant miRNAs
Accurate identification of target genes is essential to reveal the regulatory function of miRNA. Previous work on gene identification in tea plants was limited to bioinformatics prediction . The recently developed high-throughput experimental approach allows the identification of target genes for known and new miRNAs ,. A total of 514 potential targets for 13 known miRNA families and 249 new tea-specific targets for 40 novel miRNA families in tea plant were identified by degradome sequencing (Additional file 10: Table S5 and Additional file 11: Table S6). Numerous highly conserved miRNAs were identified, and did not have detectable sliced targets, including miR159, miR166, miR168, and miR394 (Additional file 1: Table S1). The novel miRNAs targeted different genes covering various predicted functions (Additional file 11: Table S6). For example, csn-smR5749 targets nine genes of the Ras-related, small GTP-binding protein family, while csn-smR2592 affects zinc finger (C3HC4-type RING finger) protein family. These two protein families have been reported to be involved in various physiological processes in plants ,. By contrast, csn-smR2845 targets the dormancy or auxin-associated protein family, which has an important function in the inhibition of growth and differentiation . Recently, an increasing number of genetic and biochemical results showed that plant miRNA-guided silencing has a widespread translational inhibitory component ,.
Neighboring nucleotides (the ninth or 11th positions) are cleavage sites for some plant miRNAs in certain cases ,. This observation accounts for the occasional positional heterogeneity of mature miRNAs  and their cleavage products ,. Upon comparing the two degradome libraries, we classified many cold-responsive miRNA targets as categories 0, 1, and 2 (Additional file 10: Table S5), thereby confirming the accuracy of our degradome analysis, because categories 0, 1, and 2 are characteristic for miRNA targets in Arabidopsis. The most obvious difference between the targets of cold-treated and cold-free miRNAs was that 254 of 371 cold-free miRNA targets belonged to categories 3 and 4 (Additional file 10: Table S5), in which cleavage abundance was lower than the median on the target transcripts. The target distribution in each library was generally similar to those recently reported in other plant species ,. Moreover, the distribution pattern and transcript abundance of the targets between the two libraries differed, suggesting that the cleavage of targets by miRNAs was affected by cold stress. We also found target transcripts that were regulated by pairs of miRNAs. The miR160/miR319 pair targeted five different targets, the miR160/miR408 pair controlled two motif elements, the miR319/miR396 pair targeted protein prenyltransferase, and the miR171/miR396 pair targeted the Ras-related, small GTP-binding family protein (Additional file 10: Table S5). This phenomenon suggests a combinatorial gene regulation pathway using a pair of miRNAs in tea plants .
Degradome sequencing confirmed that csn-miR164 targeted a LEA protein. LEA proteins represent one of the functional elements, which are important factors in maintaining the viability of an organism and its biological structures in the ametabolic low water state, such as cold, heat, and drought . Part of the mechanism could be related to protein stabilization but other mechanisms have also been suggested for LEA proteins, because LEA proteins and other hydrophilins can preserve enzyme activity in vitro after desiccation or freezing -. In addition, csn-miR319 was detected to target a zinc ion binding and trans-membrane amino acid transporter protein family. Metal ion (copper, zinc, or iron) binding is important in normal respiratory function and the process of protein damage during oxidative stress in A. thaliana. Strong correlations exist between the sets of immobilized metal affinity chromatography-interacting proteins, which are proteins predicted to contain metal-binding motifs, and protein sets that are oxidized or degraded during abiotic stress . Zinc ion binding increases the structural complexity of abscisic acid stress ripening proteins, which is associated with DNA binding ,. The WW-domain-binding protein (WBP), which was previously found only in animals, has been predicted in tea plants and targeted by csn-miR408. WW-domains are small protein modules found in various proteins that participate in cell signaling or regulation . The WBP targets proline-rich sequence-recognizing domains, which are important for selectivity in signaling . The conserved miRNA families may have functions concerning the regulation of cytothesis, abscisic acid stress, signaling, and maintaining viability of biological structures in response to cold stress in tea plants.
The degradome sequencing data shows six novel miRNAs (csn-smR35, csn-smR3146, csn-smR5749, csn-smR7277, csn-smR8111, and csn-smR9722), which differentially target some cold-responsive genes in the cold-treated library. csn-smR5749 was predicted to target a KH domain protein. The KH domain proteins are important in RNA metabolism, vegetative, and reproductive development -. These proteins also have functions in organogenesis in stem cells by interacting with other proteins . The drought-responsive protein gene family was predicted to be the target of csn-smR9722, suggesting the existence of crosstalk between drought and cold stress signaling pathways.
miRNA generally function as negative regulators of gene expression by mediating the cleavage of target mRNAs  or by repressing their translation , the cleavage of target mRNAs appear to be predominant mode of gene regulation by plant miRNAs . To verify the nature the csn-miRNA target genes and study how the csn-miRNA regulate their target gene, we performed the RACE on unigene to detect and clone the mRNA fragment corresponding precisely to the product of miRNA processing. In total, we performed RLM-5′RACE assays on six target genes i.e. the representative targets of five conserved miRNAs (Figure 6). This result showed that the degradome analysis was an efficient and powerful approach which can be successfully used to validate the sequences of miRNAs.
This study identified 106 known miRNAs and 98 potentially novel miRNAs, based on the constructed miRNAomes between two tea plant cultivars, YS and BY. A total of 18 and 14 conserved cold responsive miRNA families were identified from YS and BY, respectively, using a customized microarray. A total of 455 and 591 genes were identified as cleavage targets of miRNAs from cold treatments and control libraries, respectively, and 283 targets were present in both libraries. The RLM-5′RACE procedure was successfully used to map the cleavage sites in six target genes of C. sinensis. GO annotation revealed that highly ranked miRNA target genes were those implicated in the developmental process, regulation of transcription and stress response. These findings provide valuable information for further functional characterization of miRNAs associated with cold stress in tea plants.
Plant materials and cold stress treatment
Two tea plant cultivars, ‘Yingshuang’ (YS, a cold-tolerant tea plant cultivar) and `Baiye 1’ (BY, a cold-sensitive tea plant cultivar), were culture-grown under a 12-h light (28°C)/12-h dark (22°C) photoperiod (photo intensity 1800 Lx photos m−2.s−2) for 20 days. To analyze the miRNAs of tea leave under cold stress, tea plants were stored at 4°C for 1, 4, 8, 12, 24, and 48 h. For the controls, tea plants were stored at 28°C for 1, 4, 8, 12, 24, and 48 h. The treated and untreated young leaves were then harvested into liquid nitrogen and stored at −80°C until total RNA extraction.
Small RNA library construction and sequencing
Total RNA was isolated with TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. RQ1 RNase-Free DNase (Promega, Madison, WI, USA) was used to remove genomic DNA contamination. The quality and integrity of RNAs were examined using an Agilent 2100 Bioanalyzer. After taking equal amounts of controls and six-stage total RNAs (4°C for 1, 4, 8, 12, 24, and 48 h), these samples were pooled together for the construction of a cold-stressed small-RNA library (CL).
Total RNA was purified by electrophoretic separation on a 15% Tris-Borate-EDTA (TBE)-urea denaturing polyacrylamide gel, and small RNA regions (15 to 30 nucleotide bands) were excised and recovered. The small purified RNAs were then ligated with 5′ adaptors (Illumina, San Diego, CA, USA), and the ligation products (40 to 60 bases in length) were purified on a Novex 15% TBE-Urea gel to remove unligated adaptors. Subsequently, a 3′ adaptor (Illumina) was ligated to the 5′ ligation products. After gel purification on a Novex 10% TBE-Urea gel (Invitrogen), RNA fragments with adaptors at both ends (70 to 90 bases in length) were reverse transcribed, and the resulting cDNAs were subjected to 15 PCR cycles using the adaptor primers. The resulting small RNA library was then sequenced by Genome Analyzer GA-I (Illumina, San Diego, CA, USA).
Identification of conserved and novel miRNAs in tea plant
Clean full-length reads were obtained from raw sequences after removing all low-quality reads, contaminants, and reads smaller than 15 nt. The extracted clean reads were used to calculate length distribution, and mapped to the P. trichocarpa genome sequence, and publicly available tea EST sequences using SOAP . Sequence matching non-coding RNAs, including rRNA, tRNA, snRNA, and snoRNA, in the NCBI GenBank (http://www.ncbi.nlm.nih.gov/genbank/), repeats (http://www.girinst.org/repbase/update/index.html) and Rfam (http://www.sanger.ac.uk/software/Rfam) databases were eliminated. Unique miRNA sequences were aligned with known miRNA sequences from miRBase 19.0 (http://www.mirbase.org/) with a maximum of two mismatches. The obtained sequences were used to predict secondary structures by RNAfold software (http://www.tbi.univie.ac.at/RNA/) to identify whether these precursors form a typical miRNA stem-loop structure. The remaining unknown sRNAs were analyzed, and novel miRNAs were predicted using Mireap (http://sourceforge.net/projects/mireap/files/latest/download). The novel miRNAs were screened according to the recently published criteria of novel miRNAs . The novel miRNAs and their corresponding miRNA*s formed a stem-loop duplex with less than four mismatched bases, derived from opposite stem-arms. The novel miRNAs were typically named using a ‘miR’ prefix to denote miRNAs, and a three-letter prefix to denote the species (e.g. `csn’ representing Camellia. sinensis).
miRNA microarray chip content and hybridization of arrays
The miRNA chip included 3228 miRNA probes corresponding to the miRNA transcripts listed in Sanger miRBase release 19.0, and 283 novel miRNAs probes were found in tea plant. Microarray assay was performed using a service provider (LC Sciences, Houston, USA). Total RNAs was isolated separately from leaf samples of 4 h, 12 h, and 24 h of cold-treated and non-treated materials. Total RNA samples (4 μg to 8 μg) were then 3′-extended with a poly (A) tail using poly (A) polymerase. An oligonucleotide tag was ligated to the poly (A) tail for fluorescent dye staining. Hybridization was performed overnight on a μParaflo microfluidic chip using a micro-circulation pump (Atactic Technologies) ,. On the microfluidic chip, each detection probe consisted of a chemically modified nucleotide-coding segment complementary to the target miRNA from miRBase or other RNAs, and a spacer segment of polyethylene glycol to extend the coding segment away from the substrate. The detection probes were made by in situ synthesis using photogenerated reagent chemistry. The hybridization melting temperatures were balanced by chemical modifications of the detection probes. Hybridization was performed using 100 μL of 6xSSPE buffer (0.90 M NaCl, 60 MM Na2HPO4, 6 MM EDTA, pH 6.8) containing 25% formamide at 34°C. After RNA hybridization, tag-conjugating Cy3 dye was circulated through the microfluidic chip for staining. Fluorescence images were collected using a laser scanner (GenePix 4000B, Molecular Device), and digitized using Array-Pro image analysis software (Media Cybernetics).
Data were analyzed by first substracting the background and then normalizing the signals using a locally-weighted regression filter . An miRNA signal was accepted as detectable if it met two conditions: signal intensity higher than three times the background standard deviation, and spot CV <0.5 (CV = signal standard deviation/signal intensity). Signals from four technical replicates of each RNA derived from stressed and control plants were compared using paired and two-tailed Student’s t test. Only signals with p <0.01 were considered as significant. Clustering analysis was performed using Cluster 3.0, and the heat map was visualized using Heatmap builder and TreeView .
qRT-PCR analysis of miRNA expression
The sRNAs (<200 nt from samples treated with cold conditions (4°C) for 4, 12, and 24 h were isolated using an mirVana miRNA isolation kit (Ambion, Austin, TX, USA). Stem-loop reverse transcription PCR (RT-PCR), which can be applied in small RNA cloning and multiple assays for better specificity and efficiency, was used for reverse transcription . A stem-looped RT-PCR Primer and qRT-PCR primers (Additional file 12: Table S7) were designed. One Step Prime-Script miRNA cDNA Synthesis Kit (TAKARA) and Platinum SYBR Green qPCR SuperMix-UDG (Invitrogen) were used. The 5.8S rRNA was used as the internal control . Thereafter, qRT-PCR was performed using a total reaction volume of 20 μL, which contained 0.5 μL of diluted cDNA, 8 μM primer mix, 10.0 μL of 2× SYBR Green Mix, and 8.7 μL ddH2O. The reactions were carried out in an iCycler iQ qRT-PCR detection system (BIO-RAD) with the following amplification conditions: activation at 50°C for 2 Min; 95°C for 2 Min; followed by 39 cycles at 95°C for 15 s, and 60°C for 30 s; and final holding at 4°C. The cycle threshold (CT) was defined as the fractional cycle number at which the fluorescence passes the fixed threshold. All reactions were performed in triplicate. The qRT-PCR results were normalized as follow: sample CT values were determined and standardized based on the 5.8S gene control prime reaction, and the 2-ΔΔCT method was applied to calculate the relative changes in gene expression from qRT-PCR experiments .
Degradome library construction and target identification
To obtain the potential target mRNAs, two degradome libraries were constructed from tea leaves treated with cold (+C) and cold-free (−C) as previously described ,. In brief, poly(A)-enriched RNA molecules were isolated and ligated to an RNA oligonucleotide adaptor containing a 3′ MmeI recognition site, and the ligated products were used to generate first-strand cDNA by RT. A short PCR was then used to amplify the cDNA. To obtain sufficient quantities of DNA, the product was ligated to a double-stranded DNA adaptor, and subjected to gel purification again for PCR amplification. The final cDNA library was purified and sequenced on Illumina GAIIx following the manufacturer’s instruction.
Raw sequence reads were obtained using Illumina Pipeline v1.5 software to remove adaptor sequences and low quality sequencing reads. The extracted sequencing reads with lengths of 20 nt and 21 nt were then used to identify potentially cleaved targets by CleaveLand pipeline, as previously described ,. The degradome reads were mapped to A. thaliana and tea plant sequences of mRNA and ESTs were downloaded from NCBI (http://www.ncbi.nlm.nih.gov/). Only the perfectly matched alignments for the given read were kept and extended to 35 nt to 36 nt by adding 15 nt of upstream sequences. All resulting reads (t-signature) were reverse-complemented and aligned to the miRNA identified in our study. No more than five mismatches of the alignments were allowed. Alignments in which the 5′ degradome sequence position coincided with the 10th nucleotide of miRNA were retained and scored using a previously described method . The target was selected and categorized based on cellular component, molecular function, and biological process. All identified targets were subjected to BlastX analysis to search for similarity. GO analysis was used to uncover the miRNA-gene regulatory network on the biological process and molecular function, as previously described by Xie et al. .
Modified 5′ RNA ligase-mediated RACE for the mapping of mRNA cleavage sites
Total RNA was extracted from young leaves of different treated and CK using Trizol reagent. Poly(A)+ mRNA was purified from all kinds of pooled sample RNA using the PolyA kit (Promega, Madison, WI), according to manufacturer’s instructions. A modified procedure for RLM-5′ RACE was followed with the GeneRacer Kit (Invitrogen, CA), as described previously . The PCR amplifications were performed using the GenRacer 5′ primer and the gene-specific primer (Table 2). Nested PCR amplifications were performed using the GeneRacer 5′ nested primer and the nested gene-specific nested primer (Table 2). PCR reactions were separated by agarose gel electrophoresis, and distinct bands of the appropriate size for miRNA-mediated cleavage were purified (excised gel slices corresponded to a size range of ~200 base pairs), cloned, and sequenced.
The sRNA sequence data from this study have been submitted to Gene Expression Omnibus (GEO) under accession NO. GSE61719 at website: http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=uvmvkkcixrepvyv&acc=GSE61719.
XL YZ WF designed and conducted the experiments, YZ XZ XC CS ZZ YW MW WF analyzed the date, YZ XZ drafted the manuscript, XL supervised the project. All authors have read and approved the final version of this manuscript.
- C. sinensis :
RNA-induced silencing complex
Quantitative real-time PCR
Small nucleolar RNA
Small nuclear RNA
Minimum free energy
Minimum free energy index
- A. thaliana :
Late embryogenesis abundant
- RLM-5′ RACE:
5′ RNA ligase-mediated RACE
Expressed sequence tag
Nicotinamide adenine dinucleotide 2′-phosphate reduced tetrasodium salt
Auxin/ indole-3-acetic acid
Reverse transcription PCR
Zhang Z, Yu J, Li D, Zhang Z, Liu F, Zhou X, Wang T, Ling Y, Su Z: PMRD: plant microRNA database. Nucleic Acids Res. 2010, 38 (suppl 1): D806-D813. 10.1093/nar/gkp818.
Kurihara Y, Watanabe Y: Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci U S A. 2004, 101 (34): 12753-12758. 10.1073/pnas.0403115101.
Baumberger N, Baulcombe D: Arabidopsis ARGONAUTE1 is an RNA slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci U S A. 2005, 102 (33): 11928-11933. 10.1073/pnas.0505461102.
Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297. 10.1016/S0092-8674(04)00045-5.
Sunkar R, Chinnusamy V, Zhu J, Zhu J-K: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12 (7): 301-309. 10.1016/j.tplants.2007.05.001.
Li B, Qin Y, Duan H, Yin W, Xia X: Genome-wide characterization of new and drought stress responsive microRNAs in Populus euphratica. J Exp Bot. 2011, 62 (11): 3765-3779. 10.1093/jxb/err051.
Mahajan S, Tuteja N: Cold, salinity and drought stresses: an overview. Arch Biochem Biophys. 2005, 444 (2): 139-158. 10.1016/j.abb.2005.10.018.
Sunkar R, Zhu J-K: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Sci Signal. 2004, 16 (8): 2001-
B-h L, Henderson DA, Zhu J-K: The Arabidopsis cold-responsive transcriptome and its regulation by ICE1. The Plant Cell Online. 2005, 17 (11): 3155-3175. 10.1105/tpc.105.035568.
Zhou X, Wang G, Sutoh K, Zhu J-K, Zhang W: Identification of cold-inducible microRNAs in plants by transcriptome analysis. BBA-Gene Regul Mech. 2008, 1779 (11): 780-788.
Liu H-H, Tian X, Li Y-J, Wu C-A, Zheng C-C: Microarray-based analysis of stress-regulated microRNAs in Arabidopsis thaliana. RNA. 2008, 14 (5): 836-843. 10.1261/rna.895308.
Zhang J, Xu Y, Huan Q, Chong K: Deep sequencing of Brachypodium small RNAs at the global genome level identifies microRNAs involved in cold stress response. BMC Genomics. 2009, 10 (1): 449-10.1186/1471-2164-10-449.
Chinnusamy V, Zhu J-K, Sunkar R: Gene regulation during cold stress acclimation in plants. Plant Stress Tolerance. 2010, Springer, Berlin, 39-55. 10.1007/978-1-60761-702-0_3.
Lv D-K, Bai X, Li Y, Ding X-D, Ge Y, Cai H, Ji W, Wu N, Zhu Y-M: Profiling of cold-stress-responsive miRNAs in rice by microarrays. Gene. 2010, 459 (1): 39-47. 10.1016/j.gene.2010.03.011.
Chen L, Zhang Y, Ren Y, Xu J, Zhang Z, Wang Y: Genome-wide identification of cold-responsive and new microRNAs in Populus tomentosa by high-throughput sequencing. Biochem Biophys Res Commun. 2012, 417 (2): 892-896. 10.1016/j.bbrc.2011.12.070.
Prabu G, Mandal A: Computational identification of miRNAs and their target genes from expressed sequence tags of tea (Camellia sinensis). Genom Proteomics Bioinformatics. 2010, 8 (2): 113-10.1016/S1672-0229(10)60012-5.
Das A, Mondal TK: Computational identification of conserved microRNAs and their targets in tea (Camellia sinensis). Am J Plant Sci. 2010, 1: 77-86. 10.4236/ajps.2010.12010.
Mohanpuria P, Yadav SK: Characterization of novel small RNAs from tea (Camellia sinensis L.). Mol Biol Rep. 2012, 39 (4): 3977-3986. 10.1007/s11033-011-1178-3.
Ding D, Zhang L, Wang H, Liu Z, Zhang Z, Zheng Y: Differential expression of miRNAs in response to salt stress in maize roots. Ann Bot. 2009, 103 (1): 29-38. 10.1093/aob/mcn205.
Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ: Criteria for annotation of plant microRNAs. The Plant Cell Online. 2008, 20 (12): 3186-3190. 10.1105/tpc.108.064311.
ZHU Q-w, LUO Y-p: Identification of microRNAs and their targets in tea (Camellia sinensis). J Zhejiang Univ-SC B 2013:DOI:10.1631/jzus.B 1300006,
Anburaj J, Viswanathan C, Prabu G: Differential expression of microRNAs in dormant bud of tea [ Camellia sinensis (L.) O. Kuntze]. Plant Cell Rep 2014, DOI 10.1007/s00299-014-1589-4,
Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.
Zhang B, Pan X, Cox S, Cobb G, Anderson T: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci. 2006, 63 (2): 246-254. 10.1007/s00018-005-5467-7.
Mao W, Li Z, Xia X, Li Y, Yu J: A combined approach of high-throughput sequencing and degradome analysis reveals tissue specific expression of microRNAs and their targets in cucumber. PLoS One. 2012, 7 (3): e33040-10.1371/journal.pone.0033040.
Wang F, Li L, Liu L, Li H, Zhang Y, Yao Y, Ni Z, Gao J: High-throughput sequencing discovery of conserved and novel microRNAs in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Mol Genet Genomics. 2012, 287 (7): 555-563. 10.1007/s00438-012-0699-3.
Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of miRNA genes. PLoS One. 2007, 2 (2): e219-10.1371/journal.pone.0000219.
Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu l, Li S, Zhou H, Long C, Chen S, Hannon GJ, Qi Y: Sorting of small RNAs into Arabidopsis Argonaute complexes is directed by the 5′ terminal nucleotide. Cell. 2008, 133: 116-127. 10.1016/j.cell.2008.02.034.
Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ: Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008, 18 (10): 758-762. 10.1016/j.cub.2008.04.042.
Gao X, Gulari E, Zhou X: In situ synthesis of oligonucleotide microarrays. Biopolymers. 2004, 73 (5): 579-596. 10.1002/bip.20005.
Addo-Quaye C, Miller W, Axtell MJ: CleaveLand: a pipeline for using degradome data to find cleaved small RNA targets. Bioinformatics. 2009, 25 (1): 130-131. 10.1093/bioinformatics/btn604.
Huang J, Sun S, Xu D, Lan H, Sun H, Wang Z, Bao Y, Wang J, Tang H, Zhang H: A TFIIIA-type zinc finger protein confers multiple abiotic stress tolerances in transgenic rice (Oryza sativa L.). Plant Mol Biol. 2012, 80 (3): 337-350. 10.1007/s11103-012-9955-5.
Ma K, Xiao J, Li X, Zhang Q, Lian X: Sequence and expression analysis of the C3HC4-type RING finger gene family in rice. Gene. 2009, 444 (1): 33-45. 10.1016/j.gene.2009.05.018.
Bies-Etheve N, Gaubier-Comella P, Debures A, Lasserre E, Jobet E, Raynal M, Cooke R, Delseny M: Inventory, evolution and expression profiling diversity of the LEA (late embryogenesis abundant) protein gene family in Arabidopsis thaliana. Plant Mol Biol. 2008, 67 (1-2): 107-124. 10.1007/s11103-008-9304-x.
Shimizu-Sato S, Mori H: Control of outgrowth and dormancy in axillary buds. Plant Physiol. 2001, 127 (4): 1405-1413. 10.1104/pp.010841.
Andreasson E, Jenkins T, Brodersen P, Thorgrimsen S, Petersen NH, Zhu S, Qiu J-L, Micheelsen P, Rocher A, Petersen M: The MAP kinase substrate MKS1 is a regulator of plant defense responses. EMBO J. 2005, 24 (14): 2579-2589. 10.1038/sj.emboj.7600737.
Morikawa K, Shiina T, Murakami S, Toyoshima Y: Novel nuclear-encoded proteins interacting with a plastid sigma factor, Sig1, in Arabidopsis thaliana. FEBS Lett. 2002, 514 (2): 300-304. 10.1016/S0014-5793(02)02388-8.
Hutvagner G, Zamore PD: A microRNA in a multiple-turnover RNAi enzyme complex. Science. 2002, 297 (5589): 2056-2060. 10.1126/science.1073827.
Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science. 2002, 297 (2289): 2053-2056. 10.1126/science.1076311.
Tang G, Reinhart BJ, Bater DP, Zamore PD: A biochemical framework for RNA silencing in plants. Genes Dev. 2003, 17: 49-63. 10.1101/gad.1048103.
Yekta S, Shih I-h, Bater DP: MicroRNA-directed cleavage of HOXB8 mRNA. Science. 2004, 304 (5670): 594-596. 10.1126/science.1097434.
Kasschau KD, Xie Z, Allen E, Llave C, Chapman EJ, Krizan KA, Carrington JC: P1/HC-Pro, a viral suppressor of RNA silencing, interferes with Arabidopsis development and miRNA unction. Dev Cell. 2003, 4 (2): 205-217. 10.1016/S1534-5807(03)00025-X.
Lewis BP, Burge CB, Batel DP: Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005, 120 (1): 15-20. 10.1016/j.cell.2004.12.035.
Elbashir SM, Lendechel W, Tuschl T: RNA interference is mediated by 21-and 22-nucleotide RNAs. Genes Dev. 2001, 15: 188-200. 10.1101/gad.862301.
Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu J-K: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8 (1): 25-10.1186/1471-2229-8-25.
Pantaleo V, Szittya G, Moxon S, Miozzi L, Moulton V, Dalmay T, Burgyan J: Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis. Plant J. 2010, 62 (6): 960-976.
Chen L, Ren Y, Zhang Y, Xu J, Sun F, Zhang Z, Wang Y: Genome-wide identification and expression analysis of heat-responsive and novel microRNAs in Populus tomentosa. Gene. 2012, 504 (2): 160-165. 10.1016/j.gene.2012.05.034.
Song C, Wang C, Zhang C, Korir N, Yu H, Ma Z, Fang J: Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics. 2010, 11 (1): 431-10.1186/1471-2164-11-431.
Wang T, Chen L, Zhao M, Tian Q, Zhang W-H: Identification of drought-responsive microRNAs in Medicago truncatula by genome-wide high-throughput sequencing. BMC Genomics. 2011, 12 (1): 367-10.1186/1471-2164-12-367.
Xu M-J, Liu Q, Nisbet A, Cai X-Q, Yan C, Lin R-Q, Yuan Z-G, Song H-Q, He X-H, Zhu X-Q: Identification and characterization of microRNAs in Clonorchis sinensis of human health significance. BMC Genomics. 2010, 11 (1): 521-10.1186/1471-2164-11-521.
Shuai P, Liang D, Zhang Z, Yin W, Xia X: Identification of drought-responsive and novel Populus trichocarpa microRNAs by high-throughput sequencing and their targets using degradome analysis. BMC Genomics. 2013, 14 (1): 233-10.1186/1471-2164-14-233.
Wang L, Liu H, Li D, Chen H: Identification and characterization of maize microRNAs involved in the very early stage of seed germination. BMC Genomics. 2011, 12 (1): 154-10.1186/1471-2164-12-154.
Gao Z, Shi T, Luo X, Zhang Z, Zhuang W, Wang L: High-throughput sequencing of small RNAs and analysis of differentially expressed microRNAs associated with pistil development in Japanese apricot. BMC Genomics. 2012, 13: 371-10.1186/1471-2164-13-371.
Ng KLS, Mishra SK: De novo SVM classification of precursor microRNAs from genomic pseudo hairpins using global and intrinsic folding measures. Bioinformatics. 2007, 23 (11): 1321-1330. 10.1093/bioinformatics/btm026.
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes Dev. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.
Cheng Y, Long M: A cytosolic NADP-malic enzyme gene from rice (Oryza sativa L.) confers salt tolerance in transgenic Arabidopsis. Biotechnol Lett. 2007, 29 (7): 1129-1134. 10.1007/s10529-007-9347-0.
Vaucheret H, Mallory AC, Bartel DP: AGO1 homeostasis entails coexpression of MIR168 and AGO1 and preferential stabilization of miR168 by AGO1. Mol Cell. 2006, 22 (1): 129-136. 10.1016/j.molcel.2006.03.011.
German MA, Pillay M, Jeong D-H, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008, 26 (8): 941-946. 10.1038/nbt1417.
Zeba N, Isbat M, Kwon N-J, Lee MO, Kim SR, Hong CB: Heat-inducible C3HC4 type RING zinc finger protein gene from Capsicum annuum enhances growth of transgenic tobacco. Planta. 2009, 229 (4): 861-871. 10.1007/s00425-008-0884-0.
Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O: Widespread translational inhibition by plant miRNAs and siRNAs. Science. 2008, 320 (5880): 1185-1190. 10.1126/science.1159151.
Lanet E, Delannoy E, Sormani R, Floris M, Brodersen P, Crété P, Voinnet O, Robaglia C: Biochemical evidence for translational repression by Arabidopsis microRNAs. The Plant Cell Online. 2009, 21 (6): 1762-1768. 10.1105/tpc.108.063412.
Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14 (6): 787-799. 10.1016/j.molcel.2004.05.027.
Fattash I, Voß B, Reski R, Hess WR, Frank W: Evidence for the rapid expansion of microRNA-mediated regulation in early land plant evolution. BMC Plant Biol. 2007, 7 (1): 13-10.1186/1471-2229-7-13.
Li YF, Zheng Y, Addo-Quaye C, Zhang L, Saini A, Jagadeeswaran G, Axtell MJ, Zhang W, Sunkar R: Transcriptome-wide identification of microRNA targets in rice. Plant J. 2010, 62 (5): 742-759. 10.1111/j.1365-313X.2010.04187.x.
Tunnacliffe A, Hincha DK, Leprince O, Macherel D: LEA proteins: versatility of form and function. Dormancy and Resistance in Harsh Environments. 2010, Springer, Berlin, 91-108. 10.1007/978-3-642-12422-8_6.
Sanchez-Ballesta MT, Rodrigo MJ, Lafuente MT, Granell A, Zacarias L: Dehydrin from citrus, which confers in vitro dehydration and freezing protection activity, is constitutive and highly expressed in the flavedo of fruit but responsive to cold and water stress in leaves. J Agric Food Chem. 2004, 52 (7): 1950-1957. 10.1021/jf035216+.
Goyal K, Walton L, Tunnacliffe A: LEA proteins prevent protein aggregation due to water stress. Biochem J. 2005, 388: 151-157. 10.1042/BJ20041931.
Grelet J, Benamar A, Teyssier E, Avelange-Macherel M-H, Grunwald D, Macherel D: Identification in pea seed mitochondria of a late-embryogenesis abundant protein able to protect enzymes from drying. Plant Physiol. 2005, 137 (1): 157-167. 10.1104/pp.104.052480.
Reyes JL, Campos F, Wei H, Arora R, Yang Y, Karlson DT, Covarrubias AA: Functional dissection of hydrophilins during in vitro freeze protection. Plant Cell Environ. 2008, 31 (12): 1781-1790. 10.1111/j.1365-3040.2008.01879.x.
Nakayama K, Okawa K, Kakizaki T, Inaba T: Evaluation of the protective activities of a late embryogenesis abundant (LEA) related protein, Cor15am, during various stresses in vitro. Bioscience, Biotechnol, and Bioc. 2008, (6): 1642-1642. 10.1271/bbb.80214.
Tan Y-F, O’Toole N, Taylor NL, Millar AH: Divalent metal ions in plant mitochondria and their role in interactions with proteins and oxidative stress-induced damage to respiratory function. Plant Physiol. 2010, 152 (2): 747-761. 10.1104/pp.109.147942.
Kalifa Y, Gilad A, Konrad Z, Zaccai M, Scolnik P, Bar-Zvi D: The water-and salt-stress-regulated Asr1 (abscisic acid stress ripening) gene encodes a zinc-dependent DNA-binding protein. Biochem J. 2004, 381: 373-378. 10.1042/BJ20031800.
Goldgur Y, Rom S, Ghirlando R, Shkolnik D, Shadrin N, Konrad Z, Bar-Zvi D: Desiccation and zinc binding induce transition of tomato abscisic acid stress ripening 1, a water stress-and salt stress-regulated plant-specific protein, from unfolded to folded state. Plant Physiol. 2007, 143 (2): 617-628. 10.1104/pp.106.092965.
Lu P-J, Zhou XZ, Shen M, Lu KP: Function of WW domains as phosphoserine-or phosphothreonine-binding modules. Science. 1999, 283 (5406): 1325-1328. 10.1126/science.283.5406.1325.
Jadwin JA, Ogiue-Ikeda M, Machida K: The application of modular protein domains in proteomics. FEBS Lett. 2012, 586 (17): 2586-2596. 10.1016/j.febslet.2012.04.019.
Cheng Y, Kato N, Wang W, Li J, Chen X: Two RNA binding proteins, HEN4 and HUA1, act in the processing of AGAMOUS Pre-mRNA in Arabidopsis thaliana. Dev Cell. 2003, 4 (1): 53-66. 10.1016/S1534-5807(02)00399-4.
Mockler TC, Yu X, Shalitin D, Parikh D, Michael TP, Liou J, Huang J, Smith Z, Alonso JM, Ecker JR: Regulation of flowering time in Arabidopsis by K homology domain proteins. Proc Natl Acad Sci U S A. 2004, 101 (34): 12759-12764. 10.1073/pnas.0404552101.
José Ripoll J, Ferrándiz C, Martínez-Laborda A, Vera A: PEPPER, a novel K-homology domain gene, regulates vegetative and gynoecium development in Arabidopsis. Dev Biol. 2006, 289 (2): 346-359. 10.1016/j.ydbio.2005.10.037.
Phelps-Durr TL, Thomas J, Vahab P, Timmermans MC: Maize rough sheath2 and its Arabidopsis orthologue ASYMMETRIC LEAVES1 interact with HIRA, a predicted histone chaperone, to maintain knox gene silencing and determinacy during organogenesis. The Plant Cell Online. 2005, 17 (11): 2886-2898. 10.1105/tpc.105.035477.
Chen X: A MicroRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303 (5666): 2022-2025. 10.1126/science.1088060.
Sunkar R, Girke T, Jain PK, Zhu J: Cloning and characterization of microRNAs from rice. Plant Cell. 2005, 17 (5): 1397-1411. 10.1105/tpc.105.031682.
Hafner M, Landgraf P, Ludwig J, Rice A, Ojo T, Lin C, Holoch D, Lim C, Tuschl T: Identification of microRNAs and other small regulatory RNAs using cDNA library sequencing. Methods. 2008, 44 (1): 3-12. 10.1016/j.ymeth.2007.09.009.
Zhu Q, Hong A, Sheng N, Zhang X, Matejko A, Jun K-Y, Srivannavit O, Gulari E, Gao X, Zhou X: μParaflo™ Biochip for Nucleic Acid and Protein Analysis. Microarrays. Humana Press, New Jersey; 2007:287-312. 10.1007/978-1-59745-304-2_19.
Bolstad BM, Irizarry RA, strand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33 (20): e179-e179. 10.1093/nar/gni178.
Lu S, Sun YH, Chiang VL: Stress-responsive microRNAs in Populus. Plant J. 2008, 55 (1): 131-151. 10.1111/j.1365-313X.2008.03497.x.
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.
Allen E, Xie Z, Gustafson AM, Carrington JC: microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell. 2005, 121 (2): 207-10.1016/j.cell.2005.04.004.
Xie F, Frazier TP, Zhang B: Identification, characterization and expression analysis of MicroRNAs and their targets in the potato (Solanum tuberosum). Gene. 2011, 473 (1): 8-22. 10.1016/j.gene.2010.09.007.
Song C, Fang J, Li X, Liu H, Thomas CC: Identification and characterization of 27 conserved microRNAs in citrus. Planta. 2009, 230: 671-685. 10.1007/s00425-009-0971-x.
This study was supported by the National Natural Science Foundation of China (31370688,31470690,31370014), the earmarked fund for Modern Agro-industry Technology Research System(CARS-23), the Science and Technology Plan of Jiangsu Province (BE2009313-1, BE2011319), and the Science and Technology Plan of Suzhou (SZGD201067, WNZ1002).
The authors declare that they have no competing interests.
Electronic supplementary material
Additional file 2: Figure S1.: Mature and precursor sequences and the predicted stem-loop structures of newly identified miRNAs from C. sinensis.(PDF 797 KB)
Additional file 3: Table S2.: Novel miRNAs and potentially novel miRNAs identified from predicted RNA hairpins in C. sinensis with cold stress.(PDF 114 KB)
Additional file 4: Figure S2.: t-plot for targets of the known miRNAs found in + C and -C libraries of C. sinensis. Signature abundance throughout the length of the transcript is show. Arrows indicate signature consistent with miRNA-directed cleavage. (TIFF 1 MB)
Additional file 5: Figure S3.: t-plot for targets of the new miRNA candidates found in + C and -C libraries of C. sinensis. Signature abundance throughout the length of the transcript is show. Arrows indicate signature consistent with miRNA-directed cleavage. (TIFF 1 MB)
Additional file 6: Table S3.: Relative expression analysis of cold-responsive miRNAs in ‘Yingshuang’. The threshold of p value was set 0.01, ***indicated p value less than 0.0001, **indicated p value from 0.0001 to 0.001, *indicated p value from 0.001 to 0.01, respectively. (PDF 167 KB)
Additional file 7: Table S4.: Relative expression analysis of cold-responsive miRNAs in `Baiye 1’. The threshold of p value was set 0.01, ***indicated p value less than 0.0001, ** indicated p value from 0.0001 to 0.001, *indicated p value from 0.001 to 0.01, respectively. (PDF 252 KB)
Additional file 8: Figure S4.: Microarray analysis of the known and new miRNAs from tea plant cultivar `Yingshaung’ treated with cold and cold-free. (PDF 8 MB)
Additional file 9: Figure S5.: Microarray analysis of the known and new miRNAs from tea plant cultivar `Baiye 1’ treated with cold and cold-free. (PDF 4 MB)
Additional file 10: Table S5.: Targets for conserved miRNAs from C. sinensis treated with cold (+C) and without cold (−C). (PDF 376 KB)
Additional file 11: Table S6.: Targets for new miRNA candidates from C. sinensis treated with cold (+C) and without cold (−C). (PDF 297 KB)
Additional file 12: Table S7.: Sequences of primers used for the reverse transcription and Quantitative real-time PCR experiments. (PDF 81 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zhang, Y., Zhu, X., Chen, X. et al. Identification and characterization of cold-responsive microRNAs in tea plant (Camellia sinensis) and their targets using high-throughput sequencing and degradome analysis. BMC Plant Biol 14, 271 (2014). https://doi.org/10.1186/s12870-014-0271-x