Identification and characterization of miRNAome in root, stem, leaf and tuber developmental stages of potato (Solanum tuberosum L.) by high-throughput sequencing
- Nisha Lakhotia†1,
- Gopal Joshi†1,
- Ankur R Bhardwaj1,
- Surekha Katiyar-Agarwal2,
- Manu Agarwal1,
- Arun Jagannath1,
- Shailendra Goel1 and
- Amar Kumar1Email author
© Lakhotia et al.; licensee BioMed Central Ltd. 2014
Received: 20 August 2013
Accepted: 26 December 2013
Published: 7 January 2014
MicroRNAs (miRNAs) are ubiquitous components of endogenous plant transcriptome. miRNAs are small, single-stranded and ~21 nt long RNAs which regulate gene expression at the post-transcriptional level and are known to play essential roles in various aspects of plant development and growth. Previously, a number of miRNAs have been identified in potato through in silico analysis and deep sequencing approach. However, identification of miRNAs through deep sequencing approach was limited to a few tissue types and developmental stages. This study reports the identification and characterization of potato miRNAs in three different vegetative tissues and four stages of tuber development by high throughput sequencing.
Small RNA libraries were constructed from leaf, stem, root and four early developmental stages of tuberization and subjected to deep sequencing, followed by bioinformatics analysis. A total of 89 conserved miRNAs (belonging to 33 families), 147 potato-specific miRNAs (with star sequence) and 112 candidate potato-specific miRNAs (without star sequence) were identified. The digital expression profiling based on TPM (Transcripts Per Million) and qRT-PCR analysis of conserved and potato-specific miRNAs revealed that some of the miRNAs showed tissue specific expression (leaf, stem and root) while a few demonstrated tuberization stage-specific expressions. Targets were predicted for identified conserved and potato-specific miRNAs, and predicted targets of four conserved miRNAs, miR160, miR164, miR172 and miR171, which are ARF16 (Auxin Response Factor 16), NAM (NO APICAL MERISTEM), RAP1 (Relative to APETALA2 1) and HAM (HAIRY MERISTEM) respectively, were experimentally validated using 5′ RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends). Gene ontology (GO) analysis for potato-specific miRNAs was also performed to predict their potential biological functions.
We report a comprehensive study of potato miRNAs at genome-wide level by high-throughput sequencing and demonstrate that these miRNAs have tissue and/or developmental stage-specific expression profile. Also, predicted targets of conserved miRNAs were experimentally confirmed for the first time in potato. Our findings indicate the existence of extensive and complex small RNA population in this crop and suggest their important role in pathways involved in diverse biological processes, including tuber development.
In recent years, small non coding RNAs (small ncRNAs) have emerged as key regulators of gene expression in both plants and animals [1, 2]. In plants, there are many classes of small ncRNAs and the two major classes that have been well-studied and defined are miRNAs and siRNAs [1, 2]. Although both are functionally similar, they have different mode of biogenesis. siRNAs, 20–24 nt long, are processed from perfect double stranded RNA. They are mainly involved in heterochromatin modification, maintaining genome integrity by silencing of transgenes and transposons, defence against viruses and regulation of gene expression . There are diverse siRNA species depending on the Dicer-like (DCL) proteins involved in the production, such as natural antisense transcript-derived siRNAs (24 nt), heterochromatic siRNAs (24 nt) and trans-acting siRNAs (21 nt) . miRNAs are small, endogenous, single-stranded and ~21 nt long RNAs that regulate gene expression at the post-transcriptional level by degrading or repressing translation of targeted mRNAs, depending on the extent of complementarity between miRNA and mRNA [1, 2].
miRNAs are transcribed by RNA polymerase II into an imperfect long folded back structure, called as primary miRNA (pri-miRNA). The pri-miRNA is cleaved by Dicer-like 1 (DCL1) protein (endoribonuclease having RNase III activity), in association with other protein factors (HYL1 and others), into a precursor sequence with hairpin structure (pre-miRNA), which is further processed by DCL1 into a miRNA:miRNA* duplex . One of the strands of this duplex (guide strand or mature miRNA) gets incorporated into the silencing complex or effector complex (RISC- RNA Induced Silencing Complex) containing argonaute proteins, while the other strand (miRNA* or passenger strand) gets degraded. Within the complex, miRNA guides the target mRNA cleavage or inhibits translation with the help of one of the argonaute proteins (AGO1) . There is increasing evidence which indicates that miRNAs play vital roles in stress responses and various developmental processes, including shoot apical meristem formation, leaf morphogenesis and polarity, floral organ identity, root development, vegetative phase change and vascular development [2, 4].
A large number of miRNAs are evolutionary conserved among diverse species within the plant kingdom while several miRNAs, that are considered to be recently evolved show species-specificity and often express at lower levels relative to conserved miRNAs [5, 6]. Due to their low expression levels, most of the species-specific miRNAs remained unidentified in many plant species. However, in recent years, with the advent of high-throughput sequencing technology, both species-specific and conserved miRNAs have been identified in diverse plant species [6–15].
Potato (Solanum tuberosum L.) is the fourth most important crop worldwide, after rice, wheat and barley. The main reasons for the increasing popularity of potato are high nutritional value combined with the simplicity of its propagation by vegetative reproduction. Potato tubers are derived from the swollen tips of specialized basal lateral juvenile shoots, called stolons. Tuber development is a highly systematic and well-coordinated morpho-physiological process that involves interactions between several environmental, biochemical and genetic factors [16–18]. Since control of tuber initiation impacts many commercial aspects of crop improvement, including scheduling, size grade distribution, yield and expansion of climate zone suitable for propagation, this is a high priority area for potato research. Several transcription factors, hormones, environmental and molecular signals have been studied that regulate tuberization [16–18]. Although, the potato genome sequencing project  has revealed that both gene family expansion and recruitment of existing genes for new pathways have contributed to the evolution of tuber development in potato, this phenomena still remains a challenge for plant biologists.
To date, only miR172 and miR156 have been shown to play a role in tuberization process [20, 21] and only a few miRNAs have been shown to regulate defense genes in potato . Since the potato genome has been completely sequenced recently , it will greatly help in identifying the entire population of small RNAs and their target genes that might be involved in various developmental processes in potato. In previous studies, through an in silico approach, many conserved miRNAs were identified in potato based on miRBase repository by using potato EST, GSS and nr (non-redundant) databases (existing in the National Centre for Biotechnology Information) [23–26]. Earlier, large scale sequencing studies were limited to the identification of only conserved miRNAs in potato . A recent study reported 28 conserved and 120 potato-specific miRNA families in potato by deep sequencing of sRNA libraries from leaf and stolon tissues . However, knowledge about the potato-specific miRNA population in different vegetative tissues and developmental stages of tuberization remain largely unknown.
In this study, we report a population of potato miRNAs from vegetative organs (leaf, stem and root) and four early development stages of tuberization that might be involved in regulating developmental processes, including tuberization. By taking advantage of the recently sequenced potato genome , 33 conserved miRNA families and 147 potato-specific miRNA were identified using high-throughput sequencing and bioinformatics pipeline. We also report the expression profiles of both conserved and potato-specific miRNAs in different tissues (leaf, stem and root) and different developmental stages of tuberization based on digital and qRT-PCR expression analyses. Additionally, the putative target genes for the conserved and potato-specific miRNAs were predicted and 5′ RLM-RACE was carried out to map the cleavage site of four conserved miRNAs targets. GO analysis of potato-specific miRNAs was also performed to gain a better understanding of their functions in various developmental processes in potato.
Analysis of small RNA population
Summary of the reads obtained from all the potato small RNA libraries
Types of reads
Filter by sequence length
Filter low-complexity sequences
Filter invalid sequences
Filter by t/rRNA (matches out)
Putative Small RNA Population
Identification of conserved and potato-specific miRNAs in potato
Digital expression profiling of predicted miRNAs
The read abundance of predicted miRNAs in different libraries was normalized and expressed as transcripts per million (TPM) to measure the expression levels of these miRNAs in different samples (Additional file 1; Additional file 2; Additional file 3).
The sequencing frequencies of a majority of potato-specific miRNAs were found to be low, consistent with the hypothesis that species-specific miRNAs are often expressed at a lower level as compared to conserved miRNAs [5, 6]. Most of these miRNAs (about 65%) were detected with total TPM less than 5 (Additional file 2). However, a few miRNAs (14) were found to have total TPM more than 100. Among all the potato-specific miRNAs, miRNA53 and miRNA72 were the most abundant, with total TPM more than 15000 (Additional file 2).
Some predicted conserved and potato-specific miRNAs displayed significant variation in expression level among the different tissues (Additional file 1 and Additional file 2). For example, conserved miRNA - miR156 exhibited high expression during tuberization stages, with highest level at PT0 (Additional file 1). Among potato-specific miRNAs, miRNA72 had highest sequencing frequencies in stem and displayed significant change in expression level during tuberization stages (with highest expression at PT1) (Additional file 2). These findings suggest the probable role of these regulatory miRNAs in specific tissues or during specific developmental stages in potato. In this study, biological replicates were not used for sequencing. Hence, the digital expression data only gives an indication of variation in expression level between different tissues. This could be further confirmed by qRT-PCR as has been done for some of the miRNAs in this study (see below).
Expression profile of conserved and potato-specific miRNAs in different tissues by qRT-PCR
Expression profile of miRNAs in vegetative tissues
Among conserved miRNAs, the expression of miR164_1 and miRNA171_10 appeared to be restricted to stem and root, respectively while miR8044_1 showed abundant expression in all the examined vegetative tissues. miR394_1 exhibited the highest expression in stem and moderate expression in leaf. miR171_2 and miR157_1 showed low expression in root whereas relatively high expression was observed in leaf and stem. The expression of miR396_1 and miR172_5 was found to be most abundant in stem and moderate in root whereas the expression of miR399_1 was relatively high in root. miR172_1 had high expression level in root, low expression in stem and could rarely be detected in leaf while miR159_1 demonstrated highest expression in stem but low expression in leaf and root. The expression of miR6149_1 and miR482_1 was higher in leaf, stem and lower in root. miR171_9 showed strong expression in stem and root but weak expression in leaf and miR398_1 was found to be expressed predominantly in stem, followed by leaf and root.
Among potato-specific miRNAs, miRNA53 had higher level of expression in leaf and stem, and moderate level in root. miRNA193 appeared to be highly expressed in leaf and weakly in stem and root. The expression pattern of miRNA111, miRNA72, miRNA34 were found to be very similar - high expression in stem, followed by root and leaf while miRNA152 exhibited moderate level of expression in root and low expression in leaf, stem.
Expression profile of miRNAs during tuber developmental stages
The expression of miR164_1, miR399_1, miR157_1, miR398_1, miR171_10, miR171_9 (conserved miRNAs) and miRNA53 (potato-specific miRNAs) could barely be detected during tuberization stages. Among potato-specific miRNAs, miRNA111, miRNA72 and among conserved miRNAs, miR171_2, miR396_1, miR159_1, miR6149_1, miR8044_1 and miR482_1 had very similar pattern of expression during different tuberization stages. They showed a slight increase in expression at PT1 stage and then a gradual decrease during later stages whereas miRNA34 and miR394_1 was found to be expressed at low level in all the stages with no variation in expression levels. miR172_5 showed moderate level of expression during tuberization - similar expression level at PT1 and PT2 stage, then a slight decrease in abundance at PT3 stage. The expression of miRNA193, miRNA152 (potato-specific miRNAs) and miR172_1 was relatively high during tuberization. miRNA193 and miR172_1 showed a gradual increase in abundance from PT1–PT2, reaching the maximum expression at PT2, then a decrease thereafter whereas miRNA152 showed similar level of expression during all the stages.
It was found that digital expression data and qRT-PCR expression profiles of miRNAs in different tissues and developmental stages of tuberization do not correlate with each other. The probable reason for this variation could be the biasness of RNA ligase activity during small RNA library preparation .
Target prediction for conserved and potato-specific miRNAs
In order to explore the functional role of all the identified miRNAs in diverse biological processes, their putative target genes were predicted using the open source web server psRNATarget  with default parameters (see Methods). The transcript sequences of the potato genome (Group Phureja DM1-3 516R44) were used as a reference set. Additional file 5 and Additional file 6 provide detailed information regarding the predicted targets for conserved and potato-specific miRNAs respectively. A majority of the miRNAs were found to target more than one transcript.
Among conserved miRNA targets, most of them were found to be transcripts coding for transcription factors (Additional file 5), such as Squamosa promoter-binding protein (regulated by miR156), GRAS family transcription factors (targeted by miR171), GAMYB-like2 (targeted by miR159), APETALA2 (target of miR172), NAC domain containing protein (targeted by miR164), Auxin response factors (regulated by miR160), PHAVOLUTA-like HD-ZIPIII protein (target of miR166) and nuclear transcription factors - YA4, YA5, YA6 (targeted by miR169). In addition to transcription factors, other targets included mRNA coding for F-box family protein (miR394), disease resistance protein and fiber expressed protein (miR159), laccase (miR397), salt tolerance protein (miR157), UDP-glucoronate decarboxylase 2 (miR164), DNA binding protein (miR166 and miR396), NL25 disease resistance protein (miR482), protein phosphatase and kinase (miR390), AGO1-1 (miR168), galactose oxidase (miR6149) and proteins with unknown functions (Additional file 5). Most of the conserved miRNA targeted plant transcription factors (SBP, NAC, ARF, GRAS, YA, AP2) were found to be similar to the conserved miRNA targets predicted in Arabidopsis and other plants [2, 9, 41], underlining the role of conserved miRNAs in essential biological processes.
For potato-specific miRNAs, targets were predicted successfully for 142 out of 147 miRNAs. Target genes were not found for the remaining five miRNAs which could be due to stringent cut-off used in our prediction analysis to reduce the number of false positives and to obtain reliable results. Their targets included mRNA coding for tropinone reductase I, Cdk10/11, endo-1,4-beta-D-glucanase, XH/XS domain containing protein, ankyrin repeat containing protein, cytochrome P450, alpha amylase, glucosyltransferase and others (Additional file 6). A few potato-specific miRNAs showed lipoxygenase and protein phosphatase type 2A, that are known to be involved in potato tuberization [42, 43], as their predictive targets (Additional file 6). In addition, some potato-specific miRNAs were also found to target transcription factors including NAC domain protein, GRAS family transcription factors, MADS-box family protein and ARF domain class transcription factor (Additional file 6). These transcription factors have been known to be involved in regulating developmental processes in plants . These results imply that these potato-specific miRNAs might be involved in some specific developmental as well as essential biological processes in potato. Many of the potato-specific miRNA target transcripts that code for proteins of unknown function.
Experimental validation of conserved miRNA targets
Since the discovery of sRNAs (miRNAs & siRNAs) as regulators of gene expression in C. elegans, there has been intense interest in the sRNA biology in eukaryotic organisms. In plants, extensive studies in the past several years have clearly shown that miRNAs regulate gene expression at the post-transcriptional level by degrading or repressing translation of targeted mRNAs and play important roles during plant development and growth . Here, we report identification and characterization of a population of both conserved and potato-specific miRNAs from different vegetative tissues and various well-defined stages of tuberization. Also, we discuss the potential role of these validated miRNAs in various developmental processes in potato.
In this study, with the application of high-throughput sequencing technology, we identified 89 conserved miRNA (belonging to 33 families), 147 potato-specific miRNAs and 112 candidate potato-specific miRNAs, based on the sequenced potato genome dataset . Consistent with previous reports [5, 6], a majority of the potato-specific miRNAs predicted in our analysis, showed low sequencing frequencies as compared to conserved miRNAs in the dataset. The majority of the pre-miRNAs were found to be encoded by intergenic regions which is also consistent with earlier findings .
The analysis of sRNA population in potato revealed that 24 nt long sRNAs (which mainly contains heterochromatic siRNAs) dominated the sRNA transcriptome, in both total reads and unique sequences, as observed in many other plant species [6–15]. The high percentage of 24 nt sRNAs suggests an important role of this class in maintaining genome integrity by heterochromatic-histone modification . The high ratio of total reads/unique sequences for 21 nt length class indicated the high abundance of this size class. The majority of the identified miRNAs were 21 nt in length, which is the canonical size for miRNAs generated from DCL1 processing. Only a few miRNAs showed minor variations in size (20 or 22 nt long) which could be due to inaccuracy of DCL1 processing .
Previously, many conserved miRNAs were reported in potato by a computational approach using available EST, GSS and nr databases [23–26]. A recent study by Zhang et al. (2013) has led to the identification of 120 potato-specific miRNA families and 28 conserved miRNA families in leaf and stolon tissues by high-throughput sequencing . However, majority of potato-specific miRNAs predicted by Zhang et al. (2013) were also detected in our sequencing dataset but could not meet the criteria for miRNA annotation employed in our analysis and thus excluded from further analysis.
To observe the expression profile of identified miRNAs, qRT-PCR was performed for 15 conserved and 6 potato-specific miRNAs. Our qRT-PCR results showed that majority of potato-specific (4) and conserved miRNAs (13) expressed at relatively high level in the stem tissues. This indicates that either these miRNAs express highly in the stem tissues or most of these miRNAs are being continuously transported through the conductive vascular tissues. In case of their high expression in stem it is possible that they might play an important role in stem development and differentiation. On the other hand, potato miR172 and other miRNAs (miR399, miR395) have been shown to be involved in graft transmissible movement through the conductive vascular tissues [20, 47]. Interestingly, potato-specific miRNAs miRNA193, miRNA152 and the conserved miR172_1, miR172_5 showed significant expression during developmental stages of tuberization. This indicates a probable role for these miRNAs in potato tuberization. In our analysis, the expression profile of miR172 (known to be involved in potato tuberization) during developmental stages of tuberization was consistent with the previous report . Also, we found that miR172_1 expressed at a high level in potato roots, which was not reported previously.
To investigate the functional importance of the identified miRNAs, targets were predicted using psRNA target finder server  and the annotated potato genome transcript data . Among conserved miRNAs, miR164 was predicted to target mRNA encoding for NAC domain containing proteins. It has been also reported that miR164 is a key regulator in diverse developmental processes including floral, embryonic, vegetative organs and lateral root development [48, 49]. PHAVOULTA like HD-ZIP III transcription factor was one of the predicted targets of miR166 in our study. Previous studies showed that miR166/miR165 targets HD-ZIP III transcription factor genes and regulate a wide range of developmental processes including meristem formation, vascular development, lateral organ polarity, root and nodule development [50–52]. In our analysis, miR159 and miR319 were predicted to target mRNA coding for GAMYB like transcription factor. GAMYB transcription factors are involved in GA signalling and it has been proven that miR159-mediated regulation of GAMYB affects plant anther development and flowering time in Arabidopsis . Auxin response factor was found as a predicted target of miR160. It has been observed that miR160 targets ARF10, ARF16 and ARF17 and regulate various aspects of plant development in Arabidopsis . Squamosa-promoter binding protein was predicted as the target of miR156 and miR157. Previous studies have demonstrated that miR156 plays a vital role in controlling leaf development, apical dominance, floral transition and development [55–57] by targeting several members of the Squamosa promoter binding protein like (SPL) family of plant-specific transcription factors. Interestingly, a recent study has shown that overexpression of miR156 in potato leads to development of stolon-borne aerial tubers from distal nodes . The digital expression profile reported in this study also showed high abundance of miR156 in tuber developmental stages with highest level at PT0 (stolon development), suggesting a role of miR156 in the phase transition from basal lateral juvenile shoots to a storage organ in the process of tuber development. In our analysis, miR171 was found to target GRAS family transcription factors and HAIRY MERISTEM (a member of GRAS family transcription factor). They have been shown to be involved in regulating various aspects of plant growth and development . In potato, our qRT-PCR data showed that all the three members of miR171 expressed highly in vegetative tissues (leaf and stem-miR171_2, stem and root-miR171_9, root-miR171_10), indicating a role of miR171 in regulating potato developmental process. Previous studies have indicated that miR172 largely targets gene encoding members of APETALA2 (AP2) and AP2-like transcription factor family and plays a major role in the regulation of flowering time, vegetative phase change and floral organ identity, by cleaving and repressing the translation of its target mRNAs [59–61]. Recently, miR172 has been shown to induce tuberization by regulating long distance signals and targeting RAP1 mRNA (Relative to APETALA2 1) . Our prediction analysis showed RAP1, which is a known target of miR172 in potato , as a predicted target of miR172. This data supports the efficacy of our target prediction analysis.
Further, miRNA-mediated cleavage of target of conserved miRNAs was confirmed by 5′ RLM-RACE. Several targets of conserved miRNAs have been validated using RACE analysis in tomato, pepper and rice [8, 13, 62]. However, this is the first study in potato where 5′ RLM-RACE was performed to authenticate predicted targets of conserved miRNAs. Our RACE results confirmed that transcripts of AFR16, NAM, RAP1 and HAM were cleaved in vivo by miR160, miR164, miR172 and miR171, respectively. For potato-specific miRNAs, degradome sequencing approach is being used for mapping of the cleavage site.
Our target prediction analysis in potato and other studies on plant species, especially rice and Arabidopsis, strongly suggest that most of the conserved miRNAs target plant transcription factors (SBP, NAC, ARF, GRAS and AP2) [2, 9, 41], thereby emphasizing the role of conserved miRNA in essential biological processes in plants, including potato tuberization.
GO was performed on potato-specific miRNAs targets to gain a better understanding of their possible function . GO results revealed that majority of the targets are involved in oxidation-reduction processes, metabolic processes, regulation of transcription, transport, signal transduction and defense responses. These results suggest that they covered a broad range of physiological and metabolic processes. Thus, in contrast to conserved miRNA targets, most of the putative targets of potato-specific miRNAs were not transcription factors. A few potato-specific miRNAs were predicted to target transcripts encoding transcription factors, including NAC domain protein, GRAS family transcription factors, MADS-box family protein and ARF domain class transcription factor. These transcription factors are known to be targeted by conserved miRNAs and play a major role in various aspects of plant growth and development . Thus, it appears that the potato-specific miRNAs might play important roles in essential biological processes as well. Future studies are essential to explain the functional significance of the inverse relationship between miRNA and their target which will improve our knowledge about the regulatory role of miRNAs in potato.
With respect to potato tuberization, expression profiling by qRT-PCR and target prediction studies indicated that a few miRNAs might play a role during transition from stolon to tuber development. qRT-PCR results indicated that miR172_1, miR172_5 (conserved miRNAs) and miRNA193, miRNA152 (potato-specific miRNAs) showed significant expression during developmental stages of tuberization. Lipoxygenase (LOX) was predicted as one of the targets of miRNA150, miRNA163, miRNA37, miRNA78 and miRNA108. Protein phosphatase type 2A (PP2A) was found to be a predicted target of miRNA90. Previous studies have suggested the role of miR172, LOX and PP2A in potato tuber developmental process [20, 42, 43]. Therefore, these potential miRNAs and their target genes are being further investigated to determine their roles in potato tuber development.
In summary, we report a comprehensive study of potato miRNAs at genome-wide level using three different vegetative tissues and four early developmental stages of tuberization by high-throughput sequencing and high-end bioinformatics methods. We have identified a population of potato-specific miRNAs that could fill the gap in understanding the regulatory role of miRNAs in potato. The digital expression data was further quantified by qRT-PCR to determine up- and down-regulations of potato miRNAs in various vegetative tissues and developmental stages of tuberization. The putative target genes were predicted for conserved and potato-specific miRNAs to gain a better understanding of their function. Further, a few conserved miRNA targets were also validated by 5′ RLM-RACE and the predicted target genes for potato-specific miRNAs were subjected to GO annotation to infer the functional role of target genes.
The plantlets of Solanum tuberosum (cv. Kufri Chandramukhi) were raised in vitro on MS medium  and maintained under the 16 hr light/8 hr dark photoperiod conditions at 22°C. The leaves, stem and root were collected from 45 days old whole in vitro grown potato plantlets and immediately frozen in liquid nitrogen, subsequently stored at -80°C until used. These in vitro grown plantlets were then potted and grown under natural field conditions, during the period of November-January, in the net house at the Department of Botany, University of Delhi, India. Tissues from the four early developmental stages of tuberization were harvested from these naturalized plants, frozen in liquid nitrogen and stored at -80°C until used.
Total RNA isolation
Total RNA was isolated from the above mentioned tissues (0.2 grams each) using TRI reagent (Sigma, USA) as per the manufacturer’s instructions. The quality and quantity of RNA isolates was verified on denaturing agarose gel and by optical density (OD) measurement, respectively.
Small RNA library construction
Small RNA libraries for each tissue type were prepared independently using an Illumina Small RNA sample prep kit (Illumina, U.S.A), according to manufacturer’s instructions. Briefly, small RNAs (20–40 nt) were separated on 15% denaturing polyacrylamide gel and purified. Next, small RNAs were ligated with the 5′ and 3′adapters sequentially, eluted and gel purified. The 5′and 3' ligated small RNAs were reverse transcribed using SuperScript II Reverse Transcriptase (Invitrogen, U.S.A) and PCR amplified with Phusion DNA Polymerase. The quality and quantity of cDNA library was assessed using Agilent Bioanalyzer (Agilent Technologies, U.S.A). Small RNA library preparation and their high-throughput sequencing was performed using an Illumina Genome Analyzer IIx (Illumina, U.S.A) at DBT-funded High Throughput Sequencing Facility at University of Delhi South Campus, India.
Bioinformatic analysis of sRNA sequences
The removal of poor quality sequences and trimming of adaptor sequences from the raw sequence data was carried out using the UEA sRNA workbench 2.4 - Plant version sequence file pre-processing (http://srna-tools.cmp.uea.ac.uk/) . Sequences, smaller than 16 nt and larger than 30 nt were also removed. High quality trimmed sequences (16–30 nt in length, reads with no “N”, no more than 4 bases with quality score <10 and no more than 6 bases with quality score <13) were mapped onto plant t/rRNAs sequences from “Rfam” (excepted miRNAs), Arabidopsis tRNAs from “The Genomic tRNA Database”, and plant t/rRNA sequences from “EMBL” release 95. The matched sequences were removed from downstream analysis. UEA sRNA toolkit-Plant version filter pipeline (http://srna-tools.cmp.uea.ac.uk/) was also used to remove any t/rRNA matchouts .
Identification of conserved and potato-specific miRNAs
The unique reads were submitted to the UEA sRNA toolkit-Plant version miRCat pipeline (http://srna-tools.cmp.uea.ac.uk/)  and default criteria were used to predict potato-specific and conserved miRNAs in potato. The reads were aligned to the potato reference genome (PGSC DM assembly version 3) (http://potatogenomics.plantbiology.msu.edu/) with no mismatch. The 100 nt windows around the aligned reads were extracted from the genome and folded using RNAfold (http://rna.tbi.univie.ac.at/cgi-bin/RNAfold.cgi)  to identify precursor sequences. The resulting secondary structures were then trimmed and analysed by miRCat to confirm the characteristic miRNA hairpin structure. These potential precursors were then screened based on updated criteria for annotation of plant miRNAs developed by Meyers et al.  to recover the potential miRNA candidates - 1) The miRNA and miRNA* are derived from opposite stem-arms such that they form a duplex with two nucleotide 3′ overhangs 2) Base-pairing between the miRNA and miRNA* is extensive such that there are typically four or fewer mismatches. 3) The frequency of asymmetric bulges is one or none and size of the bulge is no more than 2 nucleotides within the miRNA/ miRNA* duplex. To identify the conserved miRNAs in potato, all the predicted miRNA sequences were mapped to known mature plant miRNAs deposited in miRBase v20 . No more than 2 mismatches were allowed during the alignment. The remaining sequences that showed 3 or more than 3 mismatches or no homology to any previously known plant miRNAs were considered as novel or potato-specific miRNAs.
Experimental validation of miRNAs by Quantitative Real-time PCR (qRT-PCR)
The expression level of miRNAs was analyzed using poly(T) adaptor RT-PCR method . Total RNA was treated with RNase-free DNase I (NEB, England) and polyadenylated using Poly(A) Tailing kit (Ambion, U.S.A) at 37°C for 60 min in a 50 μl reaction volume containing 2 μg total RNA and 0.08 units poly (A) polymerase. The polyadenylated RNA was then reverse transcribed in a 20 μl reaction mix with SuperScript II Reverse Transcriptase (Invitrogen, U.S.A) and Poly(T) adaptor. Real time PCR was carried out on CFX Connect Real-Time System (Bio-Rad, U.S.A) using hydrolysis probe and probe fast chemistry (KAPA Biosystems, U.S.A). qRT-PCR was performed containing KAPA PROBE FAST Universal 1X qPCR Master Mix (KAPA Biosystems, U.S.A), 0.5 μl cDNA template (10X), 0.5 μM each of the miRNA specific forward primer and universal reverse primer and 0.2 μM universal hydrolysis probe. The cycling conditions were as follows - denaturation at 95°C for 5 min, followed by 40 cycles of denaturation at 95°C for 20 s, annealing and extension together at 60°C for 10 s. The PCR products were further run on 2% agarose gel to verify the presence of a single band of expected size. All the reactions were run with 3 biological replicates of each sample and 2 technical replicates of each biological replicate. Two controls (no template control and no Reverse Transcription control) were included for each miRNA. For normalization, 18S rRNA was used as an endogenous control. The quantitative PCR data was analyzed using 2-∆∆Cq method where Cq is the threshold cycle . The forward primers were designed taking full mature miRNA sequence. The detail of the primers, used in this study, is provided in Additional file 8. All the oligonucleotide sequences were synthesized from IDT (Integrated DNA technologies, U.S.A).
Prediction of miRNA targets and GO analysis
Putative target genes of conserved and potato-specific miRNAs were predicted using psRNATarget: A plant small RNA target analysis server using default parameters . The parameters were as follows - Maximum expectation: 3, Length for complementarity scoring: 20, Allowed maximum energy to unpair the target site (UPE): 25, Flanking length around target site for target accessibility analysis: 17 bp, Range of central mismatch leading to translation inhibition: 9–11 nt. Target with a low expectation score and UPE is considered to be highly potential target for miRNA. The transcript sequences of the potato genome (Group Phureja DM1-3 516R44) were used as a reference set. The putative target genes of the potato-specific miRNAs were further subjected to GO analysis to better understand their function. Blast2GO was used to investigate the GO of the identified potato-specific miRNA target .
Target validation of conserved miRNAs by RLM-RACE
5′ RLM-RACE was used to determine the cleavage site of predicted targets. Briefly, polyA+ RNA was isolated using PolyATract mRNA Isolation System IV (Promega, USA). PolyA+ mRNA (25 ng) was ligated to RNA adapter provided in First Choice RLM-RACE kit (Ambion, USA) using T4 RNA ligase, followed by nested PCR. Gene-specific reverse primers and 5' RACE adapter specific primers were used for amplification of the target gene. The PCR products were eluted, cloned in pGEMT Easy vector (Promega, USA) and sequenced to map the cleavage site. The primers used for amplification were synthesized from IDT (Integrated DNA technologies, U.S.A) and is provided in Additional file 9.
Availability of supporting data
“The data set supporting the results of this article are available in the Gene Expression Omnibus repository under accession no GSE52599 (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/query/acc.cgi?acc=GSE52599)”. The mature miRNA and precursor sequences will be submitted to miRBase registry and assigned final names after final acceptance of the manuscript.
We thank DBT (Department of Biotechnology, India) for funding this research project (No.BT/PR12905/AGR/36/637/2009) to AK. We also thank University of Delhi for sanctioning R & D grant to AK. Council of Scientific and Industrial Research, India provided fellowship to NL. We thank members of our laboratory for their support throughout this project.
- Axtell JM: Classification and comparison of small RNAs from plants. Annu Rev Plant Biol. 2013, 64: 137-159. 10.1146/annurev-arplant-050312-120043.View ArticlePubMedGoogle Scholar
- Chen X: Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009, 25: 21-44. 10.1146/annurev.cellbio.042308.113417.View ArticlePubMedGoogle Scholar
- Rogers K, Chen X: Biogenesis, turnover, and mode of action of plant MicroRNAs. Plant Cell. 2013, 25: 2383-2399. 10.1105/tpc.113.113159.PubMed CentralView ArticlePubMedGoogle Scholar
- Shukla LI, Chinnusamy V, Sunkar R: The role of microRNAs and other endogenous small RNAs in plant stress responses. Biochim Biophys Acta. 2008, 1779 (11): 743-748. 10.1016/j.bbagrm.2008.04.004.View ArticlePubMedGoogle Scholar
- 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 (12): 1282-1290. 10.1038/ng1478.View ArticlePubMedGoogle Scholar
- 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): e219-10.1371/journal.pone.0000219.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Moxon S, Jing R, 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.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuo J, Zhu B, Fu D, Zhu Y, Ma Y, Chi L, Ju Z, Wang Y, Zhai B, Luo Y: Sculpting the maturation, softening and ethylene pathway: The influences of microRNAs on tomato fruits. BMC Genomics. 2012, 13: 7-10.1186/1471-2164-13-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez G, Forment J, Llave C, Pallas V, Gomez G: High-throughput sequencing, characterization and detection of new and conserved cucumber miRNAs. PLoS One. 2011, 6 (5): e19523-10.1371/journal.pone.0019523.PubMed CentralView ArticlePubMedGoogle Scholar
- Li D, Wang L, Liu X, Cui D, Chen T, Zhang H, Jiang C, Xu C, Li P, Li S, Zhao L, Chen H: Deep sequencing of maize small RNAs reveals a diverse set of MicroRNA in dry and imbibed seeds. PLoS One. 2013, 8 (1): e55107-10.1371/journal.pone.0055107.PubMed CentralView ArticlePubMedGoogle Scholar
- Chi X, Yang Q, Chen X, Wang J, Pan L, Chen M, Yang Z, He Y, Liang X, Yu S: Identification and characterization of microRNAs from peanut (Arachis hypogaea L.) by high-throughput sequencing. PLoS One. 2011, 6 (11): e27530-10.1371/journal.pone.0027530.PubMed CentralView ArticlePubMedGoogle Scholar
- Hwang D-G, Park JH, Lim JY, Kim D, Choi Y, Kim S, Reeves G, Yeom SI, Lee JS, Park M, Kim S, Choi IY, Choi D, Shin C: The hot pepper (Capsicum annuum) MicroRNA transcriptome reveals novel and conserved targets: a foundation for understanding MicroRNA functional roles in hot pepper. PLoS ONE. 2013, 8 (5): e64238-10.1371/journal.pone.0064238.PubMed CentralView ArticlePubMedGoogle Scholar
- Song C, Wang C, Zhang C, Korir NK, Yu H, Ma Z, Fang J: Deep sequencing discovery of novel and conserved microRNAs in trifoliate orange (Citrus trifoliata). BMC Genomics. 2010, 11: 431-10.1186/1471-2164-11-431.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 Res. 2008, 18 (4): 571-584. 10.1101/gr.6897308.PubMed CentralView ArticlePubMedGoogle Scholar
- Banerjee AK, Chatterjee M, Yu Y, Suh SG, Miller WA, Hannapel DJ: Dynamics of a mobile RNA of potato involved in a long-distance signaling pathway. Plant Cell. 2006, 18 (12): 3443-3457. 10.1105/tpc.106.042473.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez-Falcon M, Bou J, Prat S: Seasonal control of tuberization in potato: conserved elements with the flowering response. Annu Rev Plant Biol. 2006, 57: 151-180. 10.1146/annurev.arplant.57.032905.105224.View ArticlePubMedGoogle Scholar
- Abelenda JA, Navarro C, Prat S: From the model to the crop: genes controlling tuber formation in potato. Curr Opin Biotechnol. 2011, 22 (2): 287-292. 10.1016/j.copbio.2010.11.013.View ArticlePubMedGoogle Scholar
- Potato Genome Sequencing Consortium: Genome sequence and analysis of the tuber crop potato. Nature. 2011, 475 (7355): 189-195. 10.1038/nature10158.View ArticleGoogle Scholar
- Martin A, Adam H, Diaz-Mendoza M, Zurczak M, Gonzalez-Schain ND, Suárez-López P: Graft-transmissible induction of potato tuberization by the microRNA miR172. Development. 2009, 136 (17): 2873-2881. 10.1242/dev.031658.View ArticlePubMedGoogle Scholar
- Eviatar-Ribak T, Shalit-Kaneh A, Chappell-Maor L, Amsellem Z, Eshed Y, Lifschitz E: A cytokinin-activating enzyme promotes tuber formation in tomato. Curr Biol. 2013, 23 (12): 1057-1064. 10.1016/j.cub.2013.04.061.View ArticlePubMedGoogle Scholar
- Li F, Pignatta D, Bendix C, Brunkard JO, Cohn MM, Tung J, Sun H, Kumar P, Baker B: MicroRNA regulation of plant innate immune receptors. Proc Natl Acad Sci U S A. 2012, 109 (5): 1790-1795. 10.1073/pnas.1118282109.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang W, Luo Y, Gong X, Zeng W, Li S: Computational identification of 48 potato microRNAs and their targets. Comput Biol Chem. 2009, 33 (1): 84-93. 10.1016/j.compbiolchem.2008.07.006.View ArticlePubMedGoogle Scholar
- Xie F, Frazier TP, Zhang B: Identification, characterization and expression analysis of MicroRNAs and their targets in the potato (Solanum tuberosum). Gene. 2010, 473 (1): 8-22.View ArticlePubMedGoogle Scholar
- Yang W, Liu X, Zhang J, Feng J, Li C, Chen J: Prediction and validation of conservative microRNAs of Solanum tuberosum L. Mol Biol Rep. 2010, 37 (7): 3081-3087. 10.1007/s11033-009-9881-z.View ArticlePubMedGoogle Scholar
- Kim HJ, Baek KH, Lee BW, Choi D, Hur CG: In silico identification and characterization of microRNAs and their putative target genes in Solanaceae plants. Genome. 2011, 54 (2): 91-98. 10.1139/G10-104.View ArticlePubMedGoogle Scholar
- Tomato Genome Consortium: The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012, 485 (7400): 635-641. 10.1038/nature11119.View ArticleGoogle Scholar
- Zhang R, Marshall D, Bryan GJ, Hornyik C: Identification and Characterization of miRNA Transcriptome in Potato by High-Throughput Sequencing. PLoS One. 2013, 8 (2): e57233-10.1371/journal.pone.0057233.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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 Biol. 2007, 8 (6): R96-10.1186/gb-2007-8-6-r96.PubMed CentralView ArticlePubMedGoogle Scholar
- Qiu D, Pan X, Wilson IW, Ketchum REB, Li F, Liu M, Teng W, Zhang BH: High-throughput sequencing technology reveals that the taxoid elicitormethyl jasmonate regulates microRNA expression in Chinese yew (Taxus chinensis). Gene. 2009, 436 (1–2): 37-44.View ArticlePubMedGoogle Scholar
- Moxon S, Schwach F, MacLean D, Dalmay T, Studholme DJ, Moulton V: A toolkit for analysing large-scale plant small RNA datasets. Bioinformatics. 2008, 24 (19): 2252-2253. 10.1093/bioinformatics/btn428.View ArticlePubMedGoogle Scholar
- Gruber AR, Lorenz R, Bernhart SH, Neubock R, Hofacker IL: The Vienna RNA websuite. Nucleic Acids Res. 2008, 36: W70-W74. 10.1093/nar/gkn188.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (12): 3186-3190. 10.1105/tpc.108.064311.PubMed CentralView ArticlePubMedGoogle Scholar
- Kozomara A, Griffiths-Jones S: miRBase: integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res. 2011, 39: D152-D157. 10.1093/nar/gkq1027.PubMed CentralView ArticlePubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415. 10.1093/nar/gkg595.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhang B, Pan X, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes. Plant J. 2006, 46 (2): 243-259. 10.1111/j.1365-313X.2006.02697.x.View ArticlePubMedGoogle Scholar
- Bonnet E, Wuyts J, Rouze P, Van de Peer Y: Evidence that microRNA precursors, unlike other non-coding RNAs, have lower folding free energies than random sequences. Bioinformatics. 2004, 20 (17): 2911-2917. 10.1093/bioinformatics/bth374.View ArticlePubMedGoogle Scholar
- Sorefan K, Pais H, Hall AE, Kozomara A, Griffiths-Jones S, Moulton V, Dalmay T: Reducing ligation bias of small RNAs in libraries for next generation sequencing. Silence. 2012, 3 (1): 4-10.1186/1758-907X-3-4.PubMed CentralView ArticlePubMedGoogle Scholar
- Dai X, Zhao PX: psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011, 39: W155-W159. 10.1093/nar/gkr319.PubMed CentralView ArticlePubMedGoogle Scholar
- Song QX, Liu YF, Hu XF, Zhang WK, Ma B, Chen SY, Zhang JS: Identification of miRNAs and their target genes in developing soybean seeds by deep sequencing. BMC Plant Biol. 2011, 11: 5-10.1186/1471-2229-11-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Kolomiets MV, Hannapel DJ, Chen H, Tymeson M, Gladon RJ: Lipoxygenase is involved in the control of potato tuber development. Plant Cell. 2001, 13 (3): 613-626.PubMed CentralView ArticlePubMedGoogle Scholar
- País SM, García MN, Téllez-Iñón MT, Capiati DA: Protein phosphatases type 2A mediate tuberization siganling in Solanum tuberosum L. leaves. Planta. 2010, 232 (1): 37-49. 10.1007/s00425-010-1150-9.View ArticlePubMedGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36 (10): 3420-3435. 10.1093/nar/gkn176.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell. 1993, 75 (5): 843-854. 10.1016/0092-8674(93)90529-Y.View ArticlePubMedGoogle Scholar
- Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687. 10.1016/j.cell.2009.01.046.View ArticlePubMedGoogle Scholar
- Marín-González E, Suárez-López P: “And yet it moves”: Cell-to-cell and long-distance signalling by plant microRNAs. Plant Sci. 2012, 196: 18-30.View ArticlePubMedGoogle Scholar
- Guo HS, Xie Q, Fei JF, Chua NH: MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 2005, 17 (5): 1376-1386. 10.1105/tpc.105.030841.PubMed CentralView ArticlePubMedGoogle Scholar
- Mallory AC, Dugas DV, Bartel DP, Bartel B: MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. Curr Biol. 2004, 14 (12): 1035-1046. 10.1016/j.cub.2004.06.022.View ArticlePubMedGoogle Scholar
- Boualem A, Laporte P, Jovanovic M, Laffont C, Plet J, Combier JP, Niebel A, Crespi M, Frugier F: MicroRNA166 controls root and nodule development in Medicago truncatula. Plant J. 2008, 54 (5): 876-887. 10.1111/j.1365-313X.2008.03448.x.View ArticlePubMedGoogle Scholar
- Williams L, Grigg SP, Xie M, Christensen S, Fletcher JC: Regulation of Arabidopsis shoot apical meristem and lateral organ formation by microRNA miR166g and its AtHD-ZIP target genes. Development. 2008, 132 (16): 3657-3668.View ArticleGoogle Scholar
- Kim J, Jung JH, Reyes JL, Kim YS, Kim SY, Chung KS, Kim JA, Lee M, Lee Y, Kim VN, Chua NH, Park CM: microRNA-directed cleavage of ATHB15 mRNA regulates vascular development in Arabidopsis inflorescence stems. Plant J. 2005, 42 (1): 84-94. 10.1111/j.1365-313X.2005.02354.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Achard P, Herr A, Baulcombe DC, Harberd NP: Modulation of floral development by a gibberellin-regulated microRNA. Development. 2004, 131 (14): 3357-3365. 10.1242/dev.01206.View ArticlePubMedGoogle Scholar
- Mallory AC, Bartel DP, Bartel B: MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell. 2005, 17 (5): 1360-1375. 10.1105/tpc.105.031716.PubMed CentralView ArticlePubMedGoogle Scholar
- Xie KB, Wu CQ, Xiong LZ: Genomic organization, differential expression, and interaction of SQUAMOSA promoter-binding-like transcription factors and microRNA156 in rice. Plant Physiol. 2006, 142 (1): 280-293. 10.1104/pp.106.084475.PubMed CentralView ArticlePubMedGoogle Scholar
- Schwab R, Palatnik JF, Riester M, Schommer C, Schmid M, Weigel D: Specific effects of microRNAs on the plant transcriptome. Dev Cell. 2005, 8 (4): 517-527. 10.1016/j.devcel.2005.01.018.View ArticlePubMedGoogle Scholar
- Cardon GH, Hohmann S, Nettesheim K, Saedler H, Huijser P: Functional analysis of the Arabidopsis thaliana SBP-box gene SPL3: a novel gene involved in the floral transition. Plant J. 1997, 12 (2): 367-377. 10.1046/j.1365-313X.1997.12020367.x.View ArticlePubMedGoogle Scholar
- Hirsch S, Oldroyd GE: GRAS-domain transcription factors that regulate plant development. Plant Signal Behav. 2009, 4 (8): 698-700. 10.4161/psb.4.8.9176.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (11): 2730-2741. 10.1105/tpc.016238.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen X: A microRNA as a translational repressor of APETALA2 in Arabidopsis flower development. Science. 2004, 303 (5666): 2022-2025. 10.1126/science.1088060.View ArticlePubMedGoogle Scholar
- Lauter N, Kampani A, Carlson S, Goebel M, Moose SP: microRNA172 downregulates glossy15 to promote vegetative phase change in maize. Proc Natl Acad Sci U S A. 2005, 102 (26): 9412-9417. 10.1073/pnas.0503927102.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou M, Gu L, Li P, Song X, Wei L, Chen Z, Cao X: Degradome sequencing reveals endogenous small RNA targets in rice (Oryza sativa L. ssp. indica). Front Biol. 2010, 5 (1): 67-90. 10.1007/s11515-010-0007-8.View ArticleGoogle Scholar
- Murashige T, Skoog F: A revised medium for rapid growth and bio assays with tobacco tissue cultures. Physiol Plant. 1962, 15 (3): 473-497. 10.1111/j.1399-3054.1962.tb08052.x.View ArticleGoogle Scholar
- Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005, 39 (4): 519-525. 10.2144/000112010.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.