Analysis of small RNAs revealed differential expressions during the development of pollen and embryo sac and their targets may associate with fertility in autotetraploid rice

Background: Partial pollen and embryo sac sterilities are the two main reasons for low fertility in autotetraploid rice. Our previous study revealed that small RNAs changes may associate with pollen fertility in autotetraploid rice. However, knowledge on comparative analysis between the development of pollen and embryo sac by small RNAs in autotetraploid rice is still unknown. In the present study, WE-CLSM (whole-mount eosin B-staining confocal laser scanning microscopy) and high-throughput sequencing technology was employed to examine the cytological variations and to analyze small RNAs changes during pollen and embryo sac development in autotetraploid rice compared with its diploid counterpart. Results: A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during development of pollen and embryo sac in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of miRNAs-enriched during the development of pollen and embryo sac in autotetraploid rice revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development. However, only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression levels in pollen and embryo sac development. Of these DEM, osa-miR1436_L+3_1ss5CT and osa-miR167h-3p were associated with the female meiosis, while osa-miR159a.1 and osa-MIR159a-p5 were related with the male meiosis. 21nt-phasiRNAs were detected both during pollen and embryo sac development, while 24nt-phasiRNAs were found only in pollen development, which displayed down-regulation in autotetraploid compared to diploid rice and their spatial-temporal expression patterns were similar to osa-miR2275d. 24nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development. Conclusion: The above results not only provide the small RNAs changes during four landmark stages of pollen and embryo sac development in autotetraploid rice but also have identified specifically expressed miRNAs, especially meiosis-related miRNAs, pollen-24nt-phasiRNAs and TEs-siRNAs in autotetraploid rice. Together, these findings provide a foundation for understanding the effect of polyploidy on small RNAs expression patterns during pollen and embryo sac development that may lead to different abnormalities in autotetraploid rice. Overall design: To investigate the small RNAs associated with the pollen and embryo sac development of autotetraploid and diploid (CK) rice, anthers with four pollen development stages, including pre-meiotic interphase (PMA), meiosis (MA), single microspore stage (SCP) and bi-cellular pollen stage (BCP), and the ovary tissues were collected with their corresponding embryo sac development stages (i.e. megasporocyte formation stage (MF), megasporocyte meiosis stage (MM), functional megaspore formation stage (FMF) and eight-nucleate embryo sac developing-stage (EES)), were selected to construct 16 libraries and sequenced by Solexa high-throughput sequencing technology.