Identification and analysis of the stigma and embryo sac-preferential/specific genes in rice pistils
© The Author(s). 2017
Received: 24 August 2016
Accepted: 23 February 2017
Published: 7 March 2017
In rice, the pistil is the female reproductive organ, and it consists of two stigmas and an ovary. The stigma is capable of receiving pollen grains and guiding pollen tube growth. The ovary holds the embryo sac, which is fertilized with male gametes to produce seed. However, little is known about the gene function and regulatory networks during these processes in rice.
Here, using the RNA-Seq technique, we identified 3531 stigma-preferential genes and 703 stigma-specific genes within the rice pistils, and we verified 13 stigma-specific genes via qRT-PCR and in situ hybridization. The GO analysis showed that the transport-, localization-, membrane-, communication-, and pollination-related genes were significantly enriched in the stigma. Additionally, to identify the embryo sac-preferential/specific genes within the pistils, we compared a wild-type ovary with a mutant dst (defective stigma) ovary and found that 385 genes were down-regulated in dst. Among these genes, 122 exhibited an ovary-specific expression pattern and are thought to be embryo sac-preferential/specific genes within the pistils. Most of them were preferentially expressed, while 14 of them were specifically expressed in the pistil. Moreover, the rice homologs of some Arabidopsis embryo sac-specific genes, which played essential roles during sexual reproduction, were down-regulated in dst. Additionally, we identified 102 novel protein-coding genes, and 6 of them exhibited differences between the stigma and ovary in rice as determined using RT-PCR.
According to these rice ovary comparisons, numerous genes were preferentially or specifically expressed in the stigma, suggesting that they were involved in stigma development or pollination. The GO analysis indicated that a dry rice stigma might primarily perform its function through the cell membrane, which was different from the wet stigma of other species. Moreover, many embryo sac-preferential/specific genes within the pistils were identified and may be expressed in female rice gametophytes, implying that these genes might participate in the process of female gametophyte specialization and fertilization. Therefore, we provide the gene information for investigating the gene function and regulatory networks during female gametophyte development and fertilization. In addition, these novel genes are valuable for the supplementation and perfection of the existing transcriptome in rice, which provides an effective method of detecting novel rice genes.
KeywordsEmbryo sac female gametophyte novel gene rice RNA-Seq Stigma
In angiosperms, the life cycle consists of a haploid gametophyte stage and a diploid sporophyte stage. The male gametophyte (pollen grain) is composed of two sperm cells and one vegetative cell, whereas the female gametophyte (mature embryo sac) is composed of one egg cell, two synergid cells, one central cell with two nuclei, and three or more antipodal cells [1–4]. After the maturation of the gametophytes, pollen grains are transferred into the stigma, and then they emit pollen tubes and grow into the transmitting tissue. After that, one pollen tube grows into one of the synergid cells and releases two sperm cells. Ultimately, one sperm cell fuses with the egg cell to form the zygote, while the other fuses with the central cell to form the primary endosperm nucleus; these fusions then lead to the production of the embryo and endosperm, respectively . During this process, communication between the pollen grain and pistil occur closely and frequently to maximize reproductive success.
The pollen grain-stigma interaction occurs in the stigma, which includes two types, namely the wet stigma (e.g., in Nicotiana tabacum and Lilium longiflorum) and the dry stigma (e.g., in Arabidopsis thaliana, Zea mays and Oryza sativa). When pollen grains are deposited on papilla cells, many stigma-preferential/-specific genes begin to be expressed and play vital roles. The abnormal expression or absence of these genes or secondary metabolites would disrupt the development or function maintenance of the stigmas. In Brassica, a receptor-like kinase (RLK) called the S-receptor kinase (SRK) binds the ligand, SCR/SP11, to phosphorylate ARC1 (armadillo repeat-containing protein1), and then ARC1 mediates the degradation of the exocyst subunit genes, such as Exo70A1, which results in the recognized failure of compatible pollen grains [6–11]. Similarly, the other two stigma-specific glycoproteins, SLG (S-locus glycoprotein) and SLRI (S-locus-related 1 protein), disturb the adhesion of incompatible pollen grains via interactions with pollen coat proteins in Brassica [12–15]. In tobacco and petunias, STIG1 (stigma-specific gene 1) is involved in the regulation of pistil exudates to affect the growth of pollen tubes [16, 17]. Additionally, many secondary metabolites, such as sulfinylated azadecalin, are also essential for the germination and growth of pollen grains .
After the pollen grains germinate, the pollen tubes grow into the transmitting tissues toward the embryo sac. Signal molecules from the pistil function in prerequisite roles for pollen tube growth, especially the guidance and termination of the pollen tubes [19–21]. Peptides and glycoproteins are star molecules, which are necessary for the growth and guidance of pollen tubes . In Arabidopsis, CLV3/ESR-related peptide (CLE45) and vacuolar protein sorting 41 (AtVPS41) promote the growth of pollen tubes, while pollen tubes germinate and grow normally in ccg, myb98 and feronia, but they cannot reach the synergid cell in the first two mutants and overgrow after arriving at the synergid cell in feronia [23–27]. In Nicotiana, arabinogalactan proteins called TTS are involved in the growth of pollen tubes and may function as a Ca2+ flux capacitor in the extracellular matrix [28–30]. In Torenia fournieri, Torenia concolor, and Brassica, the LURE peptides from synergid cells are similar defensin-like CRPs, and they play an essential role in pollen tube attraction [31–33]. In monocots, ZmES4 induces rapid pollen tube rupture in maize, and OsPTB1 (POLLEN TUBE BLOCKED 1) promotes pollen tube growth to control the rice seed-setting rate [34, 35].
Rice is a model monocotyledon and one of the most important crops. OsIG1, OsMADS13, and OsLOG play essential roles during the biogenesis of the rice ovule [36–39]. To date, approximately 30 genes have been found to participate in the pollen meiosis process, and they are likely involved in the development of female gametophytes in rice [40, 41]. However, all of them are involved in the process before the cell differentiation of the female gametophyte. Although the killer-protector system at the rice S5 locus and OsDEES1 has provided insights into embryo sac development, the mechanism of cell differentiation and the unique features of different cells in the embryo sac are still unclear [42–44]. Additionally, few researchers have reported on the development mechanism of the stigmas.
In Arabidopsis, a series of female gametophyte-specific genes were identified via a comparison of the ovules with the embryo sac and without it, with the aim of dissecting the genes and the regulatory networks of certain biologic processes. Of these genes, three transcription factors, namely, MYB98, MYB64, and AGL61, have been confirmed to play essential roles in the development of female gametophytes [23, 45–47]. Similarly, in maize, the egg cell-specific gene ZmEA1 and the synergid-specific gene ZmES4 are required for completing fertilization [34, 48]. In rice using microarrays and RNA sequencing (RNA-Seq), only two female gametophyte-specific genes, ECAGL1 (LOC_Os03g18530) and a SCP-like gene (LOC_Os04g22220), in addition to some stigma-preferential genes, were identified [49–51]. However, to date, no stigma-specific gene has been reported.
In this study, to identify rice genes that are specifically expressed in the stigma and embryo sac, we adopted RNA-Seq and qPCR techniques to analyze the gene expression in the pistils, ovaries, and stigmas of the rice variety “Hwayoung” and in ovaries from the dst mutant. We obtained a series of stigma-preferential or -specific genes within the rice pistils and verified their expression via qRT-PCR and in situ hybridization. Furthermore, the GO analysis showed that the transport-, localization-, membrane-, communication-, and pollination-related genes were significantly enriched in the stigma. In addition, we found that many ovary-specific genes were down-regulated in the mutant dst compared to the HY ovary, and we believed that they were embryo sac-preferential/specific genes within the pistils. Most of them were preferentially/specifically expressed in the pistil, implying that they play an essential role during female gametophyte development and fertilization in rice. Additionally, we identified many novel protein-coding genes and verified them by using RT-PCR, displaying their value in the supplementation and perfection of the existing transcriptome in rice and providing an effective method to detect novel rice genes.
Strategy to identify stigma-specific and embryo sac-preferential/specific genes within the pistils
Based on previous studies, we found that ECAGL1 (LOC_Os03g18530) is specifically expressed in the egg cell, an SCP-like gene (LOC_Os04g22220) is a synergy cell-specific gene, and OsMADS13 (LOC_Os12g10540) is preferentially expressed in the ovule [36, 50]. In our study, ECAGL1 and LOC_Os04g22220 exhibited ovary-specific expression patterns, and their expressions in the HY ovary were found to be more than 2 times that in the dst ovary (Additional file 1: Table S1). However, OsMADS13 was abundantly expressed in the HY and dst ovaries with no significant difference (Additional file 1: Table S1). This finding provided a validation of our strategies and indicated that these embryo sac-preferential/specific genes may exhibit embryo sac-specific expression patterns.
Global analysis of RNA-Seq data
Mapping the results of the RNA-Seq reads
To assess the gene expression in different tissues, the FPKM (fragments per kilobase of gene per million) of each gene were calculated as the normalized read counts. The total length of the exon region was considered to be the length of the gene, and genes with FPKM values lower than 0.5 were not considered to be expressed genes . To evaluate the threshold, 5 genes with no expression in the rice pistil or panicle were examined [54–56]. Most genes were not detected in the 6 HY libraries, and the highest expression was 0.37 (Additional file 2: Table S2). All of the above results suggest that the threshold is useful for distinguishing among genes that are with or without expression.
Identification and analysis of stigma-/ovary-preferential genes within the pistils
To find the genes that are preferentially expressed in the ovary and stigma, we compared the HY ovary and HY stigma. The results showed that 3793 and 3531 genes were highly expressed in the ovary and the stigma, respectively (Additional file 4: Figure S2A and Additional file 5: Table S3 and Additional file 6: Table S4). To obtain further insight regarding the functional differences in tissue-preferential genes, the WEGO (Web Gene Ontology Annotation Plot) online tool was used to classify the overrepresented GO (Gene Ontology) terms .
Based on the biological process (BP), seventeen GO terms were overrepresented in the stigma. Among them, there were two GO terms associated with the location, which implied that the gradient or polarity distribution was widespread in the stigma and was necessary for sexual reproduction. Additionally, “pollination”, the four terms of response to stimulus or stress, “cell communication”, and “secondary metabolic process” were all overrepresented in the stigma. The former two terms implied that the function of the stigma was to respond to a stimulus (including pollen grains) and recognize compatible pollen grains. The latter two terms implied that communication and secondary metabolites may play essential roles in the stigma-development and pollination. Some evidence was found in other species. When pollen grains adhered to the stigma, communication developed, and recognition factors, such as the SLG and SLR1 proteins in Brassica, distinguished among the compatible and incompatible pollen grains; however, the nature of the recognition system was not clear [12, 14]. In addition, the stigma with abnormal exudates also affected the pollination, as in the STIG1 knockdown or knockout plant in tobacco and petunias [16, 17]. All of the above results implied that secretion might perform an important function in the stigma, while many genes, such as Rab GTPases, SNARE complexes and exocyst subunits, played essential roles in secretion and maintaining the normal function of the stigma [11, 59]. However, there was no GO annotation of secretion in rice, and thus we downloaded all the Arabidopsis protein sequences under the GO term “secretion” (GO: 0046903), used the BLASTP tool with an E-value of 0.00001, and then found 647 genes that were related to secretion (Additional file 7 : Table S5). Among these genes, there were 40, 28 and 55 genes coding proteins homologous to Rab GTPases, SNARE complexes and exocyst subunits, respectively [11, 59, 60]. Interestingly, the secretion-related genes, including Rab GTPases, SNARE complexes and exocyst subunits, were overrepresented in the stigma (Additional file 7 : Table S5). By contrast, the GO terms related to the biosynthetic process, metabolic process, cellular component synthesis, and organization were enriched in the ovary, implying that they played roles in the ovary and embryo sac (Fig. 4b).
Moreover, the analysis of the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways provided a confirmation of the GO analysis. The stigma-enriched pathways were primarily associated with basal and secondary metabolism, which may synthesize nutrients and signal molecules that functioned in the recognition system or guidance system of pollen grains, and the pathways were also associated with the “phagosome” and “plant hormone signal transduction” (Additional file 8: Table S6). By contrast, the enriched pathways in the ovary were primarily divided into two categories. One category was photosynthesis and the other was the cell cycle, and the cell cycle was separated into three classes (DNA metabolism, RNA metabolism, and protein synthesis) (Additional file 9: Table S7).
Identification of stigma-/ovary-specific genes within the pistil
To determine the stigma-specific genes within the pistil, we selected stigma-preferential genes, the expression of which in any replicate of the ovary was lower than 0.5 and more than 0.5 in any replicate of the stigma, and we identified 703 stigma-specific genes (Additional file 10: Table S8). Using a similar method, 1257 ovary-specific genes were identified (Additional file 11: Table S9). Based on the biological process, the stigma-specific genes were enriched in 18 GO terms, and the top two results included the “metabolic process” and “pollination”, whereas the ovary-specific genes were overrepresented in 24 GO terms. The top two findings included the “DNA metabolic process” and the “cell cycle” (Additional file 12: Figure S3A-B).
Because there is no stigma in the pistil of the dst mutant (also called the dst ovary), the stigma-specific genes may be theoretically down-regulated in the dst pistil (ovary) compared to the HY pistil. To verify the stigma-specific genes, we compared the HY pistil and the dst pistil. The results showed that 393 genes were down-regulated in the dst ovary compared with the HY pistil (Additional file 4: Figure S2B). Among these candidate genes, we determined only 37 stigma-specific genes within the pistil (Fig. 5b, Additional file 10: Table S8).
Identification of ovary-/stigma-specific genes during the rice life cycle
During the sexual reproduction process, many genes are expressed in a certain tissue or in a certain type of cells at a certain time. If the temporal and spatial expression pattern is changed, then an abortion or an abnormal seed appears. Many research studies, on genes such as ZmEA1, have provided validation for these findings .
Thus, to identify ovary/stigma-specific genes during the life cycle of rice, we downloaded the transcriptome expression data for rice genes from the RGAP (http://rice.plantbiology.msu.edu/expression.shtml). We chose expression data from tissues including a 20-day-old leaf, a 14-day-old shoot, a seedling at the four-leaf stage, an anther, a pistil before pollination, a seed at 5 days after pollination, a seed at 10 days after pollination, an embryo at 25 days after pollination and endosperm at 25 days after pollination. The threshold at which genes were considered to be “expressed” in a tissue was when their FPKM values were not zero in that tissue. We found 159 pistil-specific genes (PS) and 702 pistil-expressed genes that were not expressed in vegetative tissues (PNV) (Fig. 5c). Among the pistil-specific genes, 7 were specifically expressed in the stigma, whereas 45 genes were specifically expressed in the ovary. Additionally, we found 41 stigma-specific genes and 75 ovary-specific genes that were only expressed in reproductive tissues (Fig. 5c, Additional file 14: Table S11).
For the 655 stigma-specific genes that were not included in the PNV, we investigated all their expression in all the tissues including the leaf, shoot, seedling, panicle, anther, pistil and seed (Fig. 5c). The heat map showed that most of the genes were preferentially expressed in the anthers and panicles after they emerged from the leaf sheath (Additional file 15: Figure S4).
Identification of embryo sac-preferential/specific genes within the pistil
The embryo sac consists of four types of cells that are embedded in the ovary, and it is difficult to investigate. To obtain the candidate genes that are specifically expressed in the embryo sac, we compared the dst ovary without the embryo sac to the HY ovary, and we obtained 385 down-regulated genes in the dst ovary (Additional file 4: Figure S2C, Additional file 16: Table S12). In 2013, Anderson et al. isolated egg cells, sperm cells and vegetative cells from rice and detected their gene expression by using deep-sequencing technology. They found that 27417 genes were expressed in the egg cells (the threshold, or the point at which the gene that was believed to be expressed in the egg cell, is that when the TPM values were not zero in any biological replicate of the egg cell) . In our study, among these down-regulated genes in the dst ovary compared with the HY ovary, 371 (96.4%) genes were expressed in egg cells (Fig. 5d, Additional file 12: Table S12), and the remaining genes may be expressed in other cells of the embryo sac.
As mentioned above, the genes for the embryo sac (female gametophyte)-specific expression belonged to two types of gene sets, that is, the ovary-specific genes and down-regulated genes in the dst ovary compared with the HY ovary. Therefore, using Venny analysis, we obtained 122 ovary-specific genes that were down-regulated in the dst ovary, suggesting that they may be embryo sac-specific expression genes in the pistil (Fig. 5e, Additional file 17: Table S13). However, we did not know the exact and detailed expression pattern of the genes in the pistil. Thus, we called them genes “embryo sac-preferential/specific genes”. Based on the biological process, these genes were overrepresented in 11 terms, and the top three included the “cell cycle”, “cellular component organization”, and the “DNA metabolic process” (Additional file 12: Figure S3C).
The expression pattern of embryo sac-preferential/specific genes
Homologs of embryo sac-specific genes that function in the pistil
The female gametophyte is indispensable for sexual reproduction in plants. Many embryo sac-specific genes have been isolated from Arabidopsis, and some play essential roles . Theoretically, the rice homologs of these genes may be expressed to a high degree in the embryo sac and down-regulated in the ovary of the mutant dst compared with the HY ovary if they have similar functions. Because many genes are conserved in different types of plants, a number of known homologous genes have been identified, and they exhibit a similar expression pattern and biological function. Therefore, we identified the rice homologs of Arabidopsis functional genes and found that some were down-regulated in dst.
Female gametophyte development is regulated by a complicated gene network that primarily includes transcription factors and their targets. Thus, it is important to identify transcription factors and reveal their functions in female gametophyte development. In Arabidopsis, six known transcription factors (AtAGL23/61/80 and AtMYB64/98/119) are specifically expressed in the embryo sac, and they are involved in the development of the female gametophyte [23, 46, 47, 65, 66]. Among the twelve homologs of AtAGL80, six genes were preferentially expressed in the HY ovary and down-regulated in dst, but the other six were not expressed in all the tissues (Fig. 7c). In addition, in Arabidopsis, there were three MYB transcription factors known as AtMYB64/98/119, which were clustered together (Additional file 19: Figure S5). We identified twenty homologous MYB genes in rice, and LOC_Os11g10130 was significantly decreased in dst and specifically expressed in the HY pistil (Fig. 7a, d), suggesting that the gene was embryo sac-preferential/specific. Additionally, the expression of the other five homologs (LOC_Os01g16810, LOC_Os06g06740, LOC_Os12g13570, LOC_Os01g12860 and LOC_Os05g38460) was lower in dst than in the HY ovary (Fig. 7d). Moreover, the three MYB transcription factors (LOC_Os02g36890, LOC_Os04g38740 and LOC_Os03g38210) were down-regulated in dst compared with the HY ovary (Additional file 19: Figure S5, Additional file 17: Table S13).
Confirmation via qRT-PCR
Confirmation using mRNA in situ hybridization
Novel protein-coding genes
The novel protein-coding genes in rice
Number of transcripts
Number of genes
100 ≥ the supported reads >0
1000 ≥ the supported reads >100
2000 ≥ the supported reads >1000
the supported reads >2000
Furthermore, a t-test was performed between the HY stigma and HY ovary. The results showed that 2 novel genes were preferentially expressed in the stigma, whereas there were 3 novel genes in the ovary of those genes with 1000 or more supported reads (P < 0.05; Additional file 21: Table S15), except XLOC_015964. Moreover, to determine whether there was a difference in the expression of some genes between varieties, we tested the RNA levels from the Nipponbare pistil, HY ovaries, and HY stigmas using semi-quantitative RT-PCR (Fig. 11f). All six novel genes were present in both Nipponbare and Hwayoung, suggesting that they are widely present in different varieties of japonica. XLOC-015964 was preferentially expressed in the stigma, whereas the other three genes (XLOC-017980, XLOC-003666 and XLOC-011975) were expressed in the ovaries (Fig. 11f). The expression of the other two genes (XLOC-018698 and XLOC-008452) exhibited no significant differences between the rice ovary and stigma (Fig. 11f). These novel genes could be valuable for the supplementation and perfection of the existing rice transcriptome.
Identification of stigma-preferential/specific genes within the pistils
Despite the existence of many research reports about rice reproduction, little is known about the molecular mechanisms of the stigma-development and pollination. Thus, to dissect the molecular basis of these mechanisms in rice, we performed RNA-Seq, and we identified 3531 stigma-preferential genes. Compared with the ovary, the stigma exhibited a distinct genome-wide expression pattern (Fig. 4). The transport-, localization-, membrane-, communication-, secretion- and pollination-related genes were significantly enriched in the stigma, while genes related to the extracellular region were not (Fig. 4a, Additional file 7 : Table S5). The results showed that the rice stigma, whose transmitting tracts were made up of several compactly arranged cell layers around the vascular bundle, might perform its function primarily through the membrane, and it was different from the wet stigma in tobacco, which is composed of transmitting tissue and papilla cells with vigorous secretion activity [28, 68]. Thus, the wet stigma primarily performed its role through the extracellular matrix. In tobacco, floral-transmitting tissue-specific arabinogalactan proteins [transmitting tissue specific proteins (TTS)] were located in the extracellular matrix of the transmitting tissue. TTS attracted the pollen tubes, stimulated their growth, and were then deglycosylated to form a glycosylation gradient in the stigma and style [28, 69]. Conversely, dry stigmas disturbed the adhesion of incompatible pollen grains primarily via membrane-related receptors in the stigma. For example, Brassica recognized the incompatible pollen grains primarily through SRK and its downstream pathways . Based on the GO analysis, we also found that the genes for secondary metabolism were overrepresented in the stigma. This GO term was also enriched in the stigma-specific genes within the pistil, and many pathways related to the secondary metabolism were enriched in the KEGG analysis of stigma-preferential genes within the pistil, although few have been studied in rice (Fig. 4a, Additional file 12: Figure S3A, Additional file 8: Table S6). In Arabidopsis, a type of pistil extract called azadecalin was isolated, and it could stimulate pollen germination. Azadecalin-like molecules also stimulate pollen germination . Thus, there may be some secondary metabolites, such as azadecalin-like molecules, that play essential roles during fertilization. To date, many stigma-specific genes have been found that are involved in the fertilization process of many species. For example, a cysteine-rich protein called STIG1 (stigma-specific protein 1) was found to be specifically expressed in the stigma, and it controls exudate secretion in tobacco and petunias. PrsS (Papaver rhoeas stigma S determinant), a stigma-specific gene, is involved in the S-specific inhibition of incompatible pollen grains, which is similar to SRK in Brassica [6, 17, 70]. However, little is known about stigma development and pollination in rice. Thus, in our study, 703 stigma-specific genes were identified within the pistil, and 47 of them were preferentially expressed in the plumose stigma, whereas 7 genes were only expressed in the stigma (Fig. 5). Interestingly, among these genes, we found two cysteine-rich peptides (LOC_Os02g55810 and LOC_Os02g55510) and two S-locus-related genes (LOC_Os01g48000 and LOC_Os01g48040) (Additional file 8: Table S8). These genes may participate in the stigma development or inhibition of incompatible pollen grains. Our study sets the foundation for gene function research in relation to the stigma and pollination.
Gene conservation and divergence between rice and Arabidopsis
Thus far, many embryo sac-specific genes have been shown to be involved in the development of normal embryo sacs and the success of double-fertilization in Arabidopsis . In rice, only two embryo sac-specific genes have been identified, but their functions are still unknown . As described above, there was no embryo sac in most of the dst ovaries (Fig. 1), and the two known rice embryo sac-specific genes were down-regulated in the rice dst ovary compared with the HY ovary (Additional file 1: Table S1). Thus, the genes that were down-regulated in the rice dst ovary compared with the HY ovary may be present in the embryo sacs. Additionally, many homologous genes were down-regulated in the rice dst ovary, implying that they have a similar expression pattern in rice (Fig. 7). Among the different species, many genes with a similar expression pattern and high homology, such as OsFZP and ZmBD1, play similar roles during the transition of the spikelet meristem to a floret meristem, and OsDMC1 and AtDMC1 also perform similar functions during homologous pairing [71–75]. Hence, we concluded that some of these homologous rice genes may be indispensable for the development of the embryo sac and fertilization, and there may be similar molecular mechanisms between rice and Arabidopsis during the reproductive process.
There were also partially or highly different mechanisms between rice and Arabidopsis with the evolution of the species. For example, the knockout of AtWUSCHEL resulted in defective shoot, flower, and fruit development, while the putative rice ortholog of AtWUSCHEL, which is called OsTAB1 (TILLERS ABSENT1), only regulated the initiation of axillary meristems with no effect on the shoot apical meristems . In Arabidopsis, mutations in the FIS complex genes (MEA, FIS2, FIE and MSI1) resulted in autonomous endosperm development without fertilization, and their seed developments were arrested at the heart embryo stage [77, 78]. In rice, there were no MEA and FIS2 genes, but two FIE orthologs, FIE1 and FIE2 are present. When FIE1 and FIE2 were knockout or knockdown, plants in rice failed to show autonomous central cell proliferation without fertilization [78–80]. In correspondence with function divergences, their expression patterns were also different. In our study, both OsFIE1 (LOC_Os08g04290) and OsFIE2 (LOC_Os08g04270) were not down-regulated in the rice dst ovary compared with the HY ovary, and many other homologous genes were also not down-regulated in the rice dst ovary (Additional file 22: Figure S7A). There are five lysophosphatidyl acyltransferases (LPATs) in Arabidopsis and three in rice. Among them, AtLPAT2 is specifically expressed in the embryo sac, and it is essential for the development of the embryo sac . However, none of the three rice genes were down-regulated in dst, suggesting that they may play different roles (Additional file 22: Figure S7B). Similarly, an Arabidopsis central cell-specific gene, AtCCG (central cell guidance), acts as an essential transcription regulator for fertilization , whereas its homolog in rice is abundantly expressed in the dst ovary and wild type HY ovary with no significant difference (Additional file 22: Figure S7C). Moreover, a synergid cell-specific putative glucosyl phosphatidylinositol-anchored protein called AtLORELEI controls the release of sperm cells in Arabidopsis . Using a BLASTP search, we found six homolog genes of AtLORELEI in rice. None of them were down-regulated in the dst mutant compared with the HY ovary, implying that the genes may not have functions similar to those in Arabidopsis. Interestingly, five of the six were preferentially expressed in the rice stigma, suggesting that they may play an essential role in stigma development or during pollination (Additional file 22: Figure S7D). To confirm our results, we use in situ hybridization to detect the three gene expression of LOC_11g41900 (homologous gene of AtLPAT2), LOC02g48980 and LOC_06g19990 (homologous gene of AtLORELEI). All the three genes were present in the entire ovary, not restricted in the embryo sac (Additional file 20: Figure S6C-H). All of the above rice genes from Arabidopsis homologs may perform different functions, implying that there is a divergent mechanism of individual development between rice and Arabidopsis.
Stigma-specific genes may show weak levels of expression
In comparing the stigma and ovary in Hwayoung rice, we identified 703 genes that were specifically expressed in the stigma. Theoretically, the stigma-specific genes would be down-regulated in the mutant dst ovary (no sigma) compared with the HY pistil. Unexpectedly, there were only 37 (5.3%) stigma-specific genes that were down-regulated in the mutant dst ovary (Fig. 5b). The analysis of the expression data showed that the FPKM values of most of the stigma-specific genes were very low in the stigma, and they were not greater than 10 (Additional file 10: Table S8). Among the stigma-specific genes, with the exception of the genes that were down-regulated in the mutant dst, they exhibited weaker expression (Additional file 23: Figure S8). In Arabidopsis, the 71 genes that were down-regulated in the mutant dif were screened using microarrays . However, AtMYB64/119 and AtAGL61, which were confirmed to be female gametophyte-specific genes in the following research, were not included because of their weak expression [46, 47, 83]. This result suggested that many cell-type-/stigma-specific genes were expressed weakly and could not be screened out by comparing the wild-type with the mutant gene using high-throughput analysis. Therefore, multiple techniques are required.
Analysis of the embryo sac-preferential/specific genes was conducive to investigating the gene function and regulation network in the female rice gametophyte
In our results, the expression of 122 genes was dependent on the existence of an embryo sac and may be specifically expressed in the ovary compared with the stigma in rice (Fig. 5e). Using reverse-genetic approaches, including T-DNA insertion and small RNA silencing, the functions of the genes could be investigated further, and some had to be knocked out/down together because they may play a redundant role. In Arabidopsis, both MYB64 and MYB119, which exhibit a female gametophyte-specific expression pattern, redundantly regulated the FG5 (Female Gametophyte 5) transition. The single mutants MYB64 or MYB119 showed no apparent defect, but the double mutant, myb64myb119, failed to initiate the FG5 transition . Similarly, the five egg cell-specific genes known as the EC1s  and the five synergid cell-specific genes known as LUREs  also performed redundant functions in Arabidopsis. In our study, we found that many embryo sac-preferential/specific genes may also have redundant functions, such as the Arabidopsis EC1s’ rice homologs, the ECA1-like genes.
Furthermore, we could construct promoter::GUS/GFP vectors, transform them into rice plants, and obtain lines containing the female gametophyte-specific marker for future research on female gametophyte development. For example, the four types of cells in the female gametophyte could be isolated and collected separately using a cell-sorting approach, and more cell type-specific genes could be identified. Moreover, with the help of the promoter:: GFP lines, we could identify rice mutants that have abnormal cell fates. Additionally, the promoters of these cell-specific genes could be used to identify the cis-regulatory elements, which would help us to investigate the gene regulatory function and network involved in female gametophyte development. In conclusion, these embryo sac-preferential/specific genes could be used as starting points for research on female gametophytes in rice.
More accurate method for gene identification and validation
With technological progress, the rice genome sequence and transcriptome became more accurate and more detailed. In 2002, the draft sequences of the rice genome, including japonica and indica, were completed using whole-genome shotgun sequencing. More than 32,000 genes have been predicted via assembled sequences, covering 93% of japonica and 92% of indica rice [84, 85]. Recently, high-throughput sequencing technology has attracted the attention of scientists because it has the advantages of having a low cost, good repeatability, and high-throughput capacity. By using the high-throughput sequencing technology, the Nipponbare genome assembly was updated, and the error rate was only approximately 0.15 per 10,000 nucleotides in the original IRGSP assembly . With the help of a new reference genome, 1240 loci from the RGAP 6.1 annotation set were thrown out from the RGAP 7.0 annotation set and 4993 loci from the RAP annotation on the IRGSP build 5 genomes were excluded in the new release.
Moreover, in animals and Arabidopsis, high-throughput transcriptome sequencing has been used to explore alternative splicing and novel genes to simultaneously discover transcripts and estimate their abundance [86, 87]. In rice, high-throughput transcriptome sequencing has primarily been used to detect gene expression, but almost none has been used for detecting novel genes. In our study, 102 novel protein-coding genes were identified, and 6 of them were verified in two Hwayoung and Nipponbare cultivars (Table 2, Fig. 11). Thus, we suggest that there are more novel genes in other rice tissues besides the pistil. Thus, using high-throughput transcriptome sequencing data in different rice tissues, the novel genes were revealed and verified more accurately and effectively.
In this study, we explored numerous stigma-preferential genes within the pistil through a comparison between the ovary and stigma by using RNA-Seq, and we found that the transport-, localization-, membrane-, communication-, secretion- and pollination-related genes were significantly enriched in the stigma through GO analysis, indicating that rice stigma may perform its function primarily through the membrane and that it was different from the wet stigma. Among these genes, many were specifically expressed in the stigma and verified by qRT-PCR and in situ hybridization, implying that they may play important roles in stigma development and pollination. In addition, many embryo sac-preferential/specific genes were found through the comparison between the wild ovary and dst ovary (without embryo sac). Most of them were preferentially expressed in pistils and may only be expressed in female gametophytes, suggesting that they may have essential roles in female gametophyte development and fertilization and could be used for investigating gene functions and molecular mechanisms. Moreover, we found that many “reads” were clustered far away from the annotated genes; we identified some novel protein-coding genes in the rice pistils, suggesting that the transcriptome data would serve as a valuable resource for identifying rice genes to supplement and perfect the existing transcriptome in rice. Therefore, these studies provided an important foundation for investigating the pollination and fertilization mechanisms in rice stigmas and female gametophytes.
The Rice plants (Oryza sativa L. ssp. japonica, cultivar ‘Hwayoung’, HY; cultivar ‘Nipponbare’, NIP), for RNA-sequencing and verification by semi-quatative RT-PCR, qRT-PCR and RNA in situ hybridization, were grown in green house at Wuhan University. The HY seedlings (SL) were collected after seeds were germinated and grown for 2 weeks. After the panicles emerged from the sheath of flag leaf, the samples of un-pollnated pistils (Pi), ovaries (O), stigmas (St), and anthers (An) were dissected from Hwayoung, meanwhile, the ovaries of dst (defective stigma) mutant, obtained from RISD DB (Rice T-DNA Insertion Sequence Database), were dissected before the florets flowering. And pistils from Nipponbare (NIP pistil) were collected after the panicles emerged from the sheath of flag leaf. At last, 5-d seed of HY (S) after the florets flowering were collected.
Ovary transparency and microscopy observation
The mature florets before flowering were collected, then stained using a modified protocol of the whole-mount eosin B colorant, and observed under confocal laser scanning microscopy according to Zeng’s method . The cleared ovaries were scanned and observed under an Olympus Fluoview 1000 laser scanning confocal microscope with a 535 nm argon laser.
RNA extraction and sequencing analysis
Total RNAs for RNA-sequencing from the pistils, ovaries, stigmas of wild type plants as well as ovaries of dst were extracted using Trizol reagent (Invitrogen, USA) following the manufacturer’s protocol and digested with DNaseI, and then quantified under Nanodrop Spectrophotometer (Nanodrop Technologies, Wilmington, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Böbelingen, Germany). After that, the samples were prepared using Illumina’s kit according manufacturer’s recommendation, and sequenced on the Illumina sequencing platform (HiSeqTM 2000) by the Shanghai Oebiotech Corporation. And 125 bp paired-end reads were generated. Our files containing sequence reads and quality scores were deposited in the Short Read Archive of the National Center for Biotechnology Information (NCBI) [Accession number SRP082242].
RNA-seq data preprocessing
Raw data (raw reads) of fastq format were firstly processed using the NGS QC Toolkit , and reads containing adapter or ploy-N and low quality were removed. All the downstream analyses were based on clean data with high quality (Q30). Then the sequencing reads were mapped to the rice genome (http://rice.plantbiology.msu.edu/pub/data/Eukaryotic_Projects/o_sativa/annotation_dbs/pseudomolecules/version_7.0/all.dir/) using Tophat (http://ccb.jhu.edu/software/tophat/index.shtml), and the mapped reads were assigned to the RGAP 7.0 gene modes using Bowtie2 [52, 90, 91] with default parameters by slightly modified. The FPKM (fragments per kilobase of gene per million) and count value were calculated using eXpress .
Differential expression analysis
Differential expression analysis was performed using the DESeq R package . The FDR (False Discovery Rate) ≤0.05 and fold change ≥ 2 was used as the threshold for significantly differential expression. The hierarchical cluster analysis of samples was also performed using the DESeq R package, and PCA (principal component analysis) was performed using the plotPCA function of the DESeq2 package. The first and second principal components explained 46.2% and 15.1% of variance, respectively. Comparison of GO function categories between ovary- and stigma-preferential genes was calculated using the WEGO online webtool (http://wego.genomics.org.cn/cgi-bin/wego/index.pl) based on Pearson Chi-Square test. The GO terms were removed if any of the expected counts were less than 5 or the p-value was more than the significant level of 0.05 . The other GO and KEGG enrichment analyses of the differentially expressed genes (DEGs) were performed based on hypergeometric distribution.
Digital expression analysis
Next generation sequencing transcriptome data were downloaded from the Rice Genome Annotation Project (http://rice.plantbiology.msu.edu/expression.shtml). The average value of a gene in one tissue plus 0.000001 was used logarithmic operation. Then, the logarithmic values were used for cluster display with the help of Cluster 3.0 and Treeview software.
BLASTP searches were performed to identify homologous genes using the sequences of AtCCG, AtLORELEI, AtMYB64/98/119, AtGEX1 and secretion-related proteins (GO: 0046903) with E-value = 0.00001. The candidate homologous genes of AtMYB64/98/119 were aligned in Clustal X (version 1.83) and MEGA 6 using neighbor-joining method (the bootstrap value was 500) [94, 95]. Then, an un-rooted phylogenetic tree was drawn.
Novel protein-coding gene forecast
The clean reads were mapped onto rice genome using Tophat (http://ccb.jhu.edu/software/tophat/index.shtml). Subsequently, transcripts were assembled using Cufflinks , and the short transcripts (the length <180 bp) were removed. Then, we identified the genes that were more than 200 bp away from any IRGSP-1.0_representative gene, and calculated the potential ability of coding protein using Coding Potential Calculator (CPC, http://cpc.cbi.pku.edu.cn/) while CPC ≥ 0 was set as the threshold . Besides, the transcriptome from the RGAP 7.0 version was also compared with the coding genes (CPC ≥ 0), and the genes that did not overlapped with any gene were thought as novel genes.
Quantitative RT-PCR and semi-quantitative RT-PCR analysis
Except the RNA samples for sequencing, the other RNA was extracted with the RNApure Extraction Kit (Bioteke, Beijing, China) according to the manual of the Kit. The methods of reverse transcription and RT-PCR were performed according to Zhao et al. . The threshold cycle (CT) value was automatically calculated by the Bio-Rad CFX Manager 3.1 system software, and the delta-delta Ct method was used to calculate the relative expression levels [99, 100]. The rice eEF-1α gene was used as an internal control for quantitative RT-PCR, while Actin gene for semi-quantitative RT-PCR . The primers were listed in the Additional file 24: Table S16.
In situ hybridization
Using the DIG RNA labeling kit (Roche), the digoxigenin labeled RNA probes were generated according the manual. The protocol of whole mount in situ hybridization was similar to the Hu’s method with some modifications and the protocol of the paraffin section was similar to Brewer’s procedure [101, 102]. The primers used were listed in the Additional file 24: Table S16.
Days after pollination
- dst :
Fragments per kilobase of exon per million mapped reads
Quantitative real time polymerase chain reaction
Rice Annotation Project Database
Rice Genome Annotation Project
We would like to thank Professor Yu Zhou (Wuhan University) for suggestion about our data processing.
This research was supported by the National Basic Research Program of China (2013CB126903, 2012CB944801). The National Natural Science Foundation of China (31670312).
Availability of data and material
All supporting data are included in the article or in the additional files.
The studies were conceived and coordinated by JZ and LY. Materials were collected by LY, YZ, YH, HM, KY and YC. Data analysis, and PCR were finished by LY, and in situ hybridization was performed by YZ and TM. The manuscript writing and editing were performed by LY, YC, TM and JZ. All authors approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
This study required no ethics approval.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- McCormick S. Control of male gametophyte development. Plant Cell. 2004;16(Suppl):S142–153.PubMedPubMed CentralView ArticleGoogle Scholar
- Yadegari R, Drews GN. Female gametophyte development. Plant Cell. 2004;16(Suppl):S133–141.PubMedPubMed CentralView ArticleGoogle Scholar
- Itoh J, Nonomura K, Ikeda K, Yamaki S, Inukai Y, Yamagishi H, Kitano H, Nagato Y. Rice plant development: from zygote to spikelet. Plant Cell Physiol. 2005;46(1):23–47.PubMedView ArticleGoogle Scholar
- Zhang DB, Wilson ZA. Stamen specification and anther development in rice. Chin Sci Bull. 2009;54:2342–53.View ArticleGoogle Scholar
- Lord EM, Russell SD. The mechanisms of pollination and fertilization in plants. Annu Rev Cell Dev Biol. 2002;18:81–105.PubMedView ArticleGoogle Scholar
- Takayama S, Shimosato H, Shiba H, Funato M, Che FS, Watanabe M, Iwano M, Isogai A. Direct ligand-receptor complex interaction controls Brassica self‐incompatibility. Nature. 2001;413(6855):534–8.PubMedView ArticleGoogle Scholar
- Stone SL, Anderson EM, Mullen RT, Goring DR. ARC1 is an E3 ubiquitin ligase and promotes the ubiquitination of proteins during the rejection of self-incompatible Brassica pollen. Plant Cell. 2003;15(4):885–98.PubMedPubMed CentralView ArticleGoogle Scholar
- Murase K, Shiba H, Iwano M, Che F-S, Watanabe M, Isogai A, Takayama S. A membrane-anchored protein kinase involved in Brassica self-incompatibility signaling. Science. 2004;303(5663):1516–9.PubMedView ArticleGoogle Scholar
- Kakita M, Murase K, Iwano M, Matsumoto T, Watanabe M, Shiba H, Isogai A, Takayama S. Two distinct forms of M-locus protein kinase localize to the plasma membrane and interact directly with S-locus receptor kinase to transduce self-incompatibility signaling in Brassica rapa. Plant Cell. 2007;19:3961–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Samuel MA, Chong YT, Haasen KE, Aldea-Brydges MG, Stone SL, Goring DR. Cellular pathways regulating responses to compatible and self-incompatible pollen in Brassica and Arabidopsis stigmas intersect at Exo70A1, a putative component of the exocyst complex. Plant Cell. 2009;21(9):2655–71.PubMedPubMed CentralView ArticleGoogle Scholar
- Safavian D, Zayed Y, Indriolo E, Chapman L, Ahmed A, Goring D. RNA silencing of exocyst genes in the stigma impairs the acceptance of compatible pollen in Arabidopsis. Plant Physiol. 2015;169(4):2526–38.PubMedPubMed CentralGoogle Scholar
- Luu DT, Heizmann P, Dumas C. Pollen-stigma adhesion in kale is Not dependent on the self-(in) compatibility genotype. Plant Physiol. 1997;115:1221–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Doughty J, Dixon S, Hiscock SJ, Willis AC, Parkin IA, Dickinson HG. PCP-A1, a defensin-like Brassica pollen coat protein that binds the S locus glycoprotein, is the product of gametophytic gene expression. Plant Cell. 1998;10:1333–47.PubMedPubMed CentralView ArticleGoogle Scholar
- Luu DT, Marty-Mazars D, Trick M, Dumas C, Heizmann P. Pollen-stigma adhesion in Brassica spp involves SLG and SLR1 glycoproteins. Plant Cell. 1999;11:251–62.PubMedPubMed CentralGoogle Scholar
- Takayama S, Shiba H, Iwano M, Asano K, Hara M, Che FS, Watanabe M, Hinata K, Isogai A. Isolation and characterization of pollen coat proteins of Brassica campestris that interact with S locus-related glycoprotein 1 involved in pollen-stigma adhesion. Proc Natl Acad Sci. 2000;97(7):3765–70.PubMedPubMed CentralView ArticleGoogle Scholar
- Goldman MH, Goldberg RB, Mariani C. Female sterile tobacco plants are produced by stigma-specific cell ablation. EMBO J. 1994;13:2976–84.PubMedPubMed CentralGoogle Scholar
- Verhoeven T, Feron R, Wolters-Arts M, Edqvist J, Gerats T, Derksen J, Mariani C. STIG1 controls exudate secretion in the pistil of petunia and tobacco. Plant Physiol. 2005;138(1):153–60.PubMedPubMed CentralView ArticleGoogle Scholar
- Qin Y, Wysocki RJ, Somogyi A, et al. Sulfinylated azadecalins act as functional mimics of a pollen germination stimulant in Arabidopsis pistils. Plant J. 2011;68(5):800–15.PubMedPubMed CentralView ArticleGoogle Scholar
- Ray SM, Park SS, Ray A. Pollen tube guidance by the female gametophyte. Development. 1997;124:2489–98.PubMedGoogle Scholar
- Shimizu KK, Okada K. Attractive and repulsive interactions between female and male gametophytes in Arabidopsis pollen tube guidance. Development. 2000;127(20):4511–8.PubMedGoogle Scholar
- Higashiyama T, Takeuchi H. The mechanism and key molecules involved in pollen tube guidance. Annu Rev Plant Biol. 2015;66:393–413.PubMedView ArticleGoogle Scholar
- Qu LJ, Li L, Lan Z, Dresselhaus T. Peptide signalling during the pollen tube journey and double fertilization. J Exp Bot. 2015;66(17):5139–50.PubMedView ArticleGoogle Scholar
- Kasahara RD, Portereiko MF, Sandaklie-Nikolova L, Rabiger DS, Drews GN. MYB98 is required for pollen tube guidance and synergid cell differentiation in Arabidopsis. Plant Cell. 2005;17(11):2981–92.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen YH, Li HJ, Shi DQ, et al. The central cell plays a critical role in pollen tube guidance in Arabidopsis. Plant Cell. 2007;19(11):3563–77.PubMedPubMed CentralView ArticleGoogle Scholar
- Escobar-Restrepo JM, Huck N, Kessler S, Gagliardini V, et al. The FERONIA receptor-like kinase mediates male–female interactions during pollen tube reception. Science. 2007;317(5838):656–60.PubMedView ArticleGoogle Scholar
- Endo S, Shinohara H, Matsubayashi Y, Fukuda H. A novel pollen-pistil interaction conferring high-temperature tolerance during reproduction via CLE45 signaling. Curr Biol. 2013;23(17):1670–6.PubMedView ArticleGoogle Scholar
- Hao L, Liu J, Zhong S, Gu H, Qu LJ. AtVPS41-mediated endocytic pathway is essential for pollen tube-stigma interaction in Arabidopsis. Proc Natl Acad Sci. 2016;113(22):6307–12.PubMedPubMed CentralView ArticleGoogle Scholar
- Cheung AY, Wang H, Wu HM. A floral transmitting tissue-specific glyxoprotein attracts pollen tubes and stimulates their growth. Cell. 1995;82:383–93.PubMedView ArticleGoogle Scholar
- Lamport DTA, Várnai P. Periplasmic arabinogalactan glycoproteins act as a calcium capacitor that regulates plant growth and development. New Phytol. 2013;197(1):58–64.PubMedView ArticleGoogle Scholar
- Ma HL, Yu L, Liang RH, Zhao J. Functional studies of arabinogalactan proteins in higher plants. Scientia Sinica Vitae. 2015;45:113–23.View ArticleGoogle Scholar
- Okuda S, Tsutsui H, Shiina K, et al. Defensin-like polypeptide LUREs are pollen tube attractants secreted from synergid cells. Nature. 2009;458(7236):357–61.PubMedView ArticleGoogle Scholar
- Kanaoka MM, Kawano N, Matsubara Y, et al. Identification and characterization of TcCRP1, a pollen tube attractant from Torenia concolor. Ann Bot. 2011;108:739–47.PubMedPubMed CentralView ArticleGoogle Scholar
- Takeuchi H, Higashiyama T. A species-specific cluster of defensin-like genes encodes diffusible pollen tube attractants in Arabidopsis. PLoS Biol. 2012;10(12):e1001449.PubMedPubMed CentralView ArticleGoogle Scholar
- Amien S, Kliwer I, Kliwer I, Márton ML, et al. Defensin-like ZmES4 mediates pollen tube burst in maize via opening of the potassium channel KZM1. PLoS Biol. 2010;8(6):e1000388.PubMedPubMed CentralView ArticleGoogle Scholar
- Li S, Li W, Huang B, et al. Natural variation in PTB1 regulates rice seed setting rate by controlling pollen tube growth. Nat Commun. 2013;4:2793.PubMedGoogle Scholar
- Lopez-Dee ZP, Wittich P, Enrico Pè M, Rigola D, Del Buono I, Gorla MS, Kater MM, Colombo L. OsMADS13, a novel rice MADS-box gene expressed during ovule development. Dev Genet. 1999;25(3):237–44.PubMedView ArticleGoogle Scholar
- Dreni L, Jacchia S, Fornara F, Fornari M, Ouwerkerk PB, An G, Colombo L, Kater MM. The D-lineage MADS-box gene OsMADS13 controls ovule identity in rice. Plant J. 2007;52(4):690–9.PubMedView ArticleGoogle Scholar
- Yamaki S, Nagato Y, Kurata N, Nonomura KI. Ovule is a lateral organ finally differentiated from the terminating floral meristem in rice. Dev Biol. 2011;351(1):208–16.PubMedView ArticleGoogle Scholar
- Zhang J, Tang W, Huang Y, Niu X, Zhao Y, Han Y, Liu Y. Down-regulation of a LBD-like gene, OsIG1, leads to occurrence of unusual double ovules and developmental abnormalities of various floral organs and megagametophyte in rice. J Exp Bot. 2015;66(1):99–112.PubMedView ArticleGoogle Scholar
- Luo Q, Li Y, Shen Y, Cheng Z. Ten years of gene discovery for meiotic event control in rice. J Genet Genomics. 2014;41(3):125–37.PubMedView ArticleGoogle Scholar
- Zhang B, Wang M, Tang D, et al. XRCC3 is essential for proper double-strand break repair and homologous recombination in rice meiosis. J Exp Bot. 2015;66:erv253.Google Scholar
- Chen J, Ding J, Ouyang Y, et al. A triallelic system of S5 is a major regulator of the reproductive barrier and compatibility of indica-japonica hybrids in rice. Proc Natl Acad Sci. 2008;105(32):11436–41.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang N, Huang HJ, Ren ST, Li JJ, Sun Y, Sun DY, Zhang SQ. The rice wall-associated receptor-like kinase gene OsDEES1 plays a role in female gametophyte development. Plant Physiol. 2012;160(2):696–707.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang J, Zhao X, Cheng K, et al. A killer-protector system regulates both hybrid sterility and segregation distortion in rice. Science. 2012;337(6100):1336–40.PubMedView ArticleGoogle Scholar
- Steffen JG, Kang IH, Macfarlane J, Drews GN. Identification of genes expressed in the Arabidopsis female gametophyte. Plant J. 2007;51(2):281–92.PubMedView ArticleGoogle Scholar
- Steffen J, Kang IM, Lloyd A, Drews GN. Agl61 interacts with agl80 and is required for central cell development in Arabidopsis. Plant Physiol. 2008;148(1):259–68.PubMedPubMed CentralView ArticleGoogle Scholar
- Rabiger DS, Drews GN. Myb64 and myb119 are required for cellularization and differentiation during female gametogenesis in Arabidopsis Thaliana. PLoS Genet. 2013;9(9):e1003783.PubMedPubMed CentralView ArticleGoogle Scholar
- Márton ML, Cordts S, Broadhvest J, Dresselhaus T. Micropylar pollen tube guidance by egg apparatus 1 of maize. Science. 2005;307(5709):573–6.PubMedView ArticleGoogle Scholar
- Li M, Xu W, Yang W, Kong Z, Xue Y. Genome-wide gene expression profiling reveals conserved and novel molecular functions of the stigma in rice. Plant Physiol. 2007;144(4):1797–812.PubMedPubMed CentralView ArticleGoogle Scholar
- Ohnishi T, Takanashi H, Mogi M, et al. Distinct Gene Expression Profiles in Egg and Synergid Cells of Rice as Revealed by Cell Type-Specific Microarrays. Plant Physiol. 2011;155(2):881–91.PubMedView ArticleGoogle Scholar
- Anderson SN, Johnson CS, Jones DS, Conrad LJ, Gou X, Russell SD, Sundaresan V. Transcriptomes of isolated Oryza sativa gametes characterized by deep sequencing: evidence for distinct sex-dependent chromatin and epigenetic states before fertilization. Plant J. 2013;76(5):729–41.PubMedView ArticleGoogle Scholar
- Kawahara Y, de la Bastide M, Hamilton J, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6:4.PubMedView ArticleGoogle Scholar
- Hollender CA, Kang C, Darwish O, et al. Floral transcriptomes in woodland strawberry uncover developing receptacle and anther gene networks. Plant Physiol. 2014;165(3):1062–75.PubMedPubMed CentralView ArticleGoogle Scholar
- Ma HL, Zhao J. Genome-wide identification, classification, and expression analysis of the arabinogalactan protein gene family in rice (Oryza sativa L.). J Exp Bot. 2010;61(10):2647–68.PubMedPubMed CentralView ArticleGoogle Scholar
- Ma HL, Zhao HM, Liu Z, Zhao J. The Phytocyanin Gene Family in Rice (Oryza sativa L.): Genome-wide identification, classification and transcriptional analysis. PLoS One. 2011;6(10):e25184.PubMedPubMed CentralView ArticleGoogle Scholar
- Ma TF, Ma HL, Zhao HM, Qi HD, Zhao J. Identification, characterization, and transcription analysis of xylogen-like arabinogalactan proteins in rice (Oryza sativa L.). BMC Plant Biol. 2014;14:299.PubMedPubMed CentralView ArticleGoogle Scholar
- SEQC/MAQC-III Consortium. A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control Consortium. Nat Biotechnol. 2014;32(9):903–14.View ArticleGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006;34 suppl 2:W293–297.PubMedPubMed CentralView ArticleGoogle Scholar
- Doucet J, Lee HK, Goring DR. Pollen Acceptance or Rejection: A Tale of Two Pathways. Trends Plant Sci. 2016. doi:10.1016/j.tplants.2016.09.004.PubMedGoogle Scholar
- Diekmann Y, Seixas E, Gouw M, Tavares-Cadete F, Seabra MC, Pereira-Leal JB. Thousands of rab GTPases for the cell biologist. PLoS Comput Biol. 2011;7(10):e1002217.PubMedPubMed CentralView ArticleGoogle Scholar
- Drews GN, Wang D, Steffen JG, Schumaker KS, Yadegari R. Identification of genes expressed in the angiosperm female gametophyte. J Exp Bot. 2011;62(5):1593–9.PubMedView ArticleGoogle Scholar
- Alandete-Saez M, Ron M, Leiboff S, McCormick S. Arabidopsis Thaliana GEX1 has dual functions in gametophyte development and early embryogenesis. Plant J. 2011;68(4):620–32.PubMedView ArticleGoogle Scholar
- Sprunck S, Rademacher S, Vogler F, Gheyselinck J, Grossniklaus U, Dresselhaus T. Egg cell-secreted EC1 triggers sperm cell activation during double fertilization. Science. 2012;338(6110):1093–7.PubMedView ArticleGoogle Scholar
- Silverstein KA, Moskal WA, Wu HC, Underwood BA, Graham MA, Town CD, VandenBosch KA. Small cysteine-rich peptides resembling antimicrobial peptides have been under-predicted in plants. Plant J. 2007;51(2):262–80.PubMedView ArticleGoogle Scholar
- Arora R, Agarwal P, Ray S, Singh AK, Singh VP, Tyagi AK, Kapoor S. MADS-box gene family in rice: genome-wide identification, organization and expression profiling during reproductive development and stress. BMC Genomics. 2007;8:242.PubMedPubMed CentralView ArticleGoogle Scholar
- Colombo M, Masiero S, Vanzulli S, Lardelli P, Kater MM, Colombo L. AGL23, a type I MADS-box gene that controls female gametophyte and embryo development in Arabidopsis. Plant J. 2008;54(6):1037–48.PubMedView ArticleGoogle Scholar
- Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14(2):178–92.PubMedView ArticleGoogle Scholar
- Ciampolini F, Shivanna KR, Cresti M. Organization of the stigma and transmitting tissue of rice (Oryza sativa L.). Plant Biol. 2001;3(2):149–55.View ArticleGoogle Scholar
- Wu HM, Wang H, Cheung AY. A pollen tube growth stimulatory glycoprotein is deglycosylated by pollen tubes and displays a glycosylation gradient in the flower. Cell. 1995;82(3):395–403.PubMedView ArticleGoogle Scholar
- Foote HCC, Ride JP, Franklin-Tong VE, Walker EA, Lawrence MJ, Franklin FCH. Cloning and expression of a distinctive class of self-incompatibility (S) gene from Papaver rhoeas L. Proc Natl Acad Sci. 1994;91(6):2265–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Klimyuk VI, Jones JDG. AtDMC1, the Arabidopsis homologue of the yeast DMC1 gene: characterization, transposon-induced allelic variation and meiosis-associated expression. Plant J. 1997;11(1):1–14.PubMedView ArticleGoogle Scholar
- Colombo L, Marziani G, Masiero S, Wittich PE, Schmidt RJ, Gorla MS, Pe ME. BRANCHED SILKLESS mediates the transition from spikelet to floral meristem during Zea mays ear development. Plant J. 1998;16(3):355–63.View ArticleGoogle Scholar
- Chuck G, Muszynski M, Kellogg E, Hake S, Schmidt RJ. The control of spikelet meristem identity by the branched silkless1 gene in maize. Science. 2002;298(5596):1238–41.PubMedView ArticleGoogle Scholar
- Zhu QH, Hoque MS, Dennis ES, Upadhyaya NM. Ds tagging of BRANCHED FLORETLESS 1 (BFL1) that mediates the transition from spikelet to floret meristem in rice (Oryza sativa L). BMC Plant Biol. 2003;3:6.PubMedPubMed CentralView ArticleGoogle Scholar
- Deng ZY, Wang T. OsDMC1 is required for homologous pairing in Oryza sativa. Plant Mol Biol. 2007;65:31–42.PubMedView ArticleGoogle Scholar
- Galli M, Gallavotti A. Expanding the regulatory network for meristem size in plants. Trends Genet. 2016;32(6):372–83.PubMedView ArticleGoogle Scholar
- Yadegari R, Kinoshita T, Lotan O, et al. Mutations in the FIE and MEA genes that encode interacting polycomb proteins cause parent-of-origin effects on seed development by distinct mechanisms. Plant Cell. 2000;12(12):2367–81.PubMedPubMed CentralView ArticleGoogle Scholar
- Nallamilli BRR, Zhang J, Mujahid H, Malone BM, Bridges SM, Peng Z. Polycomb group gene OsFIE2 regulates rice (Oryza sativa) seed development and grain filling via a mechanism distinct from Arabidopsis. PLoS Genet. 2013;9(3):e1003322.PubMedPubMed CentralView ArticleGoogle Scholar
- Li S, Zhou B, Peng X, Kuang Q, Huang X, Yao J, et al. OsFIE2 plays an essential role in the regulation of rice vegetative and reproductive development. New Phytol. 2014;201(1):66–79.PubMedView ArticleGoogle Scholar
- Folsom JJ, Begcy K, Hao X, Wang D, Walia H. Rice Fertilization-Independent Endosperm1 regulates seed size under heat stress by controlling early endosperm development. Plant Physiol. 2014;165(1):238–48.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim HU, Li Y, Huang AH. Ubiquitous and endoplasmic reticulum-located lysophosphatidyl acyltransferase, LPAT2, is essential for female but not male gametophyte development in Arabidopsis. Plant Cell. 2005;17(4):1073–89.PubMedPubMed CentralView ArticleGoogle Scholar
- Capron A, Gourgues M, Neiva LS, et al. Maternal control of male-gamete delivery in Arabidopsis involves a putative GPI-anchored protein encoded by the LORELEI gene. Plant Cell. 2008;20(11):3038–49.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang D, Zhang C, Hearn DJ, Kang IH, Punwani JA, et al. Identification of transcription-factor genes expressed in the Arabidopsis female gametophyte. BMC Plant Biol. 2010;10:110.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu J, Hu S, Wang J, et al. A draft sequence of the rice genome (Oryza sativa L. ssp.indica). Science. 2002;296(5565):79–92.PubMedView ArticleGoogle Scholar
- Goff SA, Ricke D, Lan TH, et al. A draft sequence of the rice genome (Oryza sativa L. ssp. japonica). Science. 2002;296(5565):92–100.PubMedView ArticleGoogle Scholar
- Trapnell C, Williams BA, Pertea G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Loraine AE, McCormick S, Estrada A, Patel K, Qin P. RNA-seq of Arabidopsis pollen uncovers novel transcription and alternative splicing. Plant Physiol. 2013;162(2):1092–109.PubMedPubMed CentralView ArticleGoogle Scholar
- Zeng YX, Hu CY, Lu YG, Li JQ, Liu XD. Diversity of abnormal embryo sacs in indica/japonica hybrids in rice demonstrated by confocal microscopy of whole ovary. Plant Breed. 2007;126(6):574–80.View ArticleGoogle Scholar
- Patel RK, Jain M. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data. PLoS One. 2012;7(2):e30619.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim D, Salzberg SL. TopHat-Fusion: an algorithm for discovery of novel fusion transcripts. Genome Biol. 2011;12(8):R72.PubMedPubMed CentralView ArticleGoogle Scholar
- Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9(4):357–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.PubMedView ArticleGoogle Scholar
- Anders S, Huber W. Differential expression of RNA-Seq data at the gene level-the DESeq package. Heidelberg: European Molecular Biology Laboratory (EMBL); 2012.Google Scholar
- Saitou N, Nei M. The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4(4):406–25.PubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Trapnell C, Roberts A, Goff L, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.PubMedPubMed CentralView ArticleGoogle Scholar
- Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, Gao G. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35 suppl 2:W345–349.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhao H, Ma T, Wang X, Deng Y, Ma H, Zhang R, Zhao J. OsAUX1 controls lateral root initiation in rice (Oryza sativa L.). Plant Cell Environ. 2015;38(11):2208–22.PubMedView ArticleGoogle Scholar
- Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29(9):e45.PubMedPubMed CentralView ArticleGoogle Scholar
- Jiang SY, Vanitha J, Bai Y, Ramachandran S. Identification and molecular characterization of tissue-preferred rice genes and their upstream regularly sequences on a genome-wide level. BMC Plant Biol. 2014;14:331.PubMedPubMed CentralView ArticleGoogle Scholar
- Hu TX, Yu M, Zhao J. Comparative transcriptional profiling analysis of the two daughter cells from tobacco zygote reveals the transcriptome differences in the apical and basal cells. BMC Plant Biol. 2010;10(1):167.PubMedPubMed CentralView ArticleGoogle Scholar
- Brewer PB, Heisler MG, Hejátko J, Friml J, Benková E. In situ hybridization for mRNA detection in Arabidopsis tissue sections. Nat Protoc. 2006;1(3):1462–7.PubMedView ArticleGoogle Scholar