Abstract
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 (). In addition to its high nutritional value, mulberry plays a significant role in traditional medicine (), with numerous studies reporting that its active components confer various health benefits (; ). 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 (). 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) (; ). 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 (), while auxin can inhibit this process in various plant species (). In sweet cherry, abscisic acid (ABA) has been shown to accumulate significantly during fruit abscission (). 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 (). 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 (), tomato (), 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 (). 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 (). 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 (), citrus (), and mulberry (). Enzymes such as pectate lyase have also been implicated in the abscission processes across multiple species including Arabidopsis thaliana (), soybean (), banana (), and rose (). 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 (; ; ). 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 (; ).
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 (). 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 (). Transcriptomic studies have identified multiple genes associated with AZ development and formation (; Zhao et al., 2024a; ). 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 (). In contrast, mutations in JOINTLESS impairs AZ development, resulting in the failure of organ separation at the normal abscission site (). 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 (). Similarly, DEGs identified between abscised and retained young fruits were enriched in the same or related metabolic and signaling pathways ().
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 () 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 () 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 () to generate bubble plots. Correlation analysis was conducted using the Gephi software (v0.10.0) (). 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 () in R. Heatmaps were generated with the pheatmap package (). Venn diagrams were constructed using the VennDiagram package (). 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
| Genes | Forward primers (5′→3′) | Reverse primers (5′→3′) |
|---|---|---|
| PP2C51 | GTCAAGGGCTGGTCTTCCTC | TTAATGACCCTTCCACCGGC |
| NCED3 | GTCACGGTTCGGGATTTTGC | GCTCCTCCCAAGCATTCCAT |
| IAA27 | GAGCTTCTGGGAACTGGGTT | TGGCACTCGTCGTCTGATTC |
| TIR1 | GTGTTGGAGAGAGTGGTGGG | GCAGTTTCCGATGAAGACGC |
| CESA1 | GGCCTTGACACTGATGGGAA | TCCGGCTTTCTTGTGATGCT |
| XTH33 | AGCACCACCAGTACAGCATC | GAAAACTGAAGGGTACGCGC |
| EF1α | GCTGCTGAGATGAACAAGCG | GGTCTCGAACTTCCACAGGG |
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
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
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
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
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 (). 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
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 (; ). Previous studies have reported that starch and sucrose metabolic pathways are important in the fruit abscission process (). 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
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
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 (; ; ). 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 (). For example, ethylene, ABA, and jasmonic acid are generally recognized as signals that promote abscission (; ; ), while auxin and GA are considered to inhibit this process (; ).
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 (). 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 (). 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 (). Genes encoding cellulases, pectinases, PGs, expansin, xyloglucan endotransglucosylases/hydrolases, and peroxidases were found to be significantly upregulated during the abscission cycle in litchi () and tomato (). 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 (; ). 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 (). 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 (; ). 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.
Statements
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
1
AgustíM.Martínez-FuentesA.MesejoC.MarzalA.ReigC. (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
2
BastianM.HeymannS.JacomyM. (2009). Gephi: an open source software for exploring and manipulating networks. Proc. Int. AAAI Conf. Web Soc. media3, 361–362. doi: 10.1609/icwsm.v3i1.13937
3
Ben-CheikhW.Perez-BotellaJ.TadeoF. R.TalonM.Primo-MilloE. (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
4
BonghiC.RascioN.RaminaA.CasadoroG. 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
5
BrownK. M. (1997). Ethylene and abscission. Physiol. Plant100, 567–576. doi: 10.1111/j.1399-3054.1997.tb03062.x
6
ChenH.BoutrosP. 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
7
ChenL.MengY.BaiY.YuH.QianY.ZhangD.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
8
ChenJ.RenB.BianC.QinD.ZhangL.LiJ.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
9
DengJ.AhmadB.DengX.FanZ.LiuL.LuX.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
10
DengX.AhmadB.DengJ.LiuL.LuX.FanZ.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
11
DongX.LiuX.ChengL.LiR.GeS.WangS.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
12
EisenbrandG.TangW. (2010). Handbook of Chinese Medicinal Plants: Chemistry, Pharmacology, Toxicology (Weinheim, Germany: John Wiley & Sons).
13
EstornellL. H.AgustíJ.MereloP.TalónM.TadeoF. R. (2013). Elucidating mechanisms underlying organ abscission. Plant Sci.199 - 200, 48–60. doi: 10.1016/j.plantsci.2012.10.008
14
FeiY.-C.MoY.LinK.TaoL.HeX.YeF.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
15
GiuliaE.AlessandroB.MarianoD.AndreaB.BenedettoR.AngeloR. (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
16
Gómez-CadenasA.MehouachiJ.TadeoF. R.Primo-MilloE.TalonM. (2000). Hormonal regulation of fruitlet abscission induced by carbohydrate shortage in citrus. Planta.210, 636–643. doi: 10.1007/s004250050054
17
Huerga-FernándezS.DetryN.Orman-LigezaB.BouchéF.HanikenneM.PérilleuxC.et al. (2024). JOINTLESS maintains inflorescence meristem identity in tomato. Plant Cell Physiol.65, 1197–1211. doi: 10.1093/pcp/pcae046
18
JiangF.LopezA.JeonS.De FreitasS. T.YuQ.WuZ.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
19
KeneaF. T. (2021). Role of ethylene in horticulture—a review. GP Globalize Res. J. Chem.81, 81–92.
20
KoldeR.KoldeM. R. (2019). Package ‘pheatmap’1, 790. R package. Available online at: https://cran.r-project.org/web/packages/pheatmap/pheatmap.pdf.
21
KumarR.KumarV. (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
22
LeeY.HwangK. 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
23
LiJ.SuS. (2024). Abscission in plants: from mechanism to applications. Biotechnol. Adv.2, 27. doi: 10.1007/s44307-024-00033-9
24
LiC.WangY.YingP.MaW.LiJ. (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
25
LiuX.ChengL.LiR.CaiY.WangX.FuX.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
26
LoveM. I.HuberW.AndersS. (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
27
MbéguiéA. M. D.HubertO.BaurensF. C.MatsumotoT.ChilletM.Fils-LycaonB.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
28
MehouachiJ.SernaD.ZaragozaS.AgustiM.TalonM.Primo-MilloE. (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
29
MereloP.AgustíJ.ArbonaV.CostaM. L.EstornellL. H.Gómez-CadenasA.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
30
NakanoT.KimbaraJ.FujisawaM.KitagawaM.IhashiN.MaedaH.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
31
NiederhuthC. E.PatharkarO. R.WalkerJ. C. (2013). Transcriptional profiling of the Arabidopsis abscission mutant hae hsl2 by RNA-Seq. BMC Genomics14, 37. doi: 10.1186/1471-2164-14-37
32
OzturkM.KamiliA. N.AltayV.RohelaG. K. (2023). Mulberry: from botany to phytochemistry (Cham, Switzerland: Springer Nature).
33
PautotV.CrickJ.HepworthS. R. (2025). Abscission zones: cellular interfaces for the programmed separation of organs. Ann. Bot.136, mcaf034. doi: 10.1093/aob/mcaf034
34
PelegZ.BlumwaldE. (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
35
QiuZ.WenZ.HouQ.QiaoG.YangK.HongY.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
36
RoongsatthamP.MorcilloF.FooyontphanichK.JantasuriyaratC.TragoonrungS.AmblardP.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
37
RuanY.-L. (2012). Signaling role of sucrose metabolism in development. Mol. Plant5, 763–765. doi: 10.1093/mp/sss046
38
SarkhelS.ManviD.RamachandraC. (2020). Nutrition importance and health benefits of mulberry leaf extract: A review. JPP.9, 689–695. doi: 10.22271/phyto.2020.v9.i5j.12310
39
SawickiM.Aït BarkaE.ClémentC.Vaillant-GaveauN.JacquardC. (2015). Cross-talk between environmental stresses and plant metabolism during reproductive organ abscission. J. Exp. Bot.66, 1707–1719. doi: 10.1093/jxb/eru533
40
ShiY.SongB.LiangQ.SuD.LuW.LiuY.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
41
SinghA. P.PandeyS. P.RajluxmiPandeyS.NathP.SaneA. 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
42
SmékalováV.DoskočilováA.KomisG.ŠamajJ. (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
43
SriskantharajahK.El KayalW.TorkamanehD.AyyanathM. M.SaxenaP. K.SullivanA. 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
44
TaylorJ. E.WhitelawC. A. (2001). Signals in abscission. New Phytol.151, 323–340. doi: 10.1046/j.0028-646x.2001.00194.x
45
ThevenotE. 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.
46
TranbargerT. J.TadeoF. 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
47
TuckerM. L.BurkeA.MurphyC. A.ThaiV. K.EhrenfriedM. 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
48
van NockerS. (2009). Development of the abscission zone. Stewart Postharvest Rev.5, 1–6. doi: 10.2212/spr.2009.1.5
49
WeiJ.MaF.ShiS.QiX.ZhuX.YuanJ. (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
50
WickhamH. (2011). ggplot2. Wiley Interdiscip. Rev. Comput. Stat.3, 180–185. doi: 10.1002/wics.147
51
XuanD.ChunbingC.JingD.LianlianL.JuanL.XingfuZ. (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
52
YangJ.HeH.HeY.ZhengQ.LiQ.FengX.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
53
YuG.WangL.-G.HanY.HeQ.-Y. (2012). clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS.16, 284–287. doi: 10.1089/omi.2011.0118
54
YuanR.WuZ.KostenyukI. A.BurnsJ. 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
55
YuanQ.ZhaoL.ChemistryF. (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
56
ZhanW.HuangH.CaiZ.ZhaoZ.LinX.PengY.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
57
ZhangJ.WangT.ZhangF.LiuY.WangG. (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
58
ZhaoN.GengZ.ZhaoG.LiuJ.AnZ.ZhangH.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
59
ZhaoM.ShiC.-L.LiJ. (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
60
ZhuH.DardickC. D.BeersE. P.CallanhanA. M.XiaR.YuanR. (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
61
ZhuG.YeN.YangJ.PengX.ZhangJ. (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
Summary
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
Volume
16 - 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
Updates
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, jiaming922@163.com
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.