Deep sequencing on genome-wide scale reveals the unique composition and expression patterns of microRNAs in developing pollen of Oryza sativa
© Wei et al.; licensee BioMed Central Ltd. 2011
Received: 13 December 2010
Accepted: 16 June 2011
Published: 16 June 2011
Pollen development in flowering plants requires strict control of the gene expression program and genetic information stability by mechanisms possibly including the miRNA pathway. However, our understanding of the miRNA pathway in pollen development remains limited, and the dynamic profile of miRNAs in developing pollen is unknown.
Using next-generation sequencing technology, we pyrosequenced small RNA populations from rice uninucleate microspores to tricellular pollen and control sporophytic tissues at the genome-wide level. We identified 292 known miRNAs, including members of all 20 families conserved in plants, and 75 novel miRNAs. Of the 292 known miRNAs, 202 were expressed, with 103 enriched, in developing pollen. More than half of these novel miRNAs displayed pollen-or stage-specific expression. Furthermore, analyzing the 367 miRNAs and their predicted targets indicated that correlation in expression profiles of pollen-enriched known miRNAs and their targets significantly differs from that of sporophyte-enriched known miRNAs and their targets in some functional terms, while novel miRNAs appeared to negatively regulate their targets. Importantly, gene ontology abundance analysis demonstrated chromatin assembly and disassembly was important in the targets of bicellular pollen-expressed novel miRNAs. Principal component analysis revealed pollen of all three stages was discriminated from sporophytes, largely because of the novel and non-conserved known miRNAs.
Our study, for the first time, revealed the differences in composition and expression profiles of miRNAs between developing pollen and sporophytes, with novel and non-conserved known miRNAs the main contributors. Our results suggest the important roles of the miRNA pathway in pollen development.
MicroRNAs (miRNAs) and small interfering RNAs (siRNAs) are two types of small non-coding RNAs (20 to 24 nucleotides in length) identified in nearly all eukaryotes. The pool of small RNAs in plants is highly complex, consisting primarily of many low-abundant siRNAs and a small number of highly expressed 21-nucleotide sequences; most of the latter are miRNAs [1, 2]. Most miRNA loci are encoded by independent transcriptional units in intergenic regions that are transcribed by RNA polymerase II. In plants, miRNAs are processed from stem-loop regions of long primary transcripts by a Dicer-like enzyme and are loaded into silencing complexes, where they generally direct cleavage of complementary mRNAs. Although miRNAs were identified in plants just recently, studies have revealed that miRNAs play crucial roles in each major stage of plant development, often targeting the transcription factors that mediate transition from one developmental stage to the next .
Haploid pollen (also called the gametophyte) is a key regulator of sexual reproduction in flowering plants and is produced from diploid pollen mother cells via meiosis. In contrast to animals, in which products of meiosis directly develop into sperm cells, in plants, the product of meiosis undergoes a unique postmeiotic pollen development process, finally giving rise to sperm cells. During this process, haploid uninucleate microspores (UNMs) generated from meiosis first undergo asymmetric mitosis to generate bicellular pollen (BCP) consisting of a large vegetative cell and a small generative cell enclosed in the vegetative cell. The two types of cells have different fates: the vegetative cell exits the cell cycle and can develop into a polarly growing pollen tube, whereas the generative cell undergoes further mitosis to produce two sperm cells. This postmeiotic pollen development is orderly and precisely regulated ; however, the mechanism underlying the main developmental events remains largely unknown.
Recent transcriptomics studies have revealed that the number of genes expressed in pollen greatly decreases from the UNM stage to the tricellular pollen (TCP) stage, whereas stage-specific transcripts showed a 'U-type' change, with the lowest number at the BCP stage in Oryza sativa and Arabidopsis thaliana . These data suggest that fine-tuned gene expression and function are essential to guarantee pollen development. The miRNA pathways in pollen are of interest because previous transcriptomic analysis showed key transcripts involved in miRNA pathway, such as those encoding Argonaute (AGO)1, 2, 4 and 7, Dicer-like protein (DCL)1 to 3 and RNA-dependent RNA polymerase 1, 2 and 6, were turned off during pollen development . However, recent research demonstrated that several miRNA pathway genes were expressed and some enriched in pollen grains of Arabidopsis [7–10]. Many miRNAs known to function in somatic development have been identified in Arabidopsis mature pollen [9, 11]. Using 454 sequencing, Grant-Downton et al.  revealed diverse small RNAs in mature pollen of Arabidopsis, some of which were possibly specific to mature pollen. These studies indicate the existence of an miRNA pathway in pollen. However, mature pollen is terminally differentiated and at a developmental 'standstill', so the present studies could not define the importance of miRNAs in pollen development. Sequencing small RNAs from developing pollen of different stages is necessary to obtain a global picture of the temporal dynamics of small RNA diversity . Therefore, genome-wide knowledge about the composition of miRNAs and dynamic changes in the miRNAs during pollen development is important to understand the mechanism of fine-tuned pollen development.
Rice (O. sativa) is one of the most important cereal crops; it feeds half of the world's population and has been used as an excellent model system for studying monocots after Arabidopsis. Besides the importance of rice, knowledge of the molecular mechanisms underlying rice pollen development is essential to manipulate male fertility for heterosis utilization. In this study, we used Solexa high-throughout sequencing technology to sequence the small RNA population from developing rice pollen at the UNM, BCP and TCP stages, with sporophytic tissues-roots, leaves, and callus cells-used as controls. We obtained millions of high-quality readouts from each sample. Further analyses identified 292 known miRNAs, including members of all 20 families conserved in plants, and 75 novel miRNAs. Most of the known miRNAs (hereafter called kn-miRs) were expressed in developing pollen, and most of the predicted novel miRNAs (nov-miRs) were pollen specific. We also predicted 1,353 genes possibly targeted by the 367 identified miRNAs and correlated the expression profiles of miRNAs and their targets. We supply novel insights into the dynamic profiles of miRNAs in developing pollen.
Sequencing and data analysis
Statistics of small RNA sequences from the individual libraries
High quality readsa
Sequences of 18 to 30 nucleotidesa
Sequences matched to the genomeb
Repeat-associated small RNAsb
Un-annotated small RNAsb
In general, the small RNA library generated by sequencing was complex in composition. Besides the miRNAs and siRNAs, the library includes large numbers of degradation fragments derived from other coding and noncoding transcripts [13, 14]. Therefore, to annotate small RNAs, we first mapped these small RNAs of 18 to 30 nucleotides to the rice genome using SOAP . More than 80% of the small RNAs mapped perfectly to the genome (Table 1). Furthermore, we removed sequences corresponding to known non-coding RNAs (for example, rRNAs, tRNAs, small nuclear RNAs (snRNAs) and small nucleolar RNAs (snoRNAs); Table 1) and possibly degraded species of mRNAs. Finally, we obtained 1,825,534, 3,383,556 and 1,082,758 small RNAs from UNMs, BCP and TCP, respectively, and 3,866,888, 1,862,254 and 3,044,145 from callus cells, leaves and roots, respectively, which were candidates for identifying miRNAs (including the kn-miRs, repeat-associated small RNAs and un-annotated small RNAs in Table 1).
By mapping these candidates to their precursor sequences in the rice miRNA database available in miRBase [16–18], we identified 292 kn-miRs in pollen and/or sporophytic samples (Additional files 1 and 2); 202 of these kn-miRs were expressed in pollen samples, with 141 in UNMs, 152 in BCP and 143 in TCP; 281 were expressed in sporophytic tissues, with 210 in callus cells, 239 in roots and 223 in leaves. In addition, we detected 20 conserved kn-miR families in these libraries, and 16 of the families, apart from miR394, miR395, miR398 and miR408, were sequenced in developing rice pollen (Additional files 1 and 2). Thus, to a certain extent, our sequencing depth was sufficient to reflect the expression profiles of miRNAs during pollen development, and most of the conserved miRNAs were also expressed in the developing pollen.
Recent studies have revealed that some small RNAs derived from highly repeated elements bind with different AGO proteins and are involved in important biological processes, such as chromatin maintenance and transposon control [19, 20]. We used RepeatMasker software to identify small RNAs positioned at repeat loci and annotated them as repeat-associated small RNAs (Table 1). The remaining un-annotated small RNAs were further used to predict nov-miRs.
Prediction of nov-miRs
miRNAs are derived from hairpin-like precursors, originating from a single-stranded RNA transcript through sequential processing by Dicer or Dicer-like (DCL) proteins [21, 22]. miRNA precursors have a characteristic fold-back structure, which is the primary criterion to annotate nov-miRs . Therefore, we predicted nov-miRs as follows. First, by folding the flanking genome sequence of the above un-annotated small RNAs, followed by analysis of structural features, we excluded small RNAs that cannot form the characteristic fold-back structure. Second, recently evolved/evolving miRNAs have a single locus in the genome [24–26], so small RNAs with multiple loci in the rice genome were excluded. Third, to minimize noise, we also eliminated small RNAs of low abundance (with total number of reads fewer than five) and those originating from both strands, which would generate siRNA-like small RNAs. Fourth, all the remaining un-annotated small RNA sequences were subjected to 'MIREAP', which recovers most kn-miRs with only a few exceptions whose structures cannot satisfy the common features of an miRNA gene . Finally, to distinguish miRNAs from miniature inverted repeat transposable elements (MITEs), we blasted the precursor and mature sequences of the small RNAs with characteristic fold-back structure against the Oryza Repeat Database , and the homologs of repetitive sequences were discarded. From the above analyses, our predicated nov-miRs satisfied the following criteria: precursors had a characteristic fold-back structure, contained no repetitive sequences, and matched the genome only once, and most were located in the intergenic region; lengths of mature miRNA ranged from 20 to 24 nucleotides, and the number of reads was greater than five; the mature sequences could be sequenced in two or more libraries, or the miRNA* sequence could be identified in at least one library; and the targets of predicted miRNAs could be predicted using an upgraded version of miRU .
In total, we obtained 75 predicted nov-miRs (Additional file 3). Of these, 30 were expressed in UNMs, 39 in BCP, and 18 in TCP, and only 14, 18 and 12 in callus cells, leaves and roots, respectively (Additional file 4). Compared with the increased number of kn-miRs expressed in sporophytic tissues, more nov-miRs were identified in developing rice pollen, so more miRNAs still remain to be revealed in gametophytes.
Expression profiles of miRNAs during pollen development
The expression features of these nov-miRs differ from those of kn-miRs (Additional files 5 and 6). More than half of the nov-miRs were expressed only in pollen libraries and some of them only in one pollen library (clusters 1 to 3 and 6 to 8 in Figure 2b). Nov-miRs in clusters 6, 7 and 8 (Figure 2b) displayed UNM-, BCP-and TCP-specific expression, respectively. These should represent a set of pollen-or stage-specific miRNAs. Therefore, in contrast to more kn-miRs being expressed constitutively or preferentially in sporophytic tissues, more nov-miRs were expressed preferentially and even specifically in pollen.
Target prediction and function analysis
Finding regulatory mRNA targets is essential for understanding the biological functions of miRNAs. We predicted targets of all identified miRNAs using the psRNATarget web service  and found 1,068 targets for 292 kn-miRs (Additional file 8a) and 285 targets for 75 nov-miRs (Additional file 8b). In total, we predicted 1,353 genes possibly targeted by the identified 367 miRNAs, with an average of 3.7 targets per miRNA, and each kn-miR (292 versus 1,068 targets) and nov-miR (75 versus 285) appeared to have similar numbers of targets on average.
The analysis of kn-miR targets showed that a high proportion of the targets were transcription factors (Additional file 8a). Besides transcription factor targets, targets involving defense response, hormone regulation and other metabolic pathways were overrepresented (Additional file 8a). This result echoes previous research on kn-miRs and their targets in sporophytic tissues [3, 30]. Furthermore, Gene Ontology (GO) abundance analysis of targets of both pollen-and sporophyte-enriched kn-miRs revealed terms related to transcription regulation, transcription factor activity and metabolic regulation that were all statistically significant (false discovery rate P < 0.01; Additional file 9). However, hormone-signaling-related genes were overrepresented in pollen-enriched kn-miR targets, whereas lignin metabolism-related genes were overrepresented in sporophyte-enriched kn-miR targets (Additional file 9).
miRNA-induced cleavage of predicted targets
Our observation that most of nov-miRs and their targets had negatively correlated expression (Additional file 10b) implies that these nov-miRs have potential cleavage activity. To validate the cleavage events of nov-miRs, we amplified their predicted target genes using rapid amplification of 5' cDNA ends (5' RACE).
Although, to date, the canonical cleavage site of almost all of the kn-miRs was between nucleotides 10 and 11, we unexpectedly found that miR56 and miR58 cleaved targets with high frequency (10 of 10 for miR56, 5 of 16 for miR58) between nucleotides 18 and 19 and 5 and 6, respectively (Additional file 11). These cleavage events were validated by transient co-expression in N. benthamiana leaves (Additional file 11). Furthermore, cleavage of the non-conserved kn-miR osa-miR820 occurred primarily at the canonical position with a frequency of nine out of ten in a cDNA pool from seedling and inflorescence  or seven out of eight in a mixture of cDNA from panicle and embryogenic calli  but infrequently (one of ten) between nucleotides 11 and 12 from the complementary 5' end of osa-miR820 . However, we found that osa-miR820-mediated cleavage of Os03g02010 was predominantly at nucleotide 7 from the paired 3' end (5 of 14) from the mixed cDNA pool from developmental pollen, and no other degradation fragments were detected with the targeted sequence (data not show). As indicated in a study of Arabidopsis, the cleavage frequency and cleavage sites of targets for several conserved kn-miRs differed between pollen and vegetative tissues , which suggests that the differences might account for subtle regulatory mechanism.
Recently developed degradome sequencing is a useful and powerful tool to discover new targets for conserved and non-conserved miRNAs and to validate putative targets of predicted nov-miRs [33–37]. We searched degradome data for 3-week-old seedlings and young inflorescences of rice [36, 37] and starBase (sRNA target Base)  to identify the cleavage products of nov-miRs. From the degradome data, we found eight predicted targets with only one count each for eight different nov-miRs (Additional file 12), and from starBase, three predicted target fragments for three different nov-miRs were identified. Among the eleven nov-miRs, three were expressed only in sporophytic tissues and eight were enriched in developing pollen (two in UNMs, two in TCP) or in two successive stages (two in UNMs and BCP), or in both sporophytic tissues and pollen (two of five examined samples except TCP) (Additional file 12). Because of the extremely low frequency of degradation products of the targets mentioned above, detecting them by 5' RACE is almost impossible; indeed, we detected no products for them.
Using next-generation high-throughout sequencing technology, we have analyzed small RNAs in developing rice pollen from the UNM to TCP stages at the genome-wide level. This analysis revealed dynamic features of small RNA populations in the developing pollen and expression patterns of pollen-expressed miRNAs. We also revealed a set of target genes of pollen-expressed miRNAs and the possible relationship between pollen-expressed miRNAs and their targets by analysis of their expression patterns. These results provide novel insights into molecular regulation mediated by miRNAs of pollen development.
Dynamic changes in small RNA populations in developing pollen
In plants, transcriptional and post-transcriptional gene silencing is directed by genome-encoded 21-to 24-nucleotide small RNAs. miRNAs of 21 nucleotides and trans-acting siRNAs of 21 nucleotides function in post-transcriptional gene silencing by guiding mRNA degradation or translational repression [1, 39–41]. siRNAs of 24 nucleotides are implicated in DNA and histone modifications leading to transcriptional gene silencing [40–43], and 24-nucleotide miRNAs involved in DNA methylation were also reported recently . Recent studies have revealed that the size distribution of small RNAs is tissue-specific in maize . Mature pollen of Arabidopsis loses most 24-nucleotide siRNAs and gains some 21-nucleotide siRNAs, which target silencing in the pollen . Using high-throughout Solexa sequencing, we obtained more than 5 million high-quality small RNAs from pollen samples at three developmental stages and from sporophytic tissues. The size distribution of small RNAs varied with pollen development. The number of small RNAs of 24 nucleotides, which would be mainly siRNAs, peaked in UNMs and BCP, with more 20-and 21-nucleotide small RNAs (mainly miRNAs) in TCP (Figure 1). This finding suggests that the shift in small RNA populations in the developing pollen may represent a mechanism regulating the gene expression program and finally pollen development. Also, transcriptional regulation would be an important mechanism in the early phase, but post-transcription regulation would be prevalent in the late phase of pollen development in rice.
Developing pollen has a unique expression pattern of miRNAs
miRNAs are a highly diverse class of small RNAs; hundreds of mature sequences have been registered in miRBase [16–18]. However, almost all of them, if not all, were identified originally from sporophytes. We knew little about the miRNAs in pollen. miR164 and miR171 were first discovered to be expressed in pollen of Nicotiana by northern blot analysis , and miR164 was further identified in pollen of Arabidopsis . Recent studies of Arabidopsis revealed the expression of some kn-miRs in mature pollen  and the expression of critical components of small RNA pathways in developing pollen [9–11]. Our study revealed 292 kn-miRs from pollen and sporophytic tissue libraries, 202 of which were expressed in developing pollen. Among the 202 pollen-expressed kn-miRs, many (103) were even pollen enriched (see Results for details). These results indicate that most of the kn-miRs originally identified from sporophytes were also expressed or even enriched in developing pollen. In contrast, most kn-miRs detected in Arabidposis mature pollen were in low abundance , possibly because the mature pollen is in a terminal development status ready for fertilization and needs no more miRNAs for regulation. Therefore, the existence and high expression of kn-miRs in developing rice pollen demonstrates their importance during pollen development.
Impressively, among 20 miRNA families now recognized to be conserved in diverse plant species [48, 49], 16 were identified in developing pollen and sporophytic tissues in this study, and 6 (miR156, miR159, miR164, miR166, miR167 and miR396) were accumulated highly throughout all examined samples (clusters 4, 5 and 9 in Figure 2a), which implies that these conserved miRNAs are conserved among species and among sporophytes and pollen, and among developmental stages of pollen. These conserved miRNAs may be integral regulators of both sporophytic and gametophytic development, and some have house-keeping functions in plant cells.
Although most studies focus on the more conserved and nonspecific miRNA families in plants, given the fast evolution of miRNA sequences, several studies have demonstrated the presence of recently evolved plant species-specific miRNAs . Our work revealed 75 nov-miRs in the 6 libraries. Unexpectedly, more than half of these nov-miRs were expressed only in developing rice pollen (clusters 1 to 3 and 6 to 8 in Figure 2b), and many (clusters 6 to 8 in Figure 2b) were detected only in pollen samples at individual developmental stages, which demonstrates the existence of pollen-specific miRNAs at one or more stages, at least in rice. Vascular plants evolved increasingly complex sporophytes with reduced gametophyte complexity and size, and ultimately loss of independence , and the gametophyte in flowering plants has been reduced to just a few cells in pollen grains. The high proportion of nov-miRs specific to pollen may be due to the need for this evolved development of the gametophyte.
Furthermore, we found differential expression of miRNAs in developing pollen and sporophytic tissues. Principal component analysis revealed gametophytes are distinct from sporophytes with regard to the expression profiles of miRNAs, with nov-miRs and non-conserved kn-miRs as the major contributors to this distinction (Figure 3). Therefore, miRNA expression features in pollen differ significantly from those in sporophytic tissues, which suggest a unique regulation program in developing pollen. Consistent with animal sperm cells, plant gametophytes function to faithfully maintain selfish genetic information and developmental potential in establishing the subsequent generation of an individual. Several studies of mammalian organisms have demonstrated that miRNAs enriched in sperm cells are crucial for global regulation of the germ cell developmental program and for keeping selfish genetic elements under strict surveillance . Together, these data demonstrate the involvement of miRNA-mediated gene expression regulation in pollen development.
Roles of miRNAs in pollen development
Post-meiotic pollen development involves a series of fine-tuned, coordinated cellular events associated with unique transcriptomic profiles, along with differential expression of a complex set of transcription factors [5, 53], which indicates crucial roles of gene expression regulation in the development. Recent studies have shown the key function of small RNAs, particularly miRNAs and siRNAs, in regulating gene expression in both plant and animal species [1, 3, 41] and their crucial roles in regulating sperm cell development and maintaining stability of the genetic element of the cell in mammalian organisms . Many predicted targets encode transcription factors, suggesting that plant miRNAs are master regulators. Our present work found BCP expressed more miRNAs, including known and novel miRNAs (152 kn-miRs and 39 nov-miRs) than UNMs (141 kn-miRs and 30 nov-miRs) and TCP (143 kn-miRs and 18 nov-miRs), and previous research showed that BCP expressed the smallest number of transcription factors , indicating that miRNAs cooperate with transcription factors to achieve the trait of gene expression during pollen development, with BCP having the smallest number of stage-enriched transcripts during rice pollen development .
Furthermore, a growing number of miRNA-target pairs have been confirmed experimentally (reviewed by [3, 41]) and play crucial roles in plant development. For example, overexpression of miR159 repressed mRNA levels of MYB33 and MYB65 and induced male sterility as well as delayed flowering under short days [54, 55]; expression of miR167-resistant ARF6 leads to arrested ovule development and indehiscent anthers ; NAC1-miR164 pair mediates auxin signaling in lateral root emergence ; TIR1, encoding an auxin receptor, and related F-box genes are regulated by miR393 . Our study identified all of these miRNAs and found some of them were enriched in developing pollen (Additional file 1), and the corresponding targets have already been predicted. We also found that hormone-signaling-related genes were overrepresented in a set of targets of pollen-enriched kn-miRs (Additional file 8a). Hormone signaling in pollen plays crucial roles in coordinating pollen development and maturation [59–63]. Numerous studies have revealed crucial roles of PPR proteins and F-box proteins in pollen development and the importance of chromatin remodeling and DNA repair in maintaining genetic element stability [5, 64–67]. Genes associated with these functions were also predicted to be targets of pollen-enriched kn-miRs, suggesting their important roles possibly during pollen development. However, the correlation of expression profiles of pollen-enriched kn-miRs with their corresponding targets was significantly different from that of sporophyte-enriched kn-miRs with their targets (Figure 5). It is possible, therefore, that some sporophyte-originated kn-miRs have roles in different mechanisms when expressed preferentially in pollen.
In contrast to kn-miRs, more than half of the identified nov-miRs were expressed only in developing pollen, with several being stage-specific, and their targets were more versatile. Impressively, GO abundance analysis revealed that terms related to chromatin assembly and disassembly were over-represented in the set of targets of BCP-expressed nov-miRs (Figure 6). Chromatin assembly and disassembly are fundamentally important processes that are tightly linked to DNA replication and transcription, thus ensuring that cells faithfully duplicate chromosomes. This process is essential for developing pollen to coordinate germ cell development and strict maintenance of genome integrity and stability.
Several mechanisms involve chromatin assembly and disassembly, including chromatin remodeling and modification of core histones or cytosines in CpG-rich regions of the genome. Recent studies revealed histone remodeling-related genes targeted by miRNAs expressed in mature pollen of Arabidopsis, such as miR778, which targets transcripts of two histone methyltransferase genes, SUVH5 and SUVH6 [10, 68]. A gene encoding the SET domain protein OsSET1/OsiEZ1 was verified to be the target of BCP-specific miR52 and was negatively regulated by this miRNA (Figure 7). SET-domain proteins, chiefly responsible for lysine methylation of histone H3, such as SDG4/ASHR3, are involved in pollen development . Furthermore, our 5' RACE analysis, in combination with several previous reports [31, 32, 44], confirmed that osa-miR820, which is enriched in pollen, targeted a DNA cytosine methyltransferase gene (Os03g02010). Our study also revealed that osa-miR827, which is highly expressed in pollen, targeted Os04g11510, which encodes a methyl-CpG binding domain protein (Additional file 11). Both of these proteins are important to regulate the methylation of a genome [70–73]. The epigenetic modification of histones and cytosines in the genome of mammalian germ lines is essential for sperm cell development and plays roles in parental imprinting [74, 75]. These results indicate that epigenetic regulation of chromatin assembly and disassembly may be an important mechanism underlying pollen development and may involve parental imprinting. Further function studies of genes are needed to clarify these roles.
Using Solexa high-throughout sequencing technology, we sequenced the small RNA population of rice pollen at three sequential developmental stages from microspores to tricellular pollen, with sporophytic tissues-roots, leaves, and callus cells-as controls. We obtained millions of high-quality readouts from each sample and identified 292 kn-miRs and 75 nov-miRs. The miRNA composition and expression pattern of developing pollen were obviously different from those of sporophytes, with more nov-miRs enriched/specifically expressed in pollen while more kn-miRs were enriched/specifically expressed in sporophytes. Principal component analysis revealed that pollen could be differentiated from sporophytes with regard to miRNA expression profiles, with novel and non-conserved known miRNAs the main contributors to this. Furthermore, 1,068 targets were predicted for 292 known miRNAs; although no obvious differences were found by GO abundance analysis of those targets between pollen-enriched and sporophyte-enriched kn-miRs, correlation of expression profiles of pollen-enriched kn-miRs with their targets significantly differs from that of sporophyte-enriched kn-miRs with their corresponding targets in terms of transcription, hormone signaling and chromatin remodeling. We identified 285 targets for the 75 nov-miRs, which appeared to be negatively regulated. GO terms for chromatin assembly and disassembly were statistically significantly associated with the targets of BCP-expressed nov-miRs, implying that BCP would be the key point of miRNA regulation. Our data reveal for the first time comprehensive and dynamic features of miRNAs in developing pollen.
Materials and methods
Rice cultivar Zhonghua 10 (O. sativa L. ssp. japonica) was planted in a climate chamber under a 12-hour light/12-hour dark cycle at 28°C for 2 weeks, then the roots and leaves were collected. One-month-old callus cells were generated from rice embryos on N6 solid medium containing 2, 4-Dichlorophenoxyacetic acid (2 mg/l) in the dark at 25°C. Pollen at UNM, BCP and TCP stages was obtained as described previously .
Small RNA extraction
Total RNA was extracted from sporophytic tissues and pollen at each developmental stage by use of RNAplant reagents (Tiangen Biotech, Beijing, China) according to the manufacturer's instructions. Total RNA was resolved on a denatured 15% polyacrylamide gel, then the fraction from 18 to 30 nucleotides was collected.
The sequencing was performed as described . In brief, purified small RNA molecules were ligated to a pair of Solexa adaptors at the 5' and 3' ends and amplified by use of adaptor primers for 15 cycles to produce sequencing libraries. PCR products were purified and small RNA libraries were sequenced by use of Solexa.
Computational analysis of sequencing data
Small RNA reads were produced by use of the Illumina 1 G Genome Analyzer (Illumina, San Diego, CA, USA). Our raw data have been deposited in the European Nucleotide Archive by ArrayExpress  [ArrayExpress: E-MTAB-689]. After filtering the low quality reads, and trimming adaptor sequences by a modified dynamic programming algorithm , we collected short RNAs ranging from 18 to 30 nucleotides and drew size distributions. Sequences of 18 or more nucleotides were mapped to the rice genome (TIGR version 5.0) by SOAP v1.11 . Sequences perfectly matching the genome along their entire length were considered for subsequent analyses. Sequences matching known rice rRNAs, tRNAs, snRNAs and snoRNAs in the Rfam RNA family database [79–81] and NCBI GenBank database  were discarded. With the available annotation of repeat sequences in the rice genome (TIGR Oryza Repeat Database version 3.3) and use of a parallel overlap finding algorithm, small RNAs positioned at repeat loci were identified and annotated as repeat-associated small RNAs. Degraded species of mRNAs can also be sampled out by Solexa technology, and sequences overlapping with gene regions were excluded.
Sequences mapped to miRNA precursors found by BLAST search of the miRBase miRNA Database (v16.0) [16–18] were identified as kn-miRs, with a manual check to delete the false miRNAs. miRNA precursors have a characteristic fold-back structure, which can be used to predict nov-miRs. By folding the flanking genome sequence of small RNAs, then analyzing its structural features, we could identify novel miRNA candidates by use of MIREAP . Good candidate sequences were then submitted to a target prediction web server containing plant miRNAs .
Principal component analysis of miRNAs
The normalized reads (TPM) of all miRNAs detected and predicted in this study were transformed with use of the log2 scale and analyzed by SIMCA-P+ 12.0 with PCA-X (Ctr scaling) following the PLS-DA model to obtain the principal components. miRNAs with a score contribution more than the absolute value of 2 were illustrated by comparing different samples or tissue types.
Quantitative RT-PCR of miRNAs
Total RNA was isolated from pollen at each developmental stage and from sporophytic tissues. To prepare small RNA, 0.5 M NaCl and 5% PEG8000 were used to precipitate and remove high-molecular-weight RNA. Small RNA in the resulting supernatant was pelleted with a 1/10 volume of 3 M sodium acetate and 3 volumes of ethanol . The concentration of small RNAs was determined by use of the Beckman Coulter DU730 Nucleic Acid/Protein Analyzer (Fullerton, CA, USA). cDNA was synthesized in a 10-μl reverse transcription (RT) reaction using of 100 ng purified small RNA, 1 U ReverTra Ace reverse transcriptase (Toyobo, Osaka, Japan), 2 μl of 5 × RT buffer, 1 mM dNTPs and 2.4 U Ribnuclease inhibitor (TAKARA SHIZU CO., LTD, Tokyo, Japan), and 50 nM miRNA-specific stem-loop primers (Additional file 13) designed according to Chen et al. . U6 snRNA was chosen as a reference and was reverse-transcribed with reverse primer OsU6R [84, 85]. A pulsed RT procedure was used to increase the reaction specification . In brief, the RT reaction mixture was incubated at 16°C for 30 minutes, then 60 cycles for 30 s at 30°C, 30 s at 42°C, and 1 s at 50°C, finally at 95°C for 5 minutes. No RNA or RT primers or RT controls were set at the same time. SYBR® Green Realtime PCR Master Mix (Toyobo) was used to detect miRNA expression by a Stratagene Max 3000p Detection System (La Jolla, CA, USA). Briefly, cDNAs were diluted 5 times and 1 μl diluted product was used as a template in a 10-μl PCR reaction, which contained 5 μl 2 × SYBR Green Realtime PCR Master Mix and 0.25 μM of an miRNA-specific forward primer and universal reverse primer. The quantitative PCR was conducted in duplicate for 90 s at 95°C, then 40 cycles of 15 s at 95°C and 10 s at 60°C. For each PCR, dissociation curve analysis was carried out to discriminate the specific products from the primer dimers. The fold changes of miRNA in different samples were calculated by the 2Δ ΔCt method as described .
5' RACE of miRNA cleavage
Total RNA (1 μg) from pollen of equally mixed UNMs, BCP and TCP was used to synthesize 5'-RACE-ready cDNAs with N-15 random primer mix and BD Smart RACE cDNA Amplification Kit (Clontech, Palo Alto, CA, USA) according to the manufacturer's instruction. The first round of PCR involved 10 × UPM, outer gene-specific primers and Advantage 2 Polymerase Mix (Clontech). The product was diluted 50 times and then used as a template for the second round of PCR, which involved NUP and outer/inner gene-specific primers. Amplicons were separated on the gel, cloned into pMD 19-T vector (Takara) and sequenced. The outer and inner gene-specific primers were listed in Additional file 13.
Transient co-expression of miRNAs and their targets in N. benthamianaleaves
Rice genomic fragments forming fold-back structure for precursors of miR52, miR56, miR58 and osa-miR827a carrying the miRNA and/or miRNA* duplex were amplified with the primers listed in Additional file 13 and cloned into a binary vector (pTCK303) driven by the maize ubiquitin promoter. Their corresponding target fragments obtained with primers listed in Additional file 13 were inserted into the same vector. The plasmids carrying precursors of miRNAs and targets were infiltrated with separate or mixed cells as described [37, 87]. The total RNA of infiltrated leaves was isolated after 4-day growth and used to synthesize the first cDNA as templates of semi-quantitative real-time PCR as mentioned above.
rapid amplification of cDNA ends
small interfering RNA
small nuclear RNA
small nucleolar RNA
transcripts per million
This work was supported by grants from the Chinese Ministry of Science and Technology (number 2007CB947601) and the National Natural Science Foundation of China (numbers 30821007 and 31000605).
- Voinnet O: Origin, biogenesis, and activity of plant microRNAs. Cell. 2009, 136: 669-687. 10.1016/j.cell.2009.01.046.PubMedView ArticleGoogle Scholar
- Schwach F, Moxon S, Moulton V, Dalmay T: Deciphering the diversity of small RNAs in plants: the long and short of it. Brief Funct Genomic Proteomic. 2009, 8: 472-481. 10.1093/bfgp/elp024.PubMedView ArticleGoogle Scholar
- Jones-Rhoades MW, Bartel DP, Bartel B: MicroRNAS and their regulatory roles in plants. Annu Rev Plant Biol. 2006, 57: 19-53. 10.1146/annurev.arplant.57.032905.105218.PubMedView ArticleGoogle Scholar
- Borg M, Brownfield L, Twell D: Male gametophyte development: a molecular perspective. J Exp Bot. 2009, 60: 1465-1478. 10.1093/jxb/ern355.PubMedView ArticleGoogle Scholar
- Wei LQ, Xu WY, Deng ZY, Su Z, Xue Y, Wang T: Genome-scale analysis and comparison of gene expression profiles in developing and germinated pollen in Oryza sativa. BMC Genomics. 2010, 11: 338-10.1186/1471-2164-11-338.PubMedPubMed CentralView ArticleGoogle Scholar
- Pina C, Pinto F, Feijo JA, Becker D: Gene family analysis of the Arabidopsis pollen transcriptome reveals biological implications for cell growth, division control, and gene expression regulation. Plant Physiol. 2005, 138: 744-756. 10.1104/pp.104.057935.PubMedPubMed CentralView ArticleGoogle Scholar
- Slotkin RK, Vaughn M, Borges F, Tanurdzic M, Becker JD, Feijo JA, Martienssen RA: Epigenetic reprogramming and small RNA silencing of transposable elements in pollen. Cell. 2009, 136: 461-472. 10.1016/j.cell.2008.12.038.PubMedPubMed CentralView ArticleGoogle Scholar
- Borges F, Gomes G, Gardner R, Moreno N, McCormick S, Feijo JA, Becker JD: Comparative transcriptomics of Arabidopsis sperm cells. Plant Physiol. 2008, 148: 1168-1181. 10.1104/pp.108.125229.PubMedPubMed CentralView ArticleGoogle Scholar
- Grant-Downton R, Hafidh S, Twell D, Dickinson HG: Small RNA pathways are present and functional in the angiosperm male gametophyte. Mol Plant. 2009, 2: 500-512. 10.1093/mp/ssp003.PubMedView ArticleGoogle Scholar
- Grant-Downton R, Trionnaire GL, Schmid R, Rodriguez-Enriquez J, Said Hafidh2 SM, Twell D, Dickinson H: MicroRNA and tasiRNA diversity in mature pollen of Arabidopsis thaliana. BMC Genomics. 2009, 10: 643-10.1186/1471-2164-10-643.PubMedPubMed CentralView ArticleGoogle Scholar
- Chambers C, Shuai B: Profiling microRNA expression in Arabidopsis pollen using microRNA array and real-time PCR. BMC Plant Biol. 2009, 9: 87-PubMedPubMed CentralView ArticleGoogle Scholar
- Henderson IR, Zhang XY, Lu C, Johnson L, Meyers BC, Green PJ, Jacobsen SE: Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat Genet. 2006, 38: 721-725. 10.1038/ng1804.PubMedView ArticleGoogle Scholar
- Kato M, de Lencastre A, Pincus Z, Slack FJ: Dynamic expression of small non-coding RNAs, including novel microRNAs and piRNAs/21U-RNAs, during Caenorhabditis elegans development. Genome Biol. 2009, 10: R54-10.1186/gb-2009-10-5-r54.PubMedPubMed CentralView ArticleGoogle Scholar
- Reddy AM, Zheng Y, Jagadeeswaran G, Macmil SL, Graham WB, Roe BA, Desilva U, Zhang WX, Sunkar R: Cloning, characterization and expression analysis of porcine microRNAs. BMC Genomics. 2008, 10: 65-View ArticleGoogle Scholar
- Li R, Li Y, Kristiansen K, Wang J: SOAP: Short Oligonucleotide Alignment Program. Bioinformatics (Oxford). 2008, 24: 713-714. 10.1093/bioinformatics/btn025.View ArticleGoogle Scholar
- Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34: D140-144. 10.1093/nar/gkj112.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffiths-Jones S: The microRNA registry. Nucleic Acids Res. 2004, 32: D109-111. 10.1093/nar/gkh023.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffiths-Jones S, Saini HK, van Dongen S, Enright AJ: miRBase: tools for microRNA genomics. Nucleic Acids Res. 2008, 36: D154-158. 10.1093/nar/gkn221.PubMedPubMed CentralView ArticleGoogle Scholar
- Girard A, Hannon GJ: Conserved themes in small-RNA-mediated transposon control. Trends Cell Biol. 2008, 18: 136-148. 10.1016/j.tcb.2008.01.004.PubMedPubMed CentralView ArticleGoogle Scholar
- Aravin AA, Lagos-Quintana M, Yalcin A, Zavolan M, Marks D, Snyder B, Gaasterland T, Meyer J, Tuschl T: The small RNA profile during Drosophila melanogaster development. Dev Cell. 2003, 5: 337-350. 10.1016/S1534-5807(03)00228-4.PubMedView ArticleGoogle Scholar
- Kurihara Y, Watanabe Y: Arabidopsis micro-RNA biogenesis through Dicer-like 1 protein functions. Proc Natl Acad Sci USA. 2004, 101: 12753-12758. 10.1073/pnas.0403115101.PubMedPubMed CentralView ArticleGoogle Scholar
- Reinhart BJ, Weinstein EG, Rhoades MW, Bartel B, Bartel DP: MicroRNAs in plants. Genes Dev. 2002, 16: 1616-1626. 10.1101/gad.1004402.PubMedPubMed CentralView ArticleGoogle Scholar
- Meyers BC, Axtell MJ, Bartel B, Bartel DP, Baulcombe D, Bowman JL, Cao X, Carrington JC, Chen XM, Green PJ, Griffiths-Jones S, Jacobsen SE, Mallory AC, Martienssen RA, Poethig RS, Qi Y, Vaucheret H, Voinnet O, Watanabe Y, Weigel D, Zhu JK: Criteria for annotation of plant microRNAs. Plant Cell. 2008, 20: 3186-3190. 10.1105/tpc.108.064311.PubMedPubMed CentralView ArticleGoogle Scholar
- Fahlgren N, Howell MD, Kasschau KD, Chapman EJ, Sullivan CM, Cumbie JS, Givan SA, Law TF, Grant SR, Dangl JL, Carrington JC: High-throughput sequencing of Arabidopsis microRNAs: evidence for frequent birth and death of MIRNA genes. PLoS One. 2007, 2: e219-10.1371/journal.pone.0000219.PubMedPubMed CentralView ArticleGoogle Scholar
- Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK: Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008, 8: 25-10.1186/1471-2229-8-25.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: 3407-3425. 10.1101/gad.1476406.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang G, Guo G, Hu X, Zhang Y, Li Q, Li R, Zhuang R, Lu Z, He Z, Fang X, Chen L, Tian W, Tao Y, Kristiansen K, Zhang X, Li S, Yang H, Wang J, Wang J: Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome. Genome Res. 2010, 20: 646-654. 10.1101/gr.100677.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Ouyang S, Buell CR: The TIGR Plant Repeat Databases: a collective resource for the identification of repetitive sequences in plants. Nucleic Acids Res. 2004, 32: D360-D363. 10.1093/nar/gkh099.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang YJ: miRU: an automated plant miRNA target prediction server. Nucleic Acids Res. 2005, 33: W701-W704. 10.1093/nar/gki383.PubMedPubMed CentralView ArticleGoogle Scholar
- Mallory AC, Vaucheret H: Functions of microRNAs and related small RNAs in plants. Nat Genet. 2006, 38: S31-S36. 10.1038/ng1791.PubMedView ArticleGoogle Scholar
- Lu C, Jeong DH, Kulkarni K, Pillay M, Nobuta K, German R, Thatcher SR, Maher C, Zhang L, Ware D, Liu B, Cao X, Meyers BC, Green PJ: Genome-wide analysis for discovery of rice microRNAs reveals natural antisense microRNAs (nat-miRNAs). Proc Natl Acad Sci USA. 2008, 105: 4951-4956. 10.1073/pnas.0708743105.PubMedPubMed CentralView ArticleGoogle Scholar
- Luo YC, Zhou H, Li Y, Chen JY, Yang JH, Chen YQ, Qu LH: Rice embryogenic calli express a unique set of microRNAs, suggesting regulatory roles of microRNAs in plant post-embryogenic development. FEBS Lett. 2006, 580: 5111-5116. 10.1016/j.febslet.2006.08.046.PubMedView ArticleGoogle Scholar
- German MA, Pillay M, Jeong DH, Hetawal A, Luo SJ, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, German R, De Paoli E, Lu C, Schroth G, Meyers BC, Green PJ: Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008, 26: 941-946. 10.1038/nbt1417.PubMedView ArticleGoogle Scholar
- Addo-Quaye C, Snyder JA, Park YB, Li YF, Sunkar R, Axtell MJ: Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome. RNA. 2009, 15: 2112-2121. 10.1261/rna.1774909.PubMedPubMed CentralView ArticleGoogle Scholar
- Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ: Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008, 18: 758-762. 10.1016/j.cub.2008.04.042.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu L, Zhang QQ, Zhou HY, Ni FR, Wu XY, Qi YJ: Rice microRNA effector complexes and targets. Plant Cell. 2009, 21: 3421-3435. 10.1105/tpc.109.070938.PubMedPubMed CentralView ArticleGoogle Scholar
- Li YF, Zheng Y, Addo-Quaye C, Zhang L, Saini A, Jagadeeswaran G, Axtell MJ, Zhang WX, Sunkar R: Transcriptome-wide identification of microRNA targets in rice. Plant J. 2010, 62: 742-759. 10.1111/j.1365-313X.2010.04187.x.PubMedView ArticleGoogle Scholar
- Yang JH, Li JH, Shao P, Zhou H, Chen YQ, Qu LH: starBase: a database for exploring microRNA-mRNA interaction maps from Argonaute CLIP-Seq and Degradome-Seq data. Nucleic Acids Res. 2010, 39: D202-209.PubMedPubMed CentralView ArticleGoogle Scholar
- Carthew RW, Sontheimer EJ: Origins and mechanisms of miRNAs and siRNAs. Cell. 2009, 136: 642-655. 10.1016/j.cell.2009.01.035.PubMedPubMed CentralView ArticleGoogle Scholar
- Hamilton A, Voinnet O, Chappell L, Baulcombe D: Two classes of short interfering RNA in RNA silencing. EMBO J. 2002, 21: 4671-4679. 10.1093/emboj/cdf464.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen X: Small RNAs and their roles in plant development. Annu Rev Cell Dev Biol. 2009, 25: 21-44. 10.1146/annurev.cellbio.042308.113417.PubMedView ArticleGoogle Scholar
- Zilberman D, Cao X, Jacobsen SE: ARGONAUTE4 control of locus-specific siRNA accumulation and DNA and histone methylation. Science. 2003, 299: 716-719. 10.1126/science.1079695.PubMedView ArticleGoogle Scholar
- Chan SW, Zilberman D, Xie Z, Johansen LK, Carrington JC, Jacobsen SE: RNA silencing genes control de novo DNA methylation. Science. 2004, 303: 1336-10.1126/science.1095989.PubMedView ArticleGoogle Scholar
- Wu L, Zhou HY, Zhang QQ, Zhang JG, Ni FR, Liu C, Qi YJ: DNA methylation mediated by a microRNA pathway. Mol Cell. 2010, 38: 465-475. 10.1016/j.molcel.2010.03.008.PubMedView ArticleGoogle Scholar
- Wang X, Elling AA, Li X, Li N, Peng Z, He G, Sun H, Qi Y, Liu XS, Deng XW: Genome-wide and organ-specific landscapes of epigenetic modifications and their relationships to mRNA and small RNA transcriptomes in maize. Plant Cell. 2009, 21: 1053-1069. 10.1105/tpc.109.065714.PubMedPubMed CentralView ArticleGoogle Scholar
- Valoczi A, Varallyay E, Kauppinen S, Burgyan J, Havelda Z: Spatio-temporal accumulation of microRNAs is highly coordinated in developing plant tissues. Plant J. 2006, 47: 140-151. 10.1111/j.1365-313X.2006.02766.x.PubMedView ArticleGoogle Scholar
- Sieber P, Wellmer F, Gheyselinck J, Riechmann JL, Meyerowitz EM: Redundancy and specialization among plant microRNAs: role of the MIR164 family in developmental robustness. Development. 2007, 134: 1051-1060. 10.1242/dev.02817.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Cannon CH, Cobb GP, Anderson TA: Conservation and divergence of plant microRNA genes. Plant J. 2006, 46: 243-259. 10.1111/j.1365-313X.2006.02697.x.PubMedView ArticleGoogle Scholar
- Zhang BH, Pan XP, Anderson TA: Identification of 188 conserved maize microRNAs and their targets. FEBS Lett. 2006, 580: 3753-3762. 10.1016/j.febslet.2006.05.063.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: 571-584. 10.1101/gr.6897308.PubMedPubMed CentralView ArticleGoogle Scholar
- Graham LK, Wilcox LW: The origin of alternation of generations in land plants: a focus on matrotrophy and hexose transport. Philos Trans R Soc Lond B Biol Sci. 2000, 355: 757-766. 10.1098/rstb.2000.0614. discussion 766-757PubMedPubMed CentralView ArticleGoogle Scholar
- Tang F: Small RNAs in mammalian germline: Tiny for immortal. Differentiation. 2010, 79: 141-146. 10.1016/j.diff.2009.11.002.PubMedView ArticleGoogle Scholar
- Honys D, Twell D: Transcriptome analysis of haploid male gametophyte development in Arabidopsis. Genome Biol. 2004, 5: R85-10.1186/gb-2004-5-11-r85.PubMedPubMed CentralView ArticleGoogle Scholar
- Achard P, Herr A, Baulcombe DC, Harberd NP: Modulation of floral development by a gibberellin-regulated microRNA. Development. 2004, 131: 3357-3365. 10.1242/dev.01206.PubMedView ArticleGoogle Scholar
- Millar AA, Gubler F: The Arabidopsis GAMYB-like genes, MYB33 and MYB65, are microRNA-regulated genes that redundantly facilitate anther development. Plant Cell. 2005, 17: 705-721. 10.1105/tpc.104.027920.PubMedPubMed CentralView ArticleGoogle Scholar
- Wu MF, Tian Q, Reed JW: Arabidopsis microRNA167 controls patterns of ARF6 and ARF8 expression, and regulates both female and male reproduction. Development. 2006, 133: 4211-4218. 10.1242/dev.02602.PubMedView 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: 1376-1386. 10.1105/tpc.105.030841.PubMedPubMed CentralView 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. 10.1016/j.molcel.2004.05.027.PubMedView ArticleGoogle Scholar
- Cecchetti V, Altamura MM, Falasca G, Costantino P, Cardarelli M: Auxin regulates Arabidopsis anther dehiscence, pollen maturation, and filament elongation. Plant Cell. 2008, 20: 1760-1774. 10.1105/tpc.107.057570.PubMedPubMed CentralView ArticleGoogle Scholar
- de Martinis D, Cotti G, Hekker ST, Harren FJM, Mariani C: Ethylene response to pollen tube growth in Nicotiana tabacum flowers. Planta. 2002, 214: 806-812. 10.1007/s00425-001-0684-2.PubMedView ArticleGoogle Scholar
- Huang S, Cerny RE, Qi YL, Bhat D, Aydt CM, Hanson DD, Malloy KP, Ness LA: Transgenic studies on the involvement of cytokinin and gibberellin in male development. Plant Physiol. 2003, 131: 1270-1282. 10.1104/pp.102.018598.PubMedPubMed CentralView ArticleGoogle Scholar
- Park JH, Halitschke R, Kim HB, Baldwin IT, Feldmann KA, Feyereisen R: A knock-out mutation in allene oxide synthase results in male sterility and defective wound signal transduction in Arabidopsis due to a block in jasmonic acid biosynthesis. Plant J. 2002, 31: 1-12. 10.1046/j.1365-313X.2002.01328.x.PubMedView ArticleGoogle Scholar
- Cheng H, Song SS, Xiao LT, Soo HM, Cheng ZW, Xie DX, Peng JR: Gibberellin acts through jasmonate to control the expression of MYB21, MYB24, and MYB57 to promote stamen filament growth in Arabidopsis. PLoS Genet. 2009, 5: e1000440-10.1371/journal.pgen.1000440.PubMedPubMed CentralView ArticleGoogle Scholar
- Baker SM, Bronner CE, Zhang L, Plug AW, Robatzek M, Warren G, Elliott EA, Yu JA, Ashley T, Arnheim N, Flavell RA, Liskay RM: Male-mice defective in the DNA mismatch repair gene Pms2 exhibit abnormal chromosome synapsis in meiosis. Cell. 1995, 82: 309-319. 10.1016/0092-8674(95)90318-6.PubMedView ArticleGoogle Scholar
- Komori T, Ohta S, Murai N, Takakura Y, Kuraya Y, Suzuki S, Hiei Y, Imaseki H, Nitta N: Map-based cloning of a fertility restorer gene, Rf-1, in rice (Oryza sativa L.). Plant J. 2004, 37: 315-325.PubMedView ArticleGoogle Scholar
- Saze H, Kakutani T: Heritable epigenetic mutation of a transposon-flanked Arabidopsis gene due to lack of the chromatin-remodeling factor DDM1. EMBO J. 2007, 26: 3641-3652. 10.1038/sj.emboj.7601788.PubMedPubMed CentralView ArticleGoogle Scholar
- O'Toole N, Hattori M, Andres C, Iida K, Lurin C, Schmitz-Linneweber C, Sugita M, Small I: On the expansion of the pentatricopeptide repeat gene family in plants. Mol Biol Evol. 2008, 25: 1120-1128. 10.1093/molbev/msn057.PubMedView ArticleGoogle Scholar
- Le Trionnaire G, Twell D: Small RNAs in angiosperm gametophytes: from epigenetics to gamete development. Genes Dev. 2010, 24: 1081-1085. 10.1101/gad.1936110.PubMedPubMed CentralView ArticleGoogle Scholar
- Cartagena JA, Matsunaga S, Seki M, Kurihara D, Yokoyama M, Shinozaki K, Fujimoto S, Azumi Y, Uchiyama S, Fukui K: The Arabidopsis SDG4 contributes to the regulation of pollen tube growth by methylation of histone H3 lysines 4 and 36 in mature pollen. Dev Biol. 2008, 315: 355-368. 10.1016/j.ydbio.2007.12.016.PubMedView ArticleGoogle Scholar
- Ito M, Koike A, Koizumi N, Sano H: Methylated DNA-binding proteins from arabidopsis. Plant Physiol. 2003, 133: 1747-1754. 10.1104/pp.103.026708.PubMedPubMed CentralView ArticleGoogle Scholar
- Zemach A, McDaniel IE, Silva P, Zilberman D: Genome-wide evolutionary analysis of eukaryotic DNA methylation. Science. 2010, 328: 916-919. 10.1126/science.1186366.PubMedView ArticleGoogle Scholar
- Lee TF, Zhai JX, Meyers BC: Conservation and divergence in eukaryotic DNA methylation. Proc Natl Acad Sci USA. 2010, 107: 9027-9028. 10.1073/pnas.1005440107.PubMedPubMed CentralView ArticleGoogle Scholar
- Zemach A, Grafi G: Characterization of Arabidopsis thaliana methyl-CpG-binding domain (MBD) proteins. Plant J. 2003, 34: 565-572. 10.1046/j.1365-313X.2003.01756.x.PubMedView ArticleGoogle Scholar
- Rousseaux S, Caron C, Govin K, Lestrat C, Faure AK, Khochbin S: Establishment of male-specific epigenetic information. Gene. 2005, 345: 139-153. 10.1016/j.gene.2004.12.004.PubMedView ArticleGoogle Scholar
- Trasler JM: Epigenetics in spermatogenesis. Mol Cell Endocrinol. 2009, 306: 33-36. 10.1016/j.mce.2008.12.018.PubMedView ArticleGoogle Scholar
- Liu SP, Li D, Li QB, Zhao P, Xiang ZH, Xia QY: MicroRNAs of Bombyx mori identified by Solexa sequencing. BMC Genomics. 2010, 11: 148-10.1186/1471-2164-11-148.PubMedPubMed CentralView ArticleGoogle Scholar
- ArrayExpress. [http://www.ebi.ac.uk/arrayexpress/]
- Qi X, Bao FS, Xie Z: Small RNA deep sequencing reveals role for Arabidopsis thaliana RNA-dependent RNA polymerases in viral siRNA biogenesis. PLoS One. 2009, 4: e4971-10.1371/journal.pone.0004971.PubMedPubMed CentralView ArticleGoogle Scholar
- Gardner PP, Daub J, Tate JG, Nawrocki EP, Kolbe DL, Lindgreen S, Wilkinson AC, Finn RD, Griffiths-Jones S, Eddy SR, Bateman A: Rfam: updates to the RNA families database. Nucleic Acids Res. 2009, 37: D136-D140. 10.1093/nar/gkn766.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffiths-Jones S, Bateman A, Marshall M, Khanna A, Eddy SR: Rfam: an RNA family database. Nucleic Acids Res. 2003, 31: 439-441. 10.1093/nar/gkg006.PubMedPubMed CentralView ArticleGoogle Scholar
- Griffiths-Jones S, Moxon S, Marshall M, Khanna A, Sr E, Bateman A: Rfam: annotating non-coding RNAs in complete genomes. Nucleic Acids Res. 2005, 33: D121-D124.PubMedPubMed CentralView ArticleGoogle Scholar
- GenBank. [http://www.ncbi.nlm.nih.gov/genbank/]
- Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005, 33: e179-10.1093/nar/gni178.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Q, Huang Z, Xue HL, Jin CC, Ju XL, Han JDJ, Chen YG: MicroRNA miR-24 inhibits erythropoiesis by targeting activin type I receptor ALK4. Blood. 2008, 111: 588-595. 10.1182/blood-2007-05-092718.PubMedView ArticleGoogle Scholar
- Bu Y, Lu CH, Bian CJ, Wang JH, Li J, Zhang B, Li ZY, Brewer G, Zhao RC: Knockdown of Dicer in MCF-7 human breast carcinoma cells results in G1 arrest and increased sensitivity to cisplatin. Oncol Rep. 2009, 21: 13-17.PubMedGoogle Scholar
- Varkonyi-Gasic E, Wu RM, Wood M, Walton EF, Hellens RP: Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007, 3: 12-10.1186/1746-4811-3-12.PubMedPubMed CentralView ArticleGoogle Scholar
- Lacombe S, Nagasaki H, Santi C, Duval D, Piegu B, Bangratz M, Breitler J-C, Guiderdoni E, Brugidou C, Hirsch J, Cao X, Brice C, Panaud O, Karlowski WM, Sato Y, Echeverria M: Identification of precursor transcripts for 6 novel miRNAs expands the diversity on the genomic organisation and expression of miRNA genes in rice. BMC Plant Biol. 2008, 8: 123-10.1186/1471-2229-8-123.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.