- Research article
- Open Access
Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice
BMC Plant Biology volume 20, Article number: 82 (2020)
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.
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.
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.
Whole genome doubling (WGD), or polyploidy, has been known as a vital evolutionary force for genetic diversity in plants. About 35%~ 70% angiosperms have undergone at least one round of genome duplication, such as wheat, cotton and brassica . A variety of genome alterations caused by polyploidy, including creating new source of redundant genes, genome rearrangements and genome-wide subfunctionalization and neofunctionalization [2, 3]. The potential advantages of plant genome reconstruction can be utilized with the newly developed resources of genomics for crop improvement .
Autotetraploid rice, as a useful germplasm, has more biological advantages than the diploid rice that may become a new way for rice breeding [5,6,7]. However, low fertility, an unfavorable trait, is one of the major barriers in commercial production of polyploid rice [8,9,10,11]. Meiosis is an essential biological process in the life cycle of sexual reproduction in plants . Normal chromosomes segregation in meiosis is one of the most important challenges to polyploids [13, 14]. Multivalents during meiosis often generated massively in the newly formed polyploids [15,16,17,18]. Cytogenetic analysis revealed that chromosomal abnormalities, abnormal development processes of pollens and embryo sacs were the main reasons for low fertility in autotetraploid rice [16, 19, 20]. Transcriptome analysis revealed that polyploidy enhances the F1 pollen sterility multi-loci interactions which lead to meiosis abnormalities and pollen sterility in autotetraploid rice hybrids . In addition, some meiosis-related microRNAs were identified that may interrupt the chromosome behavior in pollen mother cells (PMC) of autotetraploid rice . Recently, differential expression patterns of small RNAs during embryo sac development in autotetraploid rice have been reported to associate with sterility, especially during megasporocyte .
Non-coding RNAs (ncRNAs) cannot encode proteins; however, they play a key role in regulation of transcriptional, posttranscriptional and epigenetic mechanisms. Long non-coding RNAs (lncRNAs) are one of the important types of ncRNAs, with a length of more than 200 nt, and usually exhibited tissue/developmental-stage-specific expression patterns . LncRNAs have been identified in many plants, including rice, maize, sunflower and strawberry [23,24,25], and some of them were reported to associate with reproductive processes in rice . Additionally, the lncRNAs of rice anthers before flowering, pistils before flowering, spikelets 5 days after pollination and shoots 14 days after germination have been characterized, and several important lncRNAs were exclusively identified in young panicles, which showed a profound impact on rice fertility . Recently, overexpression of LAIR increased grain yield and up-regulated the expression of several leucine-rich repeat receptor kinase genes in rice . Meanwhile, lncRNAs played a pivotal role to small RNAs population, which could act as small RNAs precursors, miRNAs targets and decoys .
Although information on small RNAs and transcriptome sequencing in autotetraploid rice are available [11, 20], the expression patterns of lncRNAs associated with meiosis in autotetraploid rice are still poorly understood. The present study was planned to investigate the expression levels of lncRNAs in anther and ovary during meiosis in an autotetraploid rice, Taichung 65-4x, and its diploid counterpart (Taichung 65) by high-throughput sequencing. We sought to distinguish the lncRNAs changes between autotetraploid and diploid rice, and to determine the relationships between gene expression profiles and spikelet sterility.
Identification and characterization of lncRNAs in rice anther and ovary
Previous study showed that the pollen fertility of Taichung 65-4x was lower than the diploid counterpart . Here, embryo sac fertility of Taichung 65-4x was observed by Whole-mount eosin B-staining confocal laser scanning microscopy (Fig. 1), which was 62.15%. Therefore, anthers and corresponding ovaries of Taichung 65 and Taichung 65-4x from pre-meiotic interphase to meiosis prophase I (early meiosis stage) were collected for transcriptome sequencing. After a genome-wide systematic filtration, a total of 4859 confident lncRNAs (FPKM > 0.5 in at least one library) were obtained, including 4018 intergenic lncRNAs, 790 antisense lncRNAs, and 51 sense lncRNAs, and used for further analysis (Additional file 1: Fig. S1A; Additional file 2: Table S1). Venn analysis revealed that most of the lncRNAs were specifically expressed in anthers/ovaries, and a total of 1282 and 2025 lncRNAs were specifically detected in the rice reproductive tissues of Taichung 65 and Taichung 65-4x, respectively (Additional file 1: Fig. S1B). Moreover, 28% lncRNAs were exclusively expressed in the anthers of Taichung 65-4x.
We further characterized the features of lncRNAs, and found that the expression levels of lncRNAs were lower than the protein-coding genes (Additional file 1: Fig. S2A). High proportions of lncRNAs were found with one exon, which differed from the protein-coding genes (Additional file 1: Fig. S2B). The transcript length distribution was not similar between lncRNAs and protein-coding genes. Most of protein-coding genes were longer than 1 kb, whereas majority of lncRNAs were less than 1 kb. The lengths of the lncRNAs, with a mean value of 938 nucleotides, was shorter than the protein-coding genes (mean = 1862 bp) (Additional file 1: Fig. S2C). In addition, shorter ORFs (open reading frames) were observed in lncRNAs (mean = 66.36 aa) than protein-coding genes (mean = 384.27 aa) in the rice reproductive tissues (Additional file 1: Fig. S2D). Moreover, only a small portion (4.4%) of rice reproductive lncRNAs was found to be conserved, including 12 in Arabidopsis, 80 in Brachypodium, 103 in Zea mays, and 161 in Sorghum bicolor. Intriguingly, 11 of 12 conserved lncRNAs in Arabidopsis were also found to be conserved in Brachypodium, Zea mays and Sorghum bicolor (Additional file 2: Table S2).
The full-length of ten lncRNAs were obtained by rapid amplification of cDNA ends (RACE) assay (Fig. 2), including TCONS_00057811, TCONS_00130461, TCONS_00068868, TCONS_00091337, TCONS_00055980, TCONS_00130471, TCONS_00111916, TCONS_00045430, TCONS_00013598 and TCONS_00121908. These lncRNAs could be amplified and existed in the rice which indicated that lncRNAs were not assembly artifacts; however, no DNA variations were found between Taichung 65 and Taichung 65-4x. TCONS_00057811 was located in intergenic region (25685721–25,686,497, − strand) on chromosome 12. The real length of TCONS_00057811 was 777 bp, and it had two exons, which was different from the gene body structure predicted by transcriptome analysis. The length of TCONS_00130461 was 2073 bp, which was also different from the transcriptome result. Not only TCONS_00057811 and TCONS_00130461, the length of the other eight lncRNAs were different from the predicted fragments (Additional file 1: Fig. S3). Furthermore, phylogenetic tree (MEGA-X) results showed that six lncRNAs were conserved among our rice dataset (re-sequencing data of more than 100 rice lines, including neo-tetraploid, autotetraploid, typical japonica/indica and wild rice lines, and unpublished data), including TCONS_00057811, TCONS_00091337, TCONS_00130461, TCONS_00130471, TCONS_00045430 and TCONS_00055980 (Additional file 1: Fig. S4).
Differentially expressed lncRNAs in anther and ovary during meiosis
To characterize the differential expression patterns of lncRNAs, we first figured out the lncRNAs associated with anther and ovary. By cluster analysis, 78 anther-preferred and 222 ovary-preferred lncRNAs, which highly/specifically expressed in anther/ovary, were identified in Taichung 65 (defined as anther- and ovary-preferred-2x), and 222 anther-preferred and 394 ovary-preferred lncRNAs were detected in Taichung 65-4x (defined as anther- and ovary-preferred-4x), respectively (Additional file 1: Fig. S5; Additional file 2: Table S1). Of the preferred lncRNAs in ovary, 133 lncRNAs were detected in both Taichung 65 and Taichung 65-4x during meiosis (co-ovary-preferred), and 34 lncRNAs were shared by the Taichung 65 and Taichung 65-4x in anther during meiosis (co-anther-preferred).
A total of 431 differentially expressed lncRNAs (DEL) were detected in anthers of Taichung 65-4x compared to the Taichung 65 (DEL-anther), of these 221 and 210 lncRNAs showed up- and down-regulation, respectively (Additional file 2: Table S3). Besides, a small portion of DEL (32 in total) was also found in ovaries of Taichung 65-4x compared to Taichung 65 (DEL-ovary), and 8 and 24 lncRNAs exhibited up- and down-regulation, respectively. 19 common DEL were detected between DEL-anther and DEL-ovary by venn analysis, and 17 out of these showed similar differential expression patterns with 3 co-up-regulated and 14 co-down-regulated DEL in autotetraploid rice (Additional file 1: Fig. S6). Furthermore, we compared the 431 DEL-anther and 32 DEL-ovary with the anther-preferred lncRNAs and ovary-preferred lncRNAs by venn analysis, respectively. 27 out of 34 co-anther-preferred lncRNAs showed differential expression patterns in Taichung 65-4x (DEL-anther: type 1). Surprisingly, these 27 DEL showed down-regulation in Taichung 65-4x compared to Taichung 65 (Additional file 1: Fig. S7A). In addition, 25 (DEL-anther: type 2) and 40 DEL (DEL-anther: type 3) were overlapped with the anther-preferred-2x and anther-preferred-4x lncRNAs, respectively. Similarly, 16 DEL were found to be associated with co-ovary-preferred lncRNAs (DEL-ovary: type 1). Interestingly, 15 out of 16 DEL were found to be down-regulated in Taichung 65-4x compared to Taichung 65, except TCONS_00019924 (Additional file 1: Fig. S7B). Six (DEL-ovary: type 2) and three DEL (DEL-ovary: type 3) were found to be overlapped with the ovary-preferred-2x and ovary-preferred-4x lncRNAs, respectively. The rest of the DEL in anther/ovary was defined as type 4.
In addition, we used BLAST search to compare the sequences and 160 lncRNAs displayed more than 90% sequence similarity with Zhang et al. . Of these, 56 showed differential expression patterns between Taichung 65-4x and Taichung 65 (Additional file 2: Table S4). Among the 56 DEL, 18 showed differential expressions in Taichung 65-4x anthers, which highly expressed in anthers , such as TCONS_00130461. Two DEL were found to be differentially expressed in Taichung 65-4x ovaries that showed pistils expression patterns in Zhang’s results, such as TCONS_00054234.
Small RNA populations and their relationships with differentially expressed lncRNAs in autotetraploid rice during meiosis
By small RNAs sequencing, 1034 miRNAs (i.e. 357 known miRNAs and 677 novel miRNAs) and 17,609 phasiRNAs (i.e. 15,992 were 21 nt in length and 1617 were 24 nt) were detected in the rice reproductive tissues (Additional file 2: Table S5 and S6). According to the model described by Boerner and McGinnis , these types of small RNAs were aligned to the 4859 confident lncRNAs to detect the small RNAs precursor in lncRNAs.
By aligning miRNA precursors to lncRNAs, 27 lncRNAs were predicted to be the precursors of 43 miRNAs (Additional file 2: Table S7). Nine of these showed differential expression patterns in autotetraploid rice anther, and one DEL-anther (Type 1) was corresponding to the differentially expressed miRNA (DEM) (i.e. TCONS_00160902 was predicted as a precursor of zma-miR2275b-5p_1ss12CG). A total of 98 lncRNAs were predicted to be the precursors of 949 phasiRNAs, including 63 lncRNAs as precursors to 433 21 nt-phasiRNAs (2.70%) and 35 lncRNAs as precursors to 516 24 nt-phasiRNAs (31.91%). Of these 98 lncRNAs, 29 showed differential expression patterns (down-regulation) in autotetraploid rice anther compared to diploid, but no one found to be differentially expressed in ovary of autotetraploid rice (Additional file 2: Table S8). Moreover, DEL-anther (Type 1) could generate most of phasiRNAs (270 phasiRNAs/13 DEL-anther) (Table 1). Interestingly, majority of the DEL (22 of 29 DEL) were found to be the precursors of 24 nt-phasiRNAs. In addition, nine down-regulated DEL were associated with 11 down-regulated phasiRNAs in Taichung 65-4x. Taken together, 39 DEL were associated with small RNAs, which could serve as precursors and showed down-regulation in autotetraploid rice (Additional file 2: Table S3).
Surprisingly, we could not identify any DEL-anther/ovary that predicted to be the target of a DEM in autotetraploid rice during meiosis by psRNATarget  (Additional file 2: Table S9). Moreover, there was no DEL-anther or DEL-ovary that could be the target mimic of the DEM in Taichung 65-4x, which decoy miRNAs and blocked the interaction between miRNAs and their authentic targets.
Repetitive elements were associated with differentially expressed lncRNAs in autotetraploid rice during meiosis
A series of transposable elements (TEs) genes in anther and ovary were observed, and the relative expression levels of TEs were more abundant in the anthers of Taichung 65 than Taichung 65-4x, especially DNA transposons (i.e. Class II TEs, such as the Helitron family) (Fig. 3A). However, no remarkable changes in TEs were observed in ovary between diploid and autotetraploid rice. Majority of the lncRNAs (73.22%) were overlapped with repeat sequences. Of these TEs-associated lncRNAs, 66.36% carried Class II and 19.08% carried retrotransposons (Class I), while the remaining 14.56% contained both of them. Helitron family was the most frequent sub-class in Class II, whereas the Gypsy family constituted a major portion of Class I (Additional file 1: Fig. S8). The relative expression levels of TEs-lncRNAs also showed similar distribution of TEs genes between Taichung 65-4x and Taichung 65 (Fig. 3B). The differential expression patterns of some TEs and TEs-lncRNAs in anther were confirmed by qRT-PCR (Additional file 1: Fig. S9).
The TEs-associated lncRNAs showed differential expressions in autotetraploid rice anther and ovary, with a percentage of 73.78% (318/431) and 81.25% (26/32) in DEL, respectively (Additional file 2: Table S10). Of the TEs-associated DEL (TEs-DEL), down-regulated TEs-DEL-anther and TEs-DEL-ovary were mostly found to be associated with TEs (Table 2). Down-regulated TEs-DEL-anther were abundantly found in Class II, except EnSpm family, whereas up-regulated TEs-DEL-anther have a great relationship with Gypsy families of Class I in Taichung 65-4x anther. In addition, Helitron and Stowaway families were more related to TEs-DEL-anther (Type 1) compared to other types, and Helitron family was also detected in TEs-DEL-ovary (Type 1) (Additional file 3: Table S11).
Target prediction of differentially expressed lncRNAs and classification of pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis-related targets
To identify the potential targets of lncRNAs could help us to uncover the functions of lncRNAs. In total, 942 potential cis-regulated target genes and 3769 potential trans-regulated target genes were detected, and their target genes also showed differential expression patterns (hereafter referred as DEG) (Additional file 3: Table S12 and S13). Among the potential targets of the DEL-anther, 25 Gene Ontology (GO) terms were significantly categorized, including flower development, DNA binding and nucleus (Additional file 3: Table S14). Beside, only eight GO terms were significantly enriched among the DEL-ovary, such as response to abiotic stimulus.
Furthermore, we evaluated relationship between the targets predicted by DEM and DEL. A total of 104 and 108 predicted targets were identified from the DEM of anther (DEM-anther) and ovary (DEM-ovary) in Taichung 65-4x compared to the diploid rice, respectively (Additional file 3: Table S15). 76 and 81 targets were specifically predicted by DEM-anther and DEM-ovary, whereas 3739 and 328 targets were specifically predicted by DEL-anther and DEL-ovary, respectively (Additional file 1: Fig. S10). GO analysis showed that plastid was specifically enriched in the targets of DEM-anther, while anatomical structure morphogenesis, regulation of cell size and cell differentiation were specifically detected in the targets of DEM-ovary. Moreover, nucleus, response to endogenous stimulus and nucleic acid binding were specifically found in the targets of DEL-anther, and response to abiotic stimulus, secondary metabolic process, and response to extracellular stimulus were specifically identified in the targets of DEL-ovary (Additional file 3: Table S16).
Moreover, we compared the targets with the previous studies to identify the most important meiosis-related targets [32,33,34,35]. By using the predicted targets of DEL-anther, 237 DEL-anther were found to be associated with 110 meiosis-related genes (DEL-PMC) (Additional file 3: Table S17). Of these targets, six rice meiosis genes, including LOC_Os01g39630  (RAD51C, annotated as DNA repair protein Rad51) predicted by TCONS_00049558, LOC_Os03g58600  (MEL1, encoded PAZ domain containing protein) predicted by TCONS_00068868, and LOC_Os01g68870  (OsMSP1, annotated as leucine-rich repeat receptor protein kinase EXS precursor), LOC_Os02g57270  (OsTBP1, annotated as MYB family transcription factor), LOC_Os08g16610  (OsRad21–3, annotated as Rad21/Rec8 like protein), and LOC_Os12g31370  (RAD51A2, annotated as DNA repair protein Rad51) predicted by TCONS_00131142, related to the male meiocytes were identified and showed down-regulation in the present study. Thirteen putative meiosis-related genes displayed down-regulation in autotetraploid rice anther, such as LOC_Os03g12730  (the homologous Arabidopsis gene is BAM1, annotated as receptor protein kinase CLAVATA1 precursor) predicted by TCONS_00068868, LOC_Os04g30790  (the homologous Arabidopsis gene is XRI1, annotated as expressed protein) and LOC_Os11g12620  (the homologous Arabidopsis gene is BAM1, annotated as receptor protein kinase CLAVATA1 precursor) predicted by TCONS_00055980, LOC_Os12g34510  (the homologous Arabidopsis gene is H2AX, encoded Core histone H2A/H2B/H3/H4 domain containing protein) predicted by TCONS_00068868 and TCONS_00121884, and LOC_Os12g41350  (the homologous Arabidopsis gene is ASY1, annotated as meiotic asynaptic mutant 1) predicted by TCONS_00057811. In addition, 13 targets showed meiosis stage-specific gene expressions [32, 33], and 81 targets exhibited differential expression patterns in autotetraploid compared to the diploid rice during meiosis . Of these 81 targets, 28 targets exhibited regulation patterns similar to Wu et al. , such as LOC_Os03g11540 (RPA1B - Putative single-stranded DNA binding complex subunit 1) and LOC_Os05g14590 (MCM6 - Putative minichromosome maintenance MCM complex subunit 6) regulated by TCONS_00055980 and TCONS_00130461, respectively.
Of the targets of DEL-ovary, 20 DEL showed relationship to the meiosis stage (DEL-EMC) (Additional file 3: Table S18), and 39 targets showed EMC specific expressions [34, 35]. Six genes, such as LOC_Os06g12740 (expressed protein) predicted by TCONS_00121908 and LOC_Os10g32740 (zinc finger family protein) predicted by TCONS_00114760, were enriched in zygotene and pachytene of EMC . 28 genes showed the highest expressions in late meiosis stage of megasporocyte , such as LOC_Os12g39420 (nucleobase-ascorbate transporter), LOC_Os06g44300 (WAX2), LOC_Os08g45170 (carboxyl-terminal peptidase) and LOC_Os05g06460 (dihydrolipoyl dehydrogenase) predicted by TCONS_00144592, TCONS_00037710, TCONS_00121908 and TCONS_00019924, respectively. Interestingly, we found that PMC and EMC meiosis-related DEG were absolutely different, and most of them displayed down-regulation in autotetraploid rice PMC (106/110) and EMC (30/39).
Additionally, 14 meiosis-related DEG displayed protein-protein interactions (PPI) (with 8 and 6 DEG belong to two sub-networks) in EMC, including LOC_Os05g06460 (targeted by TCONS_00019924) interacted with LOC_Os03g18740 (oxidoreductase, short chain dehydrogenase/reductase family, targeted by TCONS_00114760), and LOC_Os03g15960 (connectivity = 3) annotated as hsp20/alpha crystallin family protein and targeted by TCONS_00144592 (Additional file 1: Fig. S11). A huge PPI network of 47 meiosis-related DEG, predicted by DEL-PMC, was found in PMC. Two sub-networks were established, including one contained two DEG, and the remaining DEG belonged to second sub-network, which was the main network. Two meiosis DEG, LOC_Os11g12620 (receptor protein kinase CLAVATA1 precursor, regulated by TCONS_00055980) and LOC_Os12g41350 (meiotic asynaptic mutant 1, predicted by TCONS_00057811), as well as LOC_Os01g39630 (RAD51C), were found in this huge sub-network (Additional file 1: Fig. S12).
We categorized the important DELs associated with meiosis-related DEG as shown in Fig. 4, Table 3 and Table 4. The regulatory networks between DELs and DEG may important during anther and ovary development. To confirm the differential expression levels of lncRNAs and target genes, 23 DEL, including 21 DEL-Type 1 and 26 meiosis-related genes were selected to validate the expression patterns using qRT-PCR. The qRT-PCR generally confirmed the high-throughput sequencing results, proving the reliability of RNA-sequencing data used in the present study (Additional file 1: Fig. S13 and S14).
Overexpression of lncRNA57811 reduced rice seed set
To test the potential functions of the DELs, TCONS_00057811 (lncRNA57811), TCONS_00130461 (lncRNA130461) and TCONS_00055980 (lncRNA55980) were selected for further analysis, which showed anther-specific expression patterns. Of these, lncRNA57811 and lncRNA130461 were highly expressed from pre-meiotic interphase to prophase I during pollen development, whereas lncRNA55980 was highly expressed during late meiosis stage (Fig. 5). Additionally, lncRNA57811 and lncRNA130461 exhibited down-regulation at early meiosis stage, while significantly up-regulated at later meiosis stage in Taichung 65-4x compared to diploid counterpart. Then, the expression patterns of lncRNA57811 were confirmed by GUS stained assays (Fig. 6), and GUS activity was visualized only in anther at early pollen development stage i.e. from pre-meiotic stage to prophase I. The staining intensity decreased dramatically after single microspore stage.
lncRNA57811, lncRNA130461 and lncRNA55980 mutants with the background of Taichung 65 were generated by using CRISPR/Cas9 technology; however, CRISPR/Cas9 editing did not show obvious seed phenotypes in T0 generation. We further focused on the lncRNA57811 by generating its overexpressed mutants, and then planted T1 generation. The seed set of the independent CRISPR lines of lncRNA57811 in T1 generation was similar to the wild type, even though 42 bases deletion was detected on its gene body (Fig. 6, Additional file 1: Fig. S15). Interestingly, overexpression of lncRNA57811 (57811OE) showed distinctly low seed set compared to wild type in T1 generation (Fig. 7). The seed set of 57811OE-11 and 57811OE-25 lines were 38.61 and 27.55%, respectively. The pollen fertility of 57811OE was 29.70%, which was significantly lower than the wild type (95.75%). Furthermore, overexpression and CRISPR mutants of LOC_Os12g41350 (12 g) were obtained. 12 g was preferentially expressed in reproductive tissues, while poorly expressed in vegetative tissues. Overexpression of 12 g showed similar morphology to wild type (Additional file 1: Fig. S16). One homozygous mutant with an 8 bp deletion, which resulted in premature termination of the protein sequence, was identified in T0 generation. This line displayed semi-sterility with a seed set of 41.13% in T1 generation. Moreover, LOC_Os12g41350 showed down-regulation in 57811OE mutants.
LncRNAs associated with rice reproductive development shared common characteristics with other species
More than 50% of lncRNAs showed tissue specific patterns in rice anthers and strawberry pollens, which revealed that the expression patterns of lncRNAs play important role in reproductive tissues [24, 27]. In the present study, we characterized 4859 reliable lncRNAs during meiotic stages of anther and ovary in autotetraploid and diploid rice. These lncRNAs shared similar features with the other mammals and plant species, such as shorter transcript length and lower expression levels than protein-coding genes [45,46,47,48]. On the other hand, only a few rice reproductive lncRNAs with conserved sequences among the other plant species were detected, which indicated that rice reproductive lncRNAs may undergo a rapid evolvement [47, 49]. However, 11 lncRNAs were found to be conserved in Arabidopsis, Brachypodium, Zea mays and Sorghum bicolor, which might be the housekeeping lncRNAs and carry out important functions in angiosperms. Moreover, full-length of 10 lncRNAs were obtained by RACE, and six of them were found to be conserved among the neo-tetraploid, autotetraploid, diploid and wild rice, especially lncRNA57811, which indicated that these reproductive-related lncRNAs may play important functions.
LncRNAs act as small RNAs precursors in autotetraploid rice that might play potential role in PMC and EMC meiosis
The majority of lncRNAs act as precursors for small RNA, and siRNAs seem to be the major small RNA population and many lncRNAs serve as precursors for siRNAs [23, 25, 29]. According to a previous study about Arabidopsis, the siRNAs are generated by Pol IV, and their precursors are 30–40 nt , which is in contrast to the model that lncRNAs could be defined as the transcripts longer than 200 bp . So, we focused on the miRNAs and phasiRNAs, and 37 DEL served as precursors for small RNA and all of them displayed down-regulation in autotetraploid rice, which suggested that the down-regulated DEL were associated with small RNAs generation in autotetraploid rice, particularly with anther. PhasiRNAs were found to be associated with male meiocytes in rice and maize [51, 52]. Another study revealed that a SNP in PMS1T (lncRNA), nearby miR2118 recognition site, lead to differential accumulation of the 21 nt-phasiRNAs and resulted in male sterility under long-day conditions . Recently, down-regulation of 24 nt-phasiRNAs was consistent with osa-miR2275d during pollen development in autotetraploid rice . Our results demonstrated that down-regulated DEL, especially DEL-PMC (Type 1), may cause low expression of phasiRNAs, which could harm meiocytes of autotetraploid rice. LncRNAs were predicted to be miRNAs targets or target mimics in rice and maize , however, these results are hard to detect in autotetraploid rice. We speculated that lncRNAs prefer to serve as the small RNA precursors rather than the miRNAs targets or decoys in rice, even in the autotetraploid rice. Moreover, different kinds of predicted targets between DEL and DEM were detected and only few targets (23) were co-regulated by DEL and DEM in our study. These specific targets showed different categories of GO functions in autotetraploid rice anther and ovary, which suggested that lncRNAs and miRNAs take part in various functions during meiosis in autotetraploid rice.
Low expression levels of transposable elements, associated with lncRNAs, may cause great harm to male and female meiocytes in autotetraploid rice
WGD events are an obvious route to genome expansion, and the increase genome size owing to the rapid proliferation of transposable elements . Over 70% rice reproductive lncRNAs were associated with TEs in the present study, and these results were consistent with other plant species [25, 56]. Zhang et al.  indicated that hypermethylation of class II DNA transposons may suppress the expression of neighboring genes, which become a “genome shock” to adapt the genome-dosage effects in autotetraploid rice. In our previous studies, the differentially expressed siRNAs were found to be associated with TEs during male meiosis, which implied the relationship between TEs and sterility in autotetraploid rice [20, 21]. One interesting finding about Arabidopsis male meiocytes illustrated that TEs were abundant and specifically expressed in meiocytes [58, 59]. The activated TEs only detected in vegetative cells but not in sperm cells of the haploid pollen . Male meiocytes represent the final stage of the diploid phase with reactivated TEs, which may involve in chromatin structures and facilitated the precession of meiosis [58, 59]. TEs-lncRNAs detected in male sunflower meiocytes were found to be associated with the different recombination rates between domesticated and wild sunflower . Another study reported that the expression levels of TE-lncRNAs were the highest in testis , and important in the biology of normal testis . We observed low expression levels of TEs and TEs-lncRNAs, and down-regulation of TEs-DEL (Helitron and Stowaway families of Class II) in autotetraploid rice may disturb the chromatin structures or have harmful effects on meiocytes. The TEs and TEs-lncRNAs associated with meiocytes required further studies.
Meiosis-related lncRNAs may play key roles during PMC and EMC meiosis in autotetraploid rice
LncRNAs are revealed to execute direct meiotic roles during meiosis [63, 64], such as UPGRADE2, which was highly up-regulated and attained neofunctionalization in the context of apomeiosis for pollen development in Boechera species , and LDMAR, a SNP could repress the expression level of LDMAR and cause male sterility in rice . Here, 237 DEL-PMC and 20 DEL-EMC were found to be associated with 110 and 39 meiosis-related genes. Most of these meiosis-related genes showed down-regulation and have protein-protein interactions in autotetraploid rice that related to the abnormal meiosis and sterility in our previous study [8, 11], such as LOC_Os01g39630 (RAD51C) and LOC_Os02g57270 (OsTBP1), which is required for the meiosis of male gametocytes [36, 39]. We detected TCONS_00068868, the DEL-PMC (Type 1), which showed down-regulation during meiosis of autotetraploid rice anther, and its target LOC_Os03g58600 (MEL1) showed the same differential expression pattern and validated by qRT-PCR. Rice MEL1 is a key gene involved in meiosis that regulates the cell division of pre-meiotic germ cells, and the proper modification of meiotic chromosomes [37, 66]. Moreover, MEL1 was preferentially associated with 21 nt-phasiRNAs, which play important role during meiosis . We speculated that the abnormal expression of MEL1 might cause abnormal meiosis behavior in autotetraploid rice, and the down-regulation of its lncRNA regulators may triggered the abnormal expressions.
In addition, we detected ten DEL-EMC (Type 1) that showed great relationships with the 27 embryo sac meiosis-related genes [34, 35]. Though the functions of these targets remain largely unknown, they highly expressed and showed meiosis specific expressions during megasporocytes, which suggested their key roles during the meiosis procession of megasporocytes. The regulation of EMC meiosis-related targets in autotetraploid rice was confirmed by the qRT-PCR, such as LOC_Os12g39420, LOC_Os02g46460 and LOC_Os05g06460 regulated by TCONS_00144592, TCONS_00121908 and TCONS_00019924, respectively. These lncRNAs/target pairs probably played essential roles and involved in female meiocytes development of autotetraploid rice, and their abnormal expression patterns might have a strong effect on the embryo sac fertility of autotetraploid rice. Overall, the list of meiosis-DELs in rice anther/ovary will be beneficial for future studies in plant reproduction.
Overexpression of lncRNA57811 showed low seed set and pollen sterility in T1 generation, while CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology when compared to Taichung 65. lncRNA57811 showed anther specific pattern and highly expressed from pre-meiotic interphase to prophase I during pollen development, which demonstrated that lncRNA57811 play key role at early meiosis stage. The higher expression levels of lncRNA57811 were detected in Taichung 65-4x and Taichung 65 at pre meiosis, however, significantly higher expression level of lncRNA57811 was detected in Taichung 65-4x than Taichung 65 during late meiosis stage. These results indicated that the irregular expression pattern of lncRNA57811 might be the main reason for low fertility in Taichung 65-4x.
For knock-down analysis of lncRNAs, CRISPR/Cas9 editing system should not be a first choice tool. In this study, three lncRNAs were selected for generating the editing-mutant by CRISPR/Cas9; however, no obvious change was observed in morphology. This might be happened because 1) CRISPR/Cas9 usually causes frame shift and base replacement in protein sequences, but lncRNAs are the noncoding RNA and no protein. 2) lncRNAs are full of repetitive sequences, deletions/insertions mutation in one repetitive region, and their functions might be replaced by other repetitive regions. 3) the key regions of lncRNAs have not been be edited successfully so far. 4) the functions of these three lncRNAs could be replaced by other lncRNAs/genes. Recently, Wang et al.  used RNAi method to generate the knock-down mutants of lncRNA and acquired desired results. Therefore, there is a need to develop more methods to generate knock-down of lncRNAs in future. Otherwise, more specific-sgRNA targets of CRISPR/Cas9 editing system should be utilized to make sure about the whole/most length of lncRNA to be knocked-out.
LOC_Os12g41350 (meiotic asynaptic mutant 1) was the predicted cis-target of lncRNA57811, and down-regulation of LOC_Os12g41350 in Taichung 65-4x was confirmed by qRT-PCR. CRISPR/Cas9 editing of LOC_Os12g41350 showed semi-sterility in T1 generation. Additionally, LOC_Os12g41350 showed down-regulation in lncRNA57811-OE. The homologous gene of LOC_Os12g41350, namely ASY1, is an important protein required for synapsis and crossover formation during early meiosis in Arabidopsis . Cytological studies have revealed that ASY1 is localized to the regions of chromatin, associated with the axial/lateral elements of meiotic chromosomes, and required for morphogenesis of the synaptonemal complex . Mutant lines asy1 showed identical asynaptic phenotypes during male meiosis . In addition, homologous asynapsis in meiotic prophase I could induce aneuploidy in autopolyploid Arabidopsis asy1 . Moreover, a single crossover per chromosome (reduce crossover rates) could prevent multivalent formation in autotetraploid plants . We hypothesized that lncRNA57811 might be served as cis-regulatory element of the LOC_Os12g41350 in rice, which resulted in pollen sterility. However, it required further experiments to confirm the relationship between lncRNA57811 and LOC_Os12g41350, and to know whether LOC_Os12g41350 is a cis-target of lncRNA57811. The function and the regulation network of lncRNA57811 and LOC_Os12g41350 during pollen development should be further analyzed.
In this study, we utilized high-throughput sequencing to identify the long noncoding RNAs associated with meiosis in autotetraploid rice anther and ovary compared to diploid rice. A total of 4859 confident lncRNAs were identified, which shared common features with other species. Ten lncRNAs were validated by RACE assay, and six were found to be conversed in tetraploid rice. Of the DEL, 237 and 20 DEL were identified as PMC meiosis-related DEL and EMC meiosis-related DEL, respectively. Target prediction showed that these meiosis-related DEL were associated with some important meiosis genes. Interestingly, most of them were associated with the transposable elements and could act as small RNAs precursors, especially the anther/ovary-preferred lncRNAs in autotetraploid rice. Overexpression of lncRNA57811 showed low seed set and pollen sterility indicated that lncRNA played important function during pollen development in polyploid rice. Consequently, the significant differential expression profiles of long noncoding RNAs in the autotetraploid rice may play key roles in regulating the reproductive tissues, thus resulting in low fertility and poor seed set. Our findings provide a new insight on long noncoding RNA in autotetraploid rice reproduction.
Autotetraploid rice, Taichung 65-4x, and its diploid rice, Taichung 65, were used in this study. Taichung 65 is a famous japonica cultivar from Taiwan and bred by Taichung agricultural improvement farm in 1929 . Taichung 65 was kindly characterized and provided by Prof. Guiquan Zhang (College of Agriculture, South China Agricultural University), and is being used by our research group since 1989 [11, 71, 72]. Taichung 65 has been deposited in the International Rice Genebank Collection Information System (IRGCIS, http://www.irgcis.irri.org:81/grc/irgcishome.html) with the accession number of 9434 since 1962. Taichung 65-4x was developed from the chromosome doubling of Taichung 65 by colchicine and was self-crossed for more than 30 generations by our research group (Prof. Xiangdong Liu) at South China Agricultural University (Guangzhou, Guangdong, China). The voucher specimen of Taichung 65-4x has been deposited to our lab but not in any publicly available herbarium. These lines were planted at the experimental farm of South China Agricultural University under field conditions. Anthers at pollen development stages were collected from Taichung 65-4x and Taichung 65 according to Wu et al. . The ovary tissues were collected from the same spikelet (corresponding anther). All the samples were stored at − 80 °C for RNA isolation. Three biological replications were used for each tissue for high-throughput RNA-sequencing.
Transcriptome library construction and sequencing
Total RNA was extracted from the anthers and ovaries by using Trizol reagent (Invitrogen, CA, USA) following the manufacturer’s procedure. The quantity and purity of total RNA were analyzed by Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, CA, USA) with RIN number > 7.0. Approximately 10 μg of total RNA was used to deplete ribosomal RNA by using Epicentre Ribo-Zero Gold Kit (Illumina, San Diego, USA). Divalent cations under high temperature were used to break poly(A)- or poly(A) + RNA fractions is into small pieces. The final cDNA library was constructed according to the methodology of mRNA-Seq sample preparation kit (Illumina, San Diego, USA). Finally, the paired-end sequencing was performed by Illumina Hiseq2500 (LC Sciences, Hangzhou, China) with an average insert size of 300 bp (±50 bp).
Prediction and identification of lncRNAs
Prior to assembly, the low quality reads including reads containing sequencing adaptors, reads containing sequencing primer, nucleotide with quality score lower than 20 (Q < 20), were removed. After that, the clean reads were aligned to the Oryza sativa genome (ftp://ftp.jgi-psf.org/pub/compgen/phytozome/.v10.0/Osativa 204 v7.0/assembly/) using TopHat2 package . The aligned reads were then used to assemble transcripts of each sample independently, using Cufflinks program . The expression levels of the assembled transcripts were further calculated and normalized by Cufflinks using the measurement unit of fragments per kilobase of transcript per million fragments (FPKM). The prediction of lncRNAs from RNA-seq data was performed according to Sun et al. . We discarded transcripts shorter than 200 bp, exons number less than one and transcripts with FPKM < 0.5. The transcripts with five class codes information (i.e. ‘i’, ‘j’, ‘u’, ‘o’, and ‘x’) were extracted. Class codes of ‘u’, ‘i’, and ‘x’ mean that the transcripts are in the intergenic, intronic, and antisense sequence, respectively, of a known gene. Class codes of ‘o’ and ‘j’ mean that the transcripts overlapped with known exons . Then the coding potential of the remaining transcripts was evaluated using coding potential calculator (CPC) software  and Coding Noncoding Index (CNCI) software . The transcripts with CPC score < = − 1 and CNCI score < = 0 were identified as the lncRNA candidates. LncRNAs with P-value < 0.05 and |log2 (fold change ratio)| > 1 were considered as differentially expressed lncRNAs.
Small RNA library construction and sequencing
We further used about 1 μg of total RNA for the small RNA sequencing, and then prepare small RNA libraries according to protocol of TruSeq Small RNA Sample Prep Kits (Illumina, San Diego, USA). After that, we executed the single-end sequencing (36 bp) on an Illumina Hiseq2500 at the LC-BIO (Hangzhou, China) according to the manufacturer’s protocol. For the miRNAs analysis, the raw data was further processed with an in-house program, ACGT101-miR (LC Sciences, Houston, Texas, USA) . First, the reads with common RNA families (snRNA, tRNA, snoRNA, rRNA), low complexity, repeats, junk and adapter dimers were removed. Then, the unique sequences (18-25 nt in length) were BLASTed to miRBase 20.0 (ftp://mirbase.org/pub/mirbase/) to identify known miRNAs. The novel miRNAs were predicted by using RNAfold software (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi). miRNAs with P-value < 0.05 and |log2 (fold change ratio)| > 1 were considered as differentially expressed miRNAs. Targeted genes of differentially expressed miRNAs were predicted by using TargetFinder . For the phasiRNAs analysis, 21 nt- and 24 nt-phasiRNAs were systemically characterized using PhaseTank packet . PhasiRNAs with P-value < 0.05 and |log2 (fold change ratio)| > 1 were considered as differentially expressed phasiRNAs.
Characterization of lncRNAs
To investigate the conservation of lncRNAs, all the lncRNAs sequences were blasted against the genome sequences of Arabidopsis thaliana (assembly TAIR10), Zea mays (assembly B73 RefGen_v3), Brachypodium distachyon (assembly Brachypodium_distachyon_v2.0) and Sorghum bicolor (assembly Sorbi1) that download from NCBI with a cut-off E-value <1e-10 by using BLASTN. The lncRNAs with more than 20% sequence matched to other genomes were defined as conversed lncRNAs. Additionally, we also conducted the lncRNAs associated with repetitive element by Repeat Masker (http://www.repeatmasker.org/; based on the Repeat Database namely repeatmaskerlibraries-20,160,829 from http://www.girinst.org/) with the ‘slow’ option. To identify the lncRNAs that may act as precursors of miRNAs or phasiRNAs, the small RNAs were aligned to the identified lncRNAs by following the model of Boerner and McGinnis . LncRNAs as targets of miRNAs were predicted by psRNATarget , with expectation ≤3 . The target mimics were predicted according to the method of Meng et al. .
Target gene prediction and functional analysis of lncRNAs
LncRNAs may play cis- and trans-acting regulation on the genes . Two algorithms were used for cis- and trans- target prediction. The first algorithm was used by python script to detect potential cis target genes (differentially expressed genes) that are physically close to lncRNAs (differentially expressed lncRNAs). The genes transcribed within a 100 kb window upstream or downstream of lncRNAs were considered as potential cis target genes . The second algorithm was used to identify potential trans targets based on RNA sequence complementarity (mRNA-lncRNA) and RNA duplex energy prediction [80, 81]. Based on the RNA duplex energy (Free energy < − 50) assessed by RIsearch algorithms , the differentially expressed genes and lncRNAs were considered as potentially involved in trans-interaction. Then, functional analysis of the lncRNAs cis/trans targets were done by using the AgriGO . Significance was expressed as a P-value < 0.05. Protein-protein interaction networks were performed using STRING .
Quantitative real-time PCR (qRT-PCR) analysis
Reverse-transcribed of total RNAs isolated from ovaries and anthers, and qRT-PCR assay were performed according to Li et al. . The qRT-PCR products of each candidate lncRNAs were Sanger sequenced by BGI Genomics, and the amplified fragments were matched to the lncRNAs. The specificity of the amplified fragments was checked using the generated melting curve. GAPDH was used as an internal control gene. All qRT-PCR amplifications were carried out in three replications, and the results are presented as the mean ± standard deviations. The 2–ΔΔCT method was employed to calculate the relative expression level . The primers were designed by Primer 5.0 software, and are listed in Additional file 3: Table S19.
5′- and 3′-rapid amplification of cDNA ends (RACE) assay
Full-length amplification of lncRNA was performed with SMARTer RACE 5′/3′ Kit (Clontech, USA) according to the manufacturer’s instructions. PCR was performed with Universal Primer A Mix and gene-specific primers using 5′/3′-cDNA as the template (Additional file 3: Table S19). The 5′- and 3′-RACE products were gel-purified and cloned. At least 12 clones were selected for Sanger sequencing to obtain the longest fragment. Finally, long distance PCR using primers designed from the extreme 5′ and 3′ ends of the candidate lncRNA and the 5′-cDNA as a template to generating full-length lncRNA.
Selection of the target site and plasmid construction
Specific gRNA targeting lncRNA57811, lncRNA55980, and lncRNA130461 were selected using the online tool (http://skl.scau.edu.cn/targetdesign/). The gRNA expression cassettes were generated by overlapping PCR. The target sequences were ligated between the corresponding promoter and sgRNA by the first PCR. The second PCR was used to induce a BsaӀ restriction site. The amplified fragments were then assembled on the pYLCRISPR/Cas9Ubi-H binary plasmid. To construct the overexpressed lncRNA57811 vector, the cDNA of lncRNA57811 were obtained and fused into the POX vector (digestion by HindIII and BamHI) using Hieff Clone® Plus One Step Cloning Kit (Hieff, Shanghai, China). To construct lncRNA57811-Pro::GUS vector, 2.0 kb upstream of ATG sequence was amplified from Taichung 65 genomic DNA and cloned into pCAMBIA1305.1. These constructs were transferred to Agrobacterium EHA105, and then were transformed into rice callus (Taichung 65). The specific primers are listed in Additional file 3: Table S19.
Availability of data and materials
Transcriptome and small RNA sequencing data are available from the NCBI under the accession numbers PRJNA395615 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA395615) and PRJNA396435 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA396435). The accession numbers of ten lncRNAs are MN481253 (lncRNA57811), MN481254 (lncRNA130461), MN481255 (lncRNA55980), MN481256 (lncRNA91337), MN481257 (lncRNA68868), MN481258 (lncRNA111916), MN481259 (lncRNA130471), MN481260 (lncRNA45430), MN481261 (lncRNA13598) and MN481262 (lncRNA121908). The information about lncRNAs are also available in Additional file 3: Table S20. All data supporting the conclusions described here are provided in tables, figures and additional files.
differentially expressed genes
differentially expressed lncRNAs
differentially expressed miRNAs
embryo sac mother cell
pollen mother cell
Wendel JF, Jackson SA, Meyers BC, Wing RA. Evolution of plant genome architecture. Genome Biol. 2016;17:37.
Liu S, Adams KL. Dramatic change in function and expression pattern of a gene duplicated by polyploidy created a paternal effect gene in the Brassicaceae. Mol Biol Evol. 2010;27:2817–28.
Hughes TE, Langdale JA, Kelly S. The impact of widespread regulatory neofunctionalization on homeolog gene evolution following whole-genome duplication in maize. Genome Res. 2014;24:1348–55.
Mason AS, Snowdon RJ. Oilseed rape: learning about ancient and recent polyploid evolution from a recent crop species. Plant Biol. 2016;18:883–92.
Shahid MQ, Liu GF, Li JQ, Naeem M, Liu X. Heterosis and gene action study of agronomic traits in diploid and autotetraploid rice. Acta Agric Scand Sect B Soil Plant Sci. 2011;61:23–32.
Guo H, Mendrikahy JN, Xie L, Deng J, Lu Z, Wu J, et al. Transcriptome analysis of neo-tetraploid rice reveals specific differential gene expressions associated with fertility and heterosis. Sci Rep. 2017;7:40139.
Shahid MQ, Xu HM, Lin SQ, Chen ZX, Naeem M, Li YJ, et al. Genetic analysis and hybrid vigor study of grain yield and other quantitative traits in autotetraploid rice. Pak J Bot. 2012;44:237–46.
Guo H, Feng J, Lu Y, Liu X. Abnormal structure of embryo sac in autotetraploid rice. Rice Sci. 2006;13:257–64.
Shahid MQ, Li YJ, Saleem MF, Naeem M, Wei CM, Liu XD. Yield and yield components in autotetraploid and diploid rice genotypes (indica and japonica) sown in early and late seasons. Aust J Crop Sci. 2013;7:632–41.
Shahid MQ, Sun J, Wei C, Zhang P, Liu X. Studies on the abnormality of embryo sac and pollen fertility in autotetraploid rice during different growing seasons. Pak J Bot. 2010;42:7–19.
Wu J, Shahid MQ, Guo H, Yin W, Chen Z, Wang L, et al. Comparative cytological and transcriptomic analysis of pollen development in autotetraploid and diploid rice. Plant Reprod. 2014;27:181–96.
Luo Q, Li Y, Shen Y, Cheng Z. Ten years of gene discovery for meiotic event control in rice. J Genet Genomics. 2014;41:125–37.
Cifuentes M, Grandont L, Moore G, Chevre AM, Jenczewski E. Genetic regulation of meiosis in polyploid species: new insights into an old question. New Phytol. 2010;186:29–36.
Bomblies K, Madlung A. Polyploidy in the Arabidopsis genus. Chromosom Res. 2014;22:117–34.
Szadkowski E, Eber F, Huteau V, Lodé M, Huneau C, Belcram H, et al. The first meiosis of resynthesized Brassica napus, a genome blender. New Phytol. 2010;186:102–12.
He JH, Shahid MQ, Li YJ, Guo HB, Cheng XA, Liu XD, et al. Allelic interaction of F1 pollen sterility loci and abnormal chromosome behaviour caused pollen sterility in intersubspecific autotetraploid rice hybrids. J Exp Bot. 2011;62:4433–45.
Yant L, Hollister JD, Wright KM, Arnold BJ, Higgins JD, Franklin FCH, et al. Meiotic adaptation to genome duplication in Arabidopsis arenosa. Curr Biol. 2013;23:2151–6.
Lloyd A, Bomblies K. Meiosis in autopolyploid and allopolyploid Arabidopsis. Curr Opin Plant Biol. 2016;30:116–22.
Wu J, Shahid MQ, Chen L, Chen Z, Wang L, Liu X, et al. Polyploidy enhances F1 pollen sterility loci interactions that increase meiosis abnormalities and pollen sterility in autotetraploid rice. Plant Physiol. 2015;169:2700–17.
Li X, Shahid MQ, Xia J, Lu Z, Fang N, Wang L, et al. Analysis of small RNAs revealed differential expressions during pollen and embryo sac development in autotetraploid rice. BMC Genomics. 2017;18:129.
Li X, Shahid MQ, Wu J, Wang L, Liu X, Lu Y. Comparative small RNA analysis of pollen development in autotetraploid and diploid rice. Int J Mol Sci. 2016;17:499.
Taylor DH, Chu ET, Spektor R, Soloway PD. Long non-coding RNA regulation of reproduction and development. Mol Reprod Dev. 2015;82:932–56.
Li L, Eichten SR, Shimizu R, Petsch K, Yeh CT, Wu W, et al. Genome-wide discovery and characterization of maize long non-coding RNAs. Genome Biol. 2014;15:R40.
Kang C, Liu Z. Global identification and analysis of long non-coding RNAs in diploid strawberry Fragaria vesca during flower and fruit development. BMC Genomics. 2015;16:815.
Flórez-Zapata NM, Reyes-Valdés MH, Martinez O. Long non-coding RNAs are major contributors to transcriptome changes in sunflower meiocytes with different recombination rates. BMC Genomics. 2016;17:490.
Ding J, Lu Q, Ouyang Y, Mao H, Zhang P, Yao J, et al. A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice. Proc Natl Acad Sci U S A. 2012;109:2654–9.
Zhang YC, Liao JY, Li ZY, Yu Y, Zhang JP, Li QF, et al. Genome-wide screening and functional analysis identify a large number of long noncoding RNAs involved in the sexual reproduction of rice. Genome Biol. 2014;15:512.
Wang Y, Luo X, Sun F, Hu J, Zha X, Su W, et al. Overexpressing lncRNA LAIR increases grain yield and regulates neighbouring gene cluster expression in rice. Nat Commun. 2018;9:3516.
Boerner S, McGinnis KM. Computational identification and functional predictions of long noncoding RNA in Zea mays. PLoS One. 2012;7:e43047.
Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35:W345–9.
Dai X, Zhao PX. psRNATarget: a plant small RNA target analysis server. Nucleic Acids Res. 2011;39:W155–9.
Fujita M, Horiuchi Y, Ueda Y, Mizuta Y, Kubo T, Yano K, et al. Rice expression atlas in reproductive development. Plant Cell Physiol. 2010;51:2060–81.
Deveshwar P, Bovill WD, Sharma R, Able JA, Kapoor S. Analysis of anther transcriptomes to identify genes contributing to meiosis and male gametophyte development in rice. BMC Plant Biol. 2011;11:78.
Kubo T, Fujita M, Takahashi H, Nakazono M, Tsutsumi N, Kurata N. Transcriptome analysis of developing ovules in rice isolated by laser microdissection. Plant Cell Physiol. 2013;54:750–65.
Wu Y, Yang L, Cao A, Wang J. Gene expression profiles in rice developing ovules provided evidence for the role of sporophytic tissue in female gametophyte development. PLoS One. 2015;10:e141613.
Kou Y, Chang Y, Li X, Xiao J, Wang S. The rice RAD51C gene is required for the meiosis of both female and male gametocytes and the DNA repair of somatic cells. J Exp Bot. 2012;63:5323–35.
Nonomura K, Morohoshi A, Nakano M, Eiguchi M, Miyao A, Hirochika H, et al. A germ cell-specific gene of the ARGONAUTE family is essential for the progression of premeiotic mitosis and meiosis during sporogenesis in rice. Plant Cell. 2007;19:2583–94.
Nonomura K, Miyoshi K, Eiguchi M, Suzuki T, Miyao A, Hirochika H, et al. The MSP1 gene is necessary to restrict the number of cells entering into male and female sporogenesis and to initiate anther wall formation in rice. Plant Cell. 2003;15:1728–39.
Hong JP, Byun MY, Koo DH, An K, Bang JW, Chung IK, et al. Suppression of RICE TELOMERE BINDING PROTEIN 1 results in severe and gradual developmental defects accompanied by genome instability in rice. Plant Cell. 2007;19:1770–81.
Tao J, Zhang L, Chong K, Wang T. OsRAD21-3, an orthologue of yeast RAD21, is required for pollen development in Oryza sativa. Plant J. 2007;51:919–30.
Morozumi Y, Ino R, Ikawa S, Mimida N, Shimizu T, Toki S, et al. Homologous pairing activities of two rice RAD51 proteins, RAD51A1 and RAD51A2. PLoS One. 2013;8:e75451.
Hord CL, Chen C, DeYoung BJ, Clark SE, Ma H. The BAM1/BAM2 receptor-like kinases are important regulators of Arabidopsis early anther development. Plant Cell. 2006;18:1667–80.
Dean PJ, Siwiec T, Waterworth WM, Schlögelhofer P, Armstrong SJ, West CE. A novel ATM-dependent X-ray-inducible gene is essential for both plant meiosis and gametogenesis. Plant J. 2009;58:791–802.
Sanchez-Moran E, Santos JL, Jones GH, Franklin FC. ASY1 mediates AtDMC1-dependent interhomolog recombination during meiosis in Arabidopsis. Genes Dev. 2007;21:2220–33.
Pauli A, Valen E, Lin MF, Garber M, Vastenhouw NL, Levin JZ, et al. Systematic identification of long noncoding RNAs expressed during zebrafish embryogenesis. Genome Res. 2012;22:577–91.
Song X, Liu G, Huang Z, Duan W, Tan H, Li Y, et al. Temperature expression patterns of genes and their coexpression with LncRNAs revealed by RNA-Seq in non-heading Chinese cabbage. BMC Genomics. 2016;17:297.
Tian J, Song Y, Du Q, Yang X, Ci D, Chen J, et al. Population genomic analysis of gibberellin-responsive long non-coding RNAs in Populus. J Exp Bot. 2016;67:2467–82.
Gao X, Ye J, Yang C, Zhang K, Li X, Luo L, et al. Screening and evaluating of long noncoding RNAs in the puberty of goats. BMC Genomics. 2017;18:164.
Zhang G, Duan A, Zhang J, He C. Genome-wide analysis of long non-coding RNAs at the mature stage of sea buckthorn (Hippophae rhamnoides Linn) fruit. Gene. 2017;596:130–6.
Zhai J, Bischof S, Wang H, Feng S, Lee TF, Teng C, et al. A one precursor one siRNA model for pol IV-dependent siRNA biogenesis. Cell. 2015;163:445–55.
Komiya R, Ohyanagi H, Niihama M, Watanabe T, Nakano M, Kurata N, et al. Rice germline-specific Argonaute MEL1 protein binds to phasiRNAs generated from more than 700 lincRNAs. Plant J. 2014;78:385–97.
Zhai J, Zhang H, Arikit S, Huang K, Nan GL, Walbot V, et al. Spatiotemporally dynamic, cell-type-dependent premeiotic and meiotic phasiRNAs in maize anthers. Proc Natl Acad Sci U S A. 2015;112:3146–51.
Fan Y, Yang J, Mathioni SM, Yu J, Shen J, Yang X, et al. PMS1T, producing phased small-interfering RNAs, regulates photoperiod-sensitive male sterility in rice. Proc Natl Acad Sci U S A. 2016;113:15144–9.
Wu HJ, Wang ZM, Wang M, Wang XJ. Widespread long noncoding RNAs as endogenous target mimics for MicroRNAs in plants. Plant Physiol. 2013;161:1875–84.
Verde I, Abbott AG, Scalabrin S, Jung S, Shu S, Marroni F, et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity, domestication and genome evolution. Nat Genet. 2013;45:487–94.
Wang X, Ai G, Zhang C, Cui L, Wang J, Li H, et al. Expression and diversification analysis reveals transposable elements play important roles in the origin of Lycopersicon-specific lncRNAs in tomato. New Phytol. 2016;209:1442–55.
Zhang J, Liu Y, Xia EH, Yao QY, Liu XD, Gao LZ. Autotetraploid rice methylome analysis reveals methylation variation of transposable elements and their effects on gene expression. Proc Natl Acad Sci U S A. 2015;112:E7022–9.
Chen C, Farmer AD, Langley RJ, Mudge J, Crow JA, May GD, et al. Meiosis-specific gene discovery in plants: RNA-Seq applied to isolated Arabidopsis male meiocytes. BMC Plant Biol. 2010;10:280.
Yang H, Lu P, Wang Y, Ma H. The transcriptome landscape of Arabidopsis male meiocytes from high-throughput sequencing: the complexity and evolution of the meiotic process. Plant J. 2011;65:503–16.
Slotkin RK, Vaughn M, Borges F, Tanurdzić M, Becker JD, Feijó JA, et al. Epigenetic reprogramming and small RNA silencing of transposable elements in pollen. Cell. 2009;136:461–72.
Kelley D, Rinn J. Transposable elements reveal a stem cell-specific class of long noncoding RNAs. Genome Biol. 2012;13:R107.
Young JM, Whiddon JL, Yao Z, Kasinathan B, Snider L, Geng LN, et al. DUX4 binding to retroelements creates promoters that are active in FSHD muscle and testis. PLoS Genet. 2013;9:e1003947.
Hongay CF, Grisafi PL, Galitski T, Fink GR. Antisense transcription controls cell fate in Saccharomyces cerevisiae. Cell. 2006;127:735–45.
van Werven FJ, Neuert G, Hendrick N, Lardenois A, Buratowski S, van Oudenaarden A, et al. Transcription of two long noncoding RNAs mediates mating-type control of gametogenesis in budding yeast. Cell. 2012;150:1170–81.
Mau M, Corral JM, Vogel H, Melzer M, Fuchs J, Kuhlmann M, et al. The conserved chimeric transcript UPGRADE2 is associated with unreduced pollen formation and is exclusively found in apomictic Boechera species. Plant Physiol. 2013;163:1640–59.
Liu H, Nonomura K. A wide reprogramming of histone H3 modifications during male meiosis I in rice is dependent on the Argonaute protein MEL1. J Cell Sci. 2016;129:3553–61.
Armstrong SJ, Caryl AP, Jones GH, Franklin FC. Asy1, a protein required for meiotic chromosome synapsis, localizes to axis-associated chromatin in Arabidopsis and Brassica. J Cell Sci. 2002;115:3645–55.
Caryl AP, Armstrong SJ, Jones GH, Franklin FC. A homologue of the yeast HOP1 gene is inactivated in the Arabidopsis meiotic mutant asy1. Chromosoma. 2000;109:62–71.
Wei F, Zhang GS. Meiotically asynapsis-induced aneuploidy in autopolyploid Arabidopsis thaliana. J Plant Res. 2010;123:87–95.
Lin SC, Min SK. Chinese rice varieties and their family tree. Shanghai: Shanghai science and technology press; 1991. (in Chinese)
Shahid MQ, Chen FY, Li HY, Wang SZ, Chen FF, et al. Double-neutral genes, Sa-n and Sb-n, for pollen fertility in rice to overcome indica×japonica hybrid sterility. Crop Sci. 2013;53:164–76.
Zhang GQ, Lu YG. Genetic studies of the hybrid sterility in cultivated rice (Oryza sativa) I. Diallel analysis of the hybrid sterility among isogenic F1 sterile lines. Rice Sci. 1989;3:97–101.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7:562–78.
Sun L, Luo H, Bu D, Zhao G, Yu K, Zhang C, et al. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 2013;41:e166.
Fahlgren N, Carrington JC. miRNA target prediction in plants. Methods Mol Biol. 2010;592:51–7.
Guo Q, Qu X, Jin W. PhaseTank: genome-wide computational identification of phasiRNAs and their regulatory cascades. Bioinformatics. 2015;31:284–6.
Meng Y, Shao C, Wang H, Jin Y. Target mimics: an embedded layer of microRNA-involved gene regulatory networks in plants. BMC Genomics. 2012;13:197.
Liu Y, Sun Y, Li Y, Bai H, Xue F, Xu S, et al. Analyses of long non-coding RNA and mRNA profiling using RNA sequencing in chicken testis with extreme sperm motility. Sci Rep. 2017;7:9055.
Shi Y, Tu H, Chen X, Zhang Y, Chen L, Liu Z, et al. The long non-coding RNA expression profile of Coxsackievirus A16 infected RD cells identified by RNA-seq. Virol Sin. 2016;31:131–41.
Danis J, Göblös A, Bata-Csörgő Z, Kemény L, Széll M. PRINS non-coding RNA regulates nucleic acid-induced innate immune responses of human keratinocytes. Front Immunol. 2017;8:1053.
Wenzel A, Akbasli E, Gorodkin J. RIsearch: fast RNA-RNA interaction search using a simplified nearest-neighbor energy model. Bioinformatics. 2012;28:2738–46.
Du Z, Zhou X, Ling Y, Zhang Z, Su Z. agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 2010;38:W64–70.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.
Li X, Yu H, Jiao Y, Shahid MQ, Wu J, Liu X. Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice. PLoS One. 2018;13:e0201854.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method. Methods. 2001;25:402–8.
The authors thank Dr. Jinwen Wu, Zhixiong Chen, Lan Wang, and Ms. Shuhong Yu and other lab members for assistance. We are very thankful to Acad. Yaoguang Liu for providing the POX vector and the CRISPR/Cas 9 editing system.
This work was supported by the Guangzhou Science and Technology Key Program (201707020015, XD Liu), National Natural Science Foundation of China (31571625, XD Liu), the Key Realm R & D Program of Guangdong Province (2018B020202012, ZM Zhang), and China Postdoctoral Science Foundation (2018 M643101, X Li). The funders had no role in the design of the study and collection, analysis, and interpretation of data and in writing of the manuscript.
Ethics approval and consent to participate
Plant samples used in the study were not collected from national park or natural reserve. According to national and local legislation, no specific permission was required for collecting these plants. We confirm that this complies with national guidelines and no formal ethics approval was required in this particular case.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Identification of the lncRNAs in the rice reproductive tissues. Fig. S2 Characteristics of lncRNAs in rice reproductive tissues. Fig. S3 Sequence alignment results between lncRNAs obtained by RACE assay and transcripts predicted by transcriptome. Fig. S4 Phylogenic tree analysis of lncRNAs based on the rice dataset. Fig. S5 Anther and ovary-preferred lncRNAs in diploid (A) and autotetraploid rice (B). Fig. S6 Venn analysis of DEL-anther and DEL-ovary in autotetraploid rice. Fig. S7 Classification of differentially expressed lncRNAs in anther (A) and ovary (B) of autotetraploid rice. Fig. S8 Transposable elements (TEs) associated with lncRNAs in autotetraploid rice. Fig. S9 Validation of the transposable elements (TEs) associated genes and TEs-lncRNAs in autotetraploid compared to diploid rice. Fig. S10 Venn analysis of the predicted targets of differentially expressed lncRNAs (DEL) and differentially expressed miRNAs (DEM). Fig. S11 Protein-protein interactions between the EMC meiosis-related targets of differentially expressed lncRNAs. Fig. S12 Protein-protein interactions between the PMC meiosis-related targets of differentially expressed lncRNAs. Fig. S13 Validation of the differentially expressed lncRNAs (DEL) in autotetraploid and diploid rice. Fig. S14 Validation of the targets of differentially expressed lncRNAs in autotetraploid and diploid rice. Fig. S15 Agronomic traits of wild type and mutants in T1 generation. Fig. S16 Phenotypes of LOC_Os12g41350-CR and LOC_Os12g41350-OE in T1 generation.
Characteristics of all lncRNAs identified in this study. Table S2. Conservation of lncRNAs in other plant species. Table S3. Differentially expressed lncRNAs (DEL) in anther and ovary of Taichung 65-4x compared to Taichung 65. Table S4. The lncRNAs selected by aligment assay compared to a previous study. Table S5. Summary of known and predicted miRNAs detected in this study. Table S6. Overview of 21 nt- and 24 nt-phasiRNAs detected in this study. Table S7. LncRNAs corresponding to miRNA precursors. Table S8. LncRNAs corresponding to phasiRNA precursors. Table S9. A list of lncRNAs predicted to be the targets of miRNAs. Table S10. Transposable elements-associated with differentially expressed lncRNAs in autotetraploid rice.
Classification of differentially expressed lncRNAs related to transposable elements (TEs) in autotetraploid rice. Table S12. The potential target genes of differentially expressed lncRNAs detected by cis-regulation. Table S13. The potential target genes of differentially expressed lncRNAs detected by trans-regulation. Table S14. GO (Gene Ontology) analysis of predicted targets of differentially expressed lncRNAs. Table S15. Overall information about the predicted targets of differentially expressed miRNAs (DEM). Table S16. GO (Gene Ontology) analysis of predicted targets of differentially expressed lncRNAs and differentially expressed miRNAs. Table S17. The most important lncRNA-target pairs related to PMC meiosis in autotetraploid rice. Table S18. The important lncRNA-target pairs related to EMC meiosis in autotetraploid rice. Table S19. A list of primers used in the present study. Table S20. The information about ten lncRNAs detected by RACE in the present study.
About this article
Cite this article
Li, X., Shahid, M.Q., Wen, M. et al. Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice. BMC Plant Biol 20, 82 (2020). https://doi.org/10.1186/s12870-020-2290-0
- Autotetraploid rice
- Embryo sac mother cell
- lncRNAs (long non-coding RNAs)
- Male sterility
- Pollen mother cell