- Research article
- Open Access
Sculpting the maturation, softening and ethylene pathway: The influences of microRNAs on tomato fruits
BMC Genomicsvolume 13, Article number: 7 (2012)
MicroRNAs (miRNAs), a ubiquitous class of short RNAs, play vital roles in physiological and biochemical processes in plants by mediating gene silencing at post-transcriptional (PTGS) level. Tomato is a model system to study molecular basis of fleshy fruit ripening and senescence, ethylene biosynthesis and signal transduction owing to its genetic and molecular tractability. To study the functions of miRNAs in tomato fruit ripening and senescence, and their possible roles in ethylene response, the next generation sequencing method was employed to identify miRNAs in tomato fruit. Bioinformatics and molecular biology approaches were combined to profile the miRNAs expression patterns at three different fruit ripening stages and by exogenous ethylene treatment.
In addition to 7 novel miRNA families, 103 conserved miRNAs belonging to 24 families and 10 non-conserved miRNAs matching 9 families were identified in our libraries. The targets of many these miRNAs were predicted to be transcriptional factors. Other targets are known to play roles in the regulation of metabolic processes. Interestingly, some targets were predicted to be involved in fruit ripening and softening, such as Pectate Lyase, beta-galactosidase, while a few others were predicted to be involved in ethylene biosynthesis and signaling pathway, such as ACS, EIN2 and CTR1. The expression patterns of a number of such miRNAs at three ripening stages were confirmed by stem-loop RT-PCR, which showed a strong negative correlation with that of their targets. The regulation of exogenous ethylene on miRNAs expression profiles were analyzed simultaneously, and 3 down-regulated, 5 up-regulated miRNAs were found in this study.
A combination of high throughput sequencing and molecular biology approaches was used to explore the involvement of miRNAs during fruit ripening. Several miRNAs showed differential expression profiles during fruit ripening, and a number of miRNAs were influenced by ethylene treatment. The results suggest the importance of miRNAs in fruit ripening and ethylene response.
Gene expression is regulated by transcriptional and post-transcriptional pathways, which are crucial for optimizing gene output and for coordinating cellular programs . One of the recently discovered mechanisms in plants was short non-coding RNAs mediated gene silencing at post-transcriptional (PTGS) level . Short RNAs (sRNAs) are diverse and can be categorized into two major classes: short interfering RNAs (siRNAs) and microRNAs (miRNAs) [3, 4]. SiRNAs, processed from perfectly double-stranded RNA (dsRNA), posttranscriptionally silence transposons, viruses, and transgenes are important for DNA methylation [5–7]. MicroRNAs (miRNAs), a near ubiquitous class of short RNAs in plants, are orchestrated by DCL-like family from single-stranded RNA precursors that possess imperfect stem-loop foldback structures [8, 9]. Base pairing is used by mature miRNAs to guide RISCs to specific mRNAs bearing fully or partly complementary sequences . Repression of the target transcripts by miRNAs may occur through translational inhibition or slicing, with the two layers of regulation not necessarily coinciding spatially or temporally [11, 12]. Based on the sheer abundance and diversity of plant miRNAs, it is likely that most, if not all, physiological and biochemical processes in plants involve at some point the action of one or more miRNAs [9, 13, 14].
Tomato is an emblematic system to study molecular basis of fleshy fruit ripening and senescence, ethylene biosynthesis and signal transduction owing to its genetic and molecular tractability [15–17]. Recently, involvement of sRNAs in tomato fruit has been received attention. The homology search and molecular biology methods were used for characterizing novel and conserved miRNAs at the outset [18–21]. As deep sequencing technology has emerged and been employed extensively to identify miRNAs in model plants such as Arabidopsis and rice because of its high throughputs and accuracy, which make explore miRNAs in large scale possible [22–25]. The 454 pyrosequencing sequencing platform was first used to sequence tomato sRNAs from young leaves and a young green fruits of Microtom . Several conserved and non-conserved miRNAs were identified in this study, and miR156 was found to be involved in fruit ripening which raised the possibility that fruit ripening process may be under miRNA regulation 
Tomato fruit ripening and senescence are genetically regulated processes. Ripening of fleshy fruits involves evolution of ethylene, accumulation of pigments such as carotene and lycopene, development of aroma and flavor, softening of fruit tissues and increased susceptibility to pathogens [15, 16, 27]. To date, functional analysis has been carried out only for a few tomato miRNAs, most of which were validated to be involved in leaf and flower development [28–31], and a few of which were proved to be involved in the stress-response and host-pathogen interactions [32–35]. In order to study the functions of miRNAs in tomato fruit ripening and senescence, and their possible roles in ethylene pathway, the next generation sequencing method (Solexa platform) was employed to identify miRNAs in tomato fruit. Bioinformatics and molecular biology approaches were combined to profile the miRNAs expression patterns of fruits at three different ripening stages (mature green stage, breaker stage, red ripe stage). The influences of exogenous ethylene on the miRNAs expression levels were also explored.
High-throughput sequencing of short RNAs of three fruit ripening stages
To explore the possible regulatory roles of miRNAs in fleshy fruit ripening process, we select fruits at three different ripening periods (mature green stage, the breaker stage and red ripening stage) of the Solanum lycopersicum (Ailsa Craig). The sRNA fraction was subject to deep sequencing on the Illumina (Solexa)1G platform which produced 4,333,963, 5,514,197 and 5,334,158 raw reads, respectively. After removal of the adaptor sequences and filtering by sequence properties, 2,388,170, 3,599,583 and 2,274,050 redundant sRNA reads were remained for further analysis, respectively (Table 1).
Two distinguishing features of the small RNAs libraries are the populations and distributions. The composition of small RNAs often reflects roles of different categories of small RNAs in a particular tissue or species or associated biogenetic machines. Since most of the small RNAs with known functions are 20-24 nt long, we only investigated the distribution of small RNAs of 18-25 nt in the three libraries (Figure 1; Figure 2). In the three ripening stages, the 24-nt size class was the overall most abundant class of sRNAs (Figure 2), which was consistent with that of Arabidopsis  and rice , but different from that of wheat, Chinese yew and grapevine [38–40]. The higher percentage of 24-nt small RNAs in tomato fruit may reflect the complexity of the tomato genome because 24-nt siRNAs are known as predominant heterochromatin-associated sRNAs . The 21-24nt sRNAs exhibited a peak at breaker stage, and showed a sharp reduction at the red ripening stage, suggesting potential roles in ripening process.
Identification of known and novel miRNAs in tomato fruit
Conserved miRNAs were found in many plant species and have important functions in plant development and stress response . To identify the conserved miRNAs in tomato, sRNA sequences obtained by deep sequencing were compared with the currently known mature plant miRNAs in miRBase . After Blastn searches and further sequence analysis, a total of 103 conserved miRNAs (Figure 3), belonging to 24 miRNA families were identified and all of which can be consulted in two important database (http://www.mirbase.org/; http://ted.bti.cornell.edu/cgi-bin/TFGD/sRNA/miRNA.cgi).
In addition to the conserved miRNAs, there are known miRNAs that are not-conserved, but are found in only one or a few plant species such as Arabidopsis or Physcomitrella patens [8, 23]. Ten members belong to 9 non-conserved miRNAs (miR158, miR161, miR173, miR393, miR398, miR403, miR414, miR858, miR894) were present in our data sets at a low abundance, with the exception of miR403 that were represented by a relative high number (Table 2). All of the non-conserved miRNAs were validated by the construction of the small RNA cDNA library method (Figure 4) .
A number of criterions are used for evaluating whether a small RNA is a genuine miRNA, such as formation of a stable hairpin structure, starting with a 5' uridine, lower minimal free energies (MFEs) for hairpin structure of its precursors, and detection of miRNA*s [44, 45]. By using these rules, we identified 11 novel miRNAs belonging to 7 families (named miRZ1-miRZ7, Figure 5, Figure 6, additional file 1, 2). Except for miRZ6, most of the novel miRNAs were in low copies. The possible reason is that the encoding genes are not broadly conserved and expressed at much lower levels .
Prediction of the targets of identified miRNAs in tomato fruit
Target prediction for miRNAs is straightforward because it is assumed that most of them match their targets with almost perfect complementarity [9, 47]. The putative target genes for all identified conserved and non-conserved miRNAs were searched by using the web-based computer psRNA Target Server (http://plantgrn.noble.org/psRNATarget/) which identifies putative targets regulated at post-transcriptional or translational levels.
Most conserved miRNA targets that are conserved across several plant species, including Arabidopsis , rice , grape , poplar  and wheat . And most miRNA families have multiple target sites, suggesting that these miRNAs are functionally divergent. In our study, 152 targets of conserved and non-conserved miRNAs were predicted (Figure 7). The majority of such targets are various transcriptional factors including SBP (miR156), MYB (miR159, miR319, miR172), NAM (miR164), and MADS-Box (miR396) that regulate plant development  or phytohormone signal transduction . Other conserved miRNAs targets include F-box protein (miR394, miR414), ATP sulfurylase (miR395), Pectate Lyase (miR482), endo-1, 4-beta- glucanase (miR396), Laccase (miR397), all of which are involved in regulation of metabolic processes. Several targets of the miRNAs are AGO protein (miR168, miR403) which regulate their own biosynthesis pathway. In addition, targets of the conserved miRNAs include disease resistance proteins (miR472, miR482) which are related to pathogen resistance [33, 34]. Furthermore, ACS and EIN2, which are involved in ethylene biosynthesis and signal transduction , are putative targets of miR159 and miR828, respectively.
The targets of the novel miRNAs identified in our library were predicted, most of which were unknown (Additional file 1). MiRZ1 was predicted to be related to disease response, MiRZ6 was predicted to be involved in virus infection response. Interestingly, a target of miRZ7 is beta-galactosidase which is an important enzyme affects fruit softening. Another member's target of miRZ7 family is starch synthase which participates in starch biosynthesis, including transient starch .
Expression profiles of the miRNAs involved in tomato fruit ripening
It has been shown that high throughput sequencing provides an alternative way to gain the genome-wide transcript profiles of miRNAs abundance  and allows us to determine the abundance of various miRNA families and even to distinguish different members of a given family of one organism. To elucidate the potential roles of miRNAs in tomato fruit ripening process, we profiled the expression levels of known and novel miRNAs. We computed the normalised counts of known miRNA sequences, plotted them across the three ripening stages and validated some results by stem-loop RT-PCR (Figure 8).
The abundance of most known miRNAs decreased during fruit ripening, and some were found only in the mature green stage (Table 2). A few known miRNAs were expressed at high levels in fruit, such as miR157, miR162, miR164, miR166, miR168, miR172, miR396. On the other hand, other miRNAs, such as miR161, miR173, miR393, miR397, miR398 and miR414, were expressed at relatively lower levels and can only be detected in mature green stage, Particularly, the expression levels of miR159 decreased at the breaker stage and increased at red ripening stage compared with that of the mature green stage which reveals its probable specific roles in fruit ripening or ethylene pathway, as one of the predicted targets is 1-aminocyclopropane -1-carboxylate synthase which is a crucial enzyme in ethylene biosynthesis . Expression of another miRNA, miR396, increased from mature green to breaker stage, then dropped sharply in the red stage. Interestingly, targets of miR396 include MADS-Box protein and endo-1, 4-beta-glucanase, both of which are involved in fruit ripening and softening process . Moreover, 3 miRNAs (miR394, miR414, and miR482) predicted to have targets related to fruit ripening and softening and ethylene response were validated (Figure 8).
Colorless non-ripening (CNR), a member of the squamosa-promoter binding protein (SBP) family that was shown to be involved in fruit ripening , is targeted by miR156 [31, 39]. The expression levels of miR156 also decreased during fruit ripening (Table 2 and Figure 8). Two miRNAs, miR828 and miR1917, whose targets are ethylene-insensitive 2 (EIN2) and serine/threonine protein kinase (CTR1) which are vital regulators in ethylene signal transduction[17, 59], reported previously  and not shown in our sequencing results, showed reduced expression during fruit ripening (Figure 8). A novel miRNA (miRZ7) expression levels were also analyzed which was predicted to relate to fruit softening.
Overwhelming the majority of the stem-loop RT-PCR results was in agreement with sequencing results. Only two miRNAs (miR414, miRZ7) were exception which could be detected by stem-loop RT-PCR but were not detected or barely detectable by sequencing. This could be due to cross amplification of other highly homologous small RNA members or biases occurred during library generation or sequencing in some samples.
Expression profiles of the miRNAs responses to exogenous ethylene in tomato fruit
Ethylene is a regulator of various physiological and morphological responses, including seedling development, leaf and flower senescence, induction of fruit ripening, floral sex determination, resistance to pathogen infection, and adaptation to stress conditions [60–63]. Tomato has proven to be a highly successful model system for fruit development and ripening. And the role of ethylene during fruit ripening has been most thoroughly studied in tomato fruit [15, 16]. In order to find out the relationships between ethylene treatment and miRNAs expression, exogenous ethylene and 1-methylcycloproene were used to treat the mature green fruits. The expression levels of a number of known and novel miRNAs were tested by stem-loop RT-PCR. Three miRNAs (miR394, miR414, miR1917) were down regulated; five miRNAs (miR156, miR159, miR396, miR482, miRZ7) were up regulated, while one miRNA (miR828) whose target is EIN2 was not affected, by exogenous ethylene treatment (Figure 9).
High-throughput sequencing of short RNAs in tomato
High-throughput sequencing has been used to study miRNAs at the whole genome level in several model plant species, including Arabidopsis , rice , and wheat . The successful application of high-throughput sequencing technology to systemically identify plant miRNAs in recent years has greatly advanced our knowledge on the functions of miRNAs in plants.
The composition of small RNAs often reflects roles of different kinds of small RNAs in a specific tissue or species. Several plant species, such as Pinus cordata, Populus balsamifera and grapevine were shown to contain substantially more 21-nt than 24-nt sRNAs [37, 39, 65]. However, we observed an unusually high level of 24-nt sRNAs compared with the 21-nt class (Figure 1; Figure 2), as in several other plant miRNAs studies [23, 36, 37, 66]. The number of 24-nt small RNAs was almost three times that of the 21-nt class, such a high percentage of 24-nt small RNAs may reflect the complexity of the tomato genome because 24-nt siRNAs are known to be involved in heterochromatin modification, especially for genomes with high content of repetitive sequences [67, 68].
Plant MIR genes are possible to arise from gene duplication events which then evolved by random mutations into short, imperfectly paired hairpins [46, 69]. Those non-conserved miRNAs are believed to be evolutionarily recent and generally represented by single copy MIR genes . Recently, many non-conserved miRNAs were reported in several species [36, 37, 64, 66]. Some of these non-conserved miRNAs are not found in other species/families and several non-conserved miRNAs are found in different phylogenetic families. We also found nine non-conserved miRNAs in our database, although most of them were present at a very low level. However, the low number of sequence reads does not rule out that they could be expressed at a high level in specific cells. Few non-conserved miRNAs were expressed at a relatively high level, suggesting that they may represent an intermediate status between the deeply conserved and the less conserved miRNAs.
New species-specific miRNAs are considered to be young miRNAs that have evolved recently, and are often expressed at a lower level than conserved miRNAs, as was reported for Arabidopsis and wheat [23, 36, 38]. This observation is true for many new tomato miRNAs identified here. However, few new miRNAs were expressed at a high level in a tissue-specific manner (miRZ6). In some cases we observed inconsistency between the level of miRNAs identified by Solexa sequencing and experimental verification. It is possible that during library generation or sequencing some bias could occur for certain sequences in some samples.
In contrast to the 454 pyrosequencing sequencing in tomato previously , more reads were obtained in each of our three libraries, and miRNAs at different fruit ripening stages were additionally analyzed in our study. In our database, most conserved miRNAs have been detected, including 33 conserved and non-conserved miRNA families, which cover all the published miRNAs in the two databases for tomato miRNAs (http://www.mirbase.org/, http://ted.bti.cornell.edu/cgi-bin/TFGD/sRNA/miRNA.cgi).
Dalmay and colleagues studied the correlation between the short RNAs and fleshy fruit development using deep sequencing technology . However, mostly sRNAs were analyzed, which is different from our miRNAs study. But they reported several miRNAs, and validated one novel miRNA related to the glutamate accumulation which contributes to the tomato fruit taste.
Targets of known and novel miRNAs in tomato fruit
Target prediction is necessary for assessing miRNAs' putative functions. Currently, the most efficient tool available for this is the bioinformatics approach facilitated by the high degree of homology between miRNA and its target sequences in plants . Our analysis reveals that most of the predicted targets in tomato have a conserved function among a variety of plant species . Consistent with previous reports, the largest of these targets in tomato fruit were plant-specific transcription factors, such as AP2, NAC, SBP and the ARF family . The second largest targets encoded a range of different proteins implicated in a variety of metabolic processes, such as ATP sulfurylase, Pectate Lyase, endo-1, 4-beta-glucanase, Laccase. In addition, functions of several targets are largely unknown.
Although many conserved miRNA targets were predicted, only a limited number of miRNA targets were identified experimentally [31, 39]. Several miRNAs slice their targets and can be validated by RACE. Alternatively, other miRNAs inhibit target gene expression through translational arrest, such as miR156 and miR172, which have been shown to regulate their target genes (SBPs and AP2) predominantly by inhibiting their translation [73–76]. Indeed, genetic and biochemical evidence is accumulated that plant miRNA-guided silencing has a widespread translational inhibitory component [2, 77].
The novel tomato miRNAs target different genes with a wide variety of predicted functions. It may be worthy to note that those miRZ6 targets six members of disease resistance protein (CC-NBS-LRR class), which are likely to be involved in host-pathogen interactions. Surprisingly, one target of miRZ7 is beta-galactosidase which is a crucial enzyme for fruit softening, which is beneficial for clarify the cooperating roles in freshly fruit ripening and senescence. The target of another member in miRZ7 family is starch synthase which participates in regulating starch biosynthesis (Additional file 1) .
Fruit ripening and softening related miRNAs
Functional analysis has been carried out only for role of a few tomato miRNAs morphological development [28–30]. Conserved and non-conserved miRNAs regulate genes involved in fleshy fruit development were reported previously . A target gene of miR156 belong to squamosa-promoter binding protein (SBP) family called CNR was validated, which plays pivotal roles in fruit ripening . Over-expression of miR156 can cause the fruit red colour slightly lighter than the wild type, which was thought that the expression of miR156 and CNR overlaps only partially and the function of the miRNA is to suppress CNR expression in specific cell types. CNR mRNA is up-regulated in fruit at the breaker stage , coincidently, our sequence results showed a negative correlation as expected. Two important targets of miR396 are endo-1, 4-beta-glucanase and MADS-box protein, both of which are important regulators in fruit ripening and softening [78, 79]. Our sequencing results showed that the abundance of miR396 peaked at the breaker stage and decreased sharply in red ripening stage suggesting its possible function in fruit ripening. A target of miR482 is pectate lyase which is an important enzyme in fruit softening [80, 81]. Meanwhile, a novel miRNA, whose target is beta-galactosidase, which is also a key softening related enzyme [82, 83]. The expression levels of the two miRNAs were also coinciding with the assumption including the softening ripening stage (Figure 8).
Intriguingly, an in-depth study of the tomato transcriptome and proteome unraveled the regulation mechanisms during fruit ripening stages. The transcriptional regulation is mostly responsible for global reductions in plastid gene expression prior to or at the onset of fruit formation whereas post-transcriptional events become predominant at the breaker stage [84, 85]. These post-transcriptional events did not lead to significant reductions in mRNA levels, although mRNA stability was reported to be modestly affected for some genes . miRNAs and their targets have been shown to regulate in translational inhibition way in a fairly large fraction .
Ethylene biosynthesis and signal transduction related miRNAs in tomato
Recent researchers have witnessed the great progress in the area of identification of interaction between miRNAs pathways and phytohormone responses, which improves our understanding of mechanism of plant development controlled by miRNAs and hormone action to a large extent. One of the most important phytohormone is ethylene which plays vital roles during all stages of the plant life cycle, functioning through seed germination to ripening and various abiotic stress conditions . MiR159 and miR394 were reported to be associated with ethylene in rice . Most targets of miR159 are MYB, and a novel target not related to MYB in tomato recently . Another important target of miR159 is 1-aminocyclopropane-1-carboxylate synthase which plays vital roles in ethylene biosynthesis . As a typical climacteric plant, abundance ethylene was synthesized at the breaker stage approximately, and the expression level of miR159 was suppressed which match our results well. MiR394 and miR414 also target F-box family proteins which are important participants in the signal transduction pathways of different plant hormones [88, 89]. The stability of the ethylene signaling regulators EIN2 and EIN3 are modulated by the F-box proteins ETP1/2and F-box proteins EBF1/2 respectively. Ethylene Insensitive 2 (EIN2), a target of miR828, is the key positive regulator of ethylene signal transduction and other hormones such as abscisic acid, auxin, cytokinin and jasmonate and thus may represent a point of crosstalk between multiple hormone signaling pathways [90–92]. CTR1 is a negative regulator of ethylene response that likely interacts directly with receptor molecules to form a signaling complex [93, 94]. As with the ethylene receptors, all tissues evaluated express CTR1 genes and their mRNAs are differentially accumulated depending on tissue. LeCTR1 induction is associated with tissues at stages of development associated with increased ethylene, including fruit ripening and is a validated target of miR1917 . Exogenous ethylene can influence the miRNAs expression patterns, most of which shown negative correlation to their targets well except for miR828, was not affected by exogenous ethylene treatment (Figure 9). The most likely reason was that miR828 affected this target slightly and we cannot catch the discrepancy.
We have used a combination of high throughput sequencing and molecular biology approaches to explore the involvement of miRNAs during fruit ripening process. Total 33 conserved and non-conserved miRNA families and 7 novel miRNA families were indentified. The expression profiles of miRNAs at different ripening stages were analyzed and validated simultaneously, most of which were in agreement with the sequencing results. The influences of exogenous ethylene on miRNAs expression were also studied. These findings provide valuable information for further functional verification of miRNAs in tomato fruit ripening and ethylene response.
Tomato samples from three ripening stages (mature green, breaker, and red-ripening stage) of Solanum lycopersicum 'Ailsa Craig' were used to prepare for the high throughout sequencing. Total RNAs were extracted using TRIzol reagents. For each sample, the 18-30nt small RNAs were ligated with 5'- and 3'- RNA adapter by T4 RNA ligase after being purified by polyacrylamide gel electrophoresis. The RNAs were subsequently transcribed to single-stranded cDNAs using Superscript II reverse transcriptase. Thereafter the cDNAs were used as templates for double-stranded cDNA synthesis by PCR amplification using primers that anneal to adapters. The purified cDNAs were sequenced on an Illumina-Solexa 1 G Genetic Analyzer (BGI).
The adapter sequences of the Solexa sequencing results were removed. And sequences larger than 30nt and smaller than 18nt were discarded. All high quality sequences were considered as significant and further analyzed. Small RNA sequences were mapped to tomato genome (ftp://ftp.sgn.cornell.edu/tomato_genome/ITAG_devel_release/). rRNAs, tRNAs, snRNAs and snoRNAs were removed from the matched sequences through BLASTn search using NCBI Genebank database (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi/). Mismatches were not allowed in the above two approaches. The unique sequences left were aligned with known miRNAs from miRBase 10.0 (http://www.mirbase.org/) . The potential novel miRNAs were analyzed using mireap. Parameters were set based on a previous plant miRNAs study . The court number was normalized as transcript per million (TPM). Target predictions were performed using the psRNATarget (http://plantgrn.noble.org/psRNATarget/) .
Exogenous ethylene and 1-Methylcycloproene treatment
Solanum Lycopersicum (Ailsa Craig) plants were grown in soil under standard greenhouse conditions. Mature green fruits were harvested and divided into three groups, the first group was treated with 50 μl/L of ethylene for 6 h, the second group was treated with 0.5 μl/L 1- methylcyclopropene (1-MCP) for 24-h, and the third group was used as a control. After treatments, samples were frozen immediately in liquid nitrogen and stored at -80°C until RNA extraction .
Small RNA samples from above fruit samples were isolated using the miRNA isolation kit (Bioteke) according to the manufacturer's instructions. Small RNA samples were polyadenylated at 37°C for 60 min in a 50-μl reaction volume containing 0.5 μg RNA and 1.5 U poly (A) polymerase (Ambion). The reaction product was diluted to 300 μl. An equal volume of acid-phenol:chloroform was added, mixed and centrifuged. The aqueous phase was carefully removed and transferred to a new tube. Two volume of ethanol and one-tenth volume of sodium acetate (3 M, pH 5.2) were added and mixed thoroughly. The mixture was allowed to precipitate at -20°C for 2 h and then centrifuged at 12000 rpm for 20 min. The supernatant was discarded and the pellet was dissolved with 20 μl RNA-free water. To generate a small RNA cDNA (srcDNA) library, 20 μl of the tailed RNA and 1 μl of RTS primer were mixed in a 26-μl reaction volume, incubated at 65°C for 10 min, and annealed at 4°C for 20 min. Reverse transcription was carried out with reverse-transcriptase (Invitrogen) at 50°C for 60 min. Finally, the reverse transcriptase was inactivated by incubation at 70°C for 15 min. A small RNA-specific primer and a universal reverse primer were used for amplification of individual small RNAs (Additional file 3). The annealing temperature was adjusted according to the Tm of individual small RNAs. After PCR, an aliquot of the PCR products was analyzed on a 2.5% agarose gel .
Differential expression analysis of microRNAs
Total RNA was extracted from the samples using TRIzol solution (Trans, Beijing, China) and treated with RNase-free DNase I (Promaga, Beijing, China). First strand cDNA was synthesized using total RNA and reverse transcriptase (Promaga, Beijing, China). Expression levels of mature miRNAs were analyzed by Semi-quantitative RT-PCR stem-loop method [97–100]. A stem-loop containing RT primer with its 5'-end complementary to target miRNA's last 6-nt at 3'-end was designed. Reverse transcription was performed at 16°C for 30 minutes, followed by 60 cycles of pulsed RT at 30°C for 30 seconds, 42°C for 30 seconds and 50°C for 1 second. Semi-quantitative RT-PCR was performed using a forward primer containing the 5' part sequence of miRNA and a universal primer complementary to the stem-loop part of RT primer at 94°C for 2 min, followed by 21 cycles of 94°C for 15 s and 60°C for 1 min. The reaction products (about 60 bp) were analyzed by electrophoresis on a 2.5% agarose gel in 1× TAE. The primers used in this study were listed in the Additional file 4.
Mallory AC, Bouche N: MicroRNA-directed regulation: to cleave or not to cleave. Trends Plant Sci. 2008, 13 (7): 359-367. 10.1016/j.tplants.2008.03.007.
Brodersen P, Voinnet O: The diversity of RNA silencing pathways in plants. Trends Genet. 2006, 22 (5): 268-280. 10.1016/j.tig.2006.03.003.
Kim VN: Sorting out small RNAs. Cell. 2008, 133 (1): 25-26. 10.1016/j.cell.2008.03.015.
Carthew RW, Sontheimer EJ: Origins and Mechanisms of miRNAs and siRNAs. Cell. 2009, 136 (4): 642-55. 10.1016/j.cell.2009.01.035.
Baulcombe D: RNA silencing in plants. Nature. 2004, 431: 356-363. 10.1038/nature02874.
Matzke M, Kanno T, Huettel B, Daxinger L, Matzke AJ: Targets of RNA-directed DNA methylation. Curr Opin Plant Biol. 2007, 10: 512-519. 10.1016/j.pbi.2007.06.007.
Khraiwesh B, Arif MA, Seumel GI, Ossowski S, Weigel D, Reski R, Frank W: Transcriptional control of gene expression by expression by microRNAs. Cell. 2010, 140: 111-122. 10.1016/j.cell.2009.12.023.
Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.
Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.
Wu L, Belasco JG: Let me count the ways: mechanisms of gene regulation by miRNAs and siRNAs. Mol cell. 2008, 29 (1): 1-7. 10.1016/j.molcel.2007.12.010.
Llave C, Kasschau KD, Rector MA, Carrington JC: Endogenous and silencing-associated small RNAs in plants. Plant Cell. 2002, 14: 1605-1619. 10.1105/tpc.003210.
Eulalio A, Huntzinger E, Izaurralde E: Getting to the root of miRNA- mediated gene silencing. Cell. 2008, 132 (1): 9-14. 10.1016/j.cell.2007.12.024.
Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D: Control of leaf morphogenesis by microRNAs. Nature. 2003, 425: 257-263. 10.1038/nature01958.
Sunkar R, Zhu JK: Novel and stress-regulated microRNAs and other small RNAs from Arabidopsis. Plant Cell. 2004, 16: 2001-2019. 10.1105/tpc.104.022830.
Giovannoni JJ: Genetic regulation of fruit development and ripening. Plant Cell. 2004, 16: 170-180. 10.1105/tpc.019158.
Giovannoni JJ: Fruit ripening mutants yield insights into ripening control. Curr Opin Plant Biol. 2007, 10: 283-289. 10.1016/j.pbi.2007.04.008.
Cara B, Giovannoni JJ: The molecular biology of ethylene during tomato fruit development and maturation. Plant Science. 2008, 175: 106-113. 10.1016/j.plantsci.2008.03.021.
Pilcher RL, Moxon S, Pakseresht N, Moulton V, Manning K, Seymour G, Dalmay T: Identification of novel small RNAs in tomato (Solanum lycopersicum). Planta. 2007, 226: 709-717. 10.1007/s00425-007-0518-y.
Yin ZJ, Li CH, Han XL, Shen FF: Identification of conserved microRNAs and their target genes in tomato (Lycopersicon esculentum). Gene. 2008, 414: 60-66. 10.1016/j.gene.2008.02.007.
Zhang J, Zeng R, Chen J, Liu X, Liao Q: Identification of conserved microRNAs and their targets from Solanum lycopersicum Mill. Gene. 2008, 423: 1-7. 10.1016/j.gene.2008.05.023.
Itaya A, Bundschuh R, Archual AJ, Joung JG, Fei ZJ, Dai XB, Zhao PX, Tang YH, Nelson RS, Ding B: Small RNAs in tomato fruit and leaf development. Biochim Biophys Acta. 2008, 1779: 99-107.
Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome. Science. 2005, 309 (5740): 1567-1569. 10.1126/science.1114112.
Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-throughput sequencing of Arabidopsis microRNAs: Evidence for frequent birth and death of MIRNA Genes. PLoS ONE. 2007, 2 (2): 219-232. 10.1371/journal.pone.0000219.
Johnson C, Bowman L, Adai AT, Vance V, Sundaresan V: CSRDB: a small RNA integrated database and browser resource for cereals. Nucleic Acids Research. 2007, 35: 829-833. 10.1093/nar/gkl991.
Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biology. 2008, 8: 25-41. 10.1186/1471-2229-8-25.
Moxon SJR, Szittya G, Schwach F, Rusholme Pilcher RL, Moulton V, Dalmay T: Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res. 2008, 18 (10): 1602-1609. 10.1101/gr.080127.108.
White PJ: Recent advances in fruit development and ripening: an overview. J Exp Bot. 2002, 53: 1995-2000. 10.1093/jxb/erf105.
Ori N, Cohen AR, Etzioni A, Brand A, Yanai O, Shleizer S, Menda N, Amsellem Z, Efroni I, Pekker I, Alvarez JP, Blum E, Zamir D, Eshed Y: Regulation of LANCEOLATE by miR319 is required for compound-leaf development in tomato. Nat Genet. 2007, 39: 787-791. 10.1038/ng2036.
Berger Y, Harpaz-Saad S, Brand A, Melnik H, Sirding N, Alvarez JP, Zinder M, Samach A, Eshed Y, Ori N: The NAC-domain transcription factor GOBLET specifies leaflet boundaries in compound tomato leaves. Development. 2009, 136: 823-832. 10.1242/dev.031625.
Buxdorf K, Hendelman A, Stav R, Lapidot M, Ori N, Arazi T: Identification and characterization of a novel miR159 target not related to MYB in tomato. Planta. 2010, 232: 1009-1022. 10.1007/s00425-010-1231-9.
Zhang X, Zou Z, Zhang J, Zhang Y, Han Q, Hu T, Xu X, Liu H, Li H, Ye Z: Over-expression of sly-miR156a in tomato results in multiple vegetative and reproductive trait alterations and partial phenocopy of the sft mutant. FEBS Lett. 2011, 585: 435-439. 10.1016/j.febslet.2010.12.036.
Gu M, Xu K, Chen A, Zhu Y, Tang G, Xu G: Expression analysis suggests potential roles of microRNAs for phosphate and arbuscular mycorrhizal signaling in Solanum lycopersicum. Physiologia Plantarum. 2010, 138: 226-237. 10.1111/j.1399-3054.2009.01320.x.
Naqvi AR, Haq QM, Mukherjee SK: MicroRNA profiling of tomato leaf curl new delhi virus (tolcndv) infected tomato leaves indicates that deregulation of mir159/319 and mir172 might be linked with leaf curl disease. Virol J. 2010, 7: 281-296. 10.1186/1743-422X-7-281.
Stav R, Hendelman A, Buxdorf K, Arazi T: Transgenic expression of tomato bushy stunt virus silencing suppressor P19 via the pOp/LhG4 transactivation system induces viral-like symptoms in tomato. Virus Genes. 2010, 40: 19-129.
Zhang XH, Zou Z, Gong PJ, Zhang JH, Ziaf K, Li HX, Xiao FM, Ye ZB: Over-expression of microRNA169 confers enhanced drought tolerance to tomato. Biotechnol Lett. 2011, 33: 403-409. 10.1007/s10529-010-0436-0.
Rajagopalan R, Vaucheret H, Trejo J, Bartel DP: A diverse and evolutionarily fluid set of microRNAs in Arabidopsis thaliana. Genes & Development. 2006, 20 (24): 3407-3425. 10.1101/gad.1476406.
Morin RD, Aksay G, Dolgosheina E, Ebhardt HA, Magrini V, Mardis ER, Sahinalp SC, Unrau PJ: Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa. Genome Research. 2008, 18 (4): 571-584. 10.1101/gr.6897308.
Yao YY, Guo GG, Ni ZF, Sunkar R, Du JK, Zhu JK, Sun QX: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biology. 2007, 8 (6): 96-108. 10.1186/gb-2007-8-6-r96.
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.
Qiu D, Pan X, Wilson IW, Ketchum REB, Li F, Liu M, Teng W, Zhang BH: High throughput sequencing technology reveals that the taxoid elicitor methyl jasmonate regulates microRNA expression in Chinese yew (Taxus chinensis). Gene. 2009, 436 (1-2): 37-44. 10.1016/j.gene.2009.01.006.
Schwach F, Moxon S, Moulton V, Dalmay T: Deciphering the diversity of small RNAs in plants: the long and short of it. Brief Funct Genomic Proteomic. 2009, 8: 472-481. 10.1093/bfgp/elp024.
Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36: 154-158. 10.1093/nar/gkn221.
Ro Seungil, Park Chanjae, Jin Jingling, Sanders Kenton, Yan Wei: A PCR-based Method for Detection and Quantification of Small RNAs. Biochem Biophys Res Commun. 2006, 351 (3): 756-763. 10.1016/j.bbrc.2006.10.105.
Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, Chen X, Dreyfuss G, Eddy SR, Griffiths-Jones S, Marshall M, Matzke M, Ruvkun G, Tuschl T: A uniform system for microRNA annotation. RNA. 2003, 9 (3): 277-279. 10.1261/rna.2183803.
Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant microRNAs. Plant Cell. 2008, 20: 3186-3190. 10.1105/tpc.108.064311.
Axtell MJ: Evolution of microRNAs and their targets: Are all microRNAs biologically relevant?. Biochim Biophys Acta. 2008, 1779: 725-727.
Bartel DP: MicroRNAs: target recognition and regulatory functions. Cell. 2009, 136 (2): 215-33. 10.1016/j.cell.2009.01.002.
Adai A, Johnson C, Mlotshwa S, Archer-Evans S, Manocha V, Archer-Evans S, Vance V, Sundaresan V: Computational prediction of miRNAs in Arabidopsis thaliana. Genome Res. 2005, 15: 78-91. 10.1101/gr.2908205.
Sunkar R, Girke T, Zhu JK: Identification and characterization of endogenous small interfering RNAs from rice. Nucleic Acids Res. 2005, 33: 4443-4454. 10.1093/nar/gki758.
Jin W, Li N, Zhang B, Wu F, Li W: Identification and verification of microRNA in wheat (Triticum aestivum). J Plant Res. 2008, 121: 351-355. 10.1007/s10265-007-0139-3.
Carra A, Mica E, Gambino G, Pindo M, Moser C: Cloning and characterization of small non-coding RNAs from grape. Plant J. 2009, 59: 750-763. 10.1111/j.1365-313X.2009.03906.x.
Lu S, Sun YH, Chiang VL: Stress-responsive microRNAs in Populus. Plant J. 2008, 55: 131-151. 10.1111/j.1365-313X.2008.03497.x.
Rubio-Somoza I, Weigel D: MicroRNA networks and developmental plasticity in plants. Trends Plant Sci. 2011, 16 (5): 258-264. 10.1016/j.tplants.2011.03.001.
Liu Q, Chen YQ: Insights into the mechanism of plant development: Interactions of miRNAs pathway with phytohormone response. Biochemical and Biophysical Research Communications. 2009, 384: 1-5. 10.1016/j.bbrc.2009.04.028.
Schaffer AA, Petreikov M: Sucrose-to-Starch Metabolism in Tomato Fruit Undergoing Transient Starch Accumulation. Plant Physiol. 1997, 113 (3): 739-746.
Hoen PA, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RH, de Menezes RX, Boer JM, van Ommen GJ, den Dunnen JT: Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008, 36: 141-151. 10.1093/nar/gkn705.
Elitzur T, Vrebalov J, Giovannoni JJ, Goldschmidt EE, Friedman H: The regulation of MADS-box gene expression during ripening of banana and their regulatory interaction with ethylene. J Exp Bot. 2010, 61 (5): 1523-1535. 10.1093/jxb/erq017.
Manning K, Tör M, Poole M, Hong Y, Thompson AJ, King GJ, Giovannoni JJ, Seymour GB: A naturally occurring epigenetic mutation in a gene encoding an SBP-box transcription factor inhibits tomato fruit ripening. Nat Genet. 2006, 38: 948-952. 10.1038/ng1841.
Adams-Phillips L, Barry C, Kannan P, Leclercq J, Bouzayen M, Giovannoni J: Evidence that CTR1-mediated ethylene signal transduction in tomato is encoded by a multigene family whose members display distinct regulatory features. Plant Mol Biol. 2004, 54: 387-404.
Bleecker AB, Kende H: Ethylene: A gaseous signal molecule in plants. Annu Rev Cell Dev Biol. 2000, 16: 1-18. 10.1146/annurev.cellbio.16.1.1.
Guo H, Ecker JR: The ethylene signaling pathway: New insights. Curr Opin Plant Biol. 2004, 7: 40-49. 10.1016/j.pbi.2003.11.011.
Kendrick MD, Chang C: Ethylene signaling: New levels of complexity and regulation. Curr Opin Plant Biol. 2008, 11: 479-485. 10.1016/j.pbi.2008.06.011.
Zhu Z, Guo H: Genetic basis of ethylene perception and signal transduction in Arabidopsis. J Integr Plant Biol. 2008, 50: 808-815. 10.1111/j.1744-7909.2008.00710.x.
Xin M, Wang Y, Yao Y, Xie C, Peng H, Ni Z, Sun Q: Diverse set of microRNAs are responsive to powdery mildew infection and heat stress in wheat (Triticum aestivum L.). BMC Plant Biol. 2010, 10: 123-134. 10.1186/1471-2229-10-123.
Barakat A, Wall PK, Diloreto S, Depamphilis CW, Carlson JE: Conservation and divergence of microRNAs in Populus. BMC Genomics. 2007, 8: 481-496. 10.1186/1471-2164-8-481.
Szittya G, Moxon S, Santos DM, Jing R, Fevereiro MP, Moulton V, Dalmay T: High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families. BMC Genomics. 2008, 9: 593-602. 10.1186/1471-2164-9-593.
Herr AJ: Pathways through the small RNA world of plants. FEBS Lett. 2005, 579: 5879-5888. 10.1016/j.febslet.2005.08.040.
Vazquez F: Arabidopsis endogenous small RNAs: highways and byways. Trends Plant Sci. 2006, 11: 460-468. 10.1016/j.tplants.2006.07.006.
Allen E, Xie Z, Gustafson AM, Sung GH, Spatafora JW, Carrington JC: Evolution of microRNA genes by inverted duplication of target gene sequences in Arabidopsis thaliana. Nat Genet. 2004, 36: 1282-1290. 10.1038/ng1478.
Mohorianu I, Schwach F, Jing R, Lopez-Gomollon S, Moxon S, Szittya G, Sorefan K, Moulton V, Dalmay T: Profiling of short RNAs during fleshy fruit development reveals stage-specific sRNAome expression patterns. Plant J. 2011, 67 (2): 232-246. 10.1111/j.1365-313X.2011.04586.x.
Rhoades MW, Reinhart BJ, Lim LP, Burge CB, Bartel B, Bartel DP: Prediction of plant microRNA targets. Cell. 2002, 110: 513-520. 10.1016/S0092-8674(02)00863-2.
Floyd SK, Bowman JL: Gene regulation: ancient microRNA target sequences in plants. Nature. 2004, 435: 441-445.
Aukerman MJ, Sakai H: Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes. Plant Cell. 2003, 15: 2730-2741. 10.1105/tpc.016238.
Chen X: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2003, 1: 11-17.
Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D: Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005, 8: 517-527. 10.1016/j.devcel.2005.01.018.
Gandikota M, Birkenbihl RP, Hohmann S, Cardon GH, Saedler H, Huijser P: The miRNA156/157 recognition element in the 3' UTR of the Arabidopsis SBP box gene SPL3 prevents early flowering by translational inhibition in seedlings. Plant J. 2007, 49: 683-693. 10.1111/j.1365-313X.2006.02983.x.
Lanet E, Delannoy E, Sormani R, Floris M, Brodersen P, Crete P, Voinnet O, Robaglia C: Biochemical evidence for translational repression by Arabidopsis microRNAs. Plant Cell. 2009, 21: 1762-1768. 10.1105/tpc.108.063412.
Real MD, Company P, García-Agustín P, Bennett AB, González-Bosch C: Characterization of tomato endo-beta -1, 4- glucanase Cel1 protein in fruit during ripening and after fungal infection. Planta. 2004, 220 (1): 80-86. 10.1007/s00425-004-1321-7.
Flors V, Leyva Mde L, Vicedo B, Finiti I, Real MD, García-Agustín P, Bennett AB, González-Bosch C: Absence of the endo-beta-1,4-glucanases Cel1 and Cel2 reduces susceptibility to Botrytis cinerea in tomato. Plant J. 2007, 52 (6): 1027-40. 10.1111/j.1365-313X.2007.03299.x.
Ordaz-Ortiz JJ, Marcus SE, Knox JP: Cell wall microstructure analysis implicates hemicellulose polysaccharides in cell adhesion in tomato fruit pericarp parenchyma. Mol Plant. 2009, 2 (5): 910-21. 10.1093/mp/ssp049.
Marín-Rodríguez MC, Orchard J, Seymour GB: Pectate lyases, cell wall degradation and fruit softening. J Exp Bot. 2002, 53 (377): 2115-2119. 10.1093/jxb/erf089.
Rodoni L, Casadei N, Concellón A, Chaves Alicia AR, Vicente AR: Effect of short-term ozone treatments on tomato (Solanum lycopersicum L.) fruit quality and cell wall degradation. J Agric Food Chem. 2010, 58 (1): 594-599. 10.1021/jf9029145.
Moctezuma E, Smith DL, Gross KC: Antisense suppression of a beta-galactosidase gene (TB G6) in tomato increases fruit cracking. J Exp Bot. 2003, 54 (390): 2025-2033. 10.1093/jxb/erg214.
Marano MR, Carrillo N: Constitutive transcription and stable RNA accumulation in plastids during the conversion of chloroplasts to chromoplasts in ripening tomato fruits. Plant Physiol. 1992, 100: 1103-1113. 10.1104/pp.100.3.1103.
Kahla US, Bock R: Plastid transcriptomics and translatomics of tomato fruit development chloroplast-to- chromoplast differentiation: chromoplast gene expression largely serves the production of a single protein. Plant Cell. 2008, 20: 856-874. 10.1105/tpc.107.055202.
Bapat VA, Trivedi PK, Ghosh A, Sane VA, Ganapathi TR, Nath P: Ripening of fleshy fruit: molecular insight and the role of ethylene. Biotechnol Adv. 2010, 28 (1): 94-107. 10.1016/j.biotechadv.2009.10.002.
Liu Q, Zhang YC, Wang CY, Luo YC, Huang QJ, Chen SY, Zhou H, Qu LH, Chen YQ: Expression analysis of phytohormone-regulated microRNAs in rice, implying their regulation roles in plant hormone signaling. FEBS Lett. 2009, 583: 723-728. 10.1016/j.febslet.2009.01.020.
Wang X, Kong H, Ma H: F-box proteins regulate ethylene signaling and more. Genes Dev. 2009, 23 (4): 391-396. 10.1101/gad.1781609.
Qiao H, Chang KN, Yazaki J, Ecker JR: Interplay between ethylene, ETP1/ETP2 F-box proteins, and degradation of EIN2 triggers ethylene responses in Arabidopsis. Genes Dev. 2009, 23 (4): 512-521. 10.1101/gad.1765709.
Benavente LM, Alonso JM: Molecular mechanisms of ethylene signaling in Arabidopsis. Mol Biosyst. 2006, 2 (3-4): 165-73. 10.1039/b513874d.
Zhu Z, Guo H: Genetic basis of ethylene perception and signal transduction in Arabidopsis. J Integr Plant Biol. 2008, 50 (7): 808-815. 10.1111/j.1744-7909.2008.00710.x.
Frankowski K, Kesy J, Kotarba W, Kopcewicz J: Ethylene signal transduction pathway. Postepy Biochem. 2008, 54 (1): 99-106.
Zhong S, Lin Z, Grierson D: Tomato ethylene receptor-CTR interactions: visualization of NEVER-RIPE interactions with multiple CTRs at the endoplasmic reticulum. J Exp Bot. 2008, 59 (4): 965-972. 10.1093/jxb/ern021.
Adams-Phillips L, Barry C, Kannan P, Leclercq J, Bouzayen M, Giovannoni J: Evidence that CTR1-mediated ethylene signal transduction in tomato is encoded by a multigene family whose members display distinct regulatory features. Plant Mol Biol. 2004, 54 (3): 387-404.
Leclercq J, Adams-Phillips LC, Zegzouti H, Jones B, Latché A, Giovannoni JJ, Pech JC, Bouzayen M: LeCTR1, a tomato CTR1-like gene, demonstrates ethylene signaling ability in Arabidopsis and novel expression patterns in tomato. Plant Physiol. 2002, 130: 1132-1142. 10.1104/pp.009415.
Zegzouti H, Jones B, Frasse P, Marty C, Maitre B, Latch A, Pech JC, Bouzayen M: Ethylene-regulated gene expression in tomato fruit: characterization of novel ethylene-responsive and ripening-relatedgenes isolated by differential display. Plant J. 1999, 18 (6): 589-600. 10.1046/j.1365-313x.1999.00483.x.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: 179-187. 10.1093/nar/gni178.
Tang F, Hajkova P, Barton SC, Lao K, Surani MA: MicroRNA expression profiling of single whole embryonic stem cells. Nucleic Acids Res. 2006, 34: 9-15. 10.1093/nar/gnj009.
Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP: Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007, 3: 12-23. 10.1186/1746-4811-3-12.
Varkonyi-Gasic E, Gould N, Sandanayaka M, Sutherland P, MacDiarmid RM: Characterisation of microRNAs from apple (Malus domestica 'Royal Gala') vascular tissue and phloem sap. BMC Plant Biology. 2010, 10: 159-173. 10.1186/1471-2229-10-159.
We were grateful to Professor Tamas Dalmay from University of East Anglia (The United Kingdom) for the kind supply of miRNAs sequence information. This work was supported by the National Natural Sciences Foundation of China (30871741 and 30972037) and National Basic Research Program ("973" program, 2011CB100604) of China.
BZ Z and JH Z design the experiments, JH Z conducted the experiments, YZ M, LH C, YB L, DQ F, Y Z and BZ Z participate in bioinformatics analysis, YX W and Z J participate in RNA extraction, JH Z wrote the paper. All authors read and approved the final manuscript.
Jinhua Zuo, Benzhong Zhu contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.