Pomegranate flowers were collected from 10-year-old “Tunisiruanzi” trees grown in nursery of the Zhengzhou Fruit Research Institute located in Zhengzhou, Henan, China, and managed using conventional methods. BF and FMF were separated based on the morphology of the ovary, which in FMF is atrophied (Holland et al., 2009). The ovules of FMF cease differentiation when their BVD was 5.0–13.0 mm. A total of 18 accessions were sequenced, comprising three BF and three FMF, and including three biological replicates. To describe the accessions clearly, we used TNSI, TNSII, and TNSIII represent BF buds when their BVD was 3.0–5.0, 5.1–13.0, 13.1–25.0 mm, respectively. Similarly, ATNSI, ATNSII, and ATNSIII represented the FMF buds when their BVD was 3.0–5.0, 5.1–13.0, 13.1–25.0 mm, respectively. Following collection, all samples were immediately frozen in liquid nitrogen and stored at -80°C.

Total RNA extraction was performed using the CTAB (Cetyltrimethyl Ammonium Bromide) method. RNA purity was checked using a NanoDrop^® 2000 instrument (Thermo, Wilmington, CA, United States). RNA concentration and integrity were measured using the RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent Technologies, Santa Clara, CA, United States). Library preparation and sequencing reactions were carried out by 1Gene, Corp. (Hangzhou, China)¹ according to the relevant manufacturer’s instructions (Illumina, San Diego, CA, United States).

A total amount of 5 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using a NEBNext^® Ultra^TM RNA Library Prep Kit for Illumina^® (NEB, Ipswich, MA, United States) following the manufacturer’s recommendations, and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEB-Next^® First Strand Synthesis Reaction Buffer (5×) (NEB, Ipswich, MA, United States). First-strand cDNA was synthesized using random hexamer primer and M-MuLV reverse transcriptase (RNase H-). Second-strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H. Remaining overhangs were converted into blunt ends using exonuclease/polymerase activities. Following adenylation of the 3′-ends of the DNA fragments, a NEBNext adaptor with a hairpin loop structure was ligated to prepare the fragments for hybridization. In order to select cDNA fragments of preferentially 250–300 bp in length, the library fragments were purified using an AMPure XP system (Beckman Coulter, Beverly, MA, United States). The size-selected, adaptor-ligated cDNA was incubated with 3 μl of USER Enzyme (NEB) at 37°C for 15 min. PCR was then performed with Q5 Hot Start HiFi DNA polymerase (NEB, Ipswich, MA, United States), Universal PCR primers and Index (X) Primer. Finally, the PCR products were purified (AMPure XP system) and library quality was assessed using the Agilent Bioanalyzer 2100 system.

Clustering of the index-coded samples was performed on a cBot Cluster Generation System using a TruSeq PE Cluster Kit v4-cBot-HS (Illumina) according to the manufacturer’s instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq 2500 platform and paired-end reads were generated.

To verify the RNA-Seq analysis, we performed qRT-PCR analysis of the expression of nine genes, which included two genes (ETR, ERF1/2) related to the ethylene response and seven genes (three AG, two SPL, one ANT and one INO) related to pistil development. The actin gene was used as a reference (Supplementary Table S9). RNA was extracted as described above and cDNA was synthesized using a TIANscript RT Kit (TIANGEN, Beijing, China). PCR was performed using 2 × SYBR Green Real-time PCR Master Mix (Roche) on a real-time PCR instrument (Roche 480, Basle, Switzerland), with the following program: pre-incubation at 95°C for 5 min, then followed by 45 cycles of 95°C for 10 s, 60°C for 10 s, and 72°C for 10 s. The ddCt method was used to calculate gene expression levels (Livak and Schmittgen, 2001). Primer sequences are listed in Supplementary Table S9.

Pistil of FMF and BF were collected when their BVD were 3.0–5.0, 5.1–10.0, 10.1–13.0, 13.1–15, 15.1–18.0, 18.1–20.0, 21.1–25.0, and 25.1–35.0 mm (blooming). Samples were fixed in 2.5% glutaraldehyde (pH = 7.4) for > 1 week at 4°C. Following fixation, they were dehydrated using an ethanol series [30% ethanol, 20 min; 50% ethanol, 20 min; 70% ethanol, 20 min; 100% ethanol, 30 min (twice)]. The dehydrated samples were then dried in a critical-point drying apparatus (Quorum, England). Dried samples were mounted on stubs and sputter-coated with gold (FEI, America) and observed under a SEM (FEI Quanta 250, America) in Henan University. Ten floral buds were observed in each stage.

Flower buds at different stages were collected based on BVDs. Measurements of flower BVD, OTDs, SL were made using a Dsect mirror (DS-Fiec, Nikon) or by vernier caliper. Ten biological replicates were measured for each stage. Difference significance analysis between BF and FMF using independent-samples t-test, P-value represent significance of difference.

“Tunisiruanzi” pomegranate trees were sprayed with an aqueous solution of ethephon at concentrations of 150, 200, and 250 mg/L, with three biological replications per treatment. Control sprayings were with water alone. Each tree was sprayed with 1. Five liter aqueous solution of ethephon or water (control). Treatments were sprayed on whole trees once, including both floral buds and leaves at a BVD of 3.0–5.0 mm (April 22, 2016), at which stage it is not yet possible to distinguish BF from FMF.

The morphology of FMF is bell-shaped, whereas that of BF is vase-shaped (Figures 1 FMF, 1BF). Pistils of FMF and BF were collected when their BVDs were: 3.0–5.0 mm (Figures 1FMF1, 1BF1); 5.1–10.0 mm (Figures 1FMF2, 1BF2); 10.1–13.0 mm (Figures 1FMF3, 1BF3); 13.1–15.0 mm (Figures 1FMF4, 1BF4); 15.1–18.0 mm (Figures 1FMF5, 1BF5); 18.1–20.0 mm (Figures 1FMF6, 1BF6); 20.1–25.0 mm (Figures 1FMF7, 1BF7) and during flowering (BVD ≥ 25.1 mm) (Figures 1FMF8, 1BF8). The dynamics of ovule development of FMF and BF were observed (Figure 1).

There were no obvious differences observed between BF1 and FMF1 (BVD 3.0–5.0 mm), at this stage, the placenta had formed but no ovule primordia were observable. Stage BF2 (BVD 5.1–10.0 mm) was characterized by the formation of inner integument primordia. However, in FMF2, though the ovules were enlarged, no inner integument primordia had formed. At BF3 (BVD 10.1–13.0 mm), outer integument primordia formed, the inner integument showed an increase in the number of cell layers and became symmetrically enlarged, growing parallel to the nucellus by anticlinal division and cell elongation (Reiser and Fischer, 1993). By contrast, at FMF3 (BVD 10.1–13.0 mm), inner integument primordium formed in part of the ovules, but no inner integument extension and outer integument primordia were observed, ovules displayed wilting. Thereafter, ovules undergo wilting during stage FMF4–8. In BF, on the other hand, the integument continued to enlarge during stage BF4 (BVD 13.1–15.0 mm); and during the subsequent stages BF5–7 (BVD 15.1–25.0 mm), the outer integument exhibited a gradient of cell division, which showed maximal growth on the abaxial side and almost no growth on the adaxial side (Robinson-Beers et al., 1992). By the reason of this asymmetric growth, by stage BF8 (when the flowers were in bloom), the outer integument completely enclosed the inner integument and the nucleus. From this comparison of FMF and BF, it was apparent that early ovule development in FMF was indistinguishable from that in BF. One of the reasons for pistil abortion in FMF may be the termination of ovule development after stage FMF3 (BVD 10.1–13.0 mm).

At flowering time, FMF were smaller than BF (Figure 1A). FMF bloomed at a BVD of 24.54 ± 1.59 mm (Figures 2A–C), a SL of 7.97 ± 1.0 mm (Figure 2A), an OTD of 6.67 ± 1.33 mm (Figure 2B), and an OVD of 14.22 ± 1.85 mm (Figure 2C). In contrast, BF bloomed at a BVD of 40.18 ± 2.37 mm (Figures 2A–C), a SL of 18.34 ± 1.48 mm (Figure 2A), an OTD of 14.22 ± 1.85 mm (Figure 2B), and an OVD of 14.22 ± 1.85 mm (Figure 2C). Thus, during the flowering period, the OTD (P < 0.01), the SL (P < 0.05) and the OVD (P < 0.01) were clearly shorter in FMF than in BF (Supplementary Table S1).

Remarkably, in FMF, the increase in OTD, SL and OVD began to decline when BVD grow into 10.32 ± 0.77–15.82 ± 0.13 mm (Figures 2A,B), suggesting OTD, SL and OVD can be indicators for female sterility. The SEM results showed that the key stage at which ovule abortion occurred when BVD was 10.0–13.0 mm (Figure 1). The data further prove the key stage of pomegranate female sterility was 10.0–13.0 mm in BVD.

In order to denote the accessions clearly, we used TNSI, TNSII, and TNSIII represent the BF buds’ pistils when their BVD was 3.0–5.0, 5.1–13.0, 13.1–25.0 mm, respectively. Similarly, the designations ATNSI, ATNSII, and ATNSIII were used to represent the FMF buds’ pistils when their BVD was 3.0–5.0, 5.1–13.0, 13.1–25.0 mm, respectively. There were three replicates for each developmental stage. In total, 18 accessions were used for transcriptomic sequencing. Sequencing of 18 accessions by Illumina platform obtained a total of 841 million clean reads of sequence, with an average of 46.7 million clean reads and 92.08% of mapped rate for each accession (Supplementary Table S2). The raw data were uploaded to NCBI Sequence Read Archive (SRA), a total of 18 samples with accession numbers SRX2735567-SRX2735584 were uploaded⁵. Pearson r² correlation values for all replications were used to evaluate the consistency of the raw data. Most of them varied from 0.85 to 1.0, and data with a value <0.80 were removed prior to subsequent analysis (Supplementary Figure S1). The data showed that the sequencing quality was high enough for further analysis (Supplementary Figure S1). A total of 25,959 genes were detected, with only a small variation between accessions (Supplementary Table S2). The highest number of expressed genes was 22,431, obtained from ATNSI (Figure 3A), while the lowest number (22,144) was from ATNSIII (Figure 3A and Supplementary Table S2).

A total of 9662 unique DEGs were obtained. Between TNSI and ATNSI, the number of DEGs was 281, whereas it was 3,772 between TNSII and ATNSII. However, the number decreased to 2,446 when comparing TNSIII with ATNSIII (Figure 3B). Moreover, a total of 7,828 DEGs were detected between the FMF different development stages (ATNSI-ATNSII, ATNSI-ATNSIII, and ATNSII-ATNSIII) and 4,744 DEGs were detected between three stages of BF development (TNSI-TNSII, TNSI-TNSIII, and TNSII-TNSIII) (Figures 3C,D). The DEGs were listed in Supplementary Table S3.

An enrichment analysis was performed to determine whether the genes (DEGs) that were differentially expressed between FMF and BF were markedly associated with specific pathways or biological processes. The DEGs were characterized using the GO and KEGG databases. The GO terms comprise three categories: molecular function, cellular component, and biological process. A total of 4,781 unique genes were enriched between FMF and BF (TNSI-ATNSI, TNSII-ATNSII, and TNSIII-ATNSIII). For the TNSI-ATNSI DEGs, the classification terms that showed enrichment included “response to hormone” (GO:0009725), “cell death” (GO:0008219), “signal transduction” (GO:0007165), “regulation of programmed cell death” (GO:0043067), “response to ethylene” (GO:0009723), “ethylene biosynthetic process” (GO:0009693), and “regulation of cellular process (GO:0050794)” (Supplementary Table S3). In contrast, for the TNSII-ATNSII DEGs, the enriched classification terms were “cell proliferation” (GO:0008283), “DNA methylation” (GO:0006306, GO:0044728), “regulation of cell cycle” (GO:0051726), “cell cycle” (GO:0007049), “mitotic cell cycle process” (GO:1903047), “G2/M transition of mitotic cell cycle” (GO:0000086), “cell division” (GO:0051301), “regulation of gene expression, epigenetic” (GO:0040029), and “nuclear division” (GO:0000280) (Supplementary Table S3). For the TNSIII-ATNSIII comparison, the enriched terms belonged to categories including “cell cycle” (GO:0007049), “cell division” (GO:0051301), “mitotic cell cycle process” (GO:1903047), “regulation of growth” (GO:0040008), “flower development” (GO:0009908), “nuclear division” (GO:0000280), “reproductive structure development” (GO:0048608), and “floral organ morphogenesis” (GO:0048444) (Supplementary Table S3). From the SEM results, ovule development of FMF ceased after the formation of the inner-integument primordium. Thus, the DEGs in the enriched categories such as “cell cycle,” “cell division,” “floral development,” and “floral organ morphogenesis” may be candidate genes associated with pistil abortion. The enrichment analysis revealed a total of 114 DEGs related to flower development (Supplementary Table S3).

Spraying with ethephon can increase the proportion of BF (Chaudhari and Desai, 1993). Consistent with this, the category “response to ethylene” (GO:0009723) was enriched between BF and FMF, consistent with the role of ethylene in pomegranate female sterility.

Pathway assignment was carried out by KEGG to analyze the biological functions of DEGs. A total of 1,302 unique DEGs between FMF and BF were mapped into KEGG pathways, containing 41, 809 and 452 DEGs between TNSI-ATNSI, TNSII-ATNSII and TNSIII-ATNSIII comparison, respectively. For the TNSI-ATNSI comparison, pathways related to “plant hormone signal transduction” (ko04075), “plant–pathogen interaction” (ko04626), and “biosynthesis of secondary metabolites” (ko01110) were enriched (Supplementary Table S4). For the TNSII-ATNSII comparison, the main pathways included “biosynthesis of secondary metabolites” (ko01110), “plant hormone signal transduction” (ko04075), and “metabolic pathways” (ko01100) (Supplementary Table S4). Similarly, “plant hormone signal transduction” (ko04075), “metabolic pathways” and “biosynthesis of secondary metabolites” (ko01110) were enriched in the TNSIII-ATNSIII comparison (Supplementary Table S4). The process of plant hormone signal transduction (Supplementary Figure S3) was concluded to be a putative pathway affecting pomegranate FMFs’ ovule development.

Constructing a dendrogram, a total of 17 distinct modules were obtained, in which each tree branch formed a module, and each leaf in the branch represented one gene (Figure 4A and Supplementary Table S5). The correlation between different modules were showed in Supplementary Figure S2. Genes in different modules were characterized by GO (Supplementary Table S6). From the measurement of the various different indices for FMF and BF, we determined differences in OTD, OVD, and SL as the growth of BVD between FMF and BF (Figure 2). The associations between these parameters and the 17 distinct modules were then compiled (Figure 4B). The modules related to SL are those named MEantiquewhite2, MEdarkviolet, MEdarkturquoise and MEfirebrick 3 (r ≥ 0.6) (Figure 4B); the same correlations were seen for OTD and OVD (Figure 4B); the modules related to BVD are those named MEdarkslateblue, MEantiquewhite2, MEdarkviolet, MEdarkturquoise and MEfirebrick 3 (r ≥ 0.6) (Figure 4B). According to the GO enrichment, “single-organism process” (GO:0065007), “organic substance metabolic process” (GO:0071704), “biological regulation” (GO:0065007), “organic substance biosynthetic process” (GO:0044711), “small molecule metabolic process” (GO:0044281), “phosphate-containing compound metabolic process” (GO:0006796), “lipid metabolic process” (GO:0006629), “lipid biosynthetic process” (GO:0008610) and “organonitrogen compound metabolic process” (GO:1901564) were also enriched in these modules (Supplementary Table S6). Genes related to flower development were revealed by the association between modules and tissues (Figure 4C). Gene modules showing differential expression between FMF and BF are named MElavenderblush 3, MElightcyan1, MEdarkred, MEdarkolivegreen, and MEpalevioletred 3 (Figure 4C). These modules were enriched mainly in the categories “cellular process” (GO:0009987), “single-organism process” (GO:0044699), and “metabolic process” (GO:0008152). Processes associated with flower development were also enriched (GO:0048573, GO:0009911, GO:0010227, GO:0048573, GO:0009911, GO:0009908, GO:0009909, GO:0048574) (Supplementary Table S6). A total of 20 genes related to flower development were enriched by the module-tissues association (Supplementary Figure S4 and Table S7). Analysis of the association between modules and tissues revealed enrichment of 20 genes related to flower development (Supplementary Figure S3 and Table S7).

Abnormal ovule development was the main factor responsible for the abortion of pomegranate flowers. Ovule development comprises several processes: primordia initiation, specification of identity, pattern formation, and eventually morphogenesis and cellular differentiation (Grossniklaus and Schneitz, 1998). The pathway of regulation of Arabidopsis ovule development was showed in Figure 5, and key regulators include AGAMOUS (Yanofsky et al., 1990; Drews et al., 1991), SPL (Ito et al., 2004), INO (L’eon-Kloosterziel et al., 1994), ANT and BEL1 (Brown et al., 2010). A total of ten AGAMOUS-like DEGs (Gglean016200, Gglean018861, Gglean020089, Gglean026618, Gglean028014, Gglean028632, Gglean030294, Gglean002838, Gglean006764, and Gglean022339), five SPL homolog genes (Gglean016197, Gglean022102, Gglean028488, Gglean001912, and Gglean005812), two ANT homolog genes (Gglean003340 and Gglean011480), two BEL1 homolog genes (Gglean014222 and Gglean001981), and one INO homolog gene (Gglean016270) involved in ovule development were differential expressed between BF and FMF (Figure 5 and Supplementary Table S8). In the TNSII-ATNSII comparison, which was the key stage for pistil abortion in pomegranate FMF, three AG-like genes (Gglean026618, Gglean028632, and Gglean028014), one SPL homolog gene (Gglean005812), two ANT homolog genes (Gglean003340 and Gglean011480), and one INO homolog gene (Gglean016270) showed significantly different expression levels (Supplementary Table S8). The results indicated that these genes may be the main regulators for the formation of FMF.

We identified three DEGs involved in the ethylene signal transduction process: ETR (Gglean022853), EIN3 (Gglean14469), and ERF1/2 (Gglean022880) (Figure 6A and Supplementary Table S8). To verify the role of ethylene in FMF, we sprayed ethephon at the stage of TNSI (ATNSI) with different concentration. We found that at all concentrations used ethephon had a significant effect on the proportion of BF. With water sprayed as a control, the proportion of BF was 2.2 ± 0.00%. At 150 mg/L of ethephon, the proportion of BF was increased to 19.83 ± 0.51%, whereas at 200 and 250 mg/L ethephon the proportion was 11.5 ± 0.43 and 10.5 ± 0.39%, respectively. The results indicated that 150 mg/L was the optimal concentration for improving the proportion of BF in pomegranate (Supplementary Table S8). The findings using ethephon were consistent with the transcriptomics results indicating that DEGs related to ethylene response may be associated with pistil abortion in FMF (Figure 6B). Ethylene may act as an upstream factor for pomegranate FMF female sterility.

To confirm the results of RNA sequencing, seven DEGs related to ovule development including three AGAMOUS-like DEGs (Gglean026618, Gglean028632, and Gglean028014), one SPL homolog gene (Gglean005812), two ANT homolog genes (Gglean003340 and Gglean011480) and one INO homolog gene (Gglean016270) were selected for qRT-PCR analysis. In addition, two DEGs associated with ethylene signal transduction, ETR (Gglean022853) and ERF1/2 (Gglean022880), were also analyzed between BF and FMF at three stages during flower development (Figure 7). The results showed a similar expression trend to the transcriptome analysis. Three AGAMOUS-like homologous genes showed differing expression patterns. Gglean026618 was up-regulated in ATNSII in the ATNSII-TNSII comparison, and Gglean028632 displayed up-regulated in FMF in the comparisons ATNSI-TNSI, ATNSII-TNSII, and ATNSIII-TNSIII. In contrast, Gglean028014 was down-regulated in ATNSII relative to TNSII. Similarly, the SPL homolog gene Gglean005812 was down-regulated in ATNSII relative to TNSII. Two ANT homolog genes (Gglean003340 and Gglean011480) was down-regulated in ATNSII relative to TNSII and almost no expression in TNSIII and ATNSIII. The INO homolog gene (Gglean016270) displayed almost the same trend as the ANT homologs. ETR (Gglean022853), the repressor of ethylene signal transduction, was up-regulated in ATNSI relative to TNSI and ATNSIII relative to TNSIII. ERF1/2 (Gglean022880), the downstream gene in ethylene signal transduction pathway, displayed differential expression in both the ATNSI-TNSI and ATNSII-TNSII comparisons, being up-regulated in ATNSI relative to TNSI but down-regulated in ATNSII relative to TNSII.

Here, we set out to explore candidate genes associated with pistil abortion in pomegranate FMF. The key stage for pistil abortion in FMF was found to be at a BVD of 5.0–13.0 mm, when the development of the ovule ceased following the formation of the inner integument primordium. Pistil abortion is widespread in a series of fruit crops, including Japanese apricot, X. sorbifolia, and Prunus armeniaca L. (Lillecrapp et al., 1999; Gao et al., 2002; Shi et al., 2011). Although the morphological characteristics of pistil abortion are distinct between species, all instances of pistil abortion result from the abnormal development of ovules. Pomegranate FMF display termination of ovule development after the formation of the inner integument primordium, and this may be the main reason for pistil abortion in this species. In addition, in Japanese apricot it has been reported that pistil abortion may be associated with the catabolism of macromolecular nutrients in the flower buds (Shi et al., 2011). Based on module-trait co-expression analysis (Figure 4B), some processes related to macromolecular nutrients, including “organic substance biosynthetic process” (GO:0044711), “small molecule metabolic process” (GO:0044281), and “organ nitrogen compound metabolic process” (GO:1901564), were found to be overrepresented, suggesting that the abortion of pomegranate FMF may be related to the catabolism of nutrients (Figure 4B).

A comparative transcriptomic analysis between BF and FMF in different developmental stages was then undertaken to mine candidate genes. Significant progress in understanding the molecular mechanism of ovule development has been achieved in some species, such as Arabidopsis, cotton, and rice (Baker et al., 1997; Lee et al., 2006; Kubo et al., 2013). AG, SPL, ANT, BEL1 and INO have been found to be key regulators of ovule development (Ito et al., 2004; Sieber et al., 2004; Ó’Maoiléidigh et al., 2013; Wei et al., 2015). In the present study, seven DEGs related to ovule development, expressed differentially at the key stage of FMF pistil abortion, have been verified by means of qRT-PCR, including three AG homolog genes, one SPL homolog, two ANT homologs and one INO homolog. Previous studies have indicated that ANT and INO are key regulators of the formation of the ovule integument (Grossniklaus and Schneitz, 1998), and here both showed down-regulation in ATNSII relative to TNSII. Furthermore, the morphology of the ovules of pomegranate FMF was similar to that of Arabidopsis ant-72F5 ovules (Grossniklaus and Schneitz, 1998). Collectively, our results indicated that ANT may be a key regulator gene of pistil abortion in pomegranate FMF. SPL can promote the growth of ovule integuments (Balasubramanian and Schneitz, 2002). SPL was down-regulation in the key stage of pistil abortion (ATNSII-TNSII) in FMF, which indicated the role of SPL promoting pomegranate ovule development. Previous studies indicated SPL represses the expression of ANT and INO to control ovule development (Wei et al., 2015). However, in the present study, SPL, ANT and INO were down-regulated in ATNSII relative to TNSII, the possible explanation is that ANT and INO may be repressed by other upstream factors for pomegranate ovule abortion.

As important regulatory genes in flower development, AG homologs play important roles in the formation of the pistil and stamen (Bowman et al., 1991). In this study, we detected three AG homolog genes amongst the DEGs, of which two were up-regulated (Gglean026618 and Gglean028632) and one was down-regulated (Gglean028014) in ATNSII relative to TNSII. The two up-regulated genes in ATNSII were AGL19 homology gene (Gglean026618) and AGL8 homology gene (Gglean028632). Both of these two genes play important roles in floral transition (Alvarez-Buylla et al., 2000). AGL8 promotes reproductive transition through interaction with SVP and SOC1 factors (Balanzà et al., 2014). Similarly, AGL19 acts as a floral activator in FLC-independent vernalization pathway (Schönrock et al., 2006). Here we found the differential expression of these two genes between BF and FMF, suggesting AGL19 and AGL8 also play important roles in regulation of female sterility in pomegranate. The down-regulated AG gene was the homolog gene of AGL62, which can stimulate nucellus degeneration (Bertoni, 2016). Moreover, AGL62 homology gene Gglean028014 was up-regulated in ATNSIII relative to TNSIII. In this study, we observed ovule degradation in ATNSIII (Figure 1), suggesting that the up-regulation of AGL62 may underlie the female sterility by contributing to the degradation of ovule.

Ethylene has been verified as an upstream factor for ovule development in tobacco (De Martinis and Mariani, 1999). In this study, ethylene response signal factor including ETR and ERF 1/2 displayed differential expression between BF and FMF. ETR protein acts as a repressor for the transmission of ethylene signal (Hua and Meyerowitz, 1998). ETR (Gglean022853) was up-regulated in ATNSI relative to TNSI indicated the suppression of ethylene signal transduction may be a factor influencing pomegranate ovule development. Down-regulation of ERF1/2 (Gglean022880) in ATNSII relative to TNSII also supported this point. Moreover, the increase of the ratio of BF in pomegranate with external ethephon treatment further confirmed the function of ethylene in pomegranate. Jointly, our findings indicated that ethylene may act as one of the upstream regulators influencing the development of pomegranate ovule. Furthermore, ANT contains two AP2 domains homologous with the DNA binding domain of EREBPs (Shigyo et al., 2006), The down-regulation of ANT may be related to ethylene.

Collectively, we indicated that ANT homolog gene (Gglean003340, Gglean011480) and INO homolog gene (Gglean016270) may act as the downstream genes that lead to ovule abortion for pomegranate FMF. AG-Like and SPL homology genes may be upstream genes impact ANT and INO expression, but there may also exist other upstream genes repressing the expression of ANT and INO in pomegranate ovule abortion. Moreover, ethylene and nutrient may be the upstream regulators that regulated some genes related to ovule development, and the optimal spraying concentration of ethephon was 150 mg / L, which can be referenced for improving the yield of pomegranate.