Identification and characterization of small non-coding RNAs from Chinese fir by high throughput sequencing
- Li-Chuan Wan†1,
- Feng Wang†1, 2,
- Xiangqian Guo3,
- Shanfa Lu4,
- Zongbo Qiu1,
- Yuanyuan Zhao1, 2,
- Haiyan Zhang1Email author and
- Jinxing Lin1Email author
© Wan et al.; licensee BioMed Central Ltd. 2012
Received: 9 October 2011
Accepted: 5 July 2012
Published: 15 August 2012
Small non-coding RNAs (sRNAs) play key roles in plant development, growth and responses to biotic and abiotic stresses. At least four classes of sRNAs have been well characterized in plants, including repeat-associated siRNAs (rasiRNAs), microRNAs (miRNAs), trans-acting siRNAs (tasiRNAs) and natural antisense transcript-derived siRNAs. Chinese fir (Cunninghamia lanceolata) is one of the most important coniferous evergreen tree species in China. No sRNA from Chinese fir has been described to date.
To obtain sRNAs in Chinese fir, we sequenced a sRNA library generated from seeds, seedlings, leaves, stems and calli, using Illumina high throughput sequencing technology. A comprehensive set of sRNAs were acquired, including conserved and novel miRNAs, rasiRNAs and tasiRNAs. With BLASTN and MIREAP we identified a total of 115 conserved miRNAs comprising 40 miRNA families and one novel miRNA with precursor sequence. The expressions of 16 conserved and one novel miRNAs and one tasiRNA were detected by RT-PCR. Utilizing real time RT-PCR, we revealed that four conserved and one novel miRNAs displayed developmental stage-specific expression patterns in Chinese fir. In addition, 209 unigenes were predicted to be targets of 30 Chinese fir miRNA families, of which five target genes were experimentally verified by 5' RACE, including a squamosa promoter-binding protein gene, a pentatricopeptide (PPR) repeat-containing protein gene, a BolA-like family protein gene, AGO1 and a gene of unknown function. We also demonstrated that the DCL3-dependent rasiRNA biogenesis pathway, which had been considered absent in conifers, existed in Chinese fir. Furthermore, the miR390-TAS3-ARF regulatory pathway was elucidated.
We unveiled a complex population of sRNAs in Chinese fir through high throughput sequencing. This provides an insight into the composition and function of sRNAs in Chinese fir and sheds new light on land plant sRNA evolution.
KeywordsChinese fir miRNA rasiRNA tasiRNA Cunninghamia lanceolata
Identification and characterization of diverse classes of small non-coding RNAs (sRNAs) in eukaryotes has been a major research focus in recent years [1, 2]. At least four classes of sRNAs have been well characterized in plants, including heterochromatic and repeat-associated small interfering RNAs (rasiRNAs) , microRNAs (miRNAs) [4, 5], trans-acting siRNAs (tasiRNAs)  and natural antisense transcript-derived siRNAs (nat-siRNAs) . In plants, the majority of endogenous sRNAs are 24-nt rasiRNAs . They repress transposable elements and maintain genome integrity though guiding DNA methylation and histone modification [9–11]. To date, the DCL3-dependent rasiRNA biogenesis pathway has been reported in angiosperms and mosses , but considered absent from conifers . MiRNAs are produced from primary miRNAs (pri-miRNAs) and precursor miRNAs (pre-miRNAs) by DICER LIKE 1 (DCL1) cleavage [14, 15]. Mature miRNAs guide the RNA-induced silencing complexes (RISCs) to degrade target mRNA transcripts  or inhibit their translation . MiRNAs play key roles in plant development, growth, nutrient homeostasis and responses to biotic and abiotic stresses . In Arabidopsis, tasiRNA biogenesis initiates from miRNA-mediated cleavage of a non-coding primary transcript originated from the TAS locus [6, 19]. RNA-dependent RNA polymerase 6 (RDR6) synthesizes long double-stranded RNA (dsRNA) molecules, which are sliced by DCL4 into phased 21-nt tasiRNAs . In Arabidopsis, the miR390-TAS3-AUXIN RESPONSE FACTOR (ARF) pathway executes important functions in leaf development and lateral root formation . Finally, nat-siRNAs are produced through natural antisense transcription and are fashioned by DCL1 or DCL2 . The functional roles for eukaryotic nat-siRNAs that has been described to date are in environmental stress responses and developmental processes .
Identification of conserved and species-specific miRNAs usually relies on two approaches: computational prediction and experimental sequencing. By searching genomic and/or EST databases for orthologous sequences of known miRNAs and analyzing their pre-miRNA hairpin structures, many conserved miRNAs are identified from a variety of plants, such as Arabidopsis thaliana and Brassica napus. High throughput sequencing technologies, such as massively parallel signature sequencing (MPSS), 454 and sequencing-by-synthesis (SBS), have greatly facilitated the discovery of low abundant and recently evolved miRNAs in diverse plants, e.g., Triticum aestivum, Oryza sativa, and Solanum lycopersicum.
There have been reports of investigation and characterization of miRNAs in gymnosperms. Lu et al. identified 37 miRNAs from stem xylem of Pinus taeda, of which 6 miRNAs were likely associated with the fusiform rust gall disease . By sequencing of small RNA libraries generated from a Taxus chinensis cell line, Qiu et al. found that the expression levels of 17 miRNAs have been significantly altered after treatment with methyl jasmonate . Recently, a total of 18 conserved and 53 novel miRNA families were revealed in Pinus contorta. Study on miRNAs in Picea abies indicated that 7 conserved and 9 novel miRNAs participated in the temperature-dependent epigenetic memory and climatic adaptation . These reports showed that like in angiosperms, miRNA-guided post transcriptional gene regulation mechanism is important for the development, growth, stress responses and a myriad of other physiological processes in gymnosperms.
Chinese fir (Cunninghamia lanceolata Lamb. Hook) is one of the most important coniferous evergreen tree species in terms of both industrial and commercial wood supplies in China . However, to date only 407 Chinese fir EST sequences are available in the public databases. No study on Chinese fir sRNAs has been reported so far. To gain mRNA transcriptome sequences of Chinese fir, we recently conducted a high-throughput sequencing of mRNAs isolated from a mixture of tissues. In this study, we used the sequencing-by-synthesis (SBS) technology to sequence a Chinese fir sRNA library and obtained a comprehensive set of sRNAs. Furthermore, we studied expression patterns of conserved and novel miRNAs by qRT-PCR. Potential targets were predicted for most miRNAs, of which 5 target genes have been experimentally verified. Intriguingly, the DCL3-dependent rasiRNA biogenesis pathway which had been thought to be absent in conifers was found in Chinese fir. These results suggest that regulative miRNAs exist in the economically important conifer, Chinese fir, and shed new light on the sRNA evolution from mosses to flowering plants.
Chinese fir has a complex small RNA population
Statistics of sRNAs
Total raw reads
3' adaptor null reads
Insert null reads
5' adaptor contaminant reads
Small insert reads (<18 nt)
Poly(A) sequence reads
Total clean sRNA reads
Unique sequence reads (18–30 nt)
Singleton sequence reads
Unique sequence reads (>2 reads)
Annotations of sRNAs perfectly matching A. thaliana and P. trichocarpa genomes
Annotations of sRNAs perfectly matching Chinese fir mRNA transcriptome
Known miRNA (exact match)
Total clean reads
Although some sRNAs were abundant and presented hundreds of thousands times in our database, the majority of sRNAs were sequenced only a few times. For example, 2,075,775 (74%) sRNAs were sequenced only once, which suggests that Chinese fir contains a large and complex sRNA population. The sRNA singleton rate of Chinese fir (74%) was similar to that of A. thaliana O. sativa, and P. trichocarpa, which is 65% , 82% , and 73% , respectively.
Discovery of conserved miRNAs
To identify conserved miRNAs in Chinese fir, we conducted a local BLASTN search using unique sRNAs (2,815,874) against known plant miRNAs in the miRBase (version 17.0), which contains 3993 miRNAs from 47 plant species (13 April 2011). A total of 115 conserved miRNAs were identified in Chinese fir, belonging to 40 miRNA families (Additional file 2). Ninety-one (79%) conserved miRNAs start with a 5' terminal uridine residue, a conserved feature of miRNAs recognized by the AGO1 protein . In addition, 584 unique sRNAs (1,402,527, 11%) perfectly matched A. thaliana miRNA precursors and 98 sRNAs were identical to A. thaliana mature miRNAs (232), whereas 375 sRNAs (1,383,562, 10%) mapped to P. trichocarpa miRNA precursors and 121 sRNAs exactly matched P. trichocarpa mature miRNAs (234) (Table 2). These figures suggest that approximately half of mature miRNAs are conserved between C. lanceolata and the two model angiosperms.
In the present study, we have tried to find the precursor sequences for the conserved Chinese fir miRNAs. Precursor sequences of 4 miRNAs, including cln-miR162d, cln-miR164b, cln-miR166o and cln-miR1310, have been identified from the Chinese fir mRNA transcriptome database (Additional file 3). Their hairpin structures predicted by MFOLD were shown in Additional file 4. The length of Chinese fir miRNA precursors ranged from 79 to 164 nt, with a majority of which (75%) ranging from 67 to 150 nt. This result is similar with that found in Arabidopsis and rice . The minimal folding free energy indices (MFEIs) of Chinese fir miRNA precursors varied from 0.64 to 1.15, with an average of 0.83, in agreement with that of other plant miRNAs, such as Arabidopsis, rice, Glycine max Medicago truncatula Saccharum officinarum Sorghum bicolor and Zea mays. With complete genome sequences and larger EST databases, more precursor sequences will be found in Chinese fir. In addition, we identified the miRNA* sequences of cln-miR162d, cln-miR164b and cln-miR1310 in at least one of the three small RNA libraries constructed from total RNAs of Chinese fir cambium, providing further evidence that the three miRNAs are canonical miRNAs.
After Blastp searches and further sequence analyses, a variety of unigene sequences in the Chinese fir mRNA transcriptome database were identified homologous to known proteins associated with miRNA biogenesis and action (Additional file 6). For example, 37 and 67 unigene sequences were found highly homologous to Arabidopsis DCL1 and AGO1 protein genes, respectively. Together with the finding of plenty of mature miRNAs and their precursors, these data provide clear evidence that the conserved pathway for miRNA generation and action exists in Chinese fir.
Identification of novel miRNAs
The primary criterion for plant miRNA designation is the precise excision of an ~21-nucleotide miRNA/miRNA* duplex from the stem of a single stranded, stem-loop precursor . Using the algorithm MIREAP, we found one putative novel miRNA, designated as cln-miRn1, in Chinese fir. The mature and precursor sequences of cln-miRn1 were further validated by subcloning (Figure 3A and B). The precursor sequence of cln-miRn1 can form canonical hairpin structure as predicted by MFOLD and has an MFEI value of 1.18, which is similar to that of other plant miRNAs (Additional file 2). The mature sequence of cln-miRn1 also begins with a 5' terminal uridine residue. The Illumina sequencing redundancy of cln-miRn1 is 17, which is much lower than that of many conserved miRNAs. The low redundancy of cln-miRn1 indicates that it may play important roles under certain conditions or in particular developmental stages in Chinese fir. Although the miR* sequence of cln-miRn1 was not identified in the sRNA database of the present study, we found it was in a small RNA library generated from total RNA of Chinese fir cambium. Using BLASTN, we did not find homologous sequence of unigene1397, the primary sequence of cln-miRn1, in the public databases. These results indicate that cln-miRn1 may be a novel and Chinese fir-specific miRNA.
Expression profiles of conserved and novel miRNAs
To provide clues about the physiological functions of sRNAs in Chinese fir development and growth, we utilized semi-quantitative RT-PCR and real time RT-PCR to examine the expression profiles of 7 conserved miRNAs (cln-miR156a, cln-miR157a, cln-miR164a, cln-miR165a, cln-miR166a, cln-miR390 and cln-miR408a), the novel miRNA, cln-miRn1 and cln-tasiR2142 in samples of diverse developmental stages, including seeds, seedlings, leaves and stems. These sRNAs were selected for their key roles in plant development and growth [4, 21, 40–42].
Prediction of miRNA targets
Conserved miRNA targets and their putative functions
Pentatricopeptide repeat-containing protein
BolA-like family protein
Pentatricopeptide repeat-containing protein
Basic blue copper protein
MYB transcription factor
MYB transcription factor
A total of 188 putative targets were not conserved in other plant species. Among them, 55 (26%) targets currently have no functional annotation. We were unable to predict targets for 10 conserved miRNA families. It may be attributable to insufficient Chinese fir mRNA sequences. Four unigene sequences were predicted to be targets of the novel miRNA, cln-miRn1. One of the 4 targets is similar to glycosyltransferase from Lycium barbarum, one is homologous to Cl-channel clc-7 from P. trichocarpa, whereas the other two have no functional annotation. This implies that cln-miRn1 may be involved in multiple physiological and metabolic processes in Chinese fir. More studies need to be performed to elucidate the functions of cln-miRn1 in the growth and development in Chinese fir.
Validation of miRNA-guided cleavage of mRNAs
Plant mature miRNAs can guide RISC complexes to cleave target mRNAs through nucleotide complementarity. The cleavage site usually corresponds to the tenth nucleotide from the 5' end of the miRNA [44, 45]. To verify that miRNAs can regulate their target mRNA expression in Chinese fir, we performed a modified RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE) experiment, using total RNA extracted from 2-month-old seedling stems (see Methods).
Identification of DCL enzymes and phylogenetic analysis
DCL3-dependent rasiRNA generation pathway exists in Chinese fir
It has been well-characterized that DCL3 pays key roles in the biogenesis of 24-nt rasiRNAs, the major endogenous sRNAs in both Arabidopsis and moss Physcomitrella patens. However, reports on small RNAs from gymnosperms, e.g., T. chinensis and P. contorta, showed that they synthesize a diverse population of small RNAs that are 21-nt long, while fail to produce significant amounts of 24-nt sRNAs. Consistently, DCL3, the enzyme responsible for the maturation of 24-nt rasiRNAs in angiosperms, appeared to be absent from conifers . In this study, we found that the length distribution of Chinese fir sRNAs was similar to that of many angiosperms (e.g., A. thaliana and trifoliate orange) in having a major peak at 24 nt. In addition, 3 unigene sequences (unigene34518, unigene35698 and unigene37455) were identified highly homologous to O. sativa P. trichocarpa and P. patens DCL3 proteins. Moreover, we revealed other unigene sequences homologous to known proteins necessary for biogenesis and action of rasiRNAs from the Chinese fir mRNA transcriptome database, including RNA polymerase IV (21 unigene sequences), RDR2 (10 unigene sequences), HEN1 (1 unigene sequence), AGO4 (7 unigene sequences) and AGO6 (4 unigene sequences) (Additional file 10). Besides,12,017 unique sRNAs (736,626) were able to map to the sense or antisense strand of known repeat-associated RNA sequences using BLASTN analysis against the Repeat database (http://www.repeatmasker.org/cgi-bin/WEBRepeatMasker) (Table 3), indicating that these sRNAs may originate from pericentromeric regions, transposons or repetitive elements. Taken together, these results suggest that DCL3-dependent rasiRNA generation pathway may exist in Chinese fir, a gymnosperm species.
TasiRNAs belong to a plant-specific class of endogenous sRNAs. Previous studies showed that one of the conserved tasiRNAs, TAS3, plays key roles in mosses and flowering plants . However, TAS3-related genes have not been reported in gymnosperms so far. In the present investigation, we found not only one MIR390 family member, but also three small RNA sequences identical to Arabidopsis TAS3 siRNAs, including siR1778, siR1769 and tasiR2142, in the Chinese fir sRNA library (Additional file 2). To confirm the expression of cln-miR390a and cln-tasiR2142 in Chinese fir, we performed relative real-time and semi-quantitative RT-PCR analyses. The transcription levels of cln-tasiR2142 were low and did not show significant differences in the four samples examined (Figure 3). Cln-miR390a was weakly expressed in seeds, seedling and stems and very weakly in leaves (Figures 3 and 4). In addition, we identified unigene sequences homologous to known proteins necessary for tasiRNA biogenesis and action, e.g., RDR6 (11 unigene sequences), DCL4 (29 unigene sequences) and AGO7 (15 unigene sequences) from our mRNA transcriptome database (Additional file 11). Using the Web-based target prediction program psRNATarget and Chinese fir mRNA transcriptome database and 61,864 contigs from P. taeda (DFCI Pine Gene Index release 7.0) as the target databases, we found 2 contigs (TC109284 and TC85135) annotated as ARFs in P. taeda to be targets of cln-tasiR2142. Thus, our study provides strong evidence that this crucial and conserved miR390-TAS3-ARF regulatory pathway may exist and play important roles in Chinese fir.
sRNA length distribution divergence among gymnosperms
Previous studies showed the sRNA length distribution patterns of many angiosperms, e.g., Arabidopsis thaliana, Citrus trifoliata and Medicago truncatula, and mosses  have a major peak at 24 nt. However, the length distribution patterns of sRNAs from two gymnosperms, Taxus chinensis and Pinus contorta, constitute a major 21-nt peak and no obvious 24-nt peak is observed. A survey of sRNAs in 24 vascular plant species revealed that 7 conifer species fail to produce significant amounts of 24-nt sRNAs . Therefore, it seems that the sRNA length distribution patterns are significantly divergent between angiosperms and gymnosperms.
In the present study, we found that Chinese fir sRNA length distribution pattern was similar to that of many angiosperms, having a 24-nt major peak. We speculate that the divergence in sRNA length distribution pattern between Chinese fir and other gymnosperms may be due to two reasons. One reason is the cell or tissue type difference of the initial sample used to generate the small RNA library. Multiple tissues and organs in different developmental stages, including seeds, seedlings, leaves, stems and calli, were used as the initial sample in Chinese fir small RNA library construction, whereas only one tissue, a stem cell line of T. chinensis or young needles of Pinus contorta or young needle tips or leaf buds of 10 gymnosperms , were utilized as the initial samples in their experiments. Young needles, leaf buds or stems mainly consist of terminally differentiated cells, while our sample contains not only terminally differentiated cells, but also differentiating cells, such as cells in seeds and seedlings. Therefore, the divergence of sRNA length distribution patterns between Chinese fir and other gymnosperms may be partially caused by the sRNA population difference between differentiating and differentiated cells. The other reason that cannot be excluded is the species differences between Chinese fir and other gymnosperms.
Conserved miRNAs in Chinese fir
In comparison with thousands of miRNAs identified from angiosperms, the reported miRNAs from gymnosperms are limited. Only miRNAs of two conifer species, Pinus taeda and Picea abies, are listed in the miRBase (version 17.0) with a total of 40 and 37 miRNAs respectively. No research has been conducted to identify miRNAs from Chinese fir, an economically important gymnosperm in China. In this report, we conducted the first screen for Chinese fir miRNAs by deep sequencing. A total of 115 conserved miRNAs belonging to 40 miRNA families were identified in Chinese fir. For the 21 common conserved miRNA families among three model angiosperms, A. thaliana O. sativa and P. trichocarpa, 20 miRNA families were found in Chinese fir, except the MIR397 family. These results imply that the ancient miRNA biogenesis and action system is well developed in the common ancestors of gymnosperms and angiosperms.
It is well-characterized that miR156 regulates the timing of the juvenile-to-adult transition by repressing the expression of SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL) family of transcription factor genes in Arabidopsis. Ath-miR156 is also involved in regulating the onset of flowering , temporal control of trichome distribution  and securing male fertility in Arabidopsis. In the present study, we identified 18 miRNA members of the cln-MIR156/157 family by Illumina sequencing, of which cln-miR157a is the most abundant miRNA in the sRNA library. Using semi-quantitative and quantitative RT-PCR, we found that the transcript levels of cln-miR156a and cln-miR157a were down-regulated along the juvenile-to-adult transition (Figure 3 and Figure 4). Moreover, we predicted 40 mRNAs to be targets of cln-MIR156 family miRNAs (Table 4 and Additional file 7). The analysis of cln-miR156 targets showed that 2 targets were SPL family of transcription factor genes, including unigene2030 and unigene2872. In addition, unigene2872, encoding a protein homologous to the Arabidopsis transcription factor SPL7, was experimentally verified to be the target of cln-miR157a by 5' RACE. These results suggest that cln-miR156a and cln-miR157a may regulate juvenile-to-adult transition by down-regulating the expression of SPL transcription factor genes in Chinese fir.
In rice, PETER PAN SYNDROME (PPS), an Arabidopsis CONSTITUTIVE PHOTOMORPHOGENIC1 (COP1) ortholog, regulates the juvenile-to-adult and vegetative-to-reproductive phase changes by controlling the expression of osa-miR156 . In the present study, unigene9100, an Arabidopsis COP1 ortholog, was predicted to be the target of cln-miR156, indicating that COP1-miR156 pathway may be self-regulated during the phase change in Chinese fir development. Furthermore, mRNAs involving in defense responses (e.g., unigene57989) and other metabolic pathways (e.g., unigene8339) were also predicted to be targets of cln-miR156. These results indicate that cln-miR156/157 family miRNAs may be involved in multiple physiological processes in Chinese fir development.
The functions of miR168 have been intensively studied in previous research. Várallyay et al. reported that plant viruses repress AGO1 accumulation by modulating endogenous miR168 level to alleviate the anti-viral function of AGO1 protein . Recently, Li et al. revealed a transcriptional regulatory mechanism by which miR168 controls AGO1 homeostasis during ABA treatment and abiotic stress responses in Arabidopsis thaliana. In the present study, we found not only 2 MIR168 family members, but also unigene6526, an ortholog of Arabidopsis AGO1, to be the target of cln-miR168a using 5' RACE. These results imply that cln-miR168a may have important roles in viral infection, stress responses and other physiological processes by repressing the expression of AGO1 in Chinese fir.
Previous studies have established that miR164 targets NAC family transcription factor is well conserved in many plant species and is involved in the trifurcate feed-forward pathway to ensure aging-induced leaf senescence  and in the lateral root development in A. thaliana. In the present experiment, we identified 5 MIR164 family members and detected developmental stage-associated expression pattern of cln-miR164a in four Chinese fir samples by qRT-PCR. Moreover, unigene7992, a NAC family transcription factor, was verified to be the target of cln-miR164a. These results provide direct evidence that cln-miR164a may participate in multiple physiological processes during development and growth by inhibiting the expression of NAC transcription factors in Chinese fir.
DCL3-dependent rasiRNA biogenesis pathway
The DCL3-dependent rasiRNA generation pathway has been reported in both angiosperms and mosses [8, 12]. However, this pathway has been considered absent from gymnosperms because of the lack of significant amounts of 24-nt sRNAs and DCL3 in conifers .
The experiments described here demonstrated that 24-nt sRNAs were the most abundant size class of sRNAs in Chinese fir. In addition, we found that 736,626 sRNAs could map to the sense or antisense strand of known repeat-associated DNA sequences, suggesting that these sRNAs may be transcribed from pericentromeric regions, transposons or repetitive elements. Moreover, unigene sequences homologous to the proteins necessary for the biogenesis and action of the 24-nt rasiRNAs were identified from the Chinese fir mRNA transcriptome database, including Pol IV, RDR2, DCL3, HEN1, AGO4 and AGO6 (Figure 6 and Additional file 10). These results indicate that the ancient but crucial DCL3-dependent rasiRNA generation pathway may exist in Chinese fir, a gymnosperm species.
In rice, two DCL3 family members have distinct roles in the generation of 24-nucleotide small RNAs . OsDCL3a plays conserved roles in producing 24-nucleotide unphased small RNAs, which can direct cytosine DNA methylation both in cis and in trans. OsDCL3b specifically acts in 24-nucleotide phased small RNA biogenesis, which exhibits panicle- and early seed-specific expression and plays roles in stamens. In Physcomitrella patens, DCL3 is required for the accumulation of 22–24 nt siRNAs, but not 21 nt siRNAs, at Pp23SR loci . In the present investigation, we identified 3 unigene sequences, including unigene34518, unigene35698 and unigene37455, highly homologous to rice DCL3a, Physcomitrella patens DCL3, and Populus trichocarpa DCL3, respectively, suggesting that unigene34518 may be responsible for the biogenesis of 24-nucleotide unphased small RNAs, whereas unigene35698 and unigene37455 may be involved in the producing other classes of 24-nucleotide small RNAs in Chinese fir.
Using deep-sequencing technology, we obtained a comprehensive set of sRNAs in Chinese fir, including rasiRNAs, conserved and novel miRNAs and tasiRNAs. A total of 115 conserved miRNAs comprising 40 miRNA families and one novel miRNA were identified. The expressions of 16 conserved and one novel miRNAs and one tasiRNA were examined by RT-PCR. Real-time quantitative PCR results demonstrated that four conserved and one novel miRNAs displayed developmental stage-specific expression patterns. In addition, 209 unigene sequences were predicted to be the targets of 30 Chinese fir miRNA families, of which five unigene sequences were experimentally verified by 5' RACE. The DCL3-dependent rasiRNA generation pathway, which had been considered absent in conifers, was found in Chinese fir. Furthermore, the miR390-TAS3-ARF regulatory pathway was investigated. In conclusion, global identification and characterization of sRNAs from Chinese fir expand our knowledge of conifer sRNAs and provide an insight into land plant sRNA evolution.
Chinese fir plants were grown under standard greenhouse conditions. Seedlings, adult leaves and stems were harvested and stored at −80°C until use. Dry seeds were stored at 4°C. Calli derived from immature Chinese fir seeds were maintained on a Chinese fir-specific culture medium.
Small RNA isolation and Illumina sequencing
Total RNAs of seeds and calli were extracted using RNAiso-mate for plant tissue and RNAiso plus (Takara, Dalian, Liaoning, China), whereas total RNAs of seedlings, adult leaves and stems were isolated with the Concert Plant RNA Reagent (Invitrogen, Carlsbad, CA, USA), and were then treated with RNase-free DNase I (Promega, Madison, WI, USA). Equal amount of total RNAs from the 5 different samples were mixed to form a single RNA pool. Twenty micrograms of total RNAs from the pool were used and 16 to 30-nt sRNAs were purified using Novex 15% TBE-Urea gel (Invitrogen). Two adaptors were sequentially ligated to the 5' and 3' ends of purified sRNAs. The ligation products were further purified from Novex 10% TBE-Urea gel. Reverse transcriptase SuperScript II (Invitrogen) and high-fidelity DNA polymerase Phusion (New England Biolabs, Ipswich, MA, USA) were used in the following RT-PCR reaction. The amplification products were cut from Novex 6% TBE-Urea gel. The purified DNA fragments were used for sequencing on an Illumina 1 G Genome Analyzer at the Beijing Genomics Institute, Shenzhen, China.
SRNA sequence processing
The raw data were processed with the BGI sRNA analysis pipeline to filter out artifact sequences. Non-redundant sRNAs ranging from 18 to 30 nt were collected and stored in the Clean file and reads of unique sequences were recorded and submitted to NCBI [GEO: GSE24226]. Unique sRNAs in the Clean file were mapped to Chinese fir mRNA transcriptome database [GenBank:SRA053525] (Additional file 12) and A. thaliana and P. trichocarpa genomes using SOAP, according to its default settings (http://www.soap.genomics.org.cn). The A. thaliana genome sequences and their annotations were downloaded from the TIGR website (ftp://ftp.tigr.org/); The P. trichocarpa genome sequences and their annotations were downloaded from the JGI website (http://genome.jgi-psf.org/Poptr1_1/Poptr1_1.download.ftp.html). rRNA, tRNA, snoRNA and snRNA sequences were downloaded from NCBI and Rfam 9.0 (http://rfam.sanger.ac.uk/), and coordinates of genomic repeats were obtained from RepeatMasker (http://www.repeatmasker.org/PreMaskedGenomes.html). The perfectly aligned sRNA was annotated as rRNA/tRNA/snRNA/snoRNA, miRNA, repeat element, exon-sense, exon-antisense, intron-sense or intron-antisense and stored in separate files, based on the annotation of the sequence it overlapped in Chinese fir transcriptome, A. thaliana or P. trichocarpa genome. Matching sRNAs without annotation were stored in the Un-annotated file.
Identification of conserved and novel miRNAs
To identify conserved miRNAs, unique sRNAs (2,815,874) from the sRNA library and contigs (525,706) from the Chinese fir mRNA transcriptome database were utilized in local BLASTN analyses (E value was set to 0.01 and mismatches were set to less than 3) against the mature and precursor sequences of miRNAs in miRBase version 17.0 (http://www.mirbase.org/) . Overlapping contig sequences were used to form longer sequences according to their alignments to known miRNA precursor sequences in the miRBase. Sequences homologous to known miRNA precursors were annotated as candidate Chinese fir miRNA precursors. To further reveal conserved and novel miRNA precursors, the unique sRNAs (2,815,874) were aligned to the Chinese fir mRNA transcriptome database using MIREAP with default parameters (http://sourceforge.net/projects/mireap/). BLASTN searches against all nucleotide sequences in NCBI databases were performed to investigate whether these potential miRNA precursors were conserved in other plant species. Putative precursors homologous to known plant rRNAs, tRNAs or mRNAs were excluded. Mfold was used to predict the secondary structures of the candidate miRNA precursor sequences, utilizing default parameters (http://mfold.bioinfo.rpi.edu/cgi-bin/rna-form1.cgi) . The minimal folding free energy (MFE) of the sequence was set to less than or equal to −30 kcal/mol. MiRNA precursors with minimal folding free energy index (MFEI) values less than 0.60 were also discarded . Only the perfectly matched sRNA sequences and homologous sequences with precursor sequences were considered to be conserved miRNAs. Sequence with proper secondary hairpin structure, a MFEI value more than 0.60 and no homologous sequence in public databases was considered as putative novel miRNA precursor sequence.
Semi-quantitative and real-time quantitative RT-PCR
Total RNAs were extracted separately from seeds, 2-month-old seedlings, adult leaves and stems as described above. For mature miRNA expression analysis, cDNAs were synthesized from 1 μg of purified total RNAs with the NCode miRNA First-Strand cDNA Synthesis Kit (MIRC-50; Invitrogen). Forward primers were designed based on mature miRNA sequences (Additional file 13). If the Tm of a mature miRNA was <60°C, it was adjusted by adding Gs or Cs to the 5' end and/or As to the 3' end of the miRNA sequence . A primer corresponding to a 20-bp segment at the 3' end of a Chinese fir 5.8 s rRNA gene was used as a reference control. We adopted stringent annealing conditions and set the annealing temperature to 65°C for quantitative real-time RT-PCR reactions because many miRNA paralogs differ by only one nucleotide . Relative real-time quantitative PCR was performed using Toyobo’s Thunderbird SYBR qPCR mix in a 20-μl reaction volume, using 2 μl of 1:10 diluted cDNA solution as template and 0.3 μM of each primer. Triplicate reactions were conducted on a quantitative PCR machine (MX3000P, Stratagene, La Jolla, CA, USA), using the following thermal cycling conditions: 95°C for 1 min, 45 cycles of 95°C for 15 s, 65°C for 15 s, and 72°C for 5 s.
For miRNA precursor expression analysis, cDNAs were synthesized from 2 μg of purified total RNAs in 25-μl reactions containing 200 U M-MLV reverse transcriptase (Promega) and 1 μg random nonamers, according to the manufacturer’s protocol. Six pairs of primers for 5 housekeeping genes, including actin, ubiquitin, gapdh, tubulin and translation factor EF-1 alpha-like gene and two pairs of primers for cln-miR164b and cln-miRn1 precursors were designed. Relative real-time quantitative PCR was carried out as described above using the following thermal cycling conditions: 95°C for 1 min, 40 cycles of 95°C for 15 s, 55°C for 15 s and 72°C for 10 s. Among the 5 housekeeping genes, gapdh had the most similar expression profiles in the four samples and was used as a reference gene.
After PCR, a thermal denaturing cycle was carried out to determine the dissociation curve and verify the specificity of the amplification. The amplification results were analyzed using a comparative Ct method, which used an arithmetic formula, 2-ΔΔCt. Ct represents the threshold cycle. All expression levels were normalized to the arithmetic mean of the selected 5.8 s rRNA gene (for mature miRNAs) or gapdh (for precursor). The expression level in seeds or leaves was arbitrarily set to 1.
Semi-quantitative RT-PCR was conducted with the following thermal cycling parameters: 95°C for 1 min, 95°C for 15 s, 65°C for 15 s, and 72°C for 5 s for 40 cycles. The final concentration of primers was 0.3 μM. We designed 17 forward primers for 16 conserved miRNAs and tasiR2142. Amplification products were separated with 2.5% agarose gel electrophoresis. Three RT-PCR replications for each primer were conducted with independently isolated total RNAs.
To further confirm the specificity of the amplification, we analyzed PCR samples on 2.5% agarose gels with EtBr visualization of bands. Fragments were gel-purified, cloned into the pGEM-T Easy vector (Promega) and sequenced.
Target mRNA prediction
We used the Web-based psRNATarget program to identify putative targets for conserved and novel miRNAs (http://bioinfo3.noble.org/psRNATarget/). The custom plant transcript databases include 59,669 unigene sequences from the Chinese fir mRNA transcriptome database and 61,864 contigs from P. taeda (DFCI Pine Gene Index release 7.0, http://compbio.dfci.harvard.edu/cgi-bin/tgi/gimain.pl?gudb=pine). A scoring system was applied according to Zhang . Sequences with a penalizing score ≤4 were chosen as putative targets.
MiRNA-mediated cleavage of mRNA
For identification of internal cleavage sites in the target mRNAs, a modified RNA ligase-mediated rapid amplification of cDNA ends (RLM-RACE) experiment was carried out using a 5' RACE kit (Takara) [49, 63]. Total RNA was extracted from 2-month-old seedling stems and purified as described above. An RNA Oligo adapter was directly ligated to the purified RNAs (2000 ng) without calf intestinal phosphatase and tobacco acid pyrophosphatase treatment. Thirty nesting and nested gene-specific primers were synthesized and used for PCR amplifications. Ten DNA bands with expected sizes were gel purified and cloned into the pGEM-T Easy vector and sequenced. Five of the sequenced DNA bands were found to be miRNA-guided cleavage products.
Dicer alignment and phylogenetic analysis
DCL sequences from A. thaliana and P. patens were obtained from GenBank (http://www.ncbi.nlm.nih.gov/Genbank/index.html). O. sativa DCL sequences were accessed via the Institute for Genomic Research (TIGR) rice database (http://www.tigr.org/tigr-scripts/osa1_web/gbrowse/rice). Accession numbers are as follows: A. thaliana: NP_171612, NP_566199, NP_189978, NP_197532/DQ118423; O. sativa: Os03g02970, Os03g38740, Os09g14610, Os01g68120, Os10g34430, Os04g43050; P. patens: ABV31244, ABV31245, ABV31246.
Protein domains of the Chinese fir DCLs were analyzed by scanning them against the InterPro protein signature database (http://www.ebi.ac.uk/InterProScan). Domains were assigned according to pFAM predictions.
Dicer protein alignments were conducted using Bioedit with default parameters. The MEGA program was used to generate the phylogenetic tree of Dicers in A. thaliana, O. sativa, P. patens and Chinese fir, using the neighbor-joining algorithm with 1000 rounds of bootstrapping.
Hua Enhancer 1
Minimal folding free energy index
Rapid amplification of cDNA ends
repeat-associated small interfering RNA
RNA-dependent RNA polymerase
RNA-induced silencing complex
RNA ligase-mediated rapid amplification of cDNA ends
small nuclear RNA
small nucleolar RNA
trans-acting small interfering RNA
We thank Prof. Qing-Yin Zeng from Institute of Botany, Chinese Academy of Sciences and Dr. Sam Griffiths-Jones from University of Manchester for their helpful suggestion. This research was financially supported by the National High Technology Research and Development Program of China (863 Program) (2008AA10Z130), Major State Basic Research Development Program by Ministry of Science and Technology of China (2009CB118500 and 2009CB119105), the National Natural Science Foundation of China (NSFC) (30770123), and projects from Ministry of Agriculture of China (2009ZX08009-011B and 2009ZX08009-095B).
- Bartel DP: MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004, 116 (2): 281-297.PubMedView ArticleGoogle Scholar
- Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136 (4): 669-687.PubMedView ArticleGoogle Scholar
- Herr AJ, Jensen MB, Dalmay T, Baulcombe DC: RNA polymerase IV directs silencing of endogenous DNA. Science. 2005, 308 (5718): 118-120.PubMedView ArticleGoogle Scholar
- Wu G, Park MY, Conway SR, Wang JW, Weigel D, Poethig RS: The sequential action of miR156 and miR172 regulates developmental timing in Arabidopsis. Cell. 2009, 138 (4): 750-759.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu C, Jeong DH, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, et al: Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA. 2008, 105 (12): 4951-4956.PubMedPubMed CentralView ArticleGoogle Scholar
- Adenot X, Elmayan T, Lauressergues D, Boutet S, Bouché N, Gasciolli V, Vaucheret H: DRB4-dependent TAS3 trans-acting siRNAs control leaf morphology through AGO7. Curr Biol. 2006, 16 (9): 927-932.PubMedView ArticleGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005, 123 (7): 1279-1291.PubMedPubMed CentralView ArticleGoogle Scholar
- Zheng B, Wang Z, Li S, Yu B, Liu JY, Chen X: Intergenic transcription by RNA polymerase II coordinates Pol IV and Pol V in siRNA-directed transcriptional gene silencing in Arabidopsis. Genes Dev. 2009, 23 (24): 2850-2860.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith LM, Pontes O, Searle I, Yelina N, Yousafzai FK, Herr AJ, Pikaard CS, Baulcombe DC: An SNF2 protein associated with nuclear RNA silencing and the spread of a silencing signal between cells in Arabidopsis. Plant Cell. 2007, 19 (5): 1507-1521.PubMedPubMed CentralView ArticleGoogle Scholar
- Qi Y, He X, Wang XJ, Kohany O, Jurka J, Hannon GJ: Distinct catalytic and non-catalytic roles of ARGONAUTE4 in RNA-directed DNA methylation. Nature. 2006, 443 (7114): 1008-1012.PubMedView ArticleGoogle Scholar
- Henderson IR, Jacobsen SE: Tandem repeats upstream of the Arabidopsis endogene SDC recruit non-CG DNA methylation and initiate siRNA spreading. Genes Dev. 2008, 22 (12): 1597-1606.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho SH, Addo-Quaye C, Coruh C, Arif MA, Ma Z, Frank W, Axtell MJ: Physcomitrella patens DCL3 is required for 22–24 nt siRNA accumulation, suppression of retrotransposon-derived transcripts, and normal development. PLoS Genet. 2008, 4 (12): e1000314-PubMedPubMed CentralView ArticleGoogle Scholar
- Dolgosheina EV, Morin RD, Aksay G, Sahinalp SC, Magrini V, Mardis ER, Mattsson J, Unrau PJ: Conifers have a unique small RNA silencing signature. RNA. 2008, 14 (8): 1508-1515.PubMedPubMed CentralView ArticleGoogle Scholar
- Park MY, Wu G, Gonzalez-Sulser A, Vaucheret H, Poethig RS: Nuclear processing and export of microRNAs in Arabidopsis. Proc Natl Acad Sci USA. 2005, 102 (10): 3691-3696.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie Z, Allen E, Fahlgren N, Calamar A, Givan SA, Carrington JC: Expression of Arabidopsis MIRNA genes. Plant Physiol. 2005, 138 (4): 2145-2154.PubMedPubMed CentralView ArticleGoogle Scholar
- Baumberger N, Baulcombe DC: Arabidopsis ARGONAUTE1 is an RNA slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci USA. 2005, 102 (33): 11928-11933.PubMedPubMed CentralView ArticleGoogle Scholar
- Brodersen P, Sakvarelidze-Achard L, Bruun-Rasmussen M, Dunoyer P, Yamamoto YY, Sieburth L, Voinnet O: Widespread translational inhibition by plant miRNAs and siRNAs. Science. 2008, 320 (5880): 1185-1190.PubMedView ArticleGoogle Scholar
- Abdel-Ghany SE, Pilon M: MicroRNA-mediated systemic down-regulation of copper protein expression in response to low copper availability in Arabidopsis. J Biol Chem. 2008, 283 (23): 15932-15945.PubMedPubMed CentralView ArticleGoogle Scholar
- Axtell MJ, Snyder JA, Bartel DP: Common functions for diverse small RNAs of land plants. Plant Cell. 2007, 19 (6): 1750-1769.PubMedPubMed CentralView ArticleGoogle Scholar
- Schwab R, Maizel A, Ruiz-Ferrer V, Garcia D, Bayer M, Crespi M, Voinnet O, Martienssen RA: Endogenous tasiRNAs mediate non-cell autonomous effects on gene regulation in Arabidopsis thaliana. PLoS One. 2009, 4 (6): e5980-PubMedPubMed CentralView ArticleGoogle Scholar
- Montgomery TA, Howell MD, Cuperus JT, Li D, Hansen JE, Alexander AL, Chapman EJ, Fahlgren N, Allen E, Carrington JC: Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation. Cell. 2008, 133 (1): 128-141.PubMedView ArticleGoogle Scholar
- Katiyar-Agarwal S, Gao S, Vivian-Smith A, Jin H: A novel class of bacteria-induced small RNAs in Arabidopsis. Genes Dev. 2007, 21 (23): 3123-3134.PubMedPubMed CentralView ArticleGoogle Scholar
- Ron M, Alandete Saez M, Eshed Williams L, Fletcher JC, McCormick S: Proper regulation of a sperm-specific cis-nat-siRNA is essential for double fertilization in Arabidopsis. Genes Dev. 2010, 24 (10): 1010-1021.PubMedPubMed CentralView ArticleGoogle Scholar
- Adai A, Johnson C, Mlotshwa S, Archer-Evans S, Manocha V, Vance V, Sundaresan V: Computational prediction of miRNAs in Arabidopsis thaliana. Genome Res. 2005, 15 (1): 78-91.PubMedPubMed CentralView ArticleGoogle Scholar
- Xie FL, Huang SQ, Guo K, Xiang AL, Zhu YY, Nie L, Yang ZM: Computational identification of novel microRNAs and targets in Brassica napus. FEBS Lett. 2007, 581 (7): 1464-1474.PubMedView ArticleGoogle Scholar
- Yao Y, Guo G, Ni Z, Sunkar R, Du J, Zhu JK, Sun Q: Cloning and characterization of microRNAs from wheat (Triticum aestivum L.). Genome Biol. 2007, 8 (6): R96-PubMedPubMed CentralView ArticleGoogle Scholar
- Zhu QH, Spriggs A, Matthew L, Fan L, Kennedy G, Gubler F, Helliwell C: A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res. 2008, 18 (9): 1456-1465.PubMedPubMed CentralView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Lu S, Sun YH, Amerson H, Chiang VL: MicroRNAs in loblolly pine (Pinus taeda L.) and their association with fusiform rust gall development. Plant J. 2007, 51 (6): 1077-1098.PubMedView ArticleGoogle Scholar
- Qiu D, Pan X, Wilson IW, Li F, Liu M, Teng W, Zhang B: 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.PubMedView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Yakovlev IA, Fossdal CG, Johnsen Ø: MicroRNAs, the epigenetic memory and climatic adaptation in Norway spruce. New Phytol. 2010, 187 (4): 1154-1169.PubMedView ArticleGoogle Scholar
- Shi J, Zhen Y, Zheng RH: Proteome profiling of early seed development in Cunninghamia lanceolata (Lamb.) Hook. J Exp Bot. 2010, 61 (9): 2367-2381.PubMedPubMed CentralView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Barakat A, Wall PK, Diloreto S, Depamphilis CW, Carlson JE: Conservation and divergence of microRNAs in Populus. BMC Genomics. 2007, 8: 481-PubMedPubMed CentralView ArticleGoogle Scholar
- 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: e141-PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang B, Pan X, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes. Plant J. 2006, 46: 243-259.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Cox SB, Cobb GP, Anderson TA: Evidence that miRNAs are different from other RNAs. Cell Mol Life Sci. 2006, 63 (2): 246-254.PubMedView ArticleGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen X, et al: Criteria for annotation of plant MicroRNAs. Plant Cell. 2008, 20 (12): 3186-3190.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim JH, Woo HR, Kim J, Lim PO, Lee IC, Choi SH, Hwang D, Nam HG: Trifurcate feed-forward regulation of age-dependent cell death involving miR164 in Arabidopsis. Science. 2009, 323 (5917): 1053-1057.PubMedView ArticleGoogle 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.PubMedView ArticleGoogle Scholar
- Maunoury N, Vaucheret H: AGO1 and AGO2 act redundantly in miR408-mediated Plantacyanin regulation. PLoS One. 2011, 6 (12): e28729-PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Y: miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 2005, 33 (Web Server issue): W701-704.PubMedPubMed CentralView ArticleGoogle Scholar
- Llave C, Xie Z, Kasschau KD, Carrington JC: Cleavage of Scarecrow-like mRNA targets directed by a class of Arabidopsis miRNA. Science. 2002, 297: 2053-2056.PubMedView ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP: Computational identification of plant microRNAs and their targets, including a stress-induced miRNA. Mol Cell. 2004, 14: 787-799.PubMedView ArticleGoogle Scholar
- Ronemus M, Vaughn MW, Martienssen RA: MicroRNA-targeted and small interfering RNA-mediated mRNA degradation is regulated by argonaute, dicer, and RNA-dependent RNA polymerase in Arabidopsis. Plant Cell. 2006, 18 (7): 1559-1574.PubMedPubMed CentralView ArticleGoogle Scholar
- De Paola D, Cattonaro F, Pignone D, Sonnante G: The miRNAome of globe artichoke: conserved and novel micro RNAs and target analysis. BMC Genomics. 2012, 13 (1): 41-PubMedPubMed CentralView ArticleGoogle Scholar
- Xie Z, Johansen LK, Gustafson AM, Kasschau KD, Lellis AD, Zilberman D, Jacobsen SE, Carrington JC: Genetic and functional diversification of small RNA pathways in plants. PLoS Biol. 2004, 2 (5): E104-PubMedPubMed CentralView ArticleGoogle 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-PubMedPubMed CentralView ArticleGoogle Scholar
- Lelandais-Brière C, Naya L, Sallet E, Calenge F, Frugier F, Hartmann C, Gouzy J, Crespi M: Genome-wide Medicago truncatula small RNA analysis revealed novel microRNAs and isoforms differentially regulated in roots and nodules. Plant Cell. 2009, 21 (9): 2780-2796.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang JW, Park MY, Wang LJ, Koo Y, Chen XY, Weigel D, Poethig RS: MiRNA control of vegetative phase change in trees. PLoS Genet. 2010, 7 (2): e1002012-View ArticleGoogle Scholar
- Yant L, Mathieu J, Dinh TT, Ott F, Lanz C, Wollmann H, Chen X, Schmid M: Orchestration of the floral transition and floral development in Arabidopsis by the bifunctional transcription factor APETALA2. Plant Cell. 2010, 22 (7): 2156-2170.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu N, Cai WJ, Wang S, Shan CM, Wang LJ, Chen XY: Temporal control of trichome distribution by microRNA156-targeted SPL genes in Arabidopsis thaliana. Plant Cell. 2010, 22 (7): 2322-2335.PubMedPubMed CentralView ArticleGoogle Scholar
- Xing S, Salinas M, Höhmann S, Berndtgen R, Huijser P: miR156-targeted and nontargeted SBP-box transcription factors act in concert to secure male fertility in Arabidopsis. Plant Cell. 2010, 22 (12): 3935-3950.PubMedPubMed CentralView ArticleGoogle Scholar
- Tanaka N, Itoh H, Sentoku N, Kojima M, Sakakibara H, Izawa T, Itoh J, Nagato Y: The COP1 ortholog PPS regulates the juvenile-adult and vegetative-reproductive phase changes in rice. Plant Cell. 2011, 23 (6): 2143-2154.PubMedPubMed CentralView ArticleGoogle Scholar
- Várallyay E, Válóczi A, Agyi A, Burgyán J, Havelda Z: Plant virus-mediated induction of miR168 is associated with repression of ARGONAUTE1 accumulation. EMBO J. 2010, 29 (20): 3507-3519.PubMedPubMed CentralView ArticleGoogle Scholar
- Li W, Cui X, Meng Z, Huang X, Xie Q, Wu H, Jin H, Zhang D, Liang W: Transcriptional regulation of Arabidopsis MIR168a and ARGONAUTE1 homeostasis in ABA and abiotic stress responses. Plant Physiol. 2012, 158 (3): 1279-1292.PubMedPubMed CentralView ArticleGoogle 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Song X, Li P, Zhai J, Zhou M, Ma L, Liu B, Jeong DH, Nakano M, Cao S, Liu C, Chu C, Wang XJ, Green PJ, Meyers BC, Cao X: Roles of DCL4 and DCL3b in rice phased small RNA biogenesis. Plant J. 2012, 69 (3): 462-474.PubMedView ArticleGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36 (Database issue): D154-158.PubMedPubMed CentralGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31 (13): 3406-3415.PubMedPubMed CentralView ArticleGoogle Scholar
- Shi R, Chiang VL: Facile means for quantifying microRNA expression by real-time PCR. Biotechniques. 2005, 39 (4): 519-525.PubMedView ArticleGoogle Scholar
- Wan LC, Zhang H, Lu S, Zhang L, Qiu Z, Zhao Y, Zeng QY, Lin J: Transcriptome-wide identification and characterization of miRNAs from Pinus densata. BMC Genomics. 2012, 13: 132-PubMedPubMed CentralView ArticleGoogle Scholar