Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 22 September 2025

Sec. Plant Development and EvoDevo

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1605312

This article is part of the Research TopicPhytohormones as Crucial Players in Organ AbscissionView all 7 articles

Unraveling the mechanisms of fruit abscission in Morus laevigata through multi-omics approaches

Jiahu YangJiahu Yang1Sha LiSha Li1Zhennan LiZhennan Li1Yahui XuanYahui Xuan2Jiamei HeJiamei He1Mingju RuanMingju Ruan1Yuming FengYuming Feng1Zhenyang TaoZhenyang Tao1Xiaoru KangXiaoru Kang1Zhengang Li*Zhengang Li1*
  • 1Sericulture and Apiculture Research Institute, Yunnan Academy of Agricultural Sciences, Mengzi, Yunnan, China
  • 2State Key Laboratory of Resource Insects, Southwest University, Beibei, Chongqing,, China

Introduction: Morus laevigata (long-fruited mulberry) is rich in active components and possesses significant nutritive value. The fruitlet stage represents a critical period for fruit abscission, and elucidating the underlying biological mechanisms of this process can provide a theoretical foundation for breeding more stable and abscission resistant cultivars.

Methods: Fruit peduncles at the fruit set stage (April to May) were selected as experimental materials, including both abscising and non-abscising fruits. Morphological analysis of the peduncle abscission zone was conducted to examine structural differences. Additionally, transcriptomic and metabolomic analyses were performed to investigate gene expression and metabolite changes associated with fruit abscission.

Results: Morphological analysis of the peduncle abscission zone in abscising fruits revealed enlarged intercellular spaces and disorganized cell arrangements. Transcriptomic and metabolomic analyses indicated that genes and metabolites related to fruit abscission were primarily involved in plant hormone signal transduction, and starch and sucrose metabolism pathways. Auxin and abscisic acid were identified as key regulators, modulating the expression of cell wall-degrading enzymes, which facilitated cell wall loosening and degradation, ultimately leading to fruit abscission. Furthermore, alterations in energy metabolism were found to play a pivotal role in this process.

Conclusion: These findings contribute to a deeper understanding of the physiological mechanisms underlying fruit development and abscission, offering valuable insights into mulberry breeding and the sustainable advancement of modern agriculture.

1 Introduction

Morus laevigata, commonly known as the long-fruited mulberry, belongs to the genus Morus L. within the family Moraceae. It is widely cultivated in various regions worldwide, particularly in the subtropical areas of Asia (Lee and Hwang, 2017). In addition to its high nutritional value, mulberry plays a significant role in traditional medicine (Eisenbrand and Tang, 2010), with numerous studies reporting that its active components confer various health benefits (Sarkhel et al., 2020; Yuan et al., 2017). Mulberry produces a multiple fruit, which is formed from the ovaries of many flowers borne on a single inflorescence. Each small drupelet within the fruit originates from an individual flower, and the entire structure matures into a fleshy, compound fruit. The peduncle supports the inflorescence and later the developing multiple fruit, playing a vital role in nutrient transport and fruit retention (Ozturk et al., 2023). During mulberry cultivation, fruit abscission is common and severe that significantly affects fruit yield and quality, thereby reducing economic benefits. Therefore, effective regulation of fruit abscission is crucial for improving productivity and profitability in mulberry production.

Fruit abscission is a vital physiological process that enables plants to shed mature, damaged, or unnecessary organs, thereby supporting reproductive success and seed dispersal. This process is tightly regulated and involves the differentiation and activation of a specialized cell layer known as the abscission zone (AZ) (Shi et al., 2023; van Nocker, 2009). Among the key regulators, plant hormones play a central role in modulating AZ development and activity. For instance, ethylene has been reported to induce fruit abscission (Kenea, 2021), while auxin can inhibit this process in various plant species (Li and Su, 2024). In sweet cherry, abscisic acid (ABA) has been shown to accumulate significantly during fruit abscission (Qiu et al., 2021). Additionally, in mulberry, ABA-related transcription factors MaABF1 and MaABI5 have been shown to antagonistically regulate the expression of MaJOINTLESS, a key gene involved in AZ development, further highlighting the pivotal role of ABA signaling in fruit abscission (Deng et al., 2023). Overall, fruit abscission is governed by a complex hormonal network, in which the crosstalk among ethylene, auxin, ABA, and other phytohormones plays a pivotal role in coordinating the timing and progression of cell separation processes in diverse crop species, such as apple (Sriskantharajah et al., 2021), tomato (Dong et al., 2024), and cotton (Zhao et al., 2024b). Conversely, hormonal imbalances can disrupt these regulatory signals, leading to fruit drop at various developmental stages, including the young fruitlet stage (Kumar and Kumar, 2016). Activation of the AZ is also closely associated with cell wall–modifying enzymes, including cellulases, polygalacturonases (PGs), and pectin methylesterases (PMEs). These enzymes degrade components key structural components of the primary cell wall and the middle lamella, thereby weakening cell-to-cell adhesion and facilitating organ detachment (Li and Su, 2024). The essential roles of cellulases in abscission are well established, with cellulase-encoding genes identified in various plant species, such as GH9 family genes in soybean (Zhan et al., 2025), citrus (Merelo et al., 2017), and mulberry (Deng et al., 2024). Enzymes such as pectate lyase have also been implicated in the abscission processes across multiple species including Arabidopsis thaliana (Niederhuth et al., 2013), soybean (Tucker et al., 2007), banana (Mbéguié et al., 2009), and rose (Singh et al., 2011). The coordinated activity of these enzymes, particularly cellulases and pectinases, leads to the breakdown of cellulose and pectin within the AZ and is considered one of the primary drivers of organ detachment (Bonghi et al., 1992; Roongsattham et al., 2016; Taylor and Whitelaw, 2001). In addition to hormonal and enzymatic regulation, carbohydrate metabolism also contributes to the regulation of fruit abscission (Zhu et al., 2011b). Carbohydrate depletion caused by impaired photosynthesis, resulting from either defoliation or shading treatment, has been shown to induce fruit abscission in various species (Mehouachi et al., 1995; Fei et al., 2025).

Recent advances in metabolomics and transcriptomics have enabled the systematic analysis of gene expression patterns and metabolic profiles associated with fruit abscission. The AZ represents a metabolically dynamic region involving complex networks of hormones, carbohydrates, and secondary metabolites (Tranbarger and Tadeo, 2025). In blue honeysuckle (Lonicera caerulea), non-targeted metabolomics revealed that differentially accumulated metabolites (DAMs) between easy- and hard-abscission cultivars were significantly enriched in pathways such as phenylpropanoid biosynthesis, starch and sucrose metabolism, and amino sugar metabolism, along with hormone biosynthesis routes like tryptophan and carotenoid metabolism in the fruit AZ (Chen et al., 2023a). Transcriptomic studies have identified multiple genes associated with AZ development and formation (van Nocker, 2009; Zhao et al., 2024a; Pautot et al., 2025). MACROCALYX and JOINTLESS, both MADS-box transcription factors, have been shown to regulate the expression of genes related to phytohormone signaling, cell wall modification, and key developmental transcription factors specifically in the tomato pedicel AZ (Nakano et al., 2012). In contrast, mutations in JOINTLESS impairs AZ development, resulting in the failure of organ separation at the normal abscission site (Huerga-Fernández et al., 2024). In mulberry, transcriptomic analysis showed that enrichment of the DEGs between easy-to-drop and less abscission-prone fruit stalks was significantly associated with pathways including flavonoid biosynthesis, the citric acid cycle, phytohormone signaling, and amino acid biosynthesis (Deng et al., 2023). Similarly, DEGs identified between abscised and retained young fruits were enriched in the same or related metabolic and signaling pathways (Xuan et al., 2023).

Despite the progress in current research, the molecular mechanisms underlying fruit abscission in mulberry remain largely unexplored. In this study, we integrated metabolomic and transcriptomic analyses to compare peduncles from abscising and non-abscising fruits of two mulberry varieties with contrasting abscission rates. The integration of these datasets offers novel insights into the coordinated regulation of gene expression and metabolite accumulation in the mulberry AZ, providing a valuable foundation for understanding the genetic and biochemical basis of fruit shedding in this species.

2 Materials and methods

2.1 Plant materials

This study utilized fruit peduncles as experimental materials, which were collected on March 15, 2023, from the Mulberry Germplasm Resource Garden of Yunnan Province. Two M. laevigata varieties with contrasting fruit abscission rates were selected: ‘Peach Fairy,’ characterized by a low fruit abscission rate, and ‘Pink Lady,’ exhibiting a high fruit abscission rate. Three healthy, disease-free three-year-old trees of each variety were randomly selected, with each tree considered a biological replicate. Fifteen days after flowering, three types of samples were collected: non-abscising fruits from ‘Peach Fairy’ (PF), as well as non-abscising (PL-N) and abscising fruits from ‘Pink Lady’ (PL-A). Non-abscising fruits were carefully selected based on their firmness and strong attachment to the branch, showing no signs of looseness or shriveling. Fruits considered to be abscising were identified by their ease of detachment, as they fell off with slight touch and had a shriveled appearance, indicating imminent abscission. For each sample, fruits were harvested from different orientations of the same tree, and their peduncles were excised and pooled. Immediately after collection, the peduncles were flash-frozen in liquid nitrogen, transferred into test tubes, and transported to the laboratory, where they were stored at -80 °C for further analysis.

The sampling time corresponds to the peak period of fruit abscission, as determined by monitoring abscission rates. Fruit drop began approximately 10 days after flowering and reached its highest level around the sampling time. To assess abscission rates, three branches were selected on each tree, and the initial fruit number was recorded. Subsequently, the number of fruits remaining on these branches was counted at different time points to calculate the rate of fruit drop (Supplementary Table 1).

2.2 Anatomical examination of fruit peduncles

Fresh fruit peduncles were collected and fixed in a standard FAA (formalin, acetic acid, and 70% ethanol) solution for 24 h. Following fixation, the tissues were processed for paraffin embedding, with longitudinal and transverse sections cut along the longitudinal axis of the peduncles. The sections were subsequently stained using safranin O and fast Green (BaiQianDu Biotechnology, Wuhan, China). The anatomical structure of the fruit peduncles was examined under a microscope (Nikon Eclipse E100, Japan) coupled with a Nikon DS-U3 imaging system. Micrographs were also taken using a scanning electron microscope (SEM) (Regulus; Hitachi, Japan).

2.3 Broad-targeted metabolomics via ultra-performance liquid chromatography–tandem mass spectrometry analysis

Approximately 100 mg of ground sample powder was extracted with 1 mL of 70% methanol containing 0.1 mg/L lidocaine as an internal standard. The extraction was performed overnight at 4 °C, with intermittent vortexing three times to ensure thorough extraction. Following extraction, the samples were centrifuged at 10,000 g for 10 min, and the supernatant was collected. The supernatant was filtered through a 0.22 μm microporous membrane and stored in vials for subsequent analysis.

Metabolic analysis was performed on a SCIEX ExionLC™ AD system coupled with SCIEX 6500+ QTRAP mass spectrometer (Framingham, MA, USA). UPLC separation was performed using a Waters ACQUITY UPLC HSS T3 C18 column (1.8 µm, 2.1 mm × 100 mm) (Milford, MA, USA). The mobile phase consisted of two components: the aqueous phase, which was ultra-pure water with 0.04% acetic acid, and the organic phase, which was acetonitrile with 0.04% acetic acid. A gradient elution was employed, starting with a mixture of water:acetonitrile at 95:5 (v/v) for 0 min, transitioning to 5:95 (v/v) at 11.0 min, maintaining 5:95 (v/v) from 11.0 to 12.0 min, returning to 95:5 (v/v) at 12.1 min, and holding at 95:5 (v/v) until 15.0 min. The flow rate was set to 0.4 mL/min, and the column temperature was maintained at 40 °C. A 5 μL sample was injected for each analysis. The electrospray ionization (ESI) source temperature was set to 550 °C, with a mass spectrometer voltage of 5500 V. The curtain gas (CUR) pressure was 25 psi, and the collision-activated dissociation (CAD) was high. Authentic standards were purchased from BioBioPha Co., Ltd. (Kunming, China) and Sigma-Aldrich (St. Louis, MO, USA).

2.4 Plant hormone analysis via ultra-high performance liquid chromatography–multiple reaction monitoring–tandem mass spectrometry

Approximately 100 mg of ground sample powder was transferred into an Eppendorf tube. A 1000 μL aliquot of extract solvent (acetonitrile-methanol, 1:1, containing an internal standard), pre-cooled to -40 °C, was added. The samples were then homogenized, vortexed for 60 s, and sonicated in an ice-water bath for 10 min. Following this, the samples were incubated at -40 °C for 2 h. After incubation, the samples were centrifuged at 12,000 rpm for 15 min at 4 °C. A total of 900 μL of the supernatant was collected and concentrated using a rotary evaporator at 4 °C. The concentrated samples were reconstituted with 90 μL of a 50% methanol aqueous solution, vortexed for 1 min, sonicated for 120 s, and vortexed again for 1 min. The samples were then centrifuged at 12,000 rpm for 10 min at 4 °C. A 70 μL aliquot of the clear supernatant was transferred to an auto-sampler vial for analysis.

UHPLC separation was performed on an ExionLC™ AD UHPLC System (SCIEX, Framingham, MA, USA), equipped with a UPLC Kinetex C18 column (2.1 mm × 100 mm, 2.6 μm; Phenomenex, Torrance, CA, USA). Mobile phase A consisted of 0.1% formic acid in water, while mobile phase B was 0.1% formic acid in methanol. The column temperature was maintained at 25 °C, while the auto-sampler temperature was set to 4 °C. A 2 μL injection volume was used for each analysis. For mass spectrometry analysis, a 6500 QTRAP+ Mass Spectrometer (SCIEX, Framingham, MA, USA) was employed, equipped with an electrospray ionization (ESI) interface. Ion source parameters were as set as follows: ion spray voltage = ± 4500 V, ion source gas 1 = 50 psi, ion source gas 2 = 50 psi, temperature (+/-) = 450 °C/400 °C, and curtain gas = 40/35 psi.

2.5 Transcriptomic profiling using Illumina RNA-sequencing technology

Total RNA was extracted using the Trizol reagent kit (Invitrogen, Carlsbad, CA, USA) and its quality was assessed through the Agilent 2100 Bioanalyzer (Santa Clara, CA, USA) and RNase-free agarose gel electrophoresis. Eukaryotic mRNA was enriched using Oligo(dT) beads, followed by fragmentation of the mRNA into short fragments using the NEBNext Ultra RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA), and reverse transcription to synthesize cDNA. The purified double-stranded cDNA underwent end repair, A-tail addition, and Illumina adapter ligation. The adapter-ligation reaction was purified using AMPure XP Beads and size-selected through agarose gel electrophoresis. After PCR amplification, the cDNA library was constructed. Finally, the library was sequenced using the Illumina Novaseq6000 system, with sequencing performed by Gene Denovo Biotechnology Co., Ltd. (Guangzhou, China).

2.6 Bioinformatics analysis

Differential expression analysis was performed using the DESeq2 package (Love et al., 2014) in R (v4.4.1), with genes identified as differentially expressed if p-value < 0.05 and variable importance in projection (VIP) > 1. Functional enrichment analysis was conducted using the clusterProfiler package (Yu et al., 2012) in R, in which Gene Ontology (GO) enrichment was performed with the enrichGO function, while KEGG pathway enrichment was conducted using the enrichKEGG function. The significance threshold for all enrichment results was set at p-value < 0.05. Visualization was carried out using ggplot2 (Wickham, 2011) to generate bubble plots. Correlation analysis was conducted using the Gephi software (v0.10.0) (Bastian et al., 2009). Only significantly altered hormone metabolites (VIP > 1.0, p-value < 0.05) and key differentially expressed genes (DEGs) (|log2FC| > 1, p-value < 0.05) were included as network nodes. Edges represented correlations between metabolite-hormone, and metabolite-hormone interactions, with Pearson correlation |r| > 0.6 and p-value < 0.05. Partial least squares discriminant analysis (PLS-DA) was conducted using the ropls package (Thevenot, 2016) in R. Heatmaps were generated with the pheatmap package (Kolde and Kolde, 2015). Venn diagrams were constructed using the VennDiagram package (Chen and Boutros, 2011). Principal component analysis (PCA) was performed and visualized using the ggplot2 package in R.

2.7 Quantitative polymerase chain reaction analysis

The extracted RNA was electrophoresed to assess RNA quality. cDNA was synthesized using the TransScript® One-Step gDNA Removal and cDNA Synthesis SuperMix (TransGen Biotech, Beijing, China), and the qPCR system was prepared using a qPCR kit (Polymerase, MF015). qPCR amplification was performed using a real-time fluorescence quantitative PCR instrument (Hehui Biotechnology, Suzhou, China). Relative mRNA expression levels were calculated using the 2-ΔΔCt method, with EF1α as the internal reference gene for normalization. Primer sequences are listed in Table 1.

Table 1
www.frontiersin.org

Table 1. Primer sequences.

2.8 Statistical analysis

Data were analyzed using RStudio 4.4.1 and Prism 8.0.2 (GraphPad). Comparisons between two groups were performed using an unpaired two-tailed t-test. Each experiment was repeated at least three times, and data are presented as the mean ± standard deviation. A p-value < 0.05 was considered statistically significant.

3 Results

3.1 Anatomical observations of the AZ in mulberry pedicels

To reveal the differences between non-abscising and abscising fruit peduncles in mulberry, both macroscopic appearance and the microstructure of the AZ were observed. Macroscopic observations showed that the inflorescences and peduncles of PL-A samples exhibited elongated morphology and slight yellowing, suggesting potential modified nutrient allocation and hormonal regulation. (Figures 1A, B). Comparative histological structure through longitudinal (Figure 1C) and transverse (Figure 1D) sections revealed disorganized cell arrangement in PL-A peduncles. In contrast, non-abscised peduncles maintained intact cellular organization with tightly packed cells. SEM analysis of transverse sections further highlighted structural features characteristic of the AZ in the PL-A group, including a loose cellular architecture and collapsed vascular bundles (Figure 1E). In comparison, the PF group maintained an intact and compact tissue structure (Figure 1F). Moreover, cells in the PL-A peduncle were irregularly arranged, with signs of cell wall rupture and partial detachment, in stark contrast to the orderly arrangement and structural integrity of cells observed in the PF peduncle (Figures 1E, F). These structural differences suggest that the tissue structure of abscising fruit peduncles was weaker, potentially making the fruit more susceptible to abscission.

Figure 1
A series of images showing plant structures.   A: Two sets of elongated structures labeled PL-A and PF, with arrows indicating the abscission zone (AZ).   B: Cross-sections of the structures with labels for PL-A and PF.   C: Microscopic view of vascular bundles in PL-A and PF, highlighting the abscission zone and non-abscission zone (NAZ).   D: Cross-section microscopy of PL-A and PF, focusing on cellular arrangements.   E: Scanning electron microscope images of PL-A and PF cross-sections showing detailed cellular structures.   F: Detailed SEM cross-sections of PL-A and PF, marked with asterisks.

Figure 1. Microscopic structure of longitudinal and transverse sections of peduncles of abscising fruits (PL-A) and non-abscising fruit (PF). (A) Representative images of fruit-bearing peduncles from PL-A and PF. (B) Representative images of peduncles without attached fruits from PL-A and PF. (C) Microscopic structure of the longitudinal section of peduncle from PL-A and PF. AZ: abscission zone; NAZ: non- abscission zone. (D) Microscopic structure of the cross-section of peduncles from PL-A and PF. (E, F) Representative scanning electron microscope (SEM) images of the abscission zones (AZs) from PL-A and PF. Red stars denote the vascular bundles, and arrows highlight the AZ and NAZ.

3.2 Metabolomic overview revealed distinct metabolic profiles among groups

PLS-DA revealed clear separation among the three sample groups, with strong consistency and correlation observed among replicates within each group (Supplementary Figure S1A). Heatmap and volcano plot analyses further highlighted substantial differences in metabolic profiles across groups, suggesting distinct metabolic responses (Supplementary Figure S1B–E). Differential analysis identified 186, 948, and 126 DAMs between PL-A and PL-N, PL-A and PF, and PL-N and PF, respectively. The greatest metabolic differences were observed between PL-A and PF, with 882 upregulated and 65 downregulated metabolites in PL-A. (Supplementary Figure S1F). Furthermore, this comparison group also showed the largest number of uniquely accumulated metabolites (Supplementary Figure S1G).

3.3 Metabolite enrichment analysis during fruitlet abscission in mulberry

To investigate the metabolic changes during mulberry fruit abscission, we compared the categories and pathway enrichments of DAMs among groups. The results showed that between PL-A and PF, DAMs were mainly enriched in amino acids and derivatives (16.7%), flavonoids (13.5%), lipids (9.3%), and organic acids and derivatives (8.6%) (Figure 2A). Similarly, in the comparison between PL-A and PL-N, amino acids and derivatives (17%) and flavonoids (13.3%) also accounted for a large proportion (Figure 2B). These may suggest that these two classes of metabolites may play key roles in the fruit abscission process. To further analyze the functions of these DMAs, KEGG pathway enrichment analysis was performed. We found that the DAMs between PL-A and PF, significantly enriched pathways included biosynthesis of plant secondary metabolites, biosynthesis of amino acids, and butanoate metabolism (Figure 2C). In the comparison between PL-A and PL-N, in addition to the enrichment of biosynthesis of plant secondary metabolites and biosynthesis of amino acids, pathways such as 2-oxocarboxylic acid metabolism, ABC transporters, and C5-branched dibasic acid metabolism were also significantly enriched (Figure 2D). These results further suggest that amino acid metabolism, organic acid metabolism, and substance transport may play key roles in the regulation of fruit abscission.

Figure 2
A four-part figure displays data in pie charts and bubble plots. Part A is a pie chart labeled “PL-A vs. PF,” showing various chemical categories with percentages, like amino acids (16.7%) and phenolic acids (9.3%). Part B is a pie chart labeled “PL-A vs. PL-N,” indicating similar categories with different proportions, such as amino acids (17%) and phenolic acids (9.5%). Part C is a bubble plot titled “KEGG (PL-A vs. PF),” depicting pathways like valine biosynthesis with varying bubble sizes for count and color for p-value. Part D is a similar bubble plot titled “KEGG (PL-A vs. PL-N),” showing pathways like leucine degradation.

Figure 2. Metabolite enrichment analysis during the fruitlet abscission process in mulberry. (A, B) Classification pie charts of DAMs between PL-A and PF (A), and PL-A and PL-N (B). (C, D) KEGG pathway enrichment analysis of DAMs between PL-A and PF (C), and PL-A and PL-N (D).

3.4 Hormonal alterations and network connectivity during mulberry fruitlet abscission

Plant hormones are crucial regulators in fruitlet abscission, where even subtle fluctuations in hormone levels can initiate complex physiological signaling cascades that determine fruit retention or detachment. To delve deeper into the hormonal basis of mulberry fruit abscission, we performed hormone metabolomics analysis on six fruit peduncle samples (three replicates each from PF and PL-A). In total, 11 hormone metabolites were significantly differentially accumulated during abscission, with 9 upregulated and 2 downregulated (Figures 3A, B). Notably, we observed that the levels of IAA and isopentenyladenosine (IPR) in the PL-A group was significantly lower than that in the PF group (P < 0.01), while ABA level was significantly elevated in the PL-A group (P < 0.001) (Figure 3C). These results may suggest that auxin depletion and ABA accumulation serve as molecular triggers for fruitlet abscission in mulberry. Further network analysis showed that auxin, cis-zeatin, and tryptophan occupied central positions in the network, indicating their strong connectivity with other metabolites. Specifically, multiple auxin conjugates (e.g., IAA-Ala, IAA-Leu) showed strong associations with other metabolites in the network, suggesting their potential role in fruit abscission (Figure 3D). These significant hormonal changes suggest that a hormonal imbalance may occur during the process of mulberry fruitlet abscission.

Figure 3
A multi-panel figure displaying data comparisons between plant hormone levels in groups PF and PL-A. Panel A is a heatmap clustering hormone levels, with a gradient ranging from blue (low) to red (high). Panel B is a volcano plot, contrasting log2 fold changes and p-values, indicating upregulated and downregulated factors. Panel C presents box plots for IAA, IPR, and Abscisic Acid levels, showing significant differences marked with asterisks. Panel D is a complex network diagram linking various metabolites and hormones, illustrating their interactions and relationships.

Figure 3. Hormone metabolite analysis in peduncles of PF and PL-A. (A) Heatmap showing the abundance of differential hormone metabolites between the PF and PL-A groups. (B) Volcano plots displaying differential hormone metabolites between PL-A and PF, p-value < 0.05, VIP > 1. (C) Comparison of indole-3-acetic acid (IAA), isopentenyladenosine (IPR), and abscisic acid (ABA) levels between PF and PL-A. (D) Correlation network between metabolites and hormones in PL-A and PF. Green nodes represent metabolites and blue nodes represent hormone metabolites. Node size reflects the number of connections (degree), and edge thickness indicates the strength of the correlation.

3.5 Transcriptomic overview revealed distinct gene expression patterns among groups

To investigate the molecular mechanisms underlying mulberry fruitlet abscission, nine samples were collected for RNA-seq analysis. PCA results showed clear separation between the three sample groups, with relatively small within-group variation (Supplementary Figure S2A). A total of 582 uniquely expressed genes were detected in the PL-A group (Supplementary Figure S2B). Minimal transcriptomic changes were observed between PL-N and PF, with 448 upregulated and 776 downregulated DEGs (Supplementary Figures S2C, F). In contrast, a total of 3,206 DEGs were identified between PL-A and PF, including 1,446 upregulated and 1,760 downregulated genes (Supplementary Figures S2D, F). Similarly, 2,807 DEGs were identified between PL-A and PL-N, comprising 1,130 upregulated and 1,677 downregulated genes (Supplementary Figures S2E, F).

3.6 KEGG enrichment analysis of DEGs

To investigate the roles of DEGs in mulberry fruitlet abscission, KEGG pathway enrichment analysis was performed. In the comparison between PL-N and PF, significantly enriched pathways included plant-pathogen interaction, MAPK signaling pathway-plant, and carotenoid biosynthesis, with the first two pathways showing the highest enrichment levels (Figure 4A). Between PL-A and PF, as well as PL-A and PL-N, DEGs were primarily enriched in starch and sucrose metabolism, and plant hormone signal transduction pathways (Figures 4B, C), indicating their potential involvement in fruit abscission. Heatmap clustering analysis of genes enriched in the plant hormone signal transduction pathway revealed significant differences between PL-A and PF (Figure 4D). Furthermore, gene-hormone network was established to explore regulatory relationships between genes and key plant hormones. ABA occupies a central position in the regulatory network, and its interactions with multiple genes highlight its critical role in regulating the fruitlet abscission process (Figure 4E). Additionally, hormones such as gibberellin A3 (GA3) and auxin were significantly associated with multiple genes in the network, further suggesting a synergistic regulatory effect of these hormones in the abscission process.

Figure 4
Five panels showing biological data visualizations. A, B, and C display KEGG pathways with rich factors and adjusted p-values for comparisons between different plant groups, represented with dot plots. D is a clustered heatmap showing gene expression profiles across samples and groups PL-N, PF, and PL-A, with color coding for expression levels. E is a network diagram illustrating interactions between various plant-related metabolites and genes, with labeled nodes and connecting lines.

Figure 4. KEGG enrichment analysis of differentially expressed genes (DEGs). (A–C) KEGG enrichment analysis of DEGs between PL-N and PF (A), PL-A and PF (B), as well as PL-A and PL-N (C). (D) Heatmap showing the expression profiles of genes involved in the plant hormone signal transduction pathway. (E) Correlation network between key plant hormones and genes in the PL-A and PF groups. Red lines indicate positive correlations, while blue lines indicate negative correlations. Red dots represent genes, and blue dots represent hormones. Node size corresponds to the number of connections (degree), and edge thickness reflects the strength of the correlation.

3.7 DEGs involved in plant hormones, cell wall biosynthesis, and signal transduction

Previous studies have demonstrated the importance of plant hormones in fruit abscission (Yuan et al., 2005). Therefore, we performed a clustering analysis of genes involved in the biosynthesis and signal transduction of ethylene, ABA, and auxin. In the PL-A group, most of the genes related to ethylene (Figure 5A) and ABA (Figure 5B) were upregulated, while genes related to auxin (Figure 5C) were downregulated, which was consistent with the results of the hormone metabolomics analysis. Moreover, in the PL-A group, genes related to cell wall biosynthesis, such as the TBL family and CESA family, were downregulated, leading to a reduced cell wall synthesis rate and possibly weakening the strength or stability of the cell wall (Figure 5D). Meanwhile, the expression of genes related to cell wall-degrading enzymes was upregulated, causing looser cell-to-cell connections, increasing the possibility of fruit abscission (Figure 5E).

Figure 5
Clustered heatmap depicting gene expression levels across different groups (PF, PL-N, PL-A) labeled A to E. Each panel contains rows of gene identifiers and columns representing samples. The color gradient from blue to red indicates expression levels from low to high. Dendrograms illustrate hierarchical clustering, showing relationships among genes and samples.

Figure 5. DEGs involved in plant hormones, cell wall biosynthesis, and signal transduction. (A) Heatmap of genes involved in ethylene biosynthesis and signal transduction. (B) Heatmap of genes involved in ABA biosynthesis and signal transduction. (C) Heatmap of genes involved in auxin biosynthesis and signal transduction. (D) Heatmap of genes involved in cell wall biosynthesis. (E) Heatmap of genes involved in cell wall degradation.

3.8 Starch and sucrose metabolism played essential roles in fruitlet abscission

Starch and sucrose are the major carbohydrates in plants, playing crucial roles in energy supply, signal transduction, and response to environmental stress (Chen et al., 2023b; Ruan, 2012). Previous studies have reported that starch and sucrose metabolic pathways are important in the fruit abscission process (Agustí et al., 2024). Consistently, our KEGG enrichment analysis revealed significant enrichment of the starch and sucrose metabolism pathway. Notably, genes involved in this pathway exhibited distinct expression patterns between PF and PL-A. In the starch metabolism pathway, the expression of starch-degrading enzymes (such as β-amylase, BAM) was upregulated in the PL-A group, leading to the rapid degradation of starch into maltose and glucose, thus providing energy for abscission process (Figure 6). In addition, genes related to cellulose degradation, such as β-glucosidase, were also upregulated in the PL-A group, indicating enhanced cell wall breakdown within the AZ. In the trehalose metabolism pathway, trehalose-6-phosphate phosphatase (TPP) expression was upregulated in the PL-A group, which may be an adaptive response of the plant to environmental stress and resource reallocation. Trehalose, as a protective sugar, helps plants maintain cellular stability under stress conditions. During the fruit abscission process, the plant faces reallocation of energy and carbon sources, and increased expression of TPP helps release trehalose, providing an additional energy source to support plant growth and repair. These changes in gene expression revealed the potential mechanisms by which sucrose and starch metabolism influence fruit abscission (Figure 6).

Figure 6
Metabolic pathway diagram showing the conversion of UDP-glucose to D-glucose through various intermediates such as trehalose, sucrose, maltodextrin, and cellulose. Enzymes like SuSy, SPS, ISA, and β-glucosidase are involved in different steps. Heat maps indicate gene expression levels and metabolite concentrations, with color scales for upregulation and downregulation. Pathway sections labeled A to E detail specific processes and interactions.

Figure 6. Pathway analysis of gene involved in starch and sucrose metabolism. (A) Biosynthesis pathway of trehalose. (B) Biosynthesis pathway of starch. (C) Biosynthesis and hydrolysis pathway of sucrose. (D) Hydrolysis pathway of starch. (E) Hydrolysis pathway of cellulose.

3.9 q-PCR validation of genes associated with mulberry fruitlet abscission

Based on transcriptomic and hormone metabolomic analyses, ABA, auxin, as well as the synthesis and degradation of the cell wall were closely associated with the abscission of mulberry fruitlets. Therefore, six genes associated with these processes were selected for further analysis: PP2C51, NCED3, IAA27, TIR1, CESA1, and XTH33 (Figure 7). The results showed that the expression of these genes was consistent with the transcriptomic data. During fruit abscission, the expression of PP2C51, NCED3, and XTH33 was significantly upregulated (P < 0.001). Among them, PP2C51 and NCED3 are closely related to the synthesis and signal transduction of abscisic acid, while XTH33 is a gene that promotes cell wall degradation. Their upregulation indicates that increased ABA levels and cell wall loosening are key factors leading to fruitlet abscission. Furthermore, the expression of auxin synthesis and signal transduction-related genes (IAA27, TIR1) and cell wall synthesis gene (CESA1) was significantly downregulated during fruit abscission, suggesting that changes in auxin levels are also an important factor in fruitlet abscission. The q-PCR results further validated that the changes in ABA and auxin content, as well as cell wall degradation, were crucial in mulberry fruitlet abscission.

Figure 7
Bar graphs labeled A to F display relative expression levels of genes: PP2C51, NCED3, IAA27, TIR1, CESA1, and XTH33 across conditions PL-N, PL-A, and PF. Significant differences are marked by asterisks. PL-A generally shows higher expression in graphs A, B, and F, while PL-N is higher in D and E.

Figure 7. Quantitative polymerase chain reaction (q-PCR) validation of gene expression related to mulberry fruitlet abscission. Relative expression of PP2C51 (A), NCED3 (B), IAA27 (C), TIR1 (D), CESA1 (E), and XTH33 (F) in peduncles of PL-A compared with PF. ***P < 0.001 (unpaired two-tailed t-test). Data are presented as mean ± standard deviation.

4 Discussion

Mulberry fruitlet abscission is a common physiological phenomenon during fruit development, which significantly affects the yield and quality of mulberries. In-depth research into the mechanisms of mulberry fruitlet abscission is crucial for improving mulberry breeding and cultivation practices. In this study, we found that significant changes in auxin and ABA levels, downregulated expression of genes related to cell wall synthesis, and the enrichment of sucrose and starch metabolic pathways all play crucial roles in regulating fruitlet abscission.

Plant hormones play crucial roles throughout the abscission process, as they mediate the response of plant organs to environmental and developmental stress signals (Estornell et al., 2013; Peleg and Blumwald, 2011; Smékalová et al., 2014). Hormones may either promote or inhibit abscission depending on their concentration, physiological balance, and transport dynamics within plant tissues, reflecting the intricacy of this regulatory network (Sawicki et al., 2015). For example, ethylene, ABA, and jasmonic acid are generally recognized as signals that promote abscission (Brown, 1997; Giulia et al., 2013; Liu et al., 2022), while auxin and GA are considered to inhibit this process (Ben-Cheikh et al., 1997; Yang et al., 2021).

This study showed that during the process of mulberry fruitlet abscission, the significant changes in auxin and ABA levels further highlighted the central role of hormones in regulating fruitlet abscission. The accumulation of ABA in the fruitlet abscission samples significantly increased, and was positively correlated with the high expression of genes related to ABA synthesis, indicating that ABA may promote fruitlet abscission by inducing the formation of the AZ. Meanwhile, auxin plays a dual role in the regulation of fruitlet abscission. It is generally accepted that high auxin levels maintain organ attachment, while a decrease in auxin concentration is necessary to permit abscission. In our study, auxin levels were significantly reduced in the abscissing peduncle, and the expression of auxin biosynthesis-related genes also showed a downward trend. This auxin depletion may enhance the capacity of the AZ to respond to separation signals, thereby facilitating the onset of fruitlet abscission in mulberry.

Cell wall degradation is a key physiological mechanism in fruit abscission, as the structural integrity of cell wall in the AZ is directly related to the strength of attachment between the fruit and the parent plant (Merelo et al., 2017). Transcriptomic analysis in this study revealed that the expression of genes encoding cell wall-modifying enzymes, such as cellulase (CEL5 and CEL1) and expansin (EXPA6, EXPA8, and EXPA11), was upregulated in the abscissing samples, indicating that these enzymes play a key role in fruitlet abscission. During fruitlet abscission, the upregulation of cellulase expression leads to the degradation of cellulose, gradually breaking down and softening the cell wall in the fruit, weakening the connection between the fruit and the plant, and ultimately resulting in abscission (Jiang et al., 2019). The elevated expression of cell wall-degrading enzymes in the AZ has also been widely reported in other fruit abscission studies. For example, in blueberries and lychees, the activation of these enzymes represents a crucial step in the abscission process (Chen et al., 2023a). Genes encoding cellulases, pectinases, PGs, expansin, xyloglucan endotransglucosylases/hydrolases, and peroxidases were found to be significantly upregulated during the abscission cycle in litchi (Li et al., 2015) and tomato (Nakano et al., 2012). Additionally, the degradation of the cell wall may be a result of hormonal regulation. Both ABA and ethylene can promote the expression of cell wall-degrading enzymes, thereby accelerating the formation of the AZ (Jiang et al., 2019; Wei et al., 2010). Therefore, the upregulation of cell wall-degrading enzymes during mulberry fruitlet abscission may be a combined effect of ABA accumulation and ethylene signaling activation. This phenomenon further supports the diversity of plant hormones in regulating fruitlet abscission.

Metabolomic and transcriptomic analyses in this study revealed that sucrose and starch metabolism pathways play a crucial role in the mulberry fruitlet abscission process. The research shows that the energy supply to the fruit is closely related to its growth, development, and abscission, with sucrose and starch metabolism providing the necessary carbon sources and energy for the fruit (Zhang et al., 2022). The significant enrichment of sucrose and starch metabolism pathways likely provided energy to the cells to meet the demands for cell wall degradation and other metabolic processes associated with the abscission process. In normally developing fruits, sucrose metabolism provides the energy required for growth, whereas during abscission, alterations in sucrose transport and degradation may affect carbon allocation and weaken fruit attachment. Therefore, sucrose level is considered an important factor in triggering the fruitlet abscission process, as demonstrated in Citrus unshiu (Gómez-Cadenas et al., 2000; Mehouachi et al., 1995). Additionally, higher levels of ABA can inhibit starch synthesis, thereby affecting fruit growth and development (Zhu et al., 2011a). In summary, this study suggests that changes in energy metabolism pathways may support carbohydrate redistribution to fuel cell wall degradation and fruit abscission.

5 Conclusion

This study explored the molecular mechanisms and pathways involved in mulberry fruitlet abscission. Through physiological experiments, metabolomics, and transcriptomics, we observed differences between the peduncles of non-abscising and abscising mulberry fruits. The abscising samples displayed typical AZ characteristics, including large intercellular spaces and loose cellular architecture. Multi-omics research revealed a complex molecular regulatory process involving numerous genes and metabolites. These findings provide valuable insights into the molecular mechanisms and pathways of mulberry fruitlet abscission and may contribute to the development of more effective breeding strategies. Further research is needed on the biosynthesis and signaling of auxin and ABA, as well as the regulatory genes controlling major components of the cell wall during the fruitlet abscission process.

Data availability statement

The data presented in this study have been deposited in the NCBI Sequence Read Archive (SRA) under BioProject accession number PRJNA1327559.

Author contributions

JY: Conceptualization, Writing – original draft. SL: Writing – original draft. ZL (3rd Author): Writing – original draft. YX: Data curation, Writing – original draft. JH: Writing – original draft, Formal Analysis. MR: Writing – review & editing, Resources. YF: Visualization, Writing – review & editing. ZT: Writing – original draft. XK: Writing – review & editing. ZL (10th Author): Funding acquisition, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This research was supported by Yunnan Provincial Joint Special Project for Basic Agricultural Research (202401BD070001-081); Academician Zhonghuai Xiang’s Research Station (No. 202305AF150166); and Special Fund for the Construction of Modern Agricultural Industry Technology System (CARS-18).

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.1605312/full#supplementary-material

References

Agustí, M., Martínez-Fuentes, A., Mesejo, C., Marzal, A., and Reig, C. (2024). Expression of carbohydrate-related genes underlying 3,5,6-TPA-induced fruitlet abscission in citrus. Sci. Rep. 14, 26482. doi: 10.1038/s41598-024-78310-9

PubMed Abstract | Crossref Full Text | Google Scholar

Bastian, M., Heymann, S., and Jacomy, M. (2009). Gephi: an open source software for exploring and manipulating networks. Proc. Int. AAAI Conf. Web Soc. media 3, 361–362. doi: 10.1609/icwsm.v3i1.13937

Crossref Full Text | Google Scholar

Ben-Cheikh, W., Perez-Botella, J., Tadeo, F. R., Talon, M., and Primo-Millo, E. (1997). Pollination increases gibberellin levels in developing ovaries of seeded varieties of citrus. Plant Physiol. 114, 557–564. doi: 10.1104/pp.114.2.557

PubMed Abstract | Crossref Full Text | Google Scholar

Bonghi, C., Rascio, N., Ramina, A., and Casadoro, G. J. P. M. B. (1992). Cellulase and polygalacturonase involvement in the abscission of leaf and fruit explants of peach. Plant Mol. Biol. 20, 839–848. doi: 10.1007/BF00027155

PubMed Abstract | Crossref Full Text | Google Scholar

Brown, K. M. (1997). Ethylene and abscission. Physiol. Plant 100, 567–576. doi: 10.1111/j.1399-3054.1997.tb03062.x

Crossref Full Text | Google Scholar

Chen, H. and Boutros, P. C. (2011). VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC. Bioinf. 12, 1–7. doi: 10.1186/1471-2105-12-35

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, L., Meng, Y., Bai, Y., Yu, H., Qian, Y., Zhang, D., et al. (2023b). Starch and sucrose metabolism and plant hormone signaling pathways play crucial roles in aquilegia salt stress adaption. Int. J. Mol. Sci. 24, 3948. doi: 10.3390/ijms24043948

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, J., Ren, B., Bian, C., Qin, D., Zhang, L., Li, J., et al. (2023a). Transcriptomic and metabolomic analyses reveal molecular mechanisms associated with the natural abscission of blue honeysuckle (Lonicera caerulea L.) ripe fruits. Plant Physiol. Biochem. 199, 107740. doi: 10.1016/j.plaphy.2023.107740

PubMed Abstract | Crossref Full Text | Google Scholar

Deng, J., Ahmad, B., Deng, X., Fan, Z., Liu, L., Lu, X., et al. (2024). Genome-wide analysis of the mulberry (Morus abla L.) GH9 gene family and the functional characterization of MaGH9B6 during the development of the abscission zone. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1352635

PubMed Abstract | Crossref Full Text | Google Scholar

Deng, X., Ahmad, B., Deng, J., Liu, L., Lu, X., Fan, Z., et al. (2023). MaABI5 and MaABF1 transcription factors regulate the expression of MaJOINTLESS during fruit abscission in mulberry (Morus alba L.). Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1229811

PubMed Abstract | Crossref Full Text | Google Scholar

Dong, X., Liu, X., Cheng, L., Li, R., Ge, S., Wang, S., et al. (2024). SlBEL11 regulates flavonoid biosynthesis, thus fine-tuning auxin efflux to prevent premature fruit drop in tomato. J. Integr. Plant Biol. 66, 749–770. doi: 10.1111/jipb.13627

PubMed Abstract | Crossref Full Text | Google Scholar

Eisenbrand, G. and Tang, W. (2010). Handbook of Chinese Medicinal Plants: Chemistry, Pharmacology, Toxicology (Weinheim, Germany: John Wiley & Sons).

Google Scholar

Estornell, L. H., Agustí, J., Merelo, P., Talón, M., and Tadeo, F. R. (2013). Elucidating mechanisms underlying organ abscission. Plant Sci. 199 - 200, 48–60. doi: 10.1016/j.plantsci.2012.10.008

PubMed Abstract | Crossref Full Text | Google Scholar

Fei, Y.-C., Mo, Y., Lin, K., Tao, L., He, X., Ye, F., et al. (2025). Comparative transcriptome analysis reveals a possible mechanism of starvation stress-induced fruit abscission in macadamia. BMC Plant Biol. 25, 739. doi: 10.1186/s12870-025-06761-x

PubMed Abstract | Crossref Full Text | Google Scholar

Giulia, E., Alessandro, B., Mariano, D., Andrea, B., Benedetto, R., and Angelo, R. (2013). Early induction of apple fruitlet abscission is characterized by an increase of both isoprene emission and abscisic acid content. Plant Physiol. 161, 1952–19 69. doi: 10.1104/pp.112.208470

PubMed Abstract | Crossref Full Text | Google Scholar

Gómez-Cadenas, A., Mehouachi, J., Tadeo, F. R., Primo-Millo, E., and Talon, M. (2000). Hormonal regulation of fruitlet abscission induced by carbohydrate shortage in citrus. Planta. 210, 636–643. doi: 10.1007/s004250050054

PubMed Abstract | Crossref Full Text | Google Scholar

Huerga-Fernández, S., Detry, N., Orman-Ligeza, B., Bouché, F., Hanikenne, M., Périlleux, C., et al. (2024). JOINTLESS maintains inflorescence meristem identity in tomato. Plant Cell Physiol. 65, 1197–1211. doi: 10.1093/pcp/pcae046

PubMed Abstract | Crossref Full Text | Google Scholar

Jiang, F., Lopez, A., Jeon, S., De Freitas, S. T., Yu, Q., Wu, Z., et al. (2019). Disassembly of the fruit cell wall by the ripening-associated polygalacturonase and expansin influences tomato cracking. Horticult. J. 6, 17. doi: 10.1038/s41438-018-0105-3

PubMed Abstract | Crossref Full Text | Google Scholar

Kenea, F. T. (2021). Role of ethylene in horticulture—a review. GP Globalize Res. J. Chem. 81, 81–92.

Google Scholar

Kolde, R. and Kolde, M. R. (2019). Package ‘pheatmap’ 1, 790. R package. Available online at: https://cran.r-project.org/web/packages/pheatmap/pheatmap.pdf.

Google Scholar

Kumar, R. and Kumar, V. (2016). Physiological disorders in perennial woody tropical and subtropical fruit crops: A review. Indian J. Agr. Sci. 86, 703–7 17. doi: 10.56093/ijas.v86i6.58831

Crossref Full Text | Google Scholar

Lee, Y. and Hwang, K. T. (2017). Changes in physicochemical properties of mulberry fruits (Morus alba L.) during ripening. Sci. Hortic. 217, 189–196. doi: 10.1016/j.scienta.2017.01.042

Crossref Full Text | Google Scholar

Li, J. and Su, S. (2024). Abscission in plants: from mechanism to applications. Biotechnol. Adv. 2, 27. doi: 10.1007/s44307-024-00033-9

PubMed Abstract | Crossref Full Text | Google Scholar

Li, C., Wang, Y., Ying, P., Ma, W., and Li, J. (2015). Genome-wide digital transcript analysis of putative fruitlet abscission related genes regulated by ethephon in litchi. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00502

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, X., Cheng, L., Li, R., Cai, Y., Wang, X., Fu, X., et al. (2022). The HD-Zip transcription factor SlHB15A regulates abscission by modulating jasmonoyl-isoleucine biosynthesis. Plant Physiol. 189, 2396–2412. doi: 10.1093/plphys/kiac212

PubMed Abstract | Crossref Full Text | Google Scholar

Love, M. I., Huber, W., and Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21. doi: 10.1186/s13059-014-0550-8

PubMed Abstract | Crossref Full Text | Google Scholar

Mbéguié, A. M. D., Hubert, O., Baurens, F. C., Matsumoto, T., Chillet, M., Fils-Lycaon, B., et al. (2009). Expression patterns of cell wall-modifying genes from banana during fruit ripening and in relationship with finger drop. J. Exp. Bot. 60, 2021–20 34. doi: 10.1093/jxb/erp079

PubMed Abstract | Crossref Full Text | Google Scholar

Mehouachi, J., Serna, D., Zaragoza, S., Agusti, M., Talon, M., and Primo-Millo, E. (1995). Defoliation increases fruit abscission and reduces carbohydrate levels in developing fruits and woody tissues of Citrus unshiu. Plant Sci. 107, 189–197. doi: 10.1016/0168-9452(95)04111-7

Crossref Full Text | Google Scholar

Merelo, P., Agustí, J., Arbona, V., Costa, M. L., Estornell, L. H., Gómez-Cadenas, A., et al. (2017). Cell wall remodeling in abscission zone cells during ethylene-promoted fruit abscission in citrus. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00126

PubMed Abstract | Crossref Full Text | Google Scholar

Nakano, T., Kimbara, J., Fujisawa, M., Kitagawa, M., Ihashi, N., Maeda, H., et al. (2012). MACROCALYX and JOINTLESS interact in the transcriptional regulation of tomato fruit abscission zone development. Plant Physiol. 158, 439–450. doi: 10.1104/pp.111.183731

PubMed Abstract | Crossref Full Text | Google Scholar

Niederhuth, C. E., Patharkar, O. R., and Walker, J. C. (2013). Transcriptional profiling of the Arabidopsis abscission mutant hae hsl2 by RNA-Seq. BMC Genomics 14, 37. doi: 10.1186/1471-2164-14-37

PubMed Abstract | Crossref Full Text | Google Scholar

Ozturk, M., Kamili, A. N., Altay, V., and Rohela, G. K. (2023). Mulberry: from botany to phytochemistry (Cham, Switzerland: Springer Nature).

Google Scholar

Pautot, V., Crick, J., and Hepworth, S. R. (2025). Abscission zones: cellular interfaces for the programmed separation of organs. Ann. Bot. 136, mcaf034. doi: 10.1093/aob/mcaf034

PubMed Abstract | Crossref Full Text | Google Scholar

Peleg, Z. and Blumwald, E. (2011). Hormone balance and abiotic stress tolerance in crop plants. Curr. Opin. Plant Biol. 14, 290–29 5. doi: 10.1016/j.pbi.2011.02.001

PubMed Abstract | Crossref Full Text | Google Scholar

Qiu, Z., Wen, Z., Hou, Q., Qiao, G., Yang, K., Hong, Y., et al. (2021). Cross-talk between transcriptome, phytohormone and HD-ZIP gene family analysis illuminates the molecular mechanism underlying fruitlet abscission in sweet cherry (Prunus avium L). BMC Plant Biol. 21, 173. doi: 10.1186/s12870-021-02940-8

PubMed Abstract | Crossref Full Text | Google Scholar

Roongsattham, P., Morcillo, F., Fooyontphanich, K., Jantasuriyarat, C., Tragoonrung, S., Amblard, P., et al. (2016). Cellular and pectin dynamics during abscission zone development and ripe fruit abscission of the monocot oil palm. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00540

PubMed Abstract | Crossref Full Text | Google Scholar

Ruan, Y.-L. (2012). Signaling role of sucrose metabolism in development. Mol. Plant 5, 763–765. doi: 10.1093/mp/sss046

PubMed Abstract | Crossref Full Text | Google Scholar

Sarkhel, S., Manvi, D., and Ramachandra, C. (2020). Nutrition importance and health benefits of mulberry leaf extract: A review. JPP. 9, 689–695. doi: 10.22271/phyto.2020.v9.i5j.12310

Crossref Full Text | Google Scholar

Sawicki, M., Aït Barka, E., Clément, C., Vaillant-Gaveau, N., and Jacquard, C. (2015). Cross-talk between environmental stresses and plant metabolism during reproductive organ abscission. J. Exp. Bot. 66, 1707–1719. doi: 10.1093/jxb/eru533

PubMed Abstract | Crossref Full Text | Google Scholar

Shi, Y., Song, B., Liang, Q., Su, D., Lu, W., Liu, Y., et al. (2023). Molecular regulatory events of flower and fruit abscission in horticultural plants. Hortic. Plant J. 9, 867–883. doi: 10.1016/j.hpj.2023.03.008

Crossref Full Text | Google Scholar

Singh, A. P., Pandey, S. P., Rajluxmi, Pandey, S., Nath, P., and Sane, A. P. (2011). Transcriptional activation of a pectate lyase gene, RbPel1, during petal abscission in rose. Postharvest Biol. Technol. 60, 143–148. doi: 10.1016/j.postharvbio.2010.12.014

Crossref Full Text | Google Scholar

Smékalová, V., Doskočilová, A., Komis, G., and Šamaj, J. (2014). Crosstalk between secondary messengers, hormones and MAPK modules during abiotic stress signalling in plants. Biotechnol. Adv. 32, 2–11. doi: 10.1016/j.bioteChadv.2013.07.009

PubMed Abstract | Crossref Full Text | Google Scholar

Sriskantharajah, K., El Kayal, W., Torkamaneh, D., Ayyanath, M. M., Saxena, P. K., Sullivan, A. J., et al. (2021). Transcriptomics of improved fruit retention by hexanal in ‘Honeycrisp’reveals hormonal crosstalk and reduced cell wall degradation in the fruit abscission zone. Int. J. Mol. Sci. 22, 8830. doi: 10.3390/ijms22168830

PubMed Abstract | Crossref Full Text | Google Scholar

Taylor, J. E. and Whitelaw, C. A. (2001). Signals in abscission. New Phytol. 151, 323–340. doi: 10.1046/j.0028-646x.2001.00194.x

Crossref Full Text | Google Scholar

Thevenot, E. A. (2016). ropls: PCA, PLS (-DA) and OPLS (-DA) for multivariate analysis and feature selection of omics data. R Package version. 1. Available online at: https://api.semanticscholar.org/CorpusID:15214843.

Google Scholar

Tranbarger, T. J. and Tadeo, F. R. (2025). Abscission zone metabolism impacts pre-and post-harvest fruit quality: a very attaching story. Front. Plant Sci. 15. doi: 10.3389/fpls.2024.1524893

PubMed Abstract | Crossref Full Text | Google Scholar

Tucker, M. L., Burke, A., Murphy, C. A., Thai, V. K., and Ehrenfried, M. L. (2007). Gene expression profiles for cell wall-modifying proteins associated with soybean cyst nematode infection, petiole abscission, root tips, flowers, apical buds, and leaves. J. Exp. Bot. 58, 3395–3406. doi: 10.1093/jxb/erm188

PubMed Abstract | Crossref Full Text | Google Scholar

van Nocker, S. (2009). Development of the abscission zone. Stewart Postharvest Rev. 5, 1–6. doi: 10.2212/spr.2009.1.5

Crossref Full Text | Google Scholar

Wei, J., Ma, F., Shi, S., Qi, X., Zhu, X., and Yuan, J. (2010). Changes and postharvest regulation of activity and gene expression of enzymes related to cell wall degradation in ripening apple fruit. Postharvest Biol. Technol. 56, 147–154. doi: 10.1016/j.postharvbio.2009.12.003

Crossref Full Text | Google Scholar

Wickham, H. (2011). ggplot2. Wiley Interdiscip. Rev. Comput. Stat. 3, 180–185. doi: 10.1002/wics.147

Crossref Full Text | Google Scholar

Xuan, D., Chunbing, C., Jing, D., Lianlian, L., Juan, L., and Xingfu, Z. (2023). Transcriptome analysis of abscised and normal surviving young fruit of Morus alba. J. Zhejiang A&F Univ. 40, 45–54. doi: 10.11833/j.issn.2095-0756.20220205

Crossref Full Text | Google Scholar

Yang, J., He, H., He, Y., Zheng, Q., Li, Q., Feng, X., et al. (2021). TMK1-based auxin signaling regulates abscisic acid responses via phosphorylating ABI1/2 in Arabidopsis. PNAS. 118, e2102544118. doi: 10.1073/pnas.2102544118

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 16, 284–287. doi: 10.1089/omi.2011.0118

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, R., Wu, Z., Kostenyuk, I. A., and Burns, J. K. (2005). G-protein-coupled α2A-adrenoreceptor agonists differentially alter citrus leaf and fruit abscission by affecting expression of ACC synthase and ACC oxidase. J. Exp. Bot. 56, 1867–1875. doi: 10.1093/jxb/eri176

PubMed Abstract | Crossref Full Text | Google Scholar

Yuan, Q., Zhao, L., and Chemistry, F. (2017). The mulberry (Morus alba L.) fruit-A review of characteristic components and health benefits. J. Agric. Food. Chem. 65, 10383–10394. doi: 10.1021/acs.jafc.7b03614

PubMed Abstract | Crossref Full Text | Google Scholar

Zhan, W., Huang, H., Cai, Z., Zhao, Z., Lin, X., Peng, Y., et al. (2025). Genome-wide identification and characterization of soybean GH9 endo-1, 4-β-glucanases. Front. Plant Sci. 16. doi: 10.3389/fpls.2025.1597668

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Wang, T., Zhang, F., Liu, Y., and Wang, G. (2022). Comparative Analysis of the Transcriptomes of Persisting and Abscised Fruitlets: Insights into Plant Hormone and Carbohydrate Metabolism Regulated Self-Thinning of Pecan Fruitlets during the Early Stage. Curr. Issues Mol. Biol. 44, 176–193. doi: 10.3390/cimb44010013

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, N., Geng, Z., Zhao, G., Liu, J., An, Z., Zhang, H., et al. (2024b). Integrated analysis of the transcriptome and metabolome reveals the molecular mechanism regulating cotton boll abscission under low light intensity. BMC Plant Biol. 24, 182. doi: 10.1186/s12870-024-04862-7

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, M., Shi, C.-L., and Li, J. (2024a). Abscission cues generated within the abscising organ and perceived by the abscission zone in woody fruit crops. Fruit Res. 4, e014. doi: 10.48130/frures-0024-0007

Crossref Full Text | Google Scholar

Zhu, H., Dardick, C. D., Beers, E. P., Callanhan, A. M., Xia, R., and Yuan, R. (2011b). Transcriptomics of shading-induced and NAA-induced abscission in apple (Malus domestica) reveals a shared pathway involving reduced photosynthesis, alterations in carbohydrate transport and signaling and hormone crosstalk. BMC Plant Biol. 11, 138. doi: 10.1186/1471-2229-11-138

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, G., Ye, N., Yang, J., Peng, X., and Zhang, J. (2011a). Regulation of expression of starch synthesis genes by ethylene and ABA in relation to the development of rice inferior and superior spikelets. J. Exp. Bot. 62, 3907–3916. doi: 10.1093/jxb/err088

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: Morus laevigata, fruit abscission, plant carbohydrates, plant hormone, auxin, abscisic acid

Citation: Yang J, Li S, Li Z, Xuan Y, He J, Ruan M, Feng Y, Tao Z, Kang X and Li Z (2025) Unraveling the mechanisms of fruit abscission in Morus laevigata through multi-omics approaches. Front. Plant Sci. 16:1605312. doi: 10.3389/fpls.2025.1605312

Received: 03 April 2025; Accepted: 19 August 2025;
Published: 22 September 2025.

Edited by:

Lara Reale, University of Perugia, Italy

Reviewed by:

Timothy John Tranbarger, Institut de Recherche Pour le Développement (IRD), France
Lei Ming Zhao, University of Chinese Academy of Sciences, China

Copyright © 2025 Yang, Li, Li, Xuan, He, Ruan, Feng, Tao, Kang and Li. 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: Zhengang Li, amlhbWluZzkyMkAxNjMuY29t

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.