An ATP-binding cassette transporter, OsABCB24, is involved in female gametophyte development and early seed growth in rice.

Female gametogenesis has been rarely studied due to gametophyte lethality and the unavailability of related genetic resources. In this study, we identified a rice ATP-binding cassette transporter, OsABCB24, whose null function displayed a significantly reduced seed setting rate by as much as 94%-100% compared with that of the wild type (WT). The reciprocal cross of WT and mutant plants demonstrated that the female reproductive organs in mutants were functionally impaired. Confocal microscopy observations revealed that, although megasporogenesis remained unaffected in CRISPR/Cas9 osabcb24 mutants, the formation of female gametophytes was interrupted. Additionally, the structure of the syncytial nucleus was impaired during the initial stages of endosperm formation. Histochemical analysis showed that OsABCB24 was preferentially expressed at the conjunction of receptacle and ovary, spanning from the functional megaspore stage to the two-nucleate embryo sac stage. Further, OsABCB24 was identified as an endoplasmic reticulum membrane-localized protein. Notably, the overexpression of OsABCB24 triggered a 1.5- to 2-fold increase in grain production compared to the WT. Our findings showed that OsABCB24 plays a key role in both female gametophyte development and the early development of seeds.

SlCRCa is a key D-class gene controlling ovule fate determination in tomato.

Cell fate determination and primordium initiation on the placental surface are two key events for ovule formation in seed plants, which directly affect ovule density and seed yield. Despite ovules form in the marginal meristematic tissues of the carpels, angiosperm carpels evolved after the ovules. It is not clear how the development of the ovules and carpels is coordinated in angiosperms. In this study, we identify the S. lycopersicum CRABS CLAW (CRC) homologue SlCRCa as an essential determinant of ovule fate. We find that SlCRCa is not only expressed in the placental surface and ovule primordia but also functions as a D-class gene to block carpel fate and promote ovule fate in the placental surface. Loss of function of SlCRCa causes homeotic transformation of the ovules to carpels. In addition, we find low levels of the S. lycopersicum AINTEGUMENTA (ANT) homologue (SlANT2) favour the ovule initiation, whereas high levels of SlANT2 promote placental carpelization. SlCRCa forms heterodimer with tomato INNER NO OUTER (INO) and AGAMOUS (AG) orthologues, SlINO and TOMATO AGAMOUS1 (TAG1), to repress SlANT2 expression during the ovule initiation. Our study confirms that angiosperm basal ovule cells indeed retain certain carpel properties and provides mechanistic insights into the ovule initiation.

Genome-Wide Characterization and Functional Validation of the ACS Gene Family in the Chestnut Reveals Its Regulatory Role in Ovule Development.

Ovule abortion significantly contributes to a reduction in chestnut yield. Therefore, an examination of the mechanisms underlying ovule abortion is crucial for increasing chestnut yield. In our previous study, we conducted a comprehensive multiomic analysis of fertile and abortive ovules and found that ACS genes in chestnuts (CmACS) play a crucial role in ovule development. Therefore, to further study the function of ACS genes, a total of seven CmACS members were identified, their gene structures, conserved structural domains, evolutionary trees, chromosomal localization, and promoter cis-acting elements were analyzed, and their subcellular localization was predicted and verified. The spatiotemporal specificity of the expression of the seven CmACS genes was confirmed via qRT-PCR analysis. Notably, CmACS7 was exclusively expressed in the floral organs, and its expression peaked during fertilization and decreased after fertilization. The ACC levels remained consistently greater in fertile ovules than in abortive ovules. The ACSase activity of CmACS7 was identified using the genetic transformation of chestnut healing tissue. Micro Solanum lycopersicum plants overexpressing CmACS7 had a significantly greater rate of seed failure than did wild-type plants. Our results suggest that ovule fertilization activates CmACS7 and increases ACC levels, whereas an overexpression of CmACS7 leads to an increase in ACC content in the ovule prior to fertilization, which can lead to abortion. In conclusion, the present study demonstrated that chestnut ovule abortion is caused by poor fertilization and not by nutritional competition. Optimization of the pollination and fertilization of female flowers is essential for increasing chestnut yield and reducing ovule abortion.

Understanding megasporogenesis through model plants: contemporary evidence and future insights.

The megasporangium serves as a model system for understanding the concept of individual cell identity, and cell-to-cell communication in angiosperms. As development of the ovule progresses, three distinct layers, the epidermal (L1), the subepidermal or the hypodermal (L2) and the innermost layers (L3) are formed along the MMC (megaspore mother cell). The MMC, which is the primary female germline cell, is initiated as a single subepidermal cell amongst several somatic cells. MMC development is governed by various regulatory pathways involving intercellular signaling, small RNAs and DNA methylation. The programming and reprograming of a single nucellar cell to enter meiosis is governed by 'permissive' interacting processes and factors. Concomitantly, several nucellar sister cells are prevented from germline fate also by a set of 'repressive' factors. However, in certain angiosperms, anomalies in development of the female gametophyte have been observed. The sporophytic tissue surrounding the female gametophyte affects the gametophyte in multiple ways. The role of genes and transcription factors in the development of the MMC and in the regulation of various processes studied in selected model plants such as Arabidopsis is explained in detail in this paper. However, as angiosperms display enormous diversity, it is important to investigate early stages of megasporogenesis in other plant systems as well. Such studies provide valuable insights in understanding the regulation of megasporogenesis and the evolution of the female gametophyte from gymnosperms to flowering plants.

Development of pollinated and unpollinated ovules in Ginkgo biloba: unravelling the role of pollen in ovule tissue maturation.

In gymnosperms such as Ginkgo biloba, the arrival of pollen plays a key role in ovule development, before fertilization occurs. Accordingly, G. biloba female plants geographically isolated from male plants abort all their ovules after the pollination drop emission, which is the event that allows the ovule to capture pollen grains. To decipher the mechanism induced by pollination required to avoid ovule senescence and then abortion, we compared the transcriptomes of pollinated and unpollinated ovules at three time points after the end of the emission of pollination drop. Transcriptomic and in situ expression analyses revealed that several key genes involved in programmed cell death such as senescence and apoptosis, DNA replication, and cell cycle regulation were differentially expressed in unpollinated ovules compared to pollinated ovules. We provide evidence that the pollen captured by the pollination drop affects auxin local accumulation and might cause deregulation of key genes required for the ovule's programmed cell death, activating both the cell cycle regulation and DNA replication genes.

Strigolactones Might Regulate Ovule Development after Fertilization in Xanthoceras sorbifolium.

Strigolactones (SLs) were recently defined as a novel class of plant hormones that act as key regulators of diverse developmental processes and environmental responses. Much research has focused on SL biosynthesis and signaling in roots and shoots, but little is known about whether SLs are produced in early developing seeds and about their roles in ovule development after fertilization. This study revealed that the fertilized ovules and early developing pericarp in Xanthoceras sorbifolium produced minute amounts of two strigolactones: 5-deoxystrigol and strigol. Their content decreased in the plants with the addition of exogenous phosphate (Pi) compared to those without the Pi treatment. The exogenous application of an SL analog (GR24) and a specific inhibitor of SL biosynthesis (TIS108) affected early seed development and fruit set. In the Xanthoceras genome, we identified 69 potential homologs of genes involved in SL biological synthesis and signaling. Using RNA-seq to characterize the expression of these genes in the fertilized ovules, 37 genes were found to express differently in the fertilized ovules that were aborting compared to the normally developing ovules. A transcriptome analysis also revealed that in normally developing ovules after fertilization, 12 potential invertase genes were actively expressed. Hexoses (glucose and fructose) accumulated at high concentrations in normally developing ovules during syncytial endosperm development. In contrast, a low ratio of hexose and sucrose levels was detected in aborting ovules with a high strigolactone content. XsD14 virus-induced gene silencing (VIGS) increased the hexose content in fertilized ovules and induced the proliferation of endosperm free nuclei, thereby promoting early seed development and fruit set. We propose that the crosstalk between sugar and strigolactone signals may be an important part of a system that accurately regulates the abortion of ovules after fertilization. This study is useful for understanding the mechanisms underlying ovule abortion, which will serve as a guide for genetic or chemical approaches to promote seed yield in Xanthoceras.

Maternal nitric oxide homeostasis impacts female gametophyte development under optimal and stress conditions.

In adverse environments, the number of fertilizable female gametophytes (FGs) in plants is reduced, leading to increased survival of the remaining offspring. How the maternal plant perceives internal growth cues and external stress conditions to alter FG development remains largely unknown. We report that homeostasis of the stress signaling molecule nitric oxide (NO) plays a key role in controlling FG development under both optimal and stress conditions. NO homeostasis is precisely regulated by S-nitrosoglutathione reductase (GSNOR). Prior to fertilization, GSNOR protein is exclusively accumulated in sporophytic tissues and indirectly controls FG development in Arabidopsis (Arabidopsis thaliana). In GSNOR null mutants, NO species accumulated in the degenerating sporophytic nucellus, and auxin efflux into the developing FG was restricted, which inhibited FG development, resulting in reduced fertility. Importantly, restoring GSNOR expression in maternal, but not gametophytic tissues, or increasing auxin efflux substrate significantly increased the proportion of normal FGs and fertility. Furthermore, GSNOR overexpression or added auxin efflux substrate increased fertility under drought and salt stress. These data indicate that NO homeostasis is critical to normal auxin transport and maternal control of FG development, which in turn determine seed yield. Understanding this aspect of fertility control could contribute to mediating yield loss under adverse conditions.

HISTONE DEACETYLASE19 Controls Ovule Number Determination and Transmitting Tract Differentiation.

The gynoecium is critical for the reproduction of flowering plants as it contains the ovules and the tissues that foster pollen germination, growth, and guidance. These tissues, known as the reproductive tract (ReT), comprise the stigma, style, and transmitting tract (TT). The ReT and ovules originate from the carpel margin meristem (CMM) within the pistil. SHOOT MERISTEMLESS (STM) is a key transcription factor for meristem formation and maintenance. In all above-ground meristems, including the CMM, local STM downregulation is required for organ formation. However, how this downregulation is achieved in the CMM is unknown. Here, we have studied the role of HISTONE DEACETYLASE 19 (HDA19) in Arabidopsis (Arabidopsis thaliana) during ovule and ReT differentiation based on the observation that the hda19-3 mutant displays a reduced ovule number and fails to differentiate the TT properly. Fluorescence-activated cell sorting coupled with RNA-sequencing revealed that in the CMM of hda19-3 mutants, genes promoting organ development are downregulated while meristematic markers, including STM, are upregulated. HDA19 was essential to downregulate STM in the CMM, thereby allowing ovule formation and TT differentiation. STM is ectopically expressed in hda19-3 at intermediate stages of pistil development, and its downregulation by RNA interference alleviated the hda19-3 phenotype. Chromatin immunoprecipitation assays indicated that STM is a direct target of HDA19 during pistil development and that the transcription factor SEEDSTICK is also required to regulate STM via histone acetylation. Thus, we identified factors required for the downregulation of STM in the CMM, which is necessary for organogenesis and tissue differentiation.

Small RNA profiling reveals that an ovule-specific microRNA, cja-miR5179, targets a B-class MADS-box gene in Camellia japonica.

BACKGROUND AND AIMS: The functional specialization of microRNA and its target genes is often an important factor in the establishment of spatiotemporal patterns of gene expression that are essential to plant development and growth. In different plant lineages, understanding the functional conservation and divergence of microRNAs remains to be explored. METHODS: To identify small regulatory RNAs underlying floral patterning, we performed a tissue-specific profiling of small RNAs in various floral organs from single and double flower varieties (flowers characterized by multiple layers of petals) in Camellia japonica. We identified cja-miR5179, which belongs to a deeply conserved microRNA family that is conserved between angiosperms and basal plants but frequently lost in eudicots. We characterized the molecular function of cja-miR5179 and its target - a B-function MADS-box gene - through gene expression analysis and transient expression assays. KEY RESULTS: We showed that cja-miR5179 is exclusively expressed in ovule tissues at the early stage of floral development. We found that cja-miR5179 targets the coding sequences of a DEFICIENS-like B-class gene (CjDEF) mRNA, which is located in the K motif of the MADS-box domain; and the target sites of miR5179/MADS-box were consistent in Camellia and orchids. Furthermore, through a petal transient-expression assay, we showed that the BASIC PENTACYSTEINE proteins bind to the GA-rich motifs in the cja-miR5179 promoter region and suppresses its expression. CONCLUSIONS: We propose that the regulation between miR5179 and a B-class MADS-box gene in C. japonica has a deep evolutionary origin before the separation of monocots and dicots. During floral development of C. japonica, cja-miR5179 is specifically expressed in the ovule, which may be required for the inhibition of CjDEF function. This work highlights the evolutionary conservation as well as functional divergence of small RNAs in floral development.

Arabidopsis TCP4 transcription factor inhibits high temperature-induced homeotic conversion of ovules.

Abnormal high temperature (HT) caused by global warming threatens plant survival and food security, but the effects of HT on plant organ identity are elusive. Here, we show that Class II TEOSINTE BRANCHED 1/CYCLOIDEA/ PCF (TCP) transcription factors redundantly protect ovule identity under HT. The duodecuple tcp2/3/4/5/10/13/17/24/1/12/18/16 (tcpDUO) mutant displays HT-induced ovule conversion into carpelloid structures. Expression of TCP4 in tcpDUO complements the ovule identity conversion. TCP4 interacts with AGAMOUS (AG), SEPALLATA3 (SEP3), and the homeodomain transcription factor BELL1 (BEL1) to strengthen the association of BEL1 with AG-SEP3. The tcpDUO mutant synergistically interacts with bel1 and the ovule identity gene seedstick (STK) mutant stk in tcpDUO bel1 and tcpDUO stk. Our findings reveal the critical roles of Class II TCPs in maintaining ovule identity under HT and shed light on the molecular mechanisms by which ovule identity is determined by the integration of internal factors and environmental temperature.

Genetic and Phenotypic Analysis of Ovule Development in Arabidopsis.

The plant seed is a remarkable structure that represents the single most important energy source in global diets. The stages of reproductive growth preceding seed formation are particularly important since they influence the number, size, and quality of seed produced. The progenitor of the seed is the ovule, a multicellular organ that produces a female gametophyte while maintaining a range of somatic ovule cells to protect the seed and ensure it receives maternal nourishment. Ovule development has been well characterized in Arabidopsis using a range of molecular, genetic, and cytological assays. These can provide insight into the mechanistic basis for ovule development, and opportunities to explore its evolutionary conservation. In this chapter, we describe some of these methods and tools that can be used to investigate early ovule development and cell differentiation.

Genome-wide chromatin accessibility landscape and dynamics of transcription factor networks during ovule and fiber development in cotton.

BACKGROUND: The development of cotton fiber is regulated by the orchestrated binding of regulatory proteins to cis-regulatory elements associated with developmental genes. The cis-trans regulatory dynamics occurred throughout the course of cotton fiber development are elusive. Here we generated genome-wide high-resolution DNase I hypersensitive sites (DHSs) maps to understand the regulatory mechanisms of cotton ovule and fiber development. RESULTS: We generated DNase I hypersensitive site (DHS) profiles from cotton ovules at 0 and 3 days post anthesis (DPA) and fibers at 8, 12, 15, and 18 DPA. We obtained a total of 1185 million reads and identified a total of 199,351 DHSs through ~ 30% unique mapping reads. It should be noted that more than half of DNase-seq reads mapped multiple genome locations and were not analyzed in order to achieve a high specificity of peak profile and to avoid bias from repetitive genomic regions. Distinct chromatin accessibilities were observed in the ovules (0 and 3 DPA) compared to the fiber elongation stages (8, 12, 15, and 18 DPA). Besides, the chromatin accessibility during ovules was particularly elevated in genomic regions enriched with transposable elements (TEs) and genes in TE-enriched regions were involved in ovule cell division. We analyzed cis-regulatory modules and revealed the influence of hormones on fiber development from the regulatory divergence of transcription factor (TF) motifs. Finally, we constructed a reliable regulatory network of TFs related to ovule and fiber development based on chromatin accessibility and gene co-expression network. From this network, we discovered a novel TF, WRKY46, which may shape fiber development by regulating the lignin content. CONCLUSIONS: Our results not only reveal the contribution of TEs in fiber development, but also predict and validate the TFs related to fiber development, which will benefit the research of cotton fiber molecular breeding.

RLI2 regulates Arabidopsis female gametophyte and embryo development by facilitating the assembly of the translational machinery.

Eukaryotic protein translation is a complex process that requires the participation of different proteins. Defects in the translational machinery often result in embryonic lethality or severe growth defects. Here, we report that RNase L inhibitor 2/ATP-BINDING CASSETTE E2 (RLI2/ABCE2) regulates translation in Arabidopsis thaliana. Null mutation of rli2 is gametophytic and embryonic lethal, whereas knockdown of RLI2 causes pleiotropic developmental defects. RLI2 interacts with several translation-related factors. Knockdown of RLI2 affects the translational efficiency of a subset of proteins involved in translation regulation and embryo development, indicating that RLI2 has critical roles in these processes. In particular, RLI2 knockdown mutant exhibits decreased expression of genes involved in auxin signaling and female gametophyte and embryo development. Therefore, our results reveal that RLI2 facilitates assembly of the translational machinery and indirectly modulates auxin signaling to regulate plant growth and development.

The molecular chaperone mtHSC70-1 interacts with DjA30 to regulate female gametophyte development and fertility in Arabidopsis.

Arabidopsis mitochondria-targeted heat shock protein 70 (mtHSC70-1) plays important roles in the establishment of cytochrome c oxidase-dependent respiration and redox homeostasis during the vegetative growth of plants. Here, we report that knocking out the mtHSC70-1 gene led to a decrease in plant fertility; the fertility defect of the mutant was completely rescued by introducing the mtHSC70-1 gene. mtHSC70-1 mutants also showed defects in female gametophyte (FG) development, including delayed mitosis, abnormal nuclear position, and ectopic gene expression in the embryo sacs. In addition, we found that an Arabidopsis mitochondrial J-protein gene (DjA30) mutant, j30+/- , had defects in FG development and fertility similar to those of mtHSC70-1 mutant. mtHSC70-1 and DjA30 had similar expression patterns in FGs and interacted in vivo, suggesting that these two proteins might cooperate during female gametogenesis. Further, respiratory chain complex IV activity in mtHSC70-1 and DjA30 mutant embryo sacs was markedly downregulated; this led to the accumulation of mitochondrial reactive oxygen species (ROS). Scavenging excess ROS by introducing Mn-superoxide dismutase 1 or catalase 1 gene into the mtHSC70-1 mutant rescued FG development and fertility. Altogether, our results suggest that mtHSC70-1 and DjA30 are essential for the maintenance of ROS homeostasis in the embryo sacs and provide direct evidence for the roles of ROS homeostasis in embryo sac maturation and nuclear patterning, which might determine the fate of gametic and accessory cells.

Embryo sac development relies on symplastic signals from ovular integuments in Arabidopsis.

Ovules are female reproductive organs of angiosperms, consisting of sporophytic integuments surrounding female gametophytes, that is, embryo sacs. Synchronization between integument growth and embryo sac development requires intracellular communication. However, signaling routes through which cells of the two generations communicate are unclear. We report that symplastic signals through plasmodesmata (PDs) of integuments are critical for the development of female gametophytes. Genetic interferences of PD biogenesis either by functional loss of CHOLINE TRANSPORTER-LIKE1 (CTL1) or by integument-specific expression of a mutated CALLOSE SYNTHASE 3 (cals3m) compromised PD formation in integuments and reduced fertility. Close examination of pINO:cals3m or ctl1 ovules indicated that female gametophytic development was either arrested at various stages after the formation of functional megaspores. In both cases, defective ovules could not attract pollen tubes, leading to the failure of fertilization. Results presented here demonstrate a key role of the symplastic route in sporophytic control of female gametophytic development.

Control of cellularization, nuclear localization, and antipodal cell cluster development in maize embryo sacs.

The maize female gametophyte contains four cell types: two synergids, an egg cell, a central cell, and a variable number of antipodal cells. In maize, these cells are produced after three rounds of free-nuclear divisions followed by cellularization, differentiation, and proliferation of the antipodal cells. Cellularization of the eight-nucleate syncytium produces seven cells with two polar nuclei in the central cell. Nuclear localization is tightly controlled in the embryo sac. This leads to precise allocation of the nuclei into the cells upon cellularization. Nuclear positioning within the syncytium is highly correlated with their identity after cellularization. Two mutants are described with extra polar nuclei, abnormal antipodal cell morphology, and reduced antipodal cell number, as well as frequent loss of antipodal cell marker expression. Mutations in one of these genes, indeterminate gametophyte2 encoding a MICROTUBULE ASSOCIATED PROTEIN65-3 homolog, shows a requirement for MAP65-3 in cellularization of the syncytial embryo sac as well as for normal seed development. The timing of the effects of ig2 suggests that the identity of the nuclei in the syncytial female gametophyte can be changed very late before cellularization.

OAF is a DAF-like gene that controls ovule development in plants.

We previously found that the RING-type E3 ligase DEFECTIVE IN ANTHER DEHISCENCE1- (DAD1-) Activating Factor (DAF) controls anther dehiscence by activating the jasmonate biosynthetic pathway in Arabidopsis. Here, we show that in Arabidopsis, the DAF ancestor was duplicated into three genes (DAF, Ovule Activating Factor (OAF), DAFL2), which evolved divergent partial functions from their ancestor through subfunctionalization. In this case, DAF-DAD1-JA signaling regulates anther dehiscence, whereas OAF controls ovule development by negatively regulating cinnamyl alcohol dehydrogenase 9 (CAD9) activity and being negatively regulated by miR847 itself in Arabidopsis. Downregulation of OAF or upregulation of CAD9 and miR847 caused similar abortion of ovule formation due to precocious ovule lignification in transgenic Arabidopsis. Interestingly, only one DAF-like gene, PaOAF, exists in the monocot orchids, which has likely evolved through nonfunctionalization and maintains a conserved function as Arabidopsis OAF in regulating ovule development since defective ovules were observed in the virus-induced gene silencing (VIGS) PaOAF Phalaenopsis orchids. The absence of the DAF ortholog and its function in orchids is likely due to the evolution of stamens to a unique pollinium structure that lacks the feature of anther dehiscence. These findings expand the current knowledge underlying the multifunctional evolution and diverse functionalization of duplicate gene pairs within/among plants.

Transcriptomic Analysis of Hormone Signal Transduction, Carbohydrate Metabolism, Heat Shock Proteins, and SCF Complexes before and after Fertilization of Korean Pine Ovules.

The fertilization process is a critical step in plant reproduction. However, the mechanism of action and mode of regulation of the fertilization process in gymnosperms remain unclear. In this study, we investigated the molecular regulatory networks involved in the fertilization process in Korean pine ovules through anatomical observation, physiological and biochemical assays, and transcriptome sequencing technology. The morphological and physiological results indicated that fertilization proceeds through the demise of the proteinaceous vacuole, egg cell division, and pollen tube elongation. Auxin, cytokinin, soluble sugar, and soluble starch contents begin to decline upon fertilization. Transcriptomic data analysis revealed a large number of differentially expressed genes at different times before and after fertilization. These genes were primarily involved in pathways associated with plant hormone signal transduction, protein processing in the endoplasmic reticulum, fructose metabolism, and mannose metabolism. The expression levels of several key genes were further confirmed by qRT-PCR. These findings represent an important step towards understanding the mechanisms underlying morphological changes in the Korean pine ovule during fertilization, and the physiological and transcriptional analyses lay a foundation for in-depth studies of the molecular regulatory network of the Korean pine fertilization process.

Single-nucleus sequencing deciphers developmental trajectories in rice pistils.

Angiosperms possess a life cycle with an alternation of sporophyte and gametophyte generations, which happens in plant organs like pistils. Rice pistils contain ovules and receive pollen for successful fertilization to produce grains. The cellular expression profile in rice pistils is largely unknown. Here, we show a cell census of rice pistils before fertilization through the use of droplet-based single-nucleus RNA sequencing. The ab initio marker identification validated by in situ hybridization assists with cell-type annotation, revealing cell heterogeneity between ovule- and carpel-originated cells. A comparison of 1N (gametophyte) and 2N (sporophyte) nuclei identifies the developmental path of germ cells in ovules with typical resetting of pluripotency before the sporophyte-gametophyte transition, while trajectory analysis of carpel-originated cells suggests previously neglected features of epidermis specification and style function. These findings gain a systems-level view of cellular differentiation and development of rice pistils before flowering and lay a foundation for understanding female reproductive development in plants.

Ovule Transcriptome Analysis Discloses Deregulation of Genes and Pathways in Sexual and Apomictic Limonium Species (Plumbaginaceae).

The genus Limonium Mill. (sea lavenders) includes species with sexual and apomixis reproductive strategies, although the genes involved in these processes are unknown. To explore the mechanisms beyond these reproduction modes, transcriptome profiling of sexual, male sterile, and facultative apomictic species was carried out using ovules from different developmental stages. In total, 15,166 unigenes were found to be differentially expressed with apomictic vs. sexual reproduction, of which 4275 were uniquely annotated using an Arabidopsis thaliana database, with different regulations according to each stage and/or species compared. Gene ontology (GO) enrichment analysis indicated that genes related to tubulin, actin, the ubiquitin degradation process, reactive oxygen species scavenging, hormone signaling such as the ethylene signaling pathway and gibberellic acid-dependent signal, and transcription factors were found among differentially expressed genes (DEGs) between apomictic and sexual plants. We found that 24% of uniquely annotated DEGs were likely to be implicated in flower development, male sterility, pollen formation, pollen-stigma interactions, and pollen tube formation. The present study identifies candidate genes that are highly associated with distinct reproductive modes and sheds light on the molecular mechanisms of apomixis expression in Limonium sp.