Gene Expression Profiles in Rice Developing Ovules Provided Evidence for the Role of Sporophytic Tissue in Female Gametophyte Development

OsmiRNA5488 Regulates the Development of Embryo Sacs and Targets OsARF25 in Rice (Oryza sativa L.)

Small RNAs are a class of non-coding RNAs that typically range from 20 to 24 nucleotides in length. Among them, microRNAs (miRNAs) are particularly important regulators for plant development. The biological function of the conserved miRNAs has been studied extensively in plants, while that of the species-specific miRNAs has been studied in-depth. In this study, the regulatory role of a rice-specific OsmiRNA5488 (OsmiR5488) was characterized with the miR5488-overexpressed line (miR5488-OE) and miR5488-silenced line (STTM-5488). The seed-setting rate was notably reduced in miR5488-OE lines, but not in STTM-5488 lines. Cytological observation demonstrated the different types of abnormal mature embryo sacs, including the degeneration of embryo sacs and other variant types, in miR5488-OE lines. The percentage of the abnormal mature embryo sacs accounted for the reduced value of the seed-setting rate. Furthermore, OsARF25 was identified as a target of OsmiR5488 via RNA ligase-mediated 3'-amplifification of cDNA ends, dual luciferase assays, and transient expression assays. The primary root length was decreased with the increases in auxin concentrations in miR5488-OE lines compared to wild-type rice. Summarily, our results suggested that OsmiR5488 regulates the seed-setting rate and down-regulates the targeted gene OsARF25.

SMALL REPRODUCTIVE ORGANS, a SUPERMAN-like transcription factor, regulates stamen and pistil growth in rice

Organ size is determined mainly by cell division and cell expansion. Several genetic factors regulating development of plant lateral organs have been characterized, but those involved in determining reproductive organ size and separation in rice (Oryza sativa) remain unknown. We have isolated the rice gene SMALL REPRODUCTIVE ORGANS (SRO) encoding a nucleus-localized Cys2His2 (C₂ H₂ ) zinc finger protein orthologous to Arabidopsis transcription factor (TF) SUPERMAN (SUP). Combined developmental, genetic, histological and transcriptomic analyses were used to determine the function of SRO in regulating reproductive organ size. SRO affects genes involved in cell division, cell expansion and phytohormone signalling in the rice flower. SRO is specifically expressed in the first stages of stamen filament development to regulate their correct formation and separation. In addition, SRO noncell-autonomously regulates the size and functionality of male and female reproductive organs. The B-class MADS-box gene OsMADS16/SPW1 is epistatic to SRO, whereas SRO regulates reproductive organ specification and floral meristem determinacy synergistically with C-class genes OsMADS3 and OsMADS58. These findings provide insights into how an evolutionarily conserved TF has a pivotal role in reproductive organ development in core eudicots and monocots, through partially conserved expression, function and regulatory network.

A study of the heterochronic sense/antisense RNA representation in florets of sexual and apomictic Paspalum notatum

Background

Apomixis, an asexual mode of plant reproduction, is a genetically heritable trait evolutionarily related to sexuality, which enables the fixation of heterozygous genetic combinations through the development of maternal seeds. Recently, reference floral transcriptomes were generated from sexual and apomictic biotypes of Paspalum notatum, one of the most well-known plant models for the study of apomixis. However, the transcriptome dynamics, the occurrence of apomixis vs. sexual expression heterochronicity across consecutive developmental steps and the orientation of transcription (sense/antisense) remain unexplored.

Results

We produced 24 Illumina TruSeq®/ Hiseq 1500 sense/antisense floral transcriptome libraries covering four developmental stages (premeiosis, meiosis, postmeiosis, and anthesis) in biological triplicates, from an obligate apomictic and a full sexual genotype. De novo assemblies with Trinity yielded 103,699 and 100,114 transcripts for the apomictic and sexual samples respectively. A global comparative analysis involving reads from all developmental stages revealed 19,352 differentially expressed sense transcripts, of which 13,205 (68%) and 6147 (32%) were up- and down-regulated in apomictic samples with respect to the sexual ones. Interestingly, 100 differentially expressed antisense transcripts were detected, 55 (55%) of them up- and 45 (45%) down-regulated in apomictic libraries. A stage-by-stage comparative analysis showed a higher number of differentially expressed candidates due to heterochronicity discrimination: the highest number of differential sense transcripts was detected at premeiosis (23,651), followed by meiosis (22,830), postmeiosis (19,100), and anthesis (17,962), while the highest number of differential antisense transcripts were detected at anthesis (495), followed by postmeiosis (164), meiosis (120) and premeiosis (115). Members of the AP2, ARF, MYB and WRKY transcription factor families, as well as the auxin, jasmonate and cytokinin plant hormone families appeared broadly deregulated. Moreover, the chronological expression profile of several well-characterized apomixis controllers was examined in detail.

Conclusions

This work provides a quantitative sense/antisense gene expression catalogue covering several subsequent reproductive developmental stages from premeiosis to anthesis for apomictic and sexual P. notatum, with potential to reveal heterochronic expression between reproductive types and discover sense/antisense mediated regulation. We detected a contrasting transcriptional and hormonal control in apomixis and sexuality as well as specific sense/antisense modulation occurring at the onset of parthenogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07450-3.

Comparative Transcriptomic and Metabolic Analyses Reveal the Molecular Mechanism of Ovule Development in the Orchid, Cymbidium sinense

Ovule development is pivotal to plant reproduction and seed development. Cymbidium sinense (Orchidaceae) has high ornamental value due to its pleasant aroma and elegant floral morphology. The regulatory mechanism underlying ovule development in orchids, especially C. sinense, is largely unknown and information on the C. sinense genome is very scarce. In this study, a combined analysis was performed on the transcriptome and non-targeted metabolomes of 18 C. sinense 'Qi Jian Hei Mo' ovule samples. Transcriptome analysis assembled gene-related information related to six growth stages of C. sinense ovules (S1-S6, equivalent to 30, 35, 42, 46, 53, and 60 days after pollination). Illumina sequencing technology was used to obtain the complete set of transcriptome sequences of the 18 samples. A total of 81,585 unigene sequences were obtained after assembly, 24,860 (30.47%) of which were functionally annotated. Using transcriptome sequencing technology, a total of 9845 differentially expressed unigenes (DEUs) were identified in C. sinense ovules that were assigned to specific metabolic pathways according to the Kyoto Encyclopedia of Genes and Genomes (KEGG). DEUs associated with transcription factors (TFs) and phytohormones were identified and analyzed. The TFs homeobox and MADS-box were associated with C. sinense ovule development. In particular, the phytohormones associated with DEUs such as indole-3-acetic acid (IAA), cytokinin (CK), gibberellin (GA), abscisic acid (ABA), brassinosteroid (BR), and jasmonate (JA), may have important regulatory effects on C. sinense ovule development. Metabolomic analysis showed an inconsistent number of KEGG annotations of differential metabolites across comparisons (S2_vs_S4, S2_vs_S5, and S4_vs_S5 contained 23, 26, and 3 annotations, respectively) in C. sinense ovules. This study provides a valuable foundation for further understanding the regulation of orchid ovule development and formation, and establishes a theoretical background for future practical applications during orchid cultivation.

Comparative Expression Profiling Reveals Genes Involved in Megasporogenesis

Megasporogenesis is a key step during ovule development in angiosperms, but the small number and inaccessibility of these cells have hampered molecular and genome-wide studies. Thus, many questions remain regarding the molecular basis of cell specification, differentiation, and development in the female gametophyte. Here, taking advantage of the correlation between spikelet length and ovule development in rice (Oryza sativa), we studied the transcriptome dynamics of young ovules at three stages, the archesporial cell, the megaspore mother cell before meiosis, and the functional megaspore after meiosis, using expression profiling based on RNA sequencing. Our analysis showed that 5,274 genes were preferentially expressed in ovules during megasporogenesis as compared to ovules at the mature female gametophyte stage. Out of these, 958 (18.16%) genes were archesporial cell- and/or megaspore mother cell-preferential genes, and represent a significant enrichment of genes involved in hormone signal transduction and plant pathogen interaction pathways, as well as genes encoding transcription factors. The expression patterns of nine genes that were preferentially expressed in ovules of different developmental stages, including the OsERECTA2 (OsER2) receptor-like kinase gene, were confirmed by in situ hybridization. We further characterized the OsER2 loss-of-function mutant, which had an excessive number of female germline cells and an abnormal female gametophyte, suggesting that OsER2 regulates germline cell specification during megasporogenesis in rice. These results expand our understanding of the molecular control of megasporogenesis in rice and contribute to the functional studies of genes involved in megasporogenesis.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility 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.

Transcriptome analysis reveals the different compatibility between LAAA × AA and LAAA × LL in Lilium

To unveil the mechanism of the compatibility of odd-allotetraploid lily (LAAA) as female with diploid male lily, the differences of expressed unigenes in the ovaries and leaves between LAAA × AA and LAAA × LL were investigated using transcriptome analysis. The results showed the fruits of LAAA × AA well developed, while those of LAAA × LL aborted. The number of differentially expressed genes was less in the ovaries of LAAA × AA than those of LAAA × LL, but it showed opposite trend in those of leaves. The unigenes related with auxins, cytokinins, gibberellins, antioxidants, expansins, chlorophylls, carbohydrates, transport proteins were usually up-expressed in the ovaries and leaves of LAAA × AA but not in LAAA × LL; while those of abscisic acid, ethylene, jasmonic acid, and salicylic acid were increased in the ovaries or leaves of LAAA × LL but not in LAAA × AA. The up-expressed unigenes in the ovaries and leaves of LAAA × AA played positive roles in its fruit development because the products of the genes, like phytohormones and antioxidants, had functions protecting leaves from senescence or scavenging ROS, and thus LAAA was compatible with AA, while those of LAAA × LL played negative roles and caused its fruits aborted, and hence LAAA was incompatible with LL.

Identification of lncRNAs involved in rice ovule development and female gametophyte abortion by genome-wide screening and functional analysis

Background

As important female reproductive tissues, the rice (Oryza sativa L.) ovule and female gametophyte is significant in terms of their fertility. Long noncoding RNAs (lncRNAs) play important and wide-ranging roles in the growth and development of plants and have become a major research focus in recent years. Therefore, we explored the characterization and expression change of lncRNAs during ovule development and female gametophytic abortion.

Results

In our study, whole-transcriptome strand-specific RNA sequencing (ssRNA-seq) was performed in the ovules of a high-frequency female-sterile rice line (fsv1) and a wild-type rice line (Gui99) at the megaspore mother cell meiosis stage (stage 1), functional megaspore mitosis stage (stage 2) and female gametophyte mature stage (stage 3). By comparing two rice lines, we identified 152, 233, and 197 differentially expressed lncRNAs at the three ovule developmental stages. Functional analysis of the coherent target genes of these differentially expressed lncRNAs indicated that many lncRNAs participate in multiple pathways such as hormone and cellular metabolism and signal transduction. Moreover, there were many differentially expressed lncRNAs acting as the precursors of some miRNAs that are involved in the development of ovules and female gametophytes. In addition, we have found that lncRNAs can act as decoys, competing with mRNAs for binding to miRNAs to maintain the normal expression of genes related to ovule and female gametophyte development.

Conclusion

These results provide important clues for elucidating the female gametophyte abortion mechanism in rice. This study also expands our understanding about the biological functions of lncRNAs and the annotation of the rice genome.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5442-6) contains supplementary material, which is available to authorized users.

Genome-wide analysis of DNA polymorphisms, the methylome and transcriptome revealed that multiple factors are associated with low pollen fertility in autotetraploid rice

Autotetraploid rice is a useful germplasm with high biomass production; however, low fertility is the main barrier in commercial utilization. In our previous study, differential expression of meiosis-related miRNAs was found to be involved in the pollen sterility of autotetraploid rice. However, genome-wide DNA variations and methylomes associated with low fertility of autotetraploid rice are still poorly understood. Here, we measured both global DNA variations and the methylome and compared them with the transcriptome during pollen development in autotetraploid rice by high-throughput sequencing. A total of 34416 SNPs, 6993 InDels, 1003 SVs and 25 CNVs were detected, and 11367 and 41117 differentially methylated regions showed hypermethylation and hypomethylation in 02428-4x. In total, 1110 genes displayed differentially expression in 02428-4x during meiosis, of these six harbored CNVs, including four upregulated genes with gain CNVs, such as LOC_Os11g38620. We identified 122 genes by comparing with the previous data that might be associated with low fertility during pollen development in 02428-4x. Of the 122 gens, 98 were displayed methylation and differential expression, including OsMADS98, CYP703A3 and OsABCG26. The downregulation of these three genes were confirmed by qPCR during meiosis of 02428-4x, which played pivotal roles in pollen fertility. These results indicate that the low fertility of autotetraploid rice is not only caused by the differential expression of genes involved in pollen development, but also by sequence variation and differential methylation, suggesting that the reason for pollen sterility in autotetraploid rice is complex and might be affected by multiple factors.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

Identification and expression analysis of microRNAs during ovule development in rice (Oryza sativa) by deep sequencing

MicroRNA (miRNA) expression profiles during rice ovule development revealed the possible miRNA-mediated regulation between ovule sporophytic tissue and female gametophyte and the involvement of miRNAs in programmed cell death. MiRNAs are 20-24-nucleotide small RNAs that play key roles in the regulation of many growth and developmental processes in plants. Rice ovule development comprises a series of biological events, which are regulated by complex molecular mechanisms. To gain insight into miRNA-mediated regulation of rice ovule development, Illumina sequencing was used to examine the expression of miRNAs from the megaspore mother cell meiosis stage to the fertilized ovule stage. Based on the sequencing data, 486 known and 204 novel miRNAs were identified during rice ovule development. Moreover, 56, 65 and 11 differentially expressed miRNAs between adjacent developmental stages were identified. By analyzing transcriptome and degradome data, we identified 41, 65 and 12 coherent target genes for the differentially expressed miRNAs in ovule development. We found that changes in the expression of plant hormone-related miRNAs may play important roles in embryo sac development, providing evidence for cross-talk communication between sporophytic tissue and the female gametophyte. Additionally, we revealed that miRNAs may be involved in programmed cell death after fertilization. Finally, we constructed miRNA-mediated regulatory networks that are active during rice ovule development.

Identification of putative QTLs for seedling stage phosphorus starvation response in finger millet (Eleusine coracana L. Gaertn.) by association mapping and cross species synteny analysis

A germplasm assembly of 128 finger millet genotypes from 18 countries was evaluated for seedling-stage phosphorus (P) responses by growing them in P sufficient (Psuf) and P deficient (Pdef) treatments. Majority of the genotypes showed adaptive responses to low P condition. Based on phenotype behaviour using the best linear unbiased predictors for each trait, genotypes were classified into, P responsive, low P tolerant and P non-responsive types. Based on the overall phenotype performance under Pdef, 10 genotypes were identified as low P tolerants. The low P tolerant genotypes were characterised by increased shoot and root length and increased root hair induction with longer root hairs under Pdef, than under Psuf. Association mapping of P response traits using mixed linear models revealed four quantitative trait loci (QTLs). Two QTLs (qLRDW.1 and qLRDW.2) for low P response affecting root dry weight explained over 10% phenotypic variation. In silico synteny analysis across grass genomes for these QTLs identified putative candidate genes such as Ser-Thr kinase and transcription factors such as WRKY and basic helix-loop-helix (bHLH). The QTLs for response under Psuf were mapped for traits such as shoot dry weight (qHSDW.1) and root length (qHRL.1). Putative associations of these QTLs over the syntenous regions on the grass genomes revealed proximity to cytochrome P450, phosphate transporter and pectin methylesterase inhibitor (PMEI) genes. This is the first report of the extent of phenotypic variability for P response in finger millet genotypes during seedling-stage, along with the QTLs and putative candidate genes associated with P starvation tolerance.

Dynamic Regulation of Auxin Response during Rice Development Revealed by Newly Established Hormone Biosensor Markers

The hormone auxin is critical for many plant developmental processes. Unlike the model eudicot plant Arabidopsis (Arabidopsis thaliana), auxin distribution and signaling in rice tissues has not been systematically investigated due to the absence of suitable auxin response reporters. In this study we observed the conservation of auxin signaling components between Arabidopsis and model monocot crop rice (Oryza sativa), and generated complementary types of auxin biosensor constructs, one derived from the Aux/IAA-based biosensor DII-VENUS but constitutively driven by maize ubiquitin-1 promoter, and the other termed DR5-VENUS in which a synthetic auxin-responsive promoter (DR5rev ) was used to drive expression of the yellow fluorescent protein (YFP). Using the obtained transgenic lines, we observed that during the vegetative development, accumulation of DR5-VENUS signal was at young and mature leaves, tiller buds and stem base. Notably, abundant DR5-VENUS signals were observed in the cytoplasm of cortex cells surrounding lateral root primordia (LRP) in rice. In addition, auxin maxima and dynamic re-localization were seen at the initiation sites of inflorescence and spikelet primordia including branch meristems (BMs), female and male organs. The comparison of these observations among Arabidopsis, rice and maize suggests the unique role of auxin in regulating rice lateral root emergence and reproduction. Moreover, protein localization of auxin transporters PIN1 homologs and GFP tagged OsAUX1 overlapped with DR5-VENUS during spikelet development, helping validate these auxin response reporters are reliable markers in rice. This work firstly reveals the direct correspondence between auxin distribution and rice reproductive and root development at tissue and cellular level, and provides high-resolution auxin tools to probe fundamental developmental processes in rice and to establish links between auxin, development and agronomical traits like yield or root architecture.

Analysis of small RNAs revealed differential expressions during pollen and embryo sac development 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 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.

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-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.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3526-8) contains supplementary material, which is available to authorized users.