MicroRNAs (miRNAs) play key roles in plant reproduction. However, knowledge on microRNAome analysis in autotetraploid rice is rather limited. Here, high-throughput sequencing technology was employed to analyze miRNAomes during pollen development in diploid and polyploid rice. A total of 172 differentially expressed miRNAs (DEM) were detected in autotetraploid rice compared to its diploid counterpart, and 57 miRNAs were specifically expressed in autotetraploid rice. Of the 172 DEM, 115 and 61 miRNAs exhibited up- and down-regulation, respectively. Gene Ontology analysis on the targets of up-regulated DEM showed that they were enriched in transport and membrane in pre-meiotic interphase, reproduction in meiosis, and nucleotide binding in single microspore stage. osa-miR5788 and osa-miR1432-5p_R+1 were up-regulated in meiosis and their targets revealed interaction with the meiosis-related genes, suggesting that they may involve in the genes regulation associated with the chromosome behavior. Abundant 24 nt siRNAs associated with transposable elements were found in autotetraploid rice during pollen development; however, they significantly declined in diploid rice, suggesting that 24 nt siRNAs may play a role in pollen development. These findings provide a foundation for understanding the effect of polyploidy on small RNA expression patterns during pollen development that cause pollen sterility in autotetraploid rice.

Polyploid rice hybrids have a powerful biological and yield potential that may become a new way for rice breeding; however, low fertility is major hindrance in commercial utilization. Here, we developed a neo-tetraploid rice that could overcome the sterility of autotetraploid rice and produce high heterosis. Transcriptome analysis of F₁ hybrid developed by crossing neo-tetraploid with autotetraploid rice displayed 807, 663 and 866 differentially expressed genes that uniquely associated with F₁ and specific to (DEGFu-sp) anther, ovary and leaf, respectively. Of the DEGFu-sp, 1224 genes displayed nonadditive expression; 44 and 10 genes were annotated as TFs and methyltransferase or hydroxymethyltransferase, respectively. Gene ontology enrichment and co-expression analysis revealed specific differential gene expressions in the DEGFu-sp to leaf, anther and ovary, such as genes related to photosynthesis, metabolic process and transport, and co-expression network including fertility, resistance and epigenetic elements. Of the DEGFu-sp to anther, 42 meiosis stage-specific genes, eight meiosis-related genes, such as RAD51 and SMC2, were identified. We identified 38 miRNAs from DEGFu-sp to anther, and their targets were associated with pollen fertility and retrotransposon protein. Our study provides new germplasm for polyploid rice breeding, and revealed complex regulatory mechanisms that might be associated with heterosis and fertility.

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. A total of 321 and 368 differentially expressed miRNAs (DEM) were detected during pollen and embryo sac development in autotetraploid rice, respectively. Gene Ontology enrichment analysis on the targets of DEM associated with embryo sac and pollen development revealed 30 prominent functional gene classes, such as cell differentiation and signal transduction during embryo sac development, while only 7 prominent functional gene classes, such as flower development and transcription factor activity, were detected during pollen development in autotetraploid rice. The expression levels of 39 DEM, which revealed interaction with meiosis-related genes, showed opposite expression patterns during 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. 21 nt-phasiRNAs were detected during both pollen and embryo sac development, while 24 nt-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. 24 nt TEs-siRNAs were found to be up-regulated in embryo sac but down-regulated in pollen development. 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-specific-24 nt-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.

Background: Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA. Results: A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice. Conclusion: The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Common wild rice (Oryza rufipogon Griff.) is an important germplasm for rice breeding, which contains many resistance genes. Re-sequencing provides an unprecedented opportunity to explore the abundant useful genes at whole genome level. Here, we identified the nucleotide-binding site leucine-rich repeat (NBS-LRR) encoding genes by re-sequencing of two wild rice lines (i.e. Huaye 1 and Huaye 2) that were developed from common wild rice. We obtained 128 to 147 million reads with approximately 32.5-fold coverage depth, and uniquely covered more than 89.6% (> = 1 fold) of reference genomes. Two wild rice lines showed high SNP (single-nucleotide polymorphisms) variation rate in 12 chromosomes against the reference genomes of Nipponbare (japonica cultivar) and 93–11 (indica cultivar). InDels (insertion/deletion polymorphisms) count-length distribution exhibited normal distribution in the two lines, and most of the InDels were ranged from -5 to 5 bp. With reference to the Nipponbare genome sequence, we detected a total of 1,209,308 SNPs, 161,117 InDels and 4,192 SVs (structural variations) in Huaye 1, and 1,387,959 SNPs, 180,226 InDels and 5,305 SVs in Huaye 2. A total of 44.9% and 46.9% genes exhibited sequence variations in two wild rice lines compared to the Nipponbare and 93–11 reference genomes, respectively. Analysis of NBS-LRR mutant candidate genes showed that they were mainly distributed on chromosome 11, and NBS domain was more conserved than LRR domain in both wild rice lines. NBS genes depicted higher levels of genetic diversity in Huaye 1 than that found in Huaye 2. Furthermore, protein-protein interaction analysis showed that NBS genes mostly interacted with the cytochrome C protein (Os05g0420600, Os01g0885000 and BGIOSGA038922), while some NBS genes interacted with heat shock protein, DNA-binding activity, Phosphoinositide 3-kinase and a coiled coil region. We explored abundant NBS-LRR encoding genes in two common wild rice lines through genome wide re-sequencing, which proved to be a useful tool to exploit elite NBS-LRR genes in wild rice. The data here provide a foundation for future work aimed at dissecting the genetic basis of disease resistance in rice, and the two wild rice lines will be useful germplasm for the molecular improvement of cultivated rice.