Skip to main content

Main menu

  • Authors
  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
    • Focus Collections
    • Classics Collection
    • Upcoming Focus Issues
  • Submit a Manuscript
  • Advertisers
  • About
    • About the Journal
    • Editorial Board and Staff
  • Subscribers
  • Librarians
  • More
    • Alerts
    • Contact Us
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Plant Cell Teaching Tools
    • ASPB
    • Plantae

User menu

  • My alerts
  • Log in

Search

  • Advanced search
Plant Physiology
  • Other Publications
    • Plant Physiology
    • The Plant Cell
    • Plant Direct
    • The Arabidopsis Book
    • Plant Cell Teaching Tools
    • ASPB
    • Plantae
  • My alerts
  • Log in
Plant Physiology

Advanced Search

  • Authors
  • Home
  • Content
    • Current Issue
    • Archive
    • Preview Papers
    • Focus Collections
    • Classics Collection
    • Upcoming Focus Issues
  • Submit a Manuscript
  • Advertisers
  • About
    • About the Journal
    • Editorial Board and Staff
  • Subscribers
  • Librarians
  • More
    • Alerts
    • Contact Us
  • Follow plantphysiol on Twitter
  • Visit plantphysiol on Facebook
  • Visit Plantae
Research ArticleGENETICS, GENOMICS, AND MOLECULAR EVOLUTION
Open Access

Small RNA Profiling in Two Brassica napus Cultivars Identifies MicroRNAs with Oil Production- and Development-Correlated Expression and New Small RNA Classes

Ying-Tao Zhao, Meng Wang, San-Xiong Fu, Wei-Cai Yang, Cun-Kou Qi, Xiu-Jie Wang
Ying-Tao Zhao
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Meng Wang
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
San-Xiong Fu
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Wei-Cai Yang
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cun-Kou Qi
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Xiu-Jie Wang
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site

Published February 2012. DOI: https://doi.org/10.1104/pp.111.187666

  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading
  • © 2012 American Society of Plant Biologists. All rights reserved.

Abstract

MicroRNAs (miRNAs) and small interfering RNAs are important regulators of plant development and seed formation, yet their population and abundance in the oil crop Brassica napus are still not well understood, especially at different developmental stages and among cultivars with varied seed oil contents. Here, we systematically analyzed the small RNA expression profiles of Brassica napus seeds at early embryonic developmental stages in high-oil-content and low-oil-content B. napus cultivars, both cultured in two environments. A total of 50 conserved miRNAs and 9 new miRNAs were identified, together with some new miRNA targets. Expression analysis revealed some miRNAs with varied expression levels in different seed oil content cultivars or at different embryonic developmental stages. A large number of 23-nucleotide small RNAs with specific nucleotide composition preferences were also identified, which may present new classes of functional small RNAs.

Increasing numbers of small RNAs have been identified as regulatory RNAs that modulate gene expression at both the transcriptional and posttranscriptional levels. Currently, two major classes of endogenous small RNAs have been identified in plants, namely microRNAs (miRNAs) and small interfering RNAs (siRNAs; Brodersen and Voinnet, 2006; Jones-Rhoades et al., 2006; Mallory and Vaucheret, 2006; Chen, 2009; Voinnet, 2009). They are approximately 20 to 24 nucleotides long and have crucial regulatory functions in diverse biological processes. Plant miRNAs are produced from hairpin-shaped precursors and mostly 21 nucleotides long, whereas the majority of plant siRNAs are 24 nucleotides in length and generated from long double-stranded RNA duplexes or transcripts derived from inverted repeat regions (Chen, 2009; Voinnet, 2009). Endogenous plant siRNAs can be divided into several classes, including miRNA-induced transacting siRNAs (ta-siRNAs), heterochromatic siRNAs, natural antisense siRNAs, and millions of unclassified small RNA sequences (Choudhuri, 2009). Ta-siRNAs are in-phase-generated small RNAs whose production depends on the cleavage of non-protein-coding transcripts of ta-siRNA genes by miRNAs (Vazquez et al., 2004; Allen et al., 2005; Yoshikawa et al., 2005; Montgomery et al., 2008a, 2008b). Both miRNAs and ta-siRNAs can induce the cleavage of mRNAs with partially or fully complementary sequences to them. Heterochromatic siRNAs are usually produced from repeats and transposable elements (Xie et al., 2004; Lu et al., 2005; Kasschau et al., 2007; Zhang et al., 2007; Mosher et al., 2008) and inhibit gene expression at the transcriptional level by modulating DNA or histone methylation of homologous regions. Natural antisense siRNAs are generated from natural sense-antisense transcript pairs and modulate the expression of one gene in the pair (Borsani et al., 2005; Katiyar-Agarwal et al., 2006; Held et al., 2008; Zhou et al., 2009; Ron et al., 2010).

MiRNAs and siRNAs play key roles in the regulation of diverse processes related to plant development, such as seed germination (Liu et al., 2007; Reyes and Chua, 2007), organ formation (Mallory et al., 2004; Baker et al., 2005; Wang et al., 2005), developmental timing and patterning (Chuck et al., 2009; Poethig, 2009; Rubio-Somoza et al., 2009; Wollmann and Weigel, 2010), hormonal response (Achard et al., 2004; Mallory et al., 2005; Wang et al., 2005; Reyes and Chua, 2007), and stress resistance (Sunkar et al., 2007; Ruiz-Ferrer and Voinnet, 2009). Tightly regulated space- and time-specific expression patterns have been identified for some miRNAs and siRNAs. The recent application of massively parallel sequencing technology has significantly accelerated the discovery and functional study of small RNAs by producing millions of short reads in a single sequencing run (Moxon et al., 2008; Sunkar et al., 2008; Zhu et al., 2008; Hsieh et al., 2009; Ledandais-Brière et al., 2009; Li et al., 2010; Pantaleo et al., 2010). Nevertheless, except for several model species, the expression patterns and functions of most small RNAs are poorly understood in other plants.

Brassica napus is a major edible oil-producing crop. The oil accumulation of B. napus seeds usually starts around 3 weeks after flowering and peaks after another 3 weeks (Norton and Harris, 1975; Cummins et al., 1993; Eastmond and Rawsthorne, 2000). The seed oil content varies significantly among B. napus cultivars or under different growth environments. So far, there have been only a few reports on small RNA population and function analysis in B. napus (Wang et al., 2007; Xie et al., 2007; Buhtz et al., 2008; He et al., 2008; Pant et al., 2008, 2009; Huang et al., 2010; Wei et al., 2010). In addition, most reports on the oil production mechanisms in B. napus focused on the rapid oil accumulation stage; few studied the gene expression profiles of early embryos, which may have already determined the oil production abilities of B. napus seeds (Li et al., 2006; Sharma et al., 2008; Fu et al., 2009). To systematically identify miRNAs and other small RNAs that may be involved in the regulation of B. napus early embryonic development as well as its seed oil accumulation, we constructed small RNA libraries from the early siliques of high-seed-oil-content and low-seed-oil-content B. napus cultivars and profiled their small RNA expression using the Solexa sequencing technology (GenomeAnalyzer; Illumina). This work identified 50 conserved and 9 new miRNAs in B. napus, among which some miRNAs with developmentally regulated or cultivar-diverged expression patterns. We also validated several conserved targets and identified new targets for some miRNAs with expression differences in the high-seed-oil-content and low-seed-oil-content cultivars. Ta-siRNAs and a group of 23-nucleotide-long small RNAs with characteristic nucleotide composition were also identified and analyzed.

RESULTS

Deep Sequencing of B. napus Small RNAs

To investigate the small RNA expression profiles in B. napus cultivars with different oil contents as well as the impact of environmental factors on small RNA expression, a low-seed-oil-content cultivar (L cultivar; XZ088) and a high-seed-oil-content cultivar (H cultivar; XZ366) were planted in Nanjing, a major B. napus-growing riverside plain land in China, and Lhasa, an altiplanic area where the oil content of B. napus seeds usually increased about 4% as compared with those harvested in Nanjing, respectively (Supplemental Table S1). For each sample, small RNAs from developing siliques at 3, 7, 14, and 21 d after flowering (DAF) were collected and pooled together for sequencing using the Solexa sequencing technology. A total of 18,881,853 reads with lengths of 15 to 30 nucleotides were obtained from the four libraries after removing adaptors and low-quality reads, representing 1,994,899, 2,149,419, 1,232,566, and 1,811,074 nonredundant sequences from the four libraries, respectively. Since the B. napus genome sequence has not been released yet, we mapped these small RNAs to the available bacterial artificial chromosome (BAC) clone sequences of all Brassica family species. As a result, perfect matches were obtained for about 20% of the total small RNAs, representing 15% of nonredundant sequences (Supplemental Table S2).

Similar to the length distribution of small RNAs observed in other plants, the most abundant small RNAs with or without perfect genomic matches were 24 nucleotides long in all four libraries, and the second largest population was the 21-nucleotide RNAs (Supplemental Fig. S1). Unexpectedly, a relatively large number of 22- and 23-nucleotide small RNAs were obtained in our data sets, as compared with previous reports in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa; Sunkar et al., 2008; Zhu et al., 2008; Hsieh et al., 2009).

Genomic location analysis showed that 13.2% of the total small RNAs with perfect genomic matches were generated from regions encoding reported B. napus miRNAs or sequences homologous to miRNA genes in other species, and 93.1% of these sequences were 21 nucleotides long (Supplemental Table S3). Small RNAs from unannotated genomic regions counted for 58.4% of the total mapped sequences and were dominated by 24-nucleotide small RNAs, followed by 22- and 23-nucleotide RNAs (Supplemental Table S3). rRNA-derived small RNAs also counted for 19.6% of the total population and were enriched in 21-nucleotide-long sequences (Supplemental Table S3).

Last Nucleotide Preference of 23-Nucleotide Small RNAs

Previous studies have shown that in plants, most miRNA sequences start with uridine (U), whereas the majority of 24-nucleotide siRNAs have adenosine (A) as their 5′ first nucleotide (Mi et al., 2008; Rubio-Somoza et al., 2009; Voinnet, 2009; Czech and Hannon, 2011). The same trends were also observed among our cloned B. napus small RNAs (Supplemental Fig. S2). Unexpectedly, a clear preference for adenosine at the last nucleotide was observed among the 23-nucleotide small RNAs, which counted for about 45.5% of total mapped 23-nucleotide small RNAs. Detailed analysis revealed that such last-nucleotide preference did not present in the sequences of other lengths, and the difference was statistically significant (Fig. 1A).

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Nucleotide preferences of small RNAs. A, Percentage of adenosine at the last position of 21- to 24-nucleotide (nt) small RNAs. The asterisk indicates a difference of P < 0.01 (Student’s single-sample two-sided t test). B, Nucleotide preferences of the 24-nucleotide-independent 23-nucleotide small RNAs at each position. C, Percentage of adenosine at each position of the 23-nucleotide small RNAs of 162 Arabidopsis deep sequencing data sets deposited in the NCBI GEO database.

To investigate the cause of such nucleotide composition bias, we compared the genomic locations of the 23-nucleotide small RNAs with the 24-nucleotide ones in detail. Among the 519,119 23-nucleotide small RNAs with perfect genomic matches, 365,343 were substrings of 24-nucleotide small RNAs with identical 5′ or 3′ ends as the 24-nucleotide ones. Therefore, we termed this group of sequences the 24-nucleotide-dependent 23-nucleotide small RNAs. The rest of the 153,776 23-nucleotide small RNAs were 24-nucleotide independent, as they were not degradation products of the 24-nucleotide ones. Nucleotide distribution analysis revealed that the preference for adenosine at the last position was most significant among the 24-nucleotide-independent 23-nucleotide small RNAs, presented in about 51.6% of sequences (Fig. 1B), whereas the occurrence rate of adenosine at the last position of the 24-nucleotide-dependent 23-nucleotide small RNAs was about 40%, only slightly higher than that of sequences with other lengths (Supplemental Fig. S3).

To examine whether the last nucleotide preference of the 23-nucleotide small RNAs also presented among small RNAs from other species, we downloaded 162 Arabidopsis small RNA deep sequencing datasets from the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) database (Barrett et al., 2011) and analyzed their nucleotide compositions. Strong preference for adenosine was also observed among the Arabidopsis 23-nucleotide small RNAs, but not among sequences of other lengths (Fig. 1C; Supplemental Fig. S4).

Identification of Conserved and New miRNAs

We first analyzed the presence of miRNAs in the sequencing data, as miRNAs have been shown to play critical roles in many aspects of plant developmental regulation. Although hundreds of plant miRNAs have been identified, only 46 conserved miRNAs were reported in B. napus by previously computational or experimental studies, according to the miRBase release 16 record (Griffiths-Jones et al., 2008). Therefore, it would be essential to systematically identify both plant-conserved and species-specific miRNAs in B. napus. A total of 36 nonredundant reads with identical sequences to miRBase-recorded plant miRNAs and 14 additional sequences that were highly homologous to miRBase-recorded plant miRNAs had hairpin-shaped precursors and were selected as miRNAs (Supplemental Table S4). Among them, 37 had their complementary short RNAs (miRNA*s) sequenced in our datasets, and only 13 were recorded as B. napus miRNAs in the miRBase. In addition, we also predicted 9 new miRNA genes with well-formed hairpin-shaped secondary structures and cloned miRNA* sequences or confirmed targets (Supplemental Fig. S5). The distribution of small RNAs on these miRNA precursors was centered on the miRNA or miRNA* sequences, with the defined miRNAs or miRNA*s having the highest clone numbers (Supplemental Fig. S5).

Some miRNA*s Accumulated Higher Than Corresponding miRNAs

It has been widely accepted that after processing from precursors, the functional strand of the short double-stranded RNA duplex will be retained as miRNAs, whereas their complementary miRNA* sequences will be quickly degraded (Jones-Rhoades et al., 2006; Chen, 2009; Voinnet, 2009). Therefore, in most cases, the abundance of miRNAs is much higher than that of their corresponding miRNA*s. Unexpectedly, we detected that the abundance of the star sequences of two conserved miRNAs, miR160a* and miR171b*, was about seven times more than that of miRNA160a and miRNA171b (Fig. 2, A and B).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Expression of several miRNAs and their corresponding miRNA*s. A to C, Normalized sequence reads of miR160a (A), miR171b (B), and miR408 (C) and their miRNA*s in the four libraries. D, Northern blot hybridization of miR408 and miR408*.

The sequencing results showed that the expression abundance of miR408* was extremely higher than that of miR408 in all four samples (Fig. 2C). Moreover, the abundance of miR408* was higher in the H cultivar than in the L cultivar, as shown by both the sequencing results and northern-blot hybridization (Fig. 2, C and D).

Differentially Expressed miRNAs between the High-Oil-Content and Low-Oil-Content Cultivars

Among the identified miRNAs, 22 miRNA families had more than five reads and therefore were compared for their expression abundance in the two cultivars (Supplemental Tables S4 and S5; Supplemental Fig. S5). Four miRNA families, namely miR169, miR390, miR394, and miR6028, had consistently higher expression in the L cultivar when planted both in Nanjing and Lhasa (Fig. 3A; 1.5-fold change cutoff). To the contrary, miR408* and miR2111 were more abundant (1.5-fold change cutoff) in the H cultivar in both growth environments. Such expression differences were also confirmed by quantitative reverse transcription (RT)-PCR and northern-blot hybridization (Fig. 3, B and C).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Expression profile of miRNAs between the H and L cultivars. A, Comparison of miRNA abundance between the two cultivars. B and C, Quantitative RT-PCR (B) and northern-blot hybridization (C) analyses of selected differentially expressed miRNAs between the cultivars. Mixed samples containing all examined time points were used. Error bars indicate the sd of three replicates.

Most other miRNAs only had greater than 1.5-fold expression changes between the two cultivars in either Nanjing or Lhasa; thus, it is difficult to conclude whether the expression changes were due to genetic or environmental differences. For example, the sequencing results showed that the expression of miR156 was higher in the H cultivar than in the L cultivar in samples grown in Nanjing but had no difference between the two cultivars in samples collected in Lhasa. Such an expression difference was also confirmed by northern-blot hybridization (Fig. 3C). We also used quantitative RT-PCR to examine the expression of a new miRNA, miR6029. Consistent with the sequencing results, miR6029 had higher expression in the L cultivar than in the H cultivar in both growth areas, but only the expression change in the Nanjing sample was greater than 1.5-fold (Fig. 3B). These experiments confirmed the fidelity of the sequencing results and also indicated that miR156 and miR6029 may relate to B. napus oil production.

Expression Changes of miRNAs during Early Silique Development

As the small RNA sequencing data were obtained by pooling samples from four silique developmental stages together, we further examined the expression dynamics of several miRNAs at different time points. To exclude environmental effects, only samples collected from Nanjing were used in this analysis. Diversified expression patterns were observed for some miRNAs. For example, the expression of miR156 remained consistent in the L cultivar but increased significantly at 21 DAF in the H cultivar (Fig. 4A). Similarly, but reversed, the amount of miR167 increased gradually along with the examined time points and shot up significantly at 21 DAF in the L cultivar (Fig. 4B). The expression pattern of miR6029 in the L cultivar was very similar to that of miR167, but it remained consistently low in the H cultivar (Fig. 4B). On the other hand, the expression of miR390 was high at 3 DAF in both cultivars, decreased quickly at 7 DAF, and remained low until 21 DAF (Fig. 4A). The amount of both miR2111 and miR6028 increased consistently during development in the L cultivar, but miR6028 decreased gradually with time in the H cultivar, whereas miR2111 also increased in the H cultivar (Fig. 4).

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Expression changes of miRNAs at different silique developmental stages. A, Quantitative RT-PCR analysis for the relative expression of miR156, miR390, miR2111, and miR6028 during early silique development. The expression level of each high-oil-content cultivar miRNA at 3 DAF is set to 1. B, Northern-blot hybridization of miR156, miR390, miR167, and miR6029 during early silique development. Error bars indicate the sd of three replicates.

Conserved and New Targets of B. napus miRNAs

We used previously reported plant miRNA target prediction criteria (Allen et al., 2005) to search for putative targets of the conserved and new B. napus miRNAs. A total of 346 putative targets were identified for all conserved miRNAs as well as all but two new miRNAs, namely miR6031 and miR6036 (Supplemental Table S6). The targets of most conserved miRNAs were homologous to those reported miRNA targets in Arabidopsis or rice. The newly identified targets with transcription factor functions are listed in Table I.

View this table:
  • View inline
  • View popup
Table I. The newly identified miRNA targets with transcription factor functions

We next performed a 5′-RNA ligase-mediated (RLM)-RACE experiment to examine whether some miRNAs could indeed induce the cleavage of their predicted targets. The putative targets of miR156d were squamosa-promoter-binding protein-like (SPL) family genes. We identified several SPL genes containing a putative miR156d-binding site. Cleavage sites mediated by miR156d were detected in the transcripts of one SPL2 gene (TC82507) and two SPL10 genes (TC98126 and TC94184; Fig. 5A). Similarly, the miRNA::target relationships were also confirmed for miR167c and an ARF8-like gene (TC84482). We identified that transcripts of a basic-Leu zipper transcription factor family gene, a gene encoding VirE2 interacting protein1 (VIP1)-like protein, as the target of a new miRNA, miR6029 (Fig. 5A). It is worth noting that the cleavage of VIP1 transcript occurred at the ninth position of miR6029 instead of the commonly observed 10th position in all 21 RACE experiments. Similarly, higher cleavage frequency was also observed between the pairs of miR156d::TC98126 and miR156d::TC94184 at the ninth position.

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Validation and expression of selected miRNA target genes. A, 5′-RLM-RACE analysis of the cleavage on target mRNAs by corresponding miRNAs. B, Quantitative RT-PCR analysis of miR156 target gene TC98126 at different developmental stages and miR6028 target gene CD814980 in different cultivars. Error bars indicate the sd of three replicates.

We also used quantitative RT-PCR to examine the abundance of several miRNAs and their targets. Contradictory to the expression profile of miR156 (Fig. 4), the abundance of its target, TC98126, decreased significantly at 21 DAF in the H cultivar but remained almost unchanged in the L cultivar (Fig. 5B). In another example, the abundance of CD814980, a putative target of a new miRNA, miR6028, was comparable in the H and L cultivars at 3 DAF (Fig. 5B) but decreased significantly in the L cultivar at 21 DAF, which was in accordance with the increased expression of miR6028 in the L cultivar at the same time (Fig. 4).

Identification of miRNA-Derived ta-siRNA Genes

MiRNA-derived ta-siRNA is another important class of regulatory small RNAs in plants, but they have not been identified in B. napus yet. We carried out a genome-wide systemic screen for ta-siRNA genes in our data sets using methods modified from previous reports (Vazquez et al., 2004; Allen et al., 2005; Yoshikawa et al., 2005; Montgomery et al., 2008a, 2008b). Phased small RNA clusters were identified from three genes with putative miRNA targeting sites, namely TC65909, S39172450, and a gene located in BAC (region 44,041–44,663, antisense strand; Supplemental Table S7). All three sequences had high homology to the Arabidopsis ta-siRNA gene 3 (TAS3; Supplemental Fig. S6) and therefore were considered as the B. napus TAS3 genes.

Similar to the observations in Arabidopsis, the B. napus TAS3 gene TC65909 also had two miR390-binding sites located at both the upstream and downstream boundaries of the ta-siRNA production region (Fig. 6A). Consecutively produced phased small RNAs covering both the sense and antisense strands of TC65909 were identified from the majority of putative ta-siRNA production loci, with TAS3 5′ D7(+) as the most abundant sense strand-produced ta-siRNA and TAS3 5′ D6(−) as the most abundant antisense strand-produced ta-siRNA. These sense and antisense strand-generated ta-siRNAs were mainly 21 nucleotides long and could form double-stranded duplexes with two-nucleotide overhangs at their 3′ end, consistent with the previously reported features of ta-siRNAs (Supplemental Fig. S6).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

MiRNA-derived ta-siRNAs in B. napus. A, Expression abundance of ta-siRNAs generated from TAS3 gene TC65909. The locations of two miR390 target sites and TAS3 5′ D7(+) are labeled by short bars. B and C, Quantitative RT-PCR (B) and northern-blot hybridization (C) analyses of TAS3 5′ D7(+) during early silique development. D, Validation of cleavage on B. napus TC110025 (ARF3) mRNAs by TAS3 5′ D7(+) using 5′-RLM-RACE. Error bars indicate the sd of three replicates.

We next investigated the biogenesis and functions of these B. napus ta-siRNAs using TAS3 5′ D7(+), whose homolog is also the functional ta-siRNA in Arabidopsis. The abundance of TAS3 5′ D7(+) was highest at 3 DAF, then decreased significantly at 7 DAF, and remained low at the following time points (Fig. 6, B and C). Such an expression pattern was identical to that of miR390 (Fig. 4), suggesting that the expression of TAS3 5′ D7(+) was induced by miR390. Target prediction revealed that several genes of the Auxin Response Factor (ARF) family, namely ARF2 (TC81248), ARF3 (TC110025), and ARF4 (TC100380), were putative targets of TAS3 5′ D7(+). We were able to confirm the cleavage of ARF3 by TAS3 5′ D7(+) using 5′-RLM-RACE, with most identified cleavage occurring at the 10th position of the TAS3 5′ D7(+)-binding site (Fig. 6D).

DISCUSSION

MiRNAs have emerged as a new class of regulatory factors and attracted much attention during the past decade. As B. napus is an important crop for edible oil production, understanding the functions of miRNAs in regulating oil accumulation in B. napus seeds will be of great value for the cultivation of B. napus cultivars with increased oil content. In this work, we examined the expression profiles of miRNAs in high-oil-content and low-oil-content B. napus cultivars. Prior to this work, only 46 miRNAs were reported in B. napus, according to the miRBase release 16 record. Here, we enlarged the known population of B. napus miRNAs by identifying 37 unreported conserved miRNAs and 9 new miRNAs, as well as the star sequences for some miRNAs.

We identified putative target genes for 57 miRNAs using computational methods. The lack of predicted targets for other miRNAs is probably due to the unavailability of the B. napus genome sequence. The target genes of most conserved miRNAs have also been identified as miRNA targets in Arabidopsis, suggesting the functional conservation of these miRNAs. We confirmed that miRNA156 regulated SPL family genes and miR167 targeted transcripts of the ARF family in B. napus. We also identified and confirmed transcripts encoding VIP1 as the target of a new miRNA, miR6029. VIP1 is a member of basic-Leu zipper transcription factors and can regulate many genes by binding to their promoter regions (Pitzschke et al., 2009). It was not reported to be regulated by miRNA previously. Our 5′-RLM-RACE experiment obtained consistent miR6029 cleavage sites on VIP1 mRNAs in 21 individual sequencing reactions, demonstrated that VIP1 is an authentic target of miR6029. The pairing between miR6029 and the VIP1 transcript contained two mismatches and two insertions, further proving that gaps are also tolerable between plant miRNA and target pairs. In addition, all cleavage of VIP1 occurred at the ninth position of the miRNA sequence, instead of the commonly observed 10th position, probably due to the presence of an insertion at the 11th position of the pair.

It is well known that many miRNAs have time- and space-specific expression patterns (Jones-Rhoades et al., 2006; Mallory and Vaucheret, 2006; Chen, 2009; Chuck et al., 2009; Poethig, 2009; Rubio-Somoza et al., 2009; Wollmann and Weigel, 2010), yet the expression profiles of miRNAs in the early developmental stage of plant embryos are still poorly understood, especially in the oil production crop B. napus. By comparing the small RNA expression profiles of 3- to 21-DAF embryos using the massively parallel sequencing technology, we have identified some differentially expressed miRNAs between the high-oil-content and low-oil-content B. napus cultivars as well as some miRNAs with expression changes at different developmental stages within the same cultivar, and we have further confirmed their expression differences using northern-blot hybridization and quantitative RT-PCR experiments. The similar or reversed expression patterns of the selected miRNAs examined at individual time points suggested that they may have synergistic or antagonistic functions.

By comparing the high-oil-content and low-oil-content samples, we have identified and confirmed that miR156, miR167, miR390, miR2111, and two new miRNAs, miR6028 and miR6029, as well as the TAS3-produced ta-siRNAs were differentially expressed between the two cultivars in one or both growth environments. According to the varied expression levels of these small RNAs at different developmental stages and the functions of their targets, we propose that miR390 is involved in early embryo development regulation via the ta-siRNA-mediated auxin signaling pathways, whereas miR156, miR167, and miR6029 function around 21 DAF by regulating SPL, ARF, and VIP1 transcripts, which may affect the oil content of B. napus seeds in later stages through multiple pathways (Supplemental Fig. S7).

The plant genomes usually express large numbers of small RNAs, the majority of which were not classified yet. Unlike the small RNA profiles of most reported plant samples, a large number of 23-nucleotide small RNAs were identified in our data. A total of 29.6% of these 23-nucleotide small RNAs were not potential degradation products of 24-nucleotide small RNAs and exhibited a strong preference for adenosine at the 23-nucleotide position, indicating that they may be an uncharacterized class of functional small RNAs.

CONCLUSION

To the best of our knowledge, this work provides the first small RNA expression analysis in B. napus early embryos. It enlarged the population of known B. napus miRNAs and proved the presence of ta-siRNAs in B. napus. The comparison between high-oil-content and low-oil-content cultivars revealed some miRNAs that may be involved in the regulation of B. napus seed oil production. The detailed expression analysis of some miRNAs and their targets at different developmental time points provided clues for their varied roles in regulating the early development of B. napus seed, which may be conserved in other dicot plants as well. The 23-nucleotide small RNAs found in this work suggested that there may be other types of uncharacterized small RNAs remaining to be discovered.

MATERIALS AND METHODS

Silique Collection, RNA Isolation, and Small RNA Sequencing

Two Brassica napus cultivars (XZ088 and XZ366) were grown in the experimental fields in the plain land Nanjing and the altiplanic area Lhasa, respectively. Siliques of both cultivars were collected at 3, 7, 14, and 21 DAF separately. For each time point, the siliques of more than 20 plants were collected and pooled together to avoid the individual genotype differences. Total RNAs of each sample were isolated using Trizol reagent (Invitrogen). Equal amounts of total RNAs of the 3-, 7-, 14-, and 21-DAF samples of each cultivar were mixed before sending for sequencing. Small RNA sequences of XZ088 and XZ366 siliques collected in both growth environments were generated by Illumina’s Solexa sequencing technology.

Primary Analysis of the Deep Sequencing Data Sets

After trimming adapter sequences and removing sequences with low quality or length less than 15 nucleotides, the resulting small RNA sequences were compared with the Rfam database (Griffiths-Jones et al., 2005) and the NCBI nucleotide database to classify tRNA, rRNA, small nucleolar RNA, and protein-coding gene-derived small RNAs. The BAC clone sequences of B. napus, Brassica rapa, and Brassica oleracea collected in the NCBI nucleotide database by the BLASTN program (Altschul et al., 1990) were used to analyze the genomic origins of the remaining small RNAs. The BAC clone sequences of B. rapa and B. oleracea were included because the B. napus genome (AC genome) is a polyploid of the B. rapa genome (A genome) and the B. oleracea genome (C genome). All BAC clone sequences were first analyzed using the RepeatMasker software to identify transposons and repeat regions. Small RNA sequences were aligned to the BAC clones by the BLASTN program (Altschul et al., 1990), and only sequences with perfect genomic matches were retained for further analysis.

Identification of miRNAs and ta-siRNAs

Small RNAs of lengths between 20 and 24 nucleotides, with perfect matches to the nonrepeat-associated noncoding regions of BAC clones, as well as with no more than seven mapping loci on nonredundant BAC sequences, were used in miRNA prediction. According to the community-established criteria (Meyers et al., 2008), small RNA sequences satisfying all the following rules were selected as miRNAs: (1) intergenic location of miRNA candidates; (2) hairpin-shaped secondary structures of precursors; (3) stem region location of mature miRNAs on precursors; (4) fewer than seven gaps and mismatches between mature miRNAs and their pairing sequences; (5) the presence of miRNA* sequences or the confirmation of target cleavage by a small RNA. Putative precursor sequences of each small RNA were extracted by assuming that the small RNAs can be derived from either the 5′ or 3′ arm of the hairpin-shaped precursors, with a 20-nucleotide extension at either the 5′ or 3′ end of the small RNA loci and a 60- to 200-nucleotide extension (20-nucleotide increment each time) at the other end. The secondary structure of each putative small RNA precursor was evaluated by the mfold software (Markham and Zuker, 2008). Candidate miRNA precursor sequences were compared with each other using the BLASTN program to remove redundant sequences from different BAC clones. The conservation status of the identified miRNAs was analyzed by comparing their sequences with all plant miRNAs collected in miRBase (release 16).

In the process of ta-siRNA identification, all Brassica BACs, B. napus tentative consensus sequences, and UniGene sequences were screened by a 200-nucleotide width window, with a 10-nucleotide increment at each step. The number of phased small RNAs with lengths of 21 nucleotides within each window was calculated, and in-phase-produced small RNAs starting from putative miRNA target sites were identified as ta-siRNAs according to previously described criteria (Vazquez et al., 2004; Allen et al., 2005; Yoshikawa et al., 2005; Montgomery et al., 2008a, 2008b). Redundant genes were merged by the BLASTN program.

Target Gene Prediction for miRNAs and ta-siRNAs

The tentative consensus sequences from the Dana Farber Cancer Institute (DFCI) Gene Index Project of B. napus were downloaded and used for miRNA and ta-siRNA target prediction. MiRNAs and ta-siRNAs were aligned to these sequences using Smith-Waterman algorithm-based C programs (Smith et al., 1981). Putative miRNA- and siRNA-binding sites were identified using a scoring system modified from a previous report (Allen et al., 2005) as follows: (1) a penalty score of 1 for each mismatch or gap between the pairing of miRNAs and transcripts; (2) a penalty score of 0.5 for each G::U pair; (3) doubled penalty scores for mismatches, gaps, or G::U pairs located within the two to 13 positions of miRNAs. A total penalty score of no more than 5 was used as the cutoff to select miRNA and ta-siRNA targets.

Small RNA Northern-Blot Hybridization

Small RNA northern-blot hybridization was carried out following previously reported procedures (Liu et al., 2005) using 20 μg of total RNA in each experiment. The complementary sequences of small RNAs were labeled by [γ-32P]ATP and used as probes. The sequences of all probes used are listed in Supplemental Table S8.

Quantitative RT-PCR

The stem-loop RT-PCR method was used in the quantitative RT-PCR experiments of small RNAs following previously reported procedures (Chen et al., 2005; Varkonyi-Gasic et al., 2007). In the quantitative RT-PCR experiments of small RNAs, 3 μg of total RNAs treated with DNase I (New England Biolabs) was used in each reaction with U6 snRNA as the internal control. These reactions were performed using previously reported procedures (Chen et al., 2005; Varkonyi-Gasic et al., 2007). In the quantitative RT-PCR experiments of mRNAs, 1 μg of total RNAs treated with DNase I (New England Biolabs) was synthesized to cDNA by Moloney murine leukemia virus (New England Biolabs) using poly(dT) oligonucleotides with the mRNA of Bn-EF1 (for elongation factor 1-α; TC111637) as the reference gene. SYBR Green PCR Master Mix (Applied Biosystems) was used in all quantitative RT-PCR experiments. The relative fold expression changes of miRNAs and genes were both calculated using the 2 δ-δ Ct method (Livak and Schmittgen, 2001). Primers used in all quantitative RT-PCR experiments are listed in Supplemental Table S9.

5′-RLM-RACE

The RNA 5′-RLM-RACE experiments were carried out using the FirstChoice RLM-RACE Kit (Ambion) following previously reported procedures (Lu et al., 2008). All gene-specific primers used in the experiments are listed in Supplemental Table S10.

Sequence data from this article can be found in the GenBank/EMBL data libraries under accession number GSE34727.

Supplemental Data

The following materials are available in the online version of this article.

  • Supplemental Figure S1. Size distribution of small RNAs.

  • Supplemental Figure S2. Nucleotide preference at each position of small RNAs.

  • Supplemental Figure S3. Nucleotide preference at each position of the 24-nucleotide-dependent 23-nucleotide small RNAs.

  • Supplemental Figure S4. Percentage of adenosine at each position of 21-nucleotide, 22-nucleotide, and 24-nucleotide small RNAs of the 162 Arabidopsis small RNA deep sequencing data sets deposited in the NCBI GEO database.

  • Supplemental Figure S5. Secondary structures of the B. napus miRNA candidate precursors, and the expression and location of small RNAs mapped to these precursors.

  • Supplemental Figure S6. Size distribution of B. napus ta-siRNAs, and the sequence conservations of B. napus TAS3 genes.

  • Supplemental Figure S7. A proposed model of the regulation of B. napus early silique development via miRNAs, ta-siRNAs, and their target genes.

  • Supplemental Table S1. Seed oil contents of the two cultivars planted in Nanjing and Lhasa.

  • Supplemental Table S2. Summary of small RNA deep sequencing data.

  • Supplemental Table S3. The mapping loci distributions of small RNAs with different lengths.

  • Supplemental Table S4. Conserved B. napus miRNA candidates.

  • Supplemental Table S5. New B. napus miRNA candidates.

  • Supplemental Table S6. Predicted target genes of the B. napus miRNAs.

  • Supplemental Table S7. Ta-siRNA genes in B. napus.

  • Supplemental Table S8. Small RNA northern-blot hybridization probes used in this study.

  • Supplemental Table S9. Quantitative RT-PCR primers used in this study.

  • Supplemental Table S10. 5′-RLM-RACE primers used in this study.

Acknowledgments

We thank Ms. Nima Zhuoma (Tibet Academy of Agricultural and Animal Husbandry Sciences) for help with plant culturing and sample collection in Lhasa.

Footnotes

  • The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Xiu-Jie Wang (xjwang{at}genetics.ac.cn).

  • www.plantphysiol.org/cgi/doi/10.1104/pp.111.187666

  • ↵1 This work was supported by the National Natural Science Foundation of China (grant nos. 30725014 and 30921016 to X.-J.W. and grant no. 30800632 to M.W.) and the National High Technology Research and Development Program of China (grant no. 2009AA101105 to C.-K.Q.).

  • ↵2 These authors contributed equally to the article.

  • ↵3 These authors contributed equally to the article.

  • ↵[OA] Open Access articles can be viewed online without a subscription.

  • ↵[W] The online version of this article contains Web-only data.

  • Received September 22, 2011.
  • Accepted November 26, 2011.
  • Published December 2, 2011.

REFERENCES

  1. ↵
    1. Achard P,
    2. Herr A,
    3. Baulcombe DC,
    4. Harberd NP
    (2004) Modulation of floral development by a gibberellin-regulated microRNA. Development 131: 3357–3365
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Allen E,
    2. Xie Z,
    3. Gustafson AM,
    4. Carrington JC
    (2005) MicroRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121: 207–221
    OpenUrlCrossRefPubMed
  3. ↵
    1. Altschul SF,
    2. Gish W,
    3. Miller W,
    4. Myers EW,
    5. Lipman DJ
    (1990) Basic local alignment search tool. J Mol Biol 215: 403–410
    OpenUrlCrossRefPubMed
  4. ↵
    1. Baker CC,
    2. Sieber P,
    3. Wellmer F,
    4. Meyerowitz EM
    (2005) The early extra petals1 mutant uncovers a role for microRNA miR164c in regulating petal number in Arabidopsis. Curr Biol 15: 303–315
    OpenUrlCrossRefPubMed
  5. ↵
    1. Barrett T,
    2. Troup DB,
    3. Wilhite SE,
    4. Ledoux P,
    5. Evangelista C,
    6. Kim IF,
    7. Tomashevsky M,
    8. Marshall KA,
    9. Phillippy KH,
    10. Sherman PM,
    11. et al.
    (2011) NCBI GEO: archive for functional genomics data sets—10 years on. Nucleic Acids Res 39: D1005–D1010
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Borsani O,
    2. Zhu J,
    3. Verslues PE,
    4. Sunkar R,
    5. Zhu JK
    (2005) Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell 123: 1279–1291
    OpenUrlCrossRefPubMed
  7. ↵
    1. Brodersen P,
    2. Voinnet O
    (2006) The diversity of RNA silencing pathways in plants. Trends Genet 22: 268–280
    OpenUrlCrossRefPubMed
  8. ↵
    1. Buhtz A,
    2. Springer F,
    3. Chappell L,
    4. Baulcombe DC,
    5. Kehr J
    (2008) Identification and characterization of small RNAs from the phloem of Brassica napus. Plant J 53: 739–749
    OpenUrlCrossRefPubMed
  9. ↵
    1. Chen C,
    2. Ridzon DA,
    3. Broomer AJ,
    4. Zhou Z,
    5. Lee DH,
    6. Nguyen JT,
    7. Barbisin M,
    8. Xu NL,
    9. Mahuvakar VR,
    10. Andersen MR,
    11. et al.
    (2005) Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 33: e179
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Chen X
    (2009) Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol 25: 21–44
    OpenUrlCrossRefPubMed
  11. ↵
    1. Choudhuri S
    (2009) Lesser known relatives of miRNA. Biochem Biophys Res Commun 388: 177–180
    OpenUrlCrossRefPubMed
  12. ↵
    1. Chuck G,
    2. Candela H,
    3. Hake S
    (2009) Big impacts by small RNAs in plant development. Curr Opin Plant Biol 12: 81–86
    OpenUrlCrossRefPubMed
  13. ↵
    1. Cummins I,
    2. Hills MJ,
    3. Ross JH,
    4. Hobbs DH,
    5. Watson MD,
    6. Murphy DJ
    (1993) Differential, temporal and spatial expression of genes involved in storage oil and oleosin accumulation in developing rapeseed embryos: implications for the role of oleosins and the mechanisms of oil-body formation. Plant Mol Biol 23: 1015–1027
    OpenUrlCrossRefPubMed
  14. ↵
    1. Czech B,
    2. Hannon GJ
    (2011) Small RNA sorting: matchmaking for Argonautes. Nat Rev Genet 12: 19–31
    OpenUrlCrossRefPubMed
  15. ↵
    1. Eastmond PJ,
    2. Rawsthorne S
    (2000) Coordinate changes in carbon partitioning and plastidial metabolism during the development of oilseed rape embryos. Plant Physiol 122: 767–774
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Fu SX,
    2. Cheng H,
    3. Qi C
    (2009) Microarray analysis of gene expression in seeds of Brassica napus planted in Nanjing (altitude: 8.9 m), Xining (altitude: 2261.2 m) and Lhasa (altitude: 3658 m) with different oil content. Mol Biol Rep 36: 2375–2386
    OpenUrlCrossRefPubMed
  17. ↵
    1. Griffiths-Jones S,
    2. Moxon S,
    3. Marshall M,
    4. Khanna A,
    5. Eddy SR,
    6. Bateman A
    (2005) Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res 33: D121–D124
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Griffiths-Jones S,
    2. Saini HK,
    3. van Dongen S,
    4. Enright AJ
    (2008) miRBase: tools for microRNA genomics. Nucleic Acids Res 36: D154–D158
    OpenUrlAbstract/FREE Full Text
  19. ↵
    1. He XF,
    2. Fang YY,
    3. Feng L,
    4. Guo HS
    (2008) Characterization of conserved and novel microRNAs and their targets, including a TuMV-induced TIR-NBS-LRR class R gene-derived novel miRNA in Brassica. FEBS Lett 582: 2445–2452
    OpenUrlCrossRefPubMed
  20. ↵
    1. Held MA,
    2. Penning B,
    3. Brandt AS,
    4. Kessans SA,
    5. Yong W,
    6. Scofield SR,
    7. Carpita NC
    (2008) Small-interfering RNAs from natural antisense transcripts derived from a cellulose synthase gene modulate cell wall biosynthesis in barley. Proc Natl Acad Sci USA 105: 20534–20539
    OpenUrlAbstract/FREE Full Text
  21. ↵
    1. Hsieh LC,
    2. Lin SI,
    3. Shih AC,
    4. Chen JW,
    5. Lin WY,
    6. Tseng CY,
    7. Li WH,
    8. Chiou TJ
    (2009) Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol 151: 2120–2132
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Huang SQ,
    2. Xiang AL,
    3. Che LL,
    4. Chen S,
    5. Li H,
    6. Song JB,
    7. Yang ZM
    (2010) A set of miRNAs from Brassica napus in response to sulphate deficiency and cadmium stress. Plant Biotechnol J 8: 887–899
    OpenUrlCrossRefPubMed
  23. ↵
    1. Jones-Rhoades MW,
    2. Bartel DP,
    3. Bartel B
    (2006) MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol 57: 19–53
    OpenUrlCrossRefPubMed
  24. ↵
    1. Kasschau KD,
    2. Fahlgren N,
    3. Chapman EJ,
    4. Sullivan CM,
    5. Cumbie JS,
    6. Givan SA,
    7. Carrington JC
    (2007) Genome-wide profiling and analysis of Arabidopsis siRNAs. PLoS Biol 5: e57
    OpenUrlCrossRefPubMed
  25. ↵
    1. Katiyar-Agarwal S,
    2. Morgan R,
    3. Dahlbeck D,
    4. Borsani O,
    5. Villegas A Jr.,
    6. Zhu JK,
    7. Staskawicz BJ,
    8. Jin H
    (2006) A pathogen-inducible endogenous siRNA in plant immunity. Proc Natl Acad Sci USA 103: 18002–18007
    OpenUrlAbstract/FREE Full Text
  26. ↵
    1. Ledandais-Brière C,
    2. Naya L,
    3. Sallet E,
    4. Calenge F,
    5. Frugier F,
    6. Hartmann C,
    7. Gouzy J,
    8. Crespi M
    (2009) Genome-wide Medicago truncatula small RNA analysis revealed novel microRNAs and isoforms differentially regulated in roots and nodules. Plant Cell 21: 2780–2796
    OpenUrlAbstract/FREE Full Text
  27. ↵
    1. Li RJ,
    2. Wang HZ,
    3. Mao H,
    4. Lu YT,
    5. Hua W
    (2006) Identification of differentially expressed genes in seeds of two near-isogenic Brassica napus lines with different oil content. Planta 224: 952–962
    OpenUrlCrossRefPubMed
  28. ↵
    1. Li Y,
    2. Zhang Q,
    3. Zhang J,
    4. Wu L,
    5. Qi Y,
    6. Zhou JM
    (2010) Identification of microRNAs involved in pathogen-associated molecular pattern-triggered plant innate immunity. Plant Physiol 152: 2222–2231
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Liu B,
    2. Li P,
    3. Li X,
    4. Liu C,
    5. Cao S,
    6. Chu C,
    7. Cao X
    (2005) Loss of function of OsDCL1 affects microRNA accumulation and causes developmental defects in rice. Plant Physiol 139: 296–305
    OpenUrlAbstract/FREE Full Text
  30. ↵
    1. Liu PP,
    2. Montgomery TA,
    3. Fahlgren N,
    4. Kasschau KD,
    5. Nonogaki H,
    6. Carrington JC
    (2007) Repression of AUXIN RESPONSE FACTOR10 by microRNA160 is critical for seed germination and post-germination stages. Plant J 52: 133–146
    OpenUrlCrossRefPubMed
  31. ↵
    1. Livak KJ,
    2. Schmittgen TD
    (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods 25: 402–408
    OpenUrlCrossRefPubMed
  32. ↵
    1. Lu C,
    2. Jeong DH,
    3. Kulkarni K,
    4. Pillay M,
    5. Nobuta K,
    6. German R,
    7. Thatcher SR,
    8. Maher C,
    9. Zhang L,
    10. Ware D,
    11. et al.
    (2008) Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA 105: 4951–4956
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Lu C,
    2. Tej SS,
    3. Luo S,
    4. Haudenschild CD,
    5. Meyers BC,
    6. Green PJ
    (2005) Elucidation of the small RNA component of the transcriptome. Science 309: 1567–1569
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Mallory AC,
    2. Bartel DP,
    3. Bartel B
    (2005) MicroRNA-directed regulation of Arabidopsis AUXIN RESPONSE FACTOR17 is essential for proper development and modulates expression of early auxin response genes. Plant Cell 17: 1360–1375
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Mallory AC,
    2. Dugas DV,
    3. Bartel DP,
    4. Bartel B
    (2004) MicroRNA regulation of NAC-domain targets is required for proper formation and separation of adjacent embryonic, vegetative, and floral organs. Curr Biol 14: 1035–1046
    OpenUrlCrossRefPubMed
  36. ↵
    1. Mallory AC,
    2. Vaucheret H
    (2006) Functions of microRNAs and related small RNAs in plants. Nat Genet (Suppl) 38: S31–S36
    OpenUrlCrossRefPubMed
  37. ↵
    1. Markham NR,
    2. Zuker M
    (2008) UNAFold: software for nucleic acid folding and hybridization. Methods Mol Biol 453: 3–31
    OpenUrlCrossRefPubMed
  38. ↵
    1. Meyers BC,
    2. Axtell MJ,
    3. Bartel B,
    4. Bartel DP,
    5. Baulcombe D,
    6. Bowman JL,
    7. Cao X,
    8. Carrington JC,
    9. Chen X,
    10. Green PJ,
    11. et al.
    (2008) Criteria for annotation of plant microRNAs. Plant Cell 20: 3186–3190
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Mi S,
    2. Cai T,
    3. Hu Y,
    4. Chen Y,
    5. Hodges E,
    6. Ni F,
    7. Wu L,
    8. Li S,
    9. Zhou H,
    10. Long C,
    11. et al.
    (2008) Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5′ terminal nucleotide. Cell 133: 116–127
    OpenUrlCrossRefPubMed
  40. ↵
    1. Montgomery TA,
    2. Howell MD,
    3. Cuperus JT,
    4. Li D,
    5. Hansen JE,
    6. Alexander AL,
    7. Chapman EJ,
    8. Fahlgren N,
    9. Allen E,
    10. Carrington JC
    (2008a) Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation. Cell 133: 128–141
    OpenUrlCrossRefPubMed
  41. ↵
    1. Montgomery TA,
    2. Yoo SJ,
    3. Fahlgren N,
    4. Gilbert SD,
    5. Howell MD,
    6. Sullivan CM,
    7. Alexander A,
    8. Nguyen G,
    9. Allen E,
    10. Ahn JH,
    11. et al.
    (2008b) AGO1-miR173 complex initiates phased siRNA formation in plants. Proc Natl Acad Sci USA 105: 20055–20062
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Mosher RA,
    2. Schwach F,
    3. Studholme D,
    4. Baulcombe DC
    (2008) PolIVb influences RNA-directed DNA methylation independently of its role in siRNA biogenesis. Proc Natl Acad Sci USA 105: 3145–3150
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Moxon S,
    2. Jing R,
    3. Szittya G,
    4. Schwach F,
    5. Rusholme Pilcher RL,
    6. Moulton V,
    7. Dalmay T
    (2008) Deep sequencing of tomato short RNAs identifies microRNAs targeting genes involved in fruit ripening. Genome Res 18: 1602–1609
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Norton G,
    2. Harris JF
    (1975) Compositional changes in developing rape seed (Brassica napus L). Planta 123: 163–174
    OpenUrlCrossRef
  45. ↵
    1. Pant BD,
    2. Buhtz A,
    3. Kehr J,
    4. Scheible WR
    (2008) MicroRNA399 is a long-distance signal for the regulation of plant phosphate homeostasis. Plant J 53: 731–738
    OpenUrlCrossRefPubMed
  46. ↵
    1. Pant BD,
    2. Musialak-Lange M,
    3. Nuc P,
    4. May P,
    5. Buhtz A,
    6. Kehr J,
    7. Walther D,
    8. Scheible WR
    (2009) Identification of nutrient-responsive Arabidopsis and rapeseed microRNAs by comprehensive real-time polymerase chain reaction profiling and small RNA sequencing. Plant Physiol 150: 1541–1555
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Pantaleo V,
    2. Szittya G,
    3. Moxon S,
    4. Miozzi L,
    5. Moulton V,
    6. Dalmay T,
    7. Burgyan J
    (2010) Identification of grapevine microRNAs and their targets using high-throughput sequencing and degradome analysis. Plant J 62: 960–976
    OpenUrlPubMed
  48. ↵
    1. Pitzschke A,
    2. Djamei A,
    3. Teige M,
    4. Hirt H
    (2009) VIP1 response elements mediate mitogen-activated protein kinase 3-induced stress gene expression. Proc Natl Acad Sci USA 106: 18414–18419
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Poethig RS
    (2009) Small RNAs and developmental timing in plants. Curr Opin Genet Dev 19: 374–378
    OpenUrlCrossRefPubMed
  50. ↵
    1. Reyes JL,
    2. Chua NH
    (2007) ABA induction of miR159 controls transcript levels of two MYB factors during Arabidopsis seed germination. Plant J 49: 592–606
    OpenUrlCrossRefPubMed
  51. ↵
    1. Ron M,
    2. Alandete Saez M,
    3. Eshed Williams L,
    4. Fletcher JC,
    5. McCormick S
    (2010) Proper regulation of a sperm-specific cis-nat-siRNA is essential for double fertilization in Arabidopsis. Genes Dev 24: 1010–1021
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Rubio-Somoza I,
    2. Cuperus JT,
    3. Weigel D,
    4. Carrington JC
    (2009) Regulation and functional specialization of small RNA-target nodes during plant development. Curr Opin Plant Biol 12: 622–627
    OpenUrlCrossRefPubMed
  53. ↵
    1. Ruiz-Ferrer V,
    2. Voinnet O
    (2009) Roles of plant small RNAs in biotic stress responses. Annu Rev Plant Biol 60: 485–510
    OpenUrlCrossRefPubMed
  54. ↵
    1. Sharma N,
    2. Anderson M,
    3. Kumar A,
    4. Zhang Y,
    5. Giblin EM,
    6. Abrams SR,
    7. Zaharia LI,
    8. Taylor DC,
    9. Fobert PR
    (2008) Transgenic increases in seed oil content are associated with the differential expression of novel Brassica-specific transcripts. BMC Genomics 9: 619
    OpenUrlCrossRefPubMed
  55. ↵
    1. Smith TF,
    2. Waterman MS,
    3. Fitch WM
    (1981) Comparative biosequence metrics. J Mol Evol 18: 38–46
    OpenUrlCrossRefPubMed
  56. ↵
    1. Sunkar R,
    2. Chinnusamy V,
    3. Zhu J,
    4. Zhu JK
    (2007) Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci 12: 301–309
    OpenUrlCrossRefPubMed
  57. ↵
    1. Sunkar R,
    2. Zhou X,
    3. Zheng Y,
    4. Zhang W,
    5. Zhu JK
    (2008) Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol 8: 25
    OpenUrlCrossRefPubMed
  58. ↵
    1. Varkonyi-Gasic E,
    2. Wu R,
    3. Wood M,
    4. Walton EF,
    5. Hellens RP
    (2007) Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods 3: 12
    OpenUrlCrossRefPubMed
  59. ↵
    1. Vazquez F,
    2. Vaucheret H,
    3. Rajagopalan R,
    4. Lepers C,
    5. Gasciolli V,
    6. Mallory AC,
    7. Hilbert JL,
    8. Bartel DP,
    9. Crété P
    (2004) Endogenous trans-acting siRNAs regulate the accumulation of Arabidopsis mRNAs. Mol Cell 16: 69–79
    OpenUrlCrossRefPubMed
  60. ↵
    1. Voinnet O
    (2009) Origin, biogenesis, and activity of plant microRNAs. Cell 136: 669–687
    OpenUrlCrossRefPubMed
  61. ↵
    1. Wang JW,
    2. Wang LJ,
    3. Mao YB,
    4. Cai WJ,
    5. Xue HW,
    6. Chen XY
    (2005) Control of root cap formation by microRNA-targeted auxin response factors in Arabidopsis. Plant Cell 17: 2204–2216
    OpenUrlAbstract/FREE Full Text
  62. ↵
    1. Wang L,
    2. Wang MB,
    3. Tu JX,
    4. Helliwell CA,
    5. Waterhouse PM,
    6. Dennis ES,
    7. Fu TD,
    8. Fan YL
    (2007) Cloning and characterization of microRNAs from Brassica napus. FEBS Lett 581: 3848–3856
    OpenUrlCrossRefPubMed
  63. ↵
    1. Wei S,
    2. Yu B,
    3. Gruber MY,
    4. Khachatourians GG,
    5. Hegedus DD,
    6. Hannoufa A
    (2010) Enhanced seed carotenoid levels and branching in transgenic Brassica napus expressing the Arabidopsis miR156b gene. J Agric Food Chem 58: 9572–9578
    OpenUrlCrossRefPubMed
  64. ↵
    1. Wollmann H,
    2. Weigel D
    (2010) Small RNAs in flower development. Eur J Cell Biol 89: 250–257
    OpenUrlCrossRefPubMed
  65. ↵
    1. Xie FL,
    2. Huang SQ,
    3. Guo K,
    4. Xiang AL,
    5. Zhu YY,
    6. Nie L,
    7. Yang ZM
    (2007) Computational identification of novel microRNAs and targets in Brassica napus. FEBS Lett 581: 1464–1474
    OpenUrlCrossRefPubMed
  66. ↵
    1. Xie Z,
    2. Johansen LK,
    3. Gustafson AM,
    4. Kasschau KD,
    5. Lellis AD,
    6. Zilberman D,
    7. Jacobsen SE,
    8. Carrington JC
    (2004) Genetic and functional diversification of small RNA pathways in plants. PLoS Biol 2: E104
    OpenUrlCrossRefPubMed
  67. ↵
    1. Yoshikawa M,
    2. Peragine A,
    3. Park MY,
    4. Poethig RS
    (2005) A pathway for the biogenesis of trans-acting siRNAs in Arabidopsis. Genes Dev 19: 2164–2175
    OpenUrlAbstract/FREE Full Text
  68. ↵
    1. Zhang X,
    2. Henderson IR,
    3. Lu C,
    4. Green PJ,
    5. Jacobsen SE
    (2007) Role of RNA polymerase IV in plant small RNA metabolism. Proc Natl Acad Sci USA 104: 4536–4541
    OpenUrlAbstract/FREE Full Text
  69. ↵
    1. Zhou X,
    2. Sunkar R,
    3. Jin H,
    4. Zhu JK,
    5. Zhang W
    (2009) Genome-wide identification and analysis of small RNAs originated from natural antisense transcripts in Oryza sativa. Genome Res 19: 70–78
    OpenUrlAbstract/FREE Full Text
  70. ↵
    1. Zhu QH,
    2. Spriggs A,
    3. Matthew L,
    4. Fan L,
    5. Kennedy G,
    6. Gubler F,
    7. Helliwell C
    (2008) A diverse set of microRNAs and microRNA-like small RNAs in developing rice grains. Genome Res 18: 1456–1465
    OpenUrlAbstract/FREE Full Text
View Abstract
PreviousNext
Back to top

Table of Contents

Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Plant Physiology.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Small RNA Profiling in Two Brassica napus Cultivars Identifies MicroRNAs with Oil Production- and Development-Correlated Expression and New Small RNA Classes
(Your Name) has sent you a message from Plant Physiology
(Your Name) thought you would like to see the Plant Physiology web site.
Citation Tools
Small RNA Profiling in Two Brassica napus Cultivars Identifies MicroRNAs with Oil Production- and Development-Correlated Expression and New Small RNA Classes
Ying-Tao Zhao, Meng Wang, San-Xiong Fu, Wei-Cai Yang, Cun-Kou Qi, Xiu-Jie Wang
Plant Physiology Feb 2012, 158 (2) 813-823; DOI: 10.1104/pp.111.187666

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Small RNA Profiling in Two Brassica napus Cultivars Identifies MicroRNAs with Oil Production- and Development-Correlated Expression and New Small RNA Classes
Ying-Tao Zhao, Meng Wang, San-Xiong Fu, Wei-Cai Yang, Cun-Kou Qi, Xiu-Jie Wang
Plant Physiology Feb 2012, 158 (2) 813-823; DOI: 10.1104/pp.111.187666
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • RESULTS
    • DISCUSSION
    • CONCLUSION
    • MATERIALS AND METHODS
    • Acknowledgments
    • Footnotes
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

In this issue

Plant Physiology: 158 (2)
Plant Physiology
Vol. 158, Issue 2
Feb 2012
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
View this article with LENS

More in this TOC Section

  • Natural Variation for Seed Longevity and Seed Dormancy Are Negatively Correlated in Arabidopsis
  • Genomics and Localization of the Arabidopsis DHHC-Cysteine-Rich Domain S-Acyltransferase Protein Family
  • Allelic Variation in Paralogs of GDP-l-Galactose Phosphorylase Is a Major Determinant of Vitamin C Concentrations in Apple Fruit
Show more Article

Similar Articles

Our Content

  • Home
  • Current Issue
  • Plant Physiology Preview
  • Archive
  • Focus Collections
  • Classic Collections
  • The Plant Cell
  • Plant Direct
  • Plantae
  • ASPB

For Authors

  • Instructions
  • Submit a Manuscript
  • Editorial Board and Staff
  • Policies
  • Recognizing our Authors

For Reviewers

  • Instructions
  • Journal Miles
  • Policies

Other Services

  • Permissions
  • Librarian resources
  • Advertise in our journals
  • Alerts
  • RSS Feeds

Copyright © 2019 by The American Society of Plant Biologists

Powered by HighWire