Abstract
The number of seeds is an important factor limiting the palatability and processing efficiency of loquat. Loquat H30-6 has been found to bear fruits with few seeds, which is attributed to its low male fertility. To investigate the molecular mechanism underlying the low male fertility of H30-6, in this study, RNA-seq was performed using pre-meiotic, meiotic and mature anthers from H30-6 and the other two loquat varieties with normal male fertility. The results indicated that a total of 331 differentially expressed genes (DEGs) were identified between H30-6 and the other two loquats and these DEGs were mainly enriched in pathways related to jasmonic acid (JA). WGCNA revealed that module ‘green’ were highly positively correlated with seed quantity (r = 0.46), suggesting that module ‘green’ is the candidate module. Enrichment analysis showed that genes in module ‘green’ were also primarily enriched in JA biosynthetic process, regulation of JA mediated signaling pathway and response to JA, implying that JA may play a key role in the molecular mechanism underlying the low male fertility of H30-6. The detection of UPLC-MS/MS indicated that the contents of JA and its derivatives in the flower bud of H30-6 were significantly lower than those in ‘Huabai 1’ (with normal male fertility). Besides, exogenous application of MeJA improved the pollen quantity per anther and seed numbers of H30-6. Our work indicates that the suppression of JA-related gene expression and the low JA content in immature flower buds may be key factors for H30-6 male sterility. These findings provide a theoretical basis for seedless breeding of loquats and lay the foundation for the development of stamens in fruit trees.
1 Introduction
Loquat (Eriobotrya japonica Lindl.) is a popular subtropical fruit of the Rosaceae family, native to southwestern China (; Wang et al., 2017). Owing to its juiciness and nutrient-rich nature, this fruit is not only used for fresh consumption, but also in food industries for the production of juices, wines, jams, and syrups, as well as traditional Chinese medicine (; ; ). However, under natural conditions, loquat fruit bears among 1–5 large seeds, occupying 20–30% of fruit volume (; ), which affects its palatability and processing efficiency. To produce seedless loquat, breeders have made enormous efforts. Since the 1990s, breeding triploids has been one of the main methods to produce seedless loquats and several triploid seedless varieties have been released (; , ). However, to improve the fruit set of triploid loquat varieties, hormone treatments are necessary (), which is both time-consuming and labor-intensive (). Recently, diploid loquat H30-6, a descendant of white-fleshed variety ‘Huabai 1’, has been found to bear fruits with few seeds and is an important breeding material with low male fertility (Xia et al., 2022). This germplasm serves as a crucial resource for breeding seedless loquats and studying the seedless and male sterile mechanisms.
Our previous study has reported that the low male fertility of H30-6 is caused by abnormal development of pollen mother cells (Xia et al., 2022). Various factors are associated to the development of pollen, including transcription factors (Wang et al., 2023), non-coding RNAs (), reactive oxygen species () and plant hormones (; ). Among the plant hormones, jasmonic acid (JA) plays an essential role in stamen development and function. In Arabidopsis thaliana, JA biosynthesis and signaling are required in stamen development (; ; ). Numerous researches have declared that the JA biosynthetic mutants, including opr3 (), aos/dde2 (, acx1 acx5 (), lox3 lox4 (), and gh3.10-2 jar1-11 (), exhibit completely or severely male sterile, as manifested by failure of stamen filament elongation, delay of anther dehiscence, and loss of pollen vitality. Besides, the stamen development and fertility of these JA biosynthetic mutants can be rescued by exogenous application of JAs (Yan et al., 2016; ). Additionally, ZmCOI2a and ZmCOI2b, encoding the F-box protein COI1, a key component of the JA signaling pathway, redundantly regulate gametophytic male fertility in maize (). Furthermore, JA has been found to regulate stamen development in peach (), almond (Amygdalus communis L.) (), and a male-sterile somatic cybrid citrus (), suggesting that JA plays a conserved role in regulating stamen development in fruit crops. It is speculated that JA may play a role in pollen development in H30-6.
To investigate the molecular mechanism underlying the low male fertility of H30-6, in this study, transcriptomic analysis was performed with anther samples from three developmental stages (pre-meiosis, meiosis and mature pollen). Furthermore, weighted gene co-expression network analysis (WGCNA) was applied to identify the related co-expression gene modules with the pollen fertility. Results of this work provide new insights into the molecular mechanisms of the low male fertility of loquat H30-6 and lay a theoretical foundation for further seedless breeding in loquat.
2 Materials and methods
2.1 Plant samples
In this study, H30-6 with a small seed quantity and two varieties ‘Huabai 1’ (HB) and ‘Jinhua 1’ (JH) with large seed quantities were used as materials. Since HB is the parent of H30-6, it was used as a control with close genetic distance. JH, with orange-flesh, differs from H30-6 in flesh color, serving as a control with distant genetic distance. All the materials were obtained from the fruit germplasm resources of the Key Fruit Lab, of the College of Horticulture and Landscape Architecture, Southwest University. According to the classification criteria of Xia et al (Xia et al., 2022). for flower buds: flower buds of H30-6 and JH with a diameter of less than 3 mm and those of HB with a diameter of less than 2 mm can be considered as having microspore mother cells in the anthers undergoing prophase of meiosis (pre-meiotic anthers); the flower buds of H30-6 and JH with a diameter of 3 to 4 mm and those of H411 with a diameter of 2 to 3 mm can be regarded as the microspore mother cells in the anthers being in the meiosis stage (meiotic anthers); flower buds of H30-6 and JH with diameters greater than 4 mm, and those of HB with diameters greater than 3 mm, possess mature anthers. Pre-meiotic, meiotic and mature anthers were collected from H30-6, HB and JH. Three replicates were set for each sample, and each replicate included 10 inflorescences. The separated stamens were frozen in liquid nitrogen for 15 min and stored in a refrigerator at -80°C for subsequent analysis.
2.2 RNA-seq
Total RNA was extracted using the BASY spin Plant Total RNA Extraction Kit (RN90, Aidlab Biotechnologies Co.,Ltd) according to the manufacturer’s instructions. RNA purity and RNA content were examined using the Agilent 2100 Bioanalyzer. The mRNA in the total RNA was enriched and purified by Oligo (dT) magnetic bead method, and then the mRNA was interrupted to 200-300 bp fragment. The first cDNA strand was synthesized with 6-base random primers (N6) and reverse transcriptase (M-MLV). Sequencing library construction and sequencing were performed by Illumina NovaSeq X Plus (PE 150) platform (Personalbio, Shanghai, China). After removing the adaptors and low-quality reads, the clean reads were then aligned to the loquat reference genome () using HISAT2 (). To estimate the expression levels of genes, the read counts for each gene were calculated by HTSeq (), and then normalized to fragments per kilobase million (FPKM) mapped reads. Differentially expressed genes were identified using the fold change factor (|Log2(FC)| > 1) and the adjusted p-value threshold (p < 0.05).
2.3 Valication by real-time quantitative PCR
qRT-PCR was used to verify the RNA-seq results using the same samples for RNA-seq. Gene-specific primers, designed on the Primer3Plus website (https://www.primer3plus.com/), were list in Supplementary Table S1. qTOWER3 G touch (Analytic Jena, German) was used to perform qRT-PCR. EjActin, as a housekeeping gene, was used to normalize the amplification. The relative expressions of genes were calculated using the 2-ΔΔCt method ().
2.4 WGCNA
According to Langfelder and Horvath (), WGCNA was performed using R (4.3.1) with the default parameters and the soft-thresholding power (β) of 8. The top 75% of genes with the highest median absolute deviation in each sample were filtered for subsequent calculations. Based on the correlations between genes (Pearson’s correlation tests), genes were clustered into different modules. Additionally, the relationship between each module and sample was calculated (Pearson’s correlation tests) to construct topological overlap matrix. The top hub gene in each module was selected based on the highest kME (eigengene connectivity) value. Moreover, Cytoscape (3.10.3) was used for the co-expression network ().
2.5 UPLC-MS/MS analysis of JA and its derivatives
The flower buds with immature anthers of H30-6 and HB were collected to detect the contents of JA and its derivatives using ultra-performance liquid chromatography (UPLC: Agilent 1290) coupled with tandem mass spectrometer (MS) system (MS/MS: Applied Biosystems 6500 Quadrupole Trap). The samples were ground into powder in liquid nitrogen. Next, 50 mg of the powder was dissolved in 0.5 mL of a mixture of isopropyl alcohol, water, and hydrochloric acid (volume ratio of 2:1:0.002), and shaken at 100 r/min on a shaker at 4°C for 30 minutes. Then, 1 mL of dichloromethane was added to the sample extract and shaken for another 30 minutes. After centrifuging at 13,000 × g for 5 min, the lower layer solution of the extract was dried with a nitrogen blower. Subsequently, the sample was dissolved in 0.1 mL of methanol solution, centrifuged at 13,000 g for 5 min, and the supernatant was taken for UPLC-MS/MS analysis. The standard solutions of JA, MEJA, 12-OPDA, and JA-ILE with gradients of 0.1, 0.2, 0.5, 2, 5, 20, 50, and 200 ng/mL were prepared using methanol (0.1% formic acid) as the solvent.
The chromatographic conditions were as follows: 2.7 µm, 2.1 × 150 mm (Agilent Poroshell 120 SB-C18). The column temperature was set to 40°C. The mobile phase consisted of solvent A, methanol containing 0.1% formic acid, and solvent B, 0.1% formic acid aqueous solution. The detection adopted a gradient elution program: 20% A and 80% B for 3 min, 50% A and 50% B for 6 min, 90% A and 10% B for 1.6 min, 20% A and 80% B for 2.9 min. The injection volume was 2 µL, and the flow rate was 0.2 mL/min.
The mass spectrometry conditions were as follows: electrospray ionization source temperature, 400°C; ionspray voltage, 4500 V; curtain gas, ion source gas I, and gas II set at 15, 65, and 70 psi, respectively; monitoring mode, multiple reaction monitoring. In the Q-Trap6500 system, each ion pair was scanned based on the optimized declustering potential and collision energy. Analyst software (v 1.6.3) was used to process mass spectrometry data.
2.6 Observation of pollen quantity and pollen germination
According to the previous study (), the inflorescences of H30-6, containing microspore mother cells in the pre-meiotic stage, were sprayed with 2 mM methyl jasmonate (MeJA) solution or mock (water) every week until the microspore mother cells entered the post-meiosis stage (approximately one month). When blossom, pollen quantity and pollen germination of H30-6, treated with MeJA and mock, and HB were detected. Fifteen dried anthers were crushed using tweezers in 200 μL of BK medium, containing 150 g/L sucrose, 100 mg/L boric acid and 1 g/L agar, with 6 replications. Pollen quantity was counted in a blood counting chamber. Pollen germination was tested according to our previous study (Xia et al., 2022). The dried anthers were placed in 200 μL of BK medium and shaken slightly to disperse the pollen. Then, the anthers were removed. After 45 min of incubation at room temperature, 50 μL of pollen suspension was transferred into a new 1.5 mL centrifuge tube and incubated upside down for 4 h at 25°C. After being shaken, 15 μL of pollen suspension was dropped on a glass slide for the observation and imaging using an Olympus (BX35) fluorescence microscope (Tokyo, Japan). Six replicates were set for each sample. At least three photographs of different fields of view were taken per replicate and germinated and total pollen were counted.
2.7 Statistical analysis
SPSS (v22) software was used to conduct statistical analysis. To determine statistical significance among groups, one-way analysis of variance (ANOVA) was adopted when the variances were homogeneous; otherwise, the Welch’s test was selected. For differences between two groups, the independent-samples t-test was used.
3 Results
3.1 The number of seeds in the fruits of HB, JH, and H30-6
Mature fruits of HB, JH, and H30-6 were collected to count their seeds (Figure 1A). For the number of normal seeds, HB had the most seed quantities with an average of 4.25 seeds per fruit, significantly higher than the average of 2.95 seeds per fruit for JH and 1.45 seeds per fruit for H30-6 (p < 0.05) (Figure 1B). However, the numbers of abnormal seeds (longitudinal diameter < 1 cm) of H30-6 (Supplementary Figure S1) and JH were higher than that of HB (p < 0.05). Besides, there was no difference in the total number of seeds, including normal and abnormal seeds, between H30-6 and HB (Figure 1B). The number of seeds mentioned below refers to the quantity of normal seeds.
Figure 1
3.2 RNA-seq and identification of DEGs
The low male fertility of H30-6 pollen is a limiting factor for the fruit setting (Xia et al., 2022). To reveal the molecular mechanism underlying the low male fertility of H30-6, the RNA-seq was performed with pre-meiotic, meiotic and mature anthers collected from H30-6, HB and JH. For the sample names of H30-6, HB and JH, the pre-meiotic anthers were named H306-1, HB-1 and JH-1; the meiotic anthers were named H306-2, HB-2 and JH-2; similarly, the mature anthers were named H306-3, HB-3 and JH-3. Using the Illumina platform, an average of about 42 million clean reads were generated per sample. The average mapping rate was 96.03% after aligning to the loquat reference genome. A total of 45,743 genes were detected, of which 31,780 genes were expressed (FPKM > 1) in at least one sample and 12,611 were expressed in all samples. A total of 12 genes were selected randomly to perform qRT-PCR for verifying the RNA-seq results. The results showed that the relative expression trends of these 12 genes were consistent with the RNA-seq results, and the determination coefficients (R2) were all greater than 0.6, supporting the reliability of the RNA-seq results (Figure 2).
Figure 2
To evaluate the expression patterns of DEGs among samples, hierarchical clustering analysis were conducted. According to the heatmap, DEGs were clustered separately by period and clearly divided into two groups, namely mature and immature anther groups (Figure 3A). Compared with HB and JH, the number of up-regulated genes in H30-6 was lower than that of down-regulated genes at different periods (Figure 3B). Besides, the DEG number in mature anther group was higher than that in immature anther group in this comparison. For pre-meiotic anthers, 1,975 and 426 genes were up-regulated and 3,222 and 2,256 genes were down-regulated in H306-1 by comparing with HB-1 and JH-1, respectively. Regarding the meiotic anthers, there were 2,376 and 1,001 up-regulated genes and 2,928 and 1,951 down-regulated genes in H306-2 in comparison with HB-2 and JH-2, respectively. By comparison to HB-3 and JH-3, 5,188 and 5,691 genes were up-regulated, and 7,508 and 7,606 genes were down-regulated in H306-3 (Figure 3B).
Figure 3
3.3 The function analysis of DEGs in immature anthers
Since the male sterile of H30-6 is attributed to abnormal meiotic synapses (Xia et al., 2022), the pre-meiosis and meiosis phases of microspore mother cells are critical for the pollen development in H30-6. In pre-meiotic and meiotic anthers, 31 up-regulated DEGs and 280 down-regulated DEGs were shared in the comparisons between H30-6 and HB and between H30-6 and JH. (Figure 4A). To explore the function of the 311 shared DEGs, Gene Ontology term (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. For GO enrichment analysis, these genes were enriched in biological process (BP), cell component (CC), and molecular function (MF), respectively. BP terms were mainly enriched for genes related to JA, including regulation of JA mediated signaling pathway and JA biosynthetic process. CC terms were predominantly enriched for genes associated with plasma membrane and MF terms were primarily enriched for genes involved in protein serine kinase activity (Figure 4B). KEGG enrichment analysis indicated that the 311 shared DEGs were mainly enriched in zeatin biosynthesis and glycosphingolipid biosynthesis-globo and isoglobo series (Figure 4C). Zeatin is a precursors to the synthesis of cytokinin, which ooperates with JA to regulate stamen development (; Yu et al., 2019). Glycosphingolipids, a type of sphingolipid, are key components of the cell wall, and JA is involved in the metabolism of sphingolipid (). These findings indicate that JA may participate in the mechanism of the low male fertility of H30-6. Among the JA-related DEGs, compared with HB, EjAOS3 (EVM0010750) had the lowest log2(fold change) values in pre-meiotic and meiotic anthers of H30-6, which were -4.32 and -5.50, respectively, followed by EjAOC4.1 (EVM0005305), which were -2.78 and -5.03, respectively.
Figure 4
3.4 Weighted gene co-expression network analysis
To investigate the relationship between gene networks with the low male fertility of H30-6, WGCNA was performed to construct the co-expressed gene modules with a total of 22,779 expressed genes. All these genes were clustered into 8 modules (Figure 5A). The correlation between gene modules and samples revealed that module ‘green’ was highly related with HB-2 (r = 0.57) and HB-1 (r = 0.54), and subsequently with JH-1 (r = 0.096) and JH-2 (r = 0.063) (Figure 5B), indicating that module ‘green’ may be related to the meiosis behavior of their pollen mother cells. Besides, module ‘green’ were highly positively correlated with seed quantity (r = 0.46) (Figure 5B), suggesting that module ‘green’ is the candidate module associated with loquat seed number. GO term and KEGG pathway enrichment analyses were performed to explore the function of the genes in module ‘green’, a total of 398 genes. For GO enrichment analysis, BP terms were also primarily enriched for genes related to JA, including JA biosynthetic process, regulation of JA mediated signaling pathway and response to JA. CC terms were mainly enriched for genes associated with plasma membrane and MF terms were predominantly enriched for genes involved in protein binding (Figure 6A). Besides, KEGG enrichment revealed that genes in module ‘green’ were enriched in alpha-linolenic acid metabolism (Figure 6B), which is the previous step of JA biosynthesis. These results further suggest that JA play an important role in the molecular mechanism underlying the low male fertility of H30-6.
Figure 5
Figure 6
3.5 Networks of genes related to JA biosynthesis and signaling pathway
To further explore the interaction between genes related to JA in the ‘green’ module, the interaction network was plotted for the co-expressed genes connected to genes of JA biosynthesis and signaling pathways with highest intramodular connectivity (KIM > 0.3) (Figure 7A). EVM0033600 (EjNPR3) and EVM0027375 (EjLOX2.2) have the top two highest degrees, 76 and 72, respectively. NPR3 and LOX2 are important components of JA signaling and biosynthesis, respectively (Figure 7B). Expression heatmap revealed that the expressions of genes related to JA biosynthesis (LOXs, AOS, AOCs, OPRs, and OPCL1), JA signaling (JAZs, MED16, MYC2s, and NPR3), and JA response (EBP) in both pre-meiotic and meiotic anthers of HB and JH were higher than those in their mature anthers (Figures 7B, C). Notably, during pre-meiotic and meiotic period, the expression of these genes in H30-6 was the lowest among the three varieties (Figure 7C), which was consistent with their male fertilities (Xia et al., 2022). Especially in the meiotic period, the expression levels of EjAOC4.1 and EjJAZ8.2 were 31.76- and 24.17-fold lower in H30-6 than in HB, respectively. Furthermore, the relative expression levels of the 19 genes related to JA biosynthesis and signaling pathways in pre-meiotic and meiotic anthers of H30-6 and HB have been verified by qRT-PCR. The results indicated that, except for EjMYC2.2 and EjEBP, the remaining 17 genes in pre-meiotic anthers had significantly differences lower relative expression levels in H30-6 than HB (p < 0.05). Particularly, the relative expression level of EjOPR3.2 in HB was 60 times that of H30-6. In meiotic anthers, the relative expression levels of 13 genes out of the 19 genes were significantly different between H30-6 and HB (p < 0.05) (Figure 8). These results suggested that the differential expression of genes related to JA biosynthesis and signaling pathways between H30-6 and HB in pre-meiotic and meiotic anthers may be associated with their fertility.
Figure 7
Figure 8
3.6 Differences in the content of JA and its derivatives in H30-6 and HB flower buds
To explore the function of JA in the pollen development of H30-6, the contents of JA and its derivatives, including methyl jasmonate (MeJA), jasmonoyl-isoleucine (JA-Ile), and 12-Oxo-phytodienoic acid (OPDA), were detected using UPLC-MS/MS in H30-6 and HB flower buds (Figure 9A). It was found that the contents of JA and its derivatives in H30-6 were significantly lower than those in HB (p < 0.01) (Figure 9B). The contents of JA, MeJA, JA-Ile, and 12-OPDA in H30-6 were 42.79%, 14.19%, 51.43%, and 18.74% lower than those in HB, respectively. These findings imply that the contents of JA and its derivatives in flower buds may be related to the male fertilities of H30-6 and HB.
Figure 9
3.7 Effects of MeJA on the pollen quantity and activity of H30-6
To evaluate the effects of MeJA on the pollen quantity and activity of H30-6, the flower buds of H30-6 were sprayed weekly with MeJA, from anthers in the pre- to post-meiotic stage. Microscopic observation revealed that the numbers of pollen grains and pollen germinations of HB were obviously higher than that of H30-6 treated with MeJA and water (Mock), at the same field-of-view magnification (Figures 10A, B). Statistically, pollen quantity per anther of MeJA group was remarkable higher than that of Mock (water) group, but lower than that of HB (p < 0.05), and MeJA treatment increased pollen quantity per anther by 40.11% (Figure 10C). Besides, the pollen germination rate of HB was significantly higher than that of H30-6 in both MeJA and Mock groups (p < 0.05), while, there was no difference between MeJA group and Mock group (Figure 10D). These findings indicate that MeJA may improve the pollen quantity per anther of loquat.
Figure 10
3.8 Effect of MeJA on the seed quantity of H30-6
Fruits in MeJA and Mock groups were collected at maturity to count the number of seeds. Cross-sections of the fruits (Figure 10E) and statistics of seed numbers (Figure 10F) indicated that the average number of seeds per fruit in the MeJA group was 1.8, which was significantly higher than the average of 1.4 in Mock group (p < 0.05). MeJA treatment increased seed numbers by 28.57% (Figure 10F), which was lower than the increase in pollen quantity (40.11%) after MeJA treatment. It suggests that, in addition to pollen quantity, there may be other factors that contribute to the low seed yield of H30-6, such as gynoecium size and ovule number (; Yu et al., 2022). These results imply that MeJA has a positive effect on seed quantity of H30-6.
4 Discussion
4.1 Meiosis abnormality is the key factor contributing to H30-6 male sterility
Male sterility is an important breeding material to product seedless fruits in fruit trees and vegetables, as it is able to reduce costs and improve seed purity (; Zhou et al., 2023; b). In plants, male sterility usually manifests as an inability to produce functional pollen or abnormal anther development, resulting in an inability to pollinate and reproduce (Zhou et al., 2022). The development of anthers in flowering plants involves a highly delicate and complex process, and even slight abnormalities in this process can lead to male infertility (). According to phenotype, forms of male sterility can be divided into tapetal cell developmental abnormality, meiotic arrest of pollen mother cells, incomplete anther development, pollen abortion, and non-functional pollen types (). Our previous work has revealed that the abnormal meiosis of pollen mother cells (Xia et al., 2022), is the key factor contributing to male sterility in H30-6, and there are phenomena of reduced pollen quantity and low pollen germination rate.
Meiosis abnormality is a general cause of seedlessness in fruit crops. A 2.09-Mb fragment translocation was found to cause abnormalities during meiosis, resulting in the formation of a seedless watermelon (). By targeting and knocking out the genes ClREC8 and ClSPO11, which are essential for meiosis, seedless watermelons has been successfully generated using gene editing technology (; ). Abnormal meiosis led to the male sterility in the seedless hybrid of citrus (). Additionally, somatic mutations in the gene MER3 related to meiosis were found to cause male sterility in the seedless citrus mutant (). Therefore, the identification of key genes involved in meiosis will contribute to seedless breeding of loquats and other fruit crops.
4.2 The suppression of JA-related gene expression during anther development contributes to the low male fertility of H30-6
Our work revealed that the suppression of JA-related gene expression and low JA content in immature flower buds may be key factors in H30-6 male sterility, which is in accordance with the important function of JA in regulating stamen development (; ; ). JA plays a core role in coordinating pollen maturation, filament elongation, and anther dehiscence (). Numerous reports indicate that the mutation of genes in the JA pathway leads to male sterility. The homologous genes of EjLOX3, EjAOS, and EjOPR3.2 have been found to be associated with male fertility in Arabidopsis. Four independent lox3 lox4 double mutants were male sterile, displaying abnormal anther maturation and cracking defects (). The Arabidopsis mutants of aos and opr3 also exhibited severe male sterility and no seed production (; ). Since the JA pathway is conserved across different species, similar mutants in species beyond Arabidopsis have been screened, such as the rice JA-deficient mutants cpm1/osaos1 (coleoptile photomorphogenesis 1) (), cpm2/hebiba/osaoc (), and osjar1 (Xiao et al., 2014), exhibiting impaired anther dehiscence or abnormal anthesis resulting in partial male sterility (). However, similar mutants have not been identified in fruit crops. Additionally, MYC2 is a core transcription factor in JA signaling and plays an essential role in plant stamen development (). MYC2/MYC4 have been reported to enhance anther dehiscence by activating the expression of NST1 under blue light conditions (Zhang et al., 2018). NtMYC2a mediated process of pollen maturation in tobacco (Nicotiana tabacum L. cv. TN90) by regulating the starch metabolism in the pollen grains, anther walls and ovaries ().
In immature flower buds of H30-6, the expression levels of genes related to JA biosynthesis and signaling pathways as well as JA content were significantly lower than those in HB and JH. These results imply that, compared with HB and JH, natural variation may exist in the JA upstream synthesis genes (e.g., EjLOXs) or transcription factors regulating JA synthesis such as AP2/ERF (), NAC (), and WRKY () in H30-6. The potential reasons for the differences in JA signal sensitivity among H30-6, HB and JH need further exploration.
4.3 Exogenous application of MeJA increases the pollen quantity and seed numbers of H30-6
A number of studies indicated that reduced fertility caused by JA synthesis defects can be rescued by MeJA treatment (; Yuan and Zhang, 2015). Severe male sterility was observed after knock-out of AOS in Arabidopsis, which was rescued by exogenous application of MeJA (). In cotton, exogenous application of MeJA to early buds of Ghaoc2 mutant lines rescued the phenotypes of sterile pollen and indehiscent anther (). Besides, the application of MeJA significantly upregulated the pollen fertility of photoperiod-sensitive genic male sterile rice, resulting in partial fertility and even reversing the expected natural sterility (). The treatment of MeJA induced the expression of LoMYB26, promoting endothecium lignification and anther dehiscence in lily (). In this study, the treatment of MeJA improved the pollen quantity per anther (Figure 10C) and seed numbers (Figure 10F) of H30-6, but not the pollen germination rate (Figure 10D).
Although JA plays an important role in coordinating pollen maturation, filament elongation, and anther dehiscence (), its effect on pollen quantity remains unclear. We suppose that, on the one hand, MeJA treatment may promote the anther dehiscence of H30-6, leading to an increase in pollen quantity; on the other hand, MeJA treatment may have rescued the meiosis behavior of H30-6 pollen mother cells, thereby increasing the number of normal pollen grains. Multiple essential genes for meiosis have been identified in plants, including ClREC8, ClSPO11 (; ) and ClDYAD (Wu et al., 2025) in watermelon, MER3 in citrus (), and PHD transcription factor gene PP1 in maize (). More meiotic mutants have been characterized in rice, maize, tomato, wheat, Brassica, and Arabidopsis (Wang et al., 2021). Currently, there is limited research on the direct regulatory role of JA in pollen mother cell meiosis. MYC2, a core transcription factor in JA signaling, plays an essential role in plant stamen development by transcriptional regulation of downstream genes (). The MYB-MYC complex formed by the interaction of MYC2 with MYB21 and MYB24 cooperatively mediates stamen development and seed formation (). Besides, MYC transcription factors redundantly mediate JA-induced expression of MYB genes (MYB21, MYB24, MYB57, and MYB108), thereby regulating stamen development (). In rice, GAMYB has been reported to mediate meiosis in pollen mother cells by regulating the expression of bHLH142 (). It is speculated that JA may target these crucial genes for meiosis through the MYB-MYC complex to regulate H30-6 rescue of male sterility, but further verification is still needed. Previous studies have found that MeJA affects pollen germination in a dose-dependent manner. It has been reported that 0.1 and 0.25 mM MeJA enhanced while 0.5- and 1 mM MeJA decreased the pollen germination in Prunus armeniaca L (). Similarly, 0.05 mM and 0.25 mM MeJA significantly increased the total germination rate of Pinus nigra pollen, while 0.5 mM MeJA had no significant effect (). Therefore, the MeJA treatment concentration used (2 mM) in this study may not be suitable for loquat pollen germination, resulting in no significant improvement in pollen germination rate.
In fruit crops, the role of JA in stamen development has been studied. In peach, JA pathway has been reported to be involved in flower development and pollination (). JA has also been verified to promote anther development in almond (Amygdalus communis L.) (). In a male-sterile somatic cybrid citrus, JA and indole-3-acetic acid jointly influenced the expression of genes related to stamen development and interfered with the morphology of stamens and flowers (). These studies suggest conservative role of JA in regulating stamen development in fruit crops, supporting that JA takes part in the low male fertility of loquat H30-6.
5 Conclusion
To investigate the molecular mechanism underlying the low male fertility of H30-6, in this study, RNA-seq and WGCNA were performed. Both RNA-seq and WGCNA analyses revealed that suppression of JA-related gene expression during pollen development contributes to the low male fertility of H30-6. UPLC-MS/MS and exogenous application of MeJA also confirmed the importance of JA content in loquat pollen development. Our work suggests that JA is an essential factor in regulating male fertility in loquats, providing a theoretical basis for seedless breeding of loquats and laying the foundation for the development of stamens in fruit trees.
Statements
Data availability statement
RNA‐seq data of pre-meiotic, meiotic and mature anthers collected from three loquats have been deposited into the NCBI (National Center for Biotechnology Information) database with accessions number PRJNA1283964.
Author contributions
JW: Methodology, Data curation, Writing – original draft. JP: Methodology, Data curation, Writing – original draft. QX: Methodology, Writing – original draft. YD: Writing – original draft, Methodology. DJ: Supervision, Writing – review & editing. JD: Writing – review & editing, Conceptualization, Supervision, Writing – original draft. QG: Supervision, Funding acquisition, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research and/or publication of this article. This work was supported by the National Key R&D Program of China (2023YFD1600800 and 2024YFD2100601), the Chongqing Science and Technology Commission (cstc2024ycjh-bgzxm0202 and cstc2021jscx-gksbX0010), and the Characteristic Fruit Industry and Technology System Innovation Team of Chongqing Agriculture and Rural Affairs Commission [No. 2020 (3) 01 and No. CQMAITS20230505].
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Generative AI statement
The author(s) declare that no Generative AI was used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1648313/full#supplementary-material
References
1
AcostaI. F.PrzybylM. (2019). Jasmonate signaling during Arabidopsis stamen maturation. Plant Cell Physiol.60, 2648–2659. doi: 10.1093/pcp/pcz201
2
AndersS.PylP. T.HuberW. (2015). HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics31, 166–169. doi: 10.1093/bioinformatics/btu638
3
BadenesM. L.JanickJ.LinS.ZhangZ.LiangG. L.WangW. (2013). Breeding loquat. Plant Breed. Rev.37, 259–296. doi: 10.1002/9781118497869.ch5
4
BianS.TianT.DingY.YanN.WangC.FangN.et al. (2022). bHLH transcription factor NtMYC2a regulates carbohydrate metabolism during the pollen development of tobacco (Nicotiana tabacum L. cv. TN90). Plants11, 17. doi: 10.3390/plants11010017
5
BiswasK. K.NeumannR.HagaK.YatohO.IinoM. (2003). Photomorphogenesis of rice seedlings: a mutant impaired in phytochrome-mediated inhibition of coleoptile growth. Plant Cell Physiol.44, 242–254. doi: 10.1093/pcp/pcg040
6
CaiH.ZhangS.WangY.YangZ.ZhangL.ZhangJ.et al. (2025). Anther-specific expression of MsMYB35 transcription factor in alfalfa (Medicago sativa L.) and its crucial role in pollen development. Plant J.122, e70126. doi: 10.1111/tpj.70126
7
CaldelariD.WangG.FarmerE. E.DongX. (2011). Arabidopsis lox3 lox4 double mutants are male sterile and defective in global proliferative arrest. Plant Mol. Biol.75, 25–33. doi: 10.1007/s11103-010-9701-9
8
CaoL.LiC.LiH.WangZ.JiangY.GuoY.et al. (2022). Disruption of REC8 in Meiosis I led to watermelon seedless. Plant Sci.323, 111394. doi: 10.1016/j.plantsci.2022.111394
9
ChenH. Y.HsiehE. J.ChengM. C.ChenC. Y.HwangS. Y.LinT. P. (2016). ORA47 (octadecanoid-responsive AP2/ERF-domain transcription factor 47) regulates jasmonic acid and abscisic acid biosynthesis and signaling through binding to a novel cis-element. New Phytol.211, 599–613. doi: 10.1111/nph.13914
10
ChenL.LiuY. G. (2014). Male sterility and fertility restoration in crops. Annu. Rev. Plant Biol.65, 579–606. doi: 10.1146/annurev-arplant-050213-040119
11
ChengZ.SongW.ZhangX. (2022). Genic male and female sterility in veget able crops. Hortic. Res.10, uhac232. doi: 10.1093/hr/uhac232
12
CucinottaM.Di MarzoM.GuazzottiA.De FolterS.KaterM. M.ColomboL. (2020). Gynoecium size and ovule number are interconnected traits that impact seed yield. J. Exp. Bot.71, 2479–2489. doi: 10.1093/jxb/eraa050
13
CuriP. N.NogueiraP. V.AlmeidaA. B. D.CarvalhoC. D. S.PioR.PasqualM.et al. (2017). Processing potential of jellies from subtropical loquat cultivars. Food Sci. Technol.37, 70–75. doi: 10.1590/1678-457X.07216
14
DangJ.GuoQ.XiangS.HeQ.SunH.JingD.et al. (2019a). A new red-fleshed and seedless loquat cultivar with large fruit ‘Huajin Wuhe 1’. Acta Hortic. Sin.46, 2763–2764. doi: 10.16420/j.issn.0513-353x.2019-0241
15
DangJ.GuoQ.XiangS.HeQ.SunH.WuD.et al. (2019b). A new seedless loquat cultivar ‘Huayu Wuhe 1’ with large fruit. Acta Hortic. Sin.46, 2765–2766. doi: 10.16420/j.issn.0513-353x.2019-0242
16
DelfinJ. C.KannoY.SeoM.KitaokaN.MatsuuraH.TohgeT.et al. (2022). AtGH3.10 is another jasmonic acid-amido synthetase in Arabidopsis thaliana. Plant J.110, 1082–1096. doi: 10.1111/tpj.15724
17
DhimanA.SuhagR.ThakurD.GuptaV.PrabhakarP. K. (2022). Current status of loquat (Eriobotrya japonica Lindl.): Bioactive functions, preservation approaches, and processed products. Food Rev. Int.38, 286–316. doi: 10.1080/87559129.2020.1866007
18
FanY. J.DuZ. Z.HeX. Y.LiuZ. A.ZhuangJ. X.XiaoG. A.et al. (2025a). Somatic variations in the meiosis-specific gene CrMER3 confer seedlessness in a citrus bud sport. J. Integr. Plant Biol.67, 1649–1664. doi: 10.1111/jipb.13872
19
FanY. J.TanF. Q.HeX. Y.LiuZ. A.BaiH. Z.JiangN.et al. (2025). CitACOS5 is crucial for male fertility in citrus by influencing pollen exine formation. J. Exp. Bot., eraf181. doi: 10.1093/jxb/eraf181
20
FangY.GuoD.WangY.WangN.FangX.ZhangY.et al. (2024). Rice transcriptional repressor OsTIE1 controls anther dehiscence and male sterility by regulating JA biosynthesis. Plant Cell36, 1697–1717. doi: 10.1093/plcell/koae028
21
FangQ.PanX.WuZ.YuJ.LiT.XuS.et al. (2025). The LoMYB26/LoJAZ4-LoCOMT module regulates anther dehiscence via jasmonic acid-mediated endothecium lignification in lily. J. Adv. Res.005. doi: 10.1016/j.jare.2025.04.005
22
Getinbaş-GençA.VardarF. (2020). Effect of methyl jasmonate on in-vitro pollen germination and tube elongation of Pinus nigra. Protoplasma257, 1655–1665. doi: 10.1007/s00709-020-01539-4
23
GuoQ.LiX.XiangS.HeQ.SunH.WuD.et al. (2016). A new white pulp seedless loquat cultivar ‘Wuhe Guoyu’. Acta Hortic. Sin.43, 2717–2718. doi: 10.16420/j.issn.0513-353x.2016-0445
24
HeY.LiuC.ZhuL.FuM.SunY.ZengH. (2021). Jasmonic acid plays a pivotal role in pollen development and fertility regulation in different types of P(T)GMS rice lines. Int. J. Mol. Sci.22, 7926. doi: 10.3390/ijms22157926
25
HuangL. Q.ChenD. K.LiP. P.BaoH. N.LiuH. Z.YinJ.et al. (2021). Jasmonates modulate sphingolipid metabolism and accelerate cell death in the ceramide kinase mutant acd5. Plant Physiol.187, 1713–1727. doi: 10.1093/plphys/kiab362
26
HuangH.ChenY.WangS.QiT.SongS. (2022). Jasmonate action and crosstalk in flower development and fertility. J. Exp. Bot.74, 1186–1197. doi: 10.1093/jxb/erac251
27
HuangH.LiuB.LiuL.SongS. (2017). Jasmonate action in plant growth and development. J. Exp. Bot.68, 1349–1359. doi: 10.1093/jxb/erw495
28
JavedT.GaoS. J. (2023). WRKY transcription factors in plant defense. Trends Genet.39, 787–801. doi: 10.1016/j.tig.2023.07.001
29
JiangS.AnH.XuF.ZhangX. (2020). Chromosome-level genome assembly and annotation of the loquat (Eriobotrya japonica) genome. GigaScience9, giaa015. doi: 10.1093/gigascience/giaa015
30
JiangN.FengM. Q.ChengL. C.KuangL. H.LiC. C.YinZ. P.et al. (2023). Spatiotemporal profiles of gene activity in stamen delineate nucleo-cytoplasmic interaction in a male-sterile somatic cybrid citrus. Hortic. Res., uhad105. doi: 10.10.1093/hr/uhad105
31
JiangJ.FengQ.ZhaoZ.LiuQ.LiuM.WangJ.et al. (2024). Establishing a highly efficient diploid seedless watermelon production system through manipulation of the SPOROCYTELESS gene. New Phytol.244, 1128–1136. doi: 10.1111/nph.20108
32
KhanA. H.MinL.MaY.WuY.DingY.LiY.et al. (2020). High day and night temperatures distinctively disrupt fatty acid and jasmonic acid metabolism, inducing male sterility in cotton. J. Exp. Bot.71, 6128–6141. doi: 10.1093/jxb/eraa319
33
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements. Nat. Methods12, 357–360. doi: 10.1038/nmeth.3317
34
KoS. S.LiM. J.HoY. C.YuC. P.YangT. T.LinY. J.et al. (2021). Rice transcription factor GAMYB modulates bHLH142 and is homeostatically regulated by TDR during anther tapetal and pollen development. J. Exp. Bot.72, 4888–4903. doi: 10.1093/jxb/erab190
35
LangfelderP.HorvathS. (2008). WGCNA: an R package for weighted correlation network analysis. BMC Bioinf.9, 559. doi: 10.1186/1471-2105-9-559
36
LiS.HuY.YangH.TianS.WeiD.TangQ.et al. (2023). The regulatory roles of MYC TFs in plant stamen development. Plant Sci.333, 111734. doi: 10.1016/j.plantsci.2023.111734
37
LiP.TianJ.GuoC.LuoS.LiJ. (2021). Interaction of gibberellin and other hormones in almond anthers: phenotypic and physiological changes and transcriptomic reprogramming. Hortic. Res.8, 94. doi: 10.1038/s41438-021-00527-w
38
LiuJ.XiaC.DongH.LiuP.YangR.ZhangL.et al. (2022). Wheat male-sterile 2 reduces ROS levels to inhibit anther development by deactivating ROS modulator 1. Mol. Plant15, 1428–1439. doi: 10.1016/j.molp.2022.07.010
39
LuoY.LiuM.CaoJ.CaoF.ZhangL. (2022). The role of salicylic acid in plant flower development. For. Res.2, 14. doi: 10.48130/FR-2022-0014
40
MarciniakK.PrzedniczekK. (2019). Comprehensive insight into gibberellin- and jasmonate-mediated stamen development. Genes10, 811. doi: 10.3390/genes10100811
41
MesejoC.ReigC.Martínez-FuentesA.AgustíM. (2010). Parthenocarpic fruit production in loquat (Eriobotrya japonica Lindl.) by using gibberellic acid. Sci. Hortic.126, 37–41. doi: 10.1016/j.scienta.2010.06.009
42
MuradoğluF.YıldızK.BaltaF. (2010). Methyl jasmonate influences of pollen germination and pollen tube growth of apricot (Prunus Armeniaca L.). Yüzüncü Yıl Univ. J. Agric. Sci. (Turkey)20, 183–188. Available online at: https://www.academia.edu/22098010.
43
NieH.ChengC.KongJ.LiH.HuaJ. (2023). Plant non-coding RNAs function in pollen development and male sterility. Front. Plant Sci.14. doi: 10.3389/fpls.2023.1109941
44
PanL.LiY.AnX.XinZ.DongZ.JinW.et al. (2025). Maize prolonged prophase 1 is essential for meiotic progression in pollen mother cells. Crop J.13, 170–180. doi: 10.1016/j.cj.2024.12.001
45
ParkJ. H.HalitschkeR.KimH. B.BaldwinI. T.FeldmannK. A.FeyereisenR. (2002). A knock-out mutation in allene oxide synthase results in male sterility and defective wound signal transduction in Arabidopsis due to a block in jasmonic acid biosynthesis. Plant J.31, 1–12. doi: 10.1046/j.1365-313X.2002.01328.x
46
QiX.GuoS.WangD.ZhongY.ChenM.ChenC.et al. (2022). ZmCOI2a and ZmCOI2b redundantly regulate anther dehiscence and gametophytic male fertility in maize. Plant J.110, 849–862. doi: 10.1111/tpj.15708
47
QiT.HuangH.SongS.XieD. (2015). Regulation of jasmonate-mediated stamen development and seed production by a bHLH-MYB complex in Arabidopsis. Plant Cell27, 1620–1633. doi: 10.1105/tpc.15.00116
48
RamakrishnanR.BanadkaA.DubeyS.Al-KhayriJ. M.NagellaP. (2025). Advances in triploid plant production: techniques, benefits, and applications. Plant Cell Tiss. Org.160, 70. doi: 10.1007/s11240-025-03011-8
49
RaoZ.SunR.LiuS.AiW.SongL.WangX.et al. (2024). Abnormal transition from meiosis I to meiosis II induces male sterility in a seedless artificial hybrid of citrus. Mol. Breed.44, 1. doi: 10.1007/s11032-024-01521-5
50
RiemannM.HagaK.ShimizuT.OkadaK.AndoS.MochizukiS.et al. (2013). Identification of rice Allene Oxide Cyclase mutants and the function of jasmonate for defence against Magnaporthe oryzae. Plant J.74, 226–238. doi: 10.1111/tpj.12115
51
SadamatsuM.KeawladdagornV.OharaH.OhkawaK.MatsuiH. (2004). Induction of parthenocarpic fruit growth with endogenous gibberellins of loquat. Acta Hortic.653, 67–70. doi: 10.17660/ActaHortic.2004.653.8
52
SchilmillerA. L.KooA. J. K.HoweG. A. (2006). Functional diversification of acyl-coenzyme a oxidases in jasmonic acid biosynthesis and action. Plant Physiol.143, 812–824. doi: 10.1104/pp.106.092916
53
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al. (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res.13, 2498–2504. doi: 10.1101/gr.1239303
54
SherifS.El-SharkawyI.MathurJ.RavindranP.KumarP.PaliyathG.et al. (2015). A sta ble JAZ protein from peach mediates the transition from outcrossing to self-pollination. BMC Biol.13, 11. doi: 10.1186/s12915-015-0124-6
55
SongS.QiT.HuangH.XieD. (2013). Regulation of stamen development by coordinated actions of jasmonate, auxin, and gibberellin in Arabidopsis. Mol. Plant6, 1065–1073. doi: 10.1093/mp/sst054
56
StintziA.BrowseJ. (2000). The Arabidopsis male-sterile mutant, opr3, lacks the 12-oxophytodienoic acid reductase required for jasmonate synthesis. Proc. Natl. Acad. Sci. U.S.A.97, 10625–10630. doi: 10.1073/pnas.190264497
57
SuL.FangL.ZhuZ.ZhangL.SunX.WangY.et al. (2020). The transcription factor VaNAC17 from grapevine (Vitis amurensis) enhances drought tolerance by modulating jasmonic acid biosynthesis in transgenic Arabidopsis. Plant Cell Rep.39, 621–634. doi: 10.1007/s00299-020-02519-x
58
SuW.JingY.LinS.YueZ.YangX.XuJ.et al. (2021). Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe. Proc. Natl. Acad. Sci. U.S.A.118, e2101767118. doi: 10.1073/pnas.2101767118
59
TianS.GeJ.AiG.JiangJ.LiuQ.ChenX.et al. (2021). A 2.09 Mb fragment translocation on chromosome 6 causes abnormalities during meiosis and leads to less seed watermelon. Hortic. Res.8, 256. doi: 10.1038/s41438-021-00687-9
60
VilanovaS.BadenesM. L.Martínez CalvoJ.LlácerG. (2001). Analysis of loquat germplasm (Eriobotrya japonica Lindl) by RAPD molecular markers. Euphytica121, 25–29. doi: 10.1023/A:1012051207948
61
Von MalekB.van der GraaffE.SchneitzK.KellerB. (2002). The Arabidopsis male-sterile mutant dde2-2 is defective in the ALLENE OXIDE SYNTHASE gene encoding one of the key enzymes of the jasmonic acid biosynthesis pathway. Planta216, 187–192. doi: 10.1007/s00425-002-0906-2
62
WangJ.LiY.WangX.CaoK.ZhuG.FangW.et al. (2022). Betulin, synthesized by PpCYP716A1, is a key endogenous defensive metabolite of peach against aphids. J. Agric. Food Chem.70, 12865–12877. doi: 10.1021/acs.jafc.2c04422
63
WangY.ShahidM. Q.LinS.ChenC.HuC. (2017). Footprints of domestication revealed by RAD-tag resequencing in loquat: SNP data reveals a non-significant domestication bottleneck and a single domestication event. BMC Genomics18, 354. doi: 10.1186/s12864-017-3738-y
64
WangY.Van RengsW. M. J.ZaidanM. W.UnderwoodC. J. (2021). Meiosis in crops: from genes to genomes. J. Exp. Bot.72, 6091–6109. doi: 10.1093/jxb/erab217
65
WangY.ZhouH.HeY.ShenX.LinS.HuangL. (2023). MYB transcription factors and their roles in the male reproductive development of flowering plants. Plant Sci.335, 111811. doi: 10.1016/j.plantsci.2023.111811
66
WuX.ZhangY.LiW.TianC.FengQ.WangJ.et al. (2025). Unlocking the role of ClDYAD in initiating meiosis: A functional analysis in watermelon. Hortic. Plant J.11, 1260–1273. doi: 10.1016/j.hpj.2024.09.005
67
XiaQ.DangJ.WangP.LiangS.WeiX.LiX.et al. (2022). Low female gametophyte fertility contributes to the low seed formation of the diploid loquat [Eriobotrya Japonica (Thunb.) Lindl.] line H30-6. Front. Plant Sci.13. doi: 10.3389/fpls.2022.882965
68
XiaoY.ChenY.CharnikhovaT.MulderP. P. J.HeijmansJ.HoogenboomA.et al. (2014). OsJAR1 is required for JA-regulated floret opening and anther dehiscence in rice. Plant Mol. Biol.86, 19–33. doi: 10.1007/s11103-014-0212-y
69
YanJ.LiS.GuM.YaoR.LiY.ChenJ.et al. (2016). Endogenous bioactive jasmonate is composed of a set of (+)-7-iso-JA-amino acid conjugates. Plant Physiol.172, 2154–2164. doi: 10.1104/pp.16.00906
70
YuJ.JiangM.GuoC. (2019). Crop pollen development under drought: From the phenotype to the mechanism. Int. J. Mol. Sci.20, 1550. doi: 10.3390/ijms20071550
71
YuS. X.JiangY. T.LinW. H. (2022). Ovule initiation: the essential step controlling offspring number in Arabidopsis. J. Integr. Plant Biol.64, 1469–1486. doi: 10.1111/jipb.13314
72
YuanZ.ZhangD. (2015). Roles of jasmonate signalling in plant inflorescence and flower development. Curr. Opin. Plant Biol.27, 44–51. doi: 10.1016/j.pbi.2015.05.024
73
ZhangQ.XieZ.ZhangR.XuP.LiuH.YangH.et al. (2018). Blue light regulates secondary cell wall thickening via MYC2/MYC4 activation of the NST1-directed transcriptional network in Arabidopsis. Plant Cell30, 2512–2528. doi: 10.1105/tpc.18.00315
74
ZhouM.DengL.YuanG.ZhaoW.MaM.SunC.et al. (2023). Rapid generation of a tomato male sterility system and its feasible application in hybrid seed production. Theor. Appl. Genet.136, 197. doi: 10.1007/s00122-023-04428-5
75
ZhouL.WangY.XuX.YanD.YuW.MiaoY.et al. (2022). Conjunctive analyses of BSA-Seq and BSR-Seq unveil the Msβ-GAL and MsJMT as key candidate genes for cytoplasmic male sterility in alfalfa (Medicago sativa L.). Int. J. Mol. Sci.23, 7172. doi: 10.3390/ijms23137172
Summary
Keywords
loquat, low male fertility, seedless fruit, transcriptome, WGCNA, jasmonic acid
Citation
Wang J, Pan J, Xia Q, Duan Y, Jing D, Dang J and Guo Q (2025) Combined analysis of transcriptome and WGCNA reveals the mediating role of JA in the low male fertility of loquat H30-6. Front. Plant Sci. 16:1648313. doi: 10.3389/fpls.2025.1648313
Received
17 June 2025
Accepted
27 August 2025
Published
12 September 2025
Volume
16 - 2025
Edited by
Sonia Osorio, University of Malaga, Spain
Reviewed by
Ma Shiwei, Fujian Agriculture and Forestry University, China
Roohollah Shamloo-Dashtpagerdi, Higher Education Center of Eghlid, Iran
Updates
Copyright
© 2025 Wang, Pan, Xia, Duan, Jing, Dang and Guo.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
* Correspondence: Qigao Guo, mailto:qgguo@126.com; Jiangbo Dang, mailto:dangjiangbo@126.com
†These authors have contributed equally to this work
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.