Preprint
Article

Effects of Gibberellic Acid on Soluble Sugar Content, Organic Acid Composition, Endogenous Hormone Levels, and Carbon Sink Strength in Shine Muscat Grapes during Berry Development Stage

Altmetrics

Downloads

107

Views

30

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

20 February 2024

Posted:

20 February 2024

You are already at the latest version

Alerts
Abstract
This study investigated the effects of exogenous gibberellic acid (GA) treatments (elongating cluster + seedless + expanding, T1; seedless + expanding, T2; expanding, T3 and water, CK) on the content of sugars, organic acids, and endogenous hormones and sink strength.Results showed that T2 treatment displayed the highest fructose and glucose levels at 100 days after treatment (DAT), whereas its effect on tartaric acid, malic acid, and citric acid concentrations at 80 and 100 DAT was relatively weak. Under GA3 treatments, both GA3, IAA, and CTK contents increased, whereas ABA contents decreased at 1, 2, 4, 8, and 48 h. Analysis of sugar phloem unloading revealed that T2 treatment exhibited the highest values during softening and ripening stages. Our findings indicate that appropriate GA3 application can positively influence sink strength by regulating sink size and activity, including berry size enlargement, sugar phloem unloading and sugar accumulation in grape sink cells.
Keywords: 
Subject: Biology and Life Sciences  -   Horticulture

1. Introduction

Gibberellin, a pivotal phytohormone that interacts synergistically with other plant hormones, functions as a central regulator of diverse developmental processes within plants. Gibberellic acid (GA), a specific form of gibberellin, plays a critical role in promoting parthenocarpy and has garnered extensive attention for application in the cultivation of seedless fruits, augmentation of berry dimensions, and elongation of rachis across a wide spectrum of grape cultivars [1,2,3,4]. GA3, in conjunction with auxin, plays a pivotal role in stimulating cell division and expansion. This interaction regulates fruit development and subsequent enlargement following fertilization. Recently, the consumption of fruits has increased owing to their high internal and external appearance quality. This trend is associated with economic growth and changing consumer preferences. To meet market demands, producing grapes with standardized cluster length and uniform and large berry size has become crucial. Consequently, gibberellin finds widespread application in production.
Presently, a number of studies have demonstrated the efficacy of gibberellin as an effective agent in enhancing the elongation of grape inflorescence. Sun [5] reported that the application of gibberellin to Cabernet Franc grape led to a proportional increase in peduncle length with an increase in gibberellin concentration. Similarly, Wang et al.[6] observed that grape varieties such as Midknight Beauty, Zaoheibao, and Summer Black exhibited elongated fruit pedicels following gibberellin treatments at a concentration of 10 mg/kg. Yang et al.[7] revealed that the application of exogenous GA3 not only inhibited the synthesis of endogenous gibberellin but also orchestrated the modulation of gibberellin signal transduction, thereby promoting inflorescence elongation.
Prior research has highlighted the utilization of GA3 as a viable strategy in grape cultivation, particularly for seed management. Notably, the application of GA3 solution at a concentration of 100 mg L−1 prior to full bloom has demonstrated its potential to produce seedless cultivars and prompt seed abortion in seeded cultivars [1]. Han and Lee [8] observed that GA3 exhibited efficacy in not only fostering fruit enlargement but also augmenting cluster length, cluster weight, and berry weight. Korkutal et al.[4]further substantiated the multifaceted impacts of gibberellins, revealing their role in promoting stem elongation, triggering flowering induction, stimulating pollen tube growth, yielding seedless fruits, and increasing the size of seedless berries. Zhao[9] expounded upon the benefits of employing GA3 and streptomycin before the full-bloom stage, coupled with GA3 application alone post this stage. This approach was found to notably increase berry size of “Shenxiangwuhe” and “12-17,” concurrently facilitating seedlessness in the latter. For the grape variety “Zhuosexiang,” the combined administration of GA3 and forchlorfenuron (CPPU) before full bloom, followed by GA3 application post full bloom, exhibited comparable benefits, yielding a remarkable seedless rate of 71%. Notably, GA3 application emerged as the most effective approach in enhancing single berry weight and seedless rate, as underscored by Zhao et al.[9].
Sink strength, which has been elucidated as the competitive capacity of an organ to efficiently intake photoassimilates, can be considered as a product of two critical components: sink size and sink activity [10]. Sink size pertains to the physical limitation represented by the total biomass of the sink organ. Notably, the quantification of cells, both in terms of their number and size in the sink organ, can serve as a suitable metric for evaluating sink size. Sink activity is regarded as the physiological limitation that governs the import of assimilates into a sink organ. It is determined by three pivotal physiological processes. First, it involves the unloading of assimilates from the phloem, followed by the subsequent transport of sugars beyond the phloem, leading to their absorption by the sink organ. Second, the sink organ’s own respiratory consumption contributes to sink activity. Third, the accumulation of carbohydrates within sink organs also influences sink activity. Prior research has indicated that phytohormones can modulate sink strength. Notably, hormones such as GA3, cytokinin (CTK), abscisic acid (ABA), and auxin have been recognized as participants in this regulatory process [11,12]. However, the underlying mechanism through which GA specifically influences sink strength in perennial fruit crops, including grape, remains elusive.
Shine Muscat is a prominent table grape variety extensively cultivated in China. This variety is widely known for its excellent characteristics such as large size, sweet taste, good storage capacity, and aromatic qualities; however, certain issues, such as the occurrence of seeded fruit, bitter seeds and pericarp, and tightly packed bunches, remain persistent. GA is widely employed to cultivate larger, seedless berries and enhance overall fruit quality. Previous investigations have demonstrated that GA can induce inflorescence elongation, foster the development of seedless fruits, and augment the size of seedless berries in different grape varieties. Our previous research has also indicated that GA influences the aromatic components in Shine Muscat berries [13]. However, a research gap remains, particularly in understanding the interrelationship among the contents of soluble sugars, organic acids, and endogenous hormones and sink strength in Shine Muscat grapes subjected to GA3 treatments. Therefore, the present study aimed to investigate how GA3 influences the contents of soluble sugars and organic acids, preliminarily explore the interrelationship among these components, and clarify the mechanism underlying the fruit sink strength modulation by GA3. Furthermore, this research seeks to contribute substantial data to support the cultivation of table grapes and the utilization of GA in China.

2. Materials and Methods

2.1. Field conditions and materials

The experiment was conducted at the Jinniushan vineyard in Tai’an, Shandong, China (36.127°N, 117.004°E). For the study, Shine Muscat (Vitis labruscana) grapes, which were self-rooted and planted in 2012, were used. The planting arrangement involved a plant density of 3.0 per hectare (with a spacing of 1.0 m between plants and 3.0 m between rows) within a rain shelter cultivation system. For the experimentation, a total of 30 healthy and well-developed vines were meticulously selected, with each treatment group comprising three replicates. On each vine, four shoots were trained, and for each treatment, a single shoot was allocated. The selected shoots were pruned to retain one basal cluster with one inflorescence. Approximately one week before flowering, the inflorescence was pruned, leaving a uniform 5 cm tip length to ensure synchronized flowering across treatments. The vines were subjected to GA3 (ProGibb 40, Valent BioScience, Walnut Creek, CA, USA) treatment at varying concentrations and during distinct periods. The experimentation comprised four distinct treatment groups, including three GA treatment groups and one control (CK) group, as detailed in Table 1.
Note:
1. Rachis elongation treatment: Twenty days prior to anthesis, inflorescences were immersed in a solution containing 5 mg/L GA3.
2. Seedlessness treatment: A day after full bloom, inflorescences were immersed in a solution containing 25 mg/L GA3.
3. Expansion treatment: Within two weeks after bloom withering, inflorescences were immersed in a solution containing 30 mg/L GA3.
4. CK treatment: Inflorescences were treated with water.

2.2. Tissue collection

The cluster length of 30 berries for each treatment was assessed at specific time intervals: 0 h, 12 h, 24 h, 3 d, and 7 d after the rachis elongation treatment. Subsequently, for each treatment, 120 grape seeds were randomly chosen from 30 berries on a regular basis at 7:00 am every 20 days [20, 40, 60, 80, and 100 days after treatment (DAT)] post expansion treatment. The samples were segregated into two groups. The first group samples, comprising 60 seeds from each treatment, were collected and transported in an ice box to the laboratory. The seeds were then analyzed for berry weight, longitudinal diameter (LD), transverse diameter (TD), total soluble solids (TSS), and titratable acidity (TA). The second group samples, comprising 60 seeds from each treatment, were frozen in liquid nitrogen and immediately stored in a refrigerator at −80 ℃. These samples were intended for the assessment of soluble sugars (glucose, fructose, and sucrose) and organic acids (tartaric acid, citric acid, and malic acid) in the grape berries. Moreover, samples comprising 30 seeds with consistent fruit diameter, obtained from 30 berries of each treatment, were collected at specific intervals: 0, 1, 2, 4, 8, 24, and 48 h post expansion treatment. Following liquid nitrogen treatment, these samples were enveloped in tinfoil, refrigerated, and subsequently brought back to the laboratory. The samples were then stored in a low-temperature refrigerator (−80 ℃) for the determination of indole-3-acetic acid (IAA), CTK, GA, and ABA contents in grape berries within one month.

2.3. Measurements of berry growth

The weight of individual fruits was determined using an electronic balance, with an accuracy of 0.01 g. For each analysis, a total of 30 berries were randomly selected from 20 to 30 clusters and each of the three replications. LD and TD of the fruits were measured using a vernier caliper. TSS was assessed by utilizing an aliquot of grapevine juice with the aid of a digital refractometer (Automatic Refractometer SMART-1, Atago, Tokyo, Japan), which is denoted as Brix. Total acidity was quantified through an acid-base titration method, employing 0.1 mol L−1 NaOH to achieve a pH of 8.1. The values are expressed as grams of citric acid per kilogram of fresh weight.

2.4. High-performance liquid chromatography analysis of glucose, fructose, and sucrose

Sugars were extracted following the procedure: for each treatment, 1 g of homogenized grape berry material was accurately weighed and then diluted to 10.0 mL using ultrapure water (Millipore, Bedford, MA, USA). The solution was subsequently incubated for 20 min in a water bath set at 35 °C. Afterward, the supernatant was subjected to centrifugation at 21,000×g for 10 min at room temperature (BHG-Hermle Z 365, Wehingen, Germany). The extraction process was repeated three times, and the resulting supernatants were pooled. The liquid supernatant underwent filtration through a 0.22-μm, 13-mm diameter syringe filter (Shanghai Xingya Purification Material Factory, Shanghai, China). Subsequently, the filtered solution was employed for the analysis of glucose, sucrose, and fructose contents. High-performance liquid chromatography (HPLC; Waters, Milford, MA, USA) was employed for the analysis. The separation conditions used for the soluble sugar analysis were as follows: detector, differential refractive index detector (RID); column, YMC- Pack Polyamine II (4.6 mm × 250 mm); phase, acetonitrile/water = 75:25 (v/v); flow rate, 0.8 mL/min; injection amount, 10 μL; column temperature, 40 °C; analysis time, 20 min. The eluted peaks were detected utilizing an RID-10A differential refractive index detector (Shimadzu Co., Ltd., Kyoto, Japan). The quantification of glucose, fructose, and sucrose was conducted using standard curves of authentic compounds. Each treatment was replicated three times.

2.5. High-performance liquid chromatography analysis of tartaric acid, citric acid, and malic acid

Organic acid analyses were conducted utilizing a Waters series 515 chromatography unit, which was equipped with two 515 pumps and a 2487 dual UV detector (Waters Alliance 2695 HPLC) operating at a wavelength of 210 nm. The separation conditions used for organic acid analysis were as follows: column, Thermo Hypersil COLD aQ (4.6 mm × 250 mm, 5 μm); phase, 10 mmol/L NH4H2PO4 (pH = 2.3)/methanol = 98/2 (v/v); flow rate, 0.8 mL/min; injection amount, 10 μL; column temperature, 25 °C; analysis time, 20 min. Quantification of tartaric acid, citric acid, and malic acid was performed using standard curves of authentic compounds. The analysis encompassed extracts from three replicate tissue samples.

2.6. Extraction and determination of GA3, IAA, CTK, and ABA

HPLC (Nexera LC-30A, Shimadzu Co., Ltd.) was employed for analysis of GA3, IAA, CTK, and ABA. The standard hormone sample utilized was a product of Sigma Company. The chromatographic conditions were as follows: mobile phase, methanol: 0.8% glacial acetic acid solution = 55:45 (v/v); flow rate, 0.8 mL/min; injection amount, 10 μL; column temperature, 30 ℃; detection wavelength, 254 nm. The analysis comprised extracts from three replicate tissue samples.

2.7. Analysis of sugar phloem unloading

The sugar phloem unloading was assessed using the fruit cup method. During the softening stage (August 19) and ripening stage (September 28), one sunward-facing grain in grape berries was randomly selected for each treatment. A cross of about 1 cm was delicately marked at the grape’s umbilical region using a scalpel. This process ensured that the flesh remained unharmed and the grape flesh tissue was not cut. Subsequently, the peel was gently separated from the pedicel using tweezers and carefully cut with anatomical scissors. The prepared Mes buffer, comprising 5 mmol/L Mes (2-(n-morph) ethanesulfonic acid), 2 mol/L CaCl2, 100 mmol/L mannitol (D-mannitol), and 0.2% (w/v) polyvinylpyrrolidone (PVPP), was then used to rinse the peel. The peeled grape fruit was meticulously placed into a 20 mm infusion tube filled with Mes buffer. The opening was sealed with a film, and the fruit cup was secured on the cluster. The outer layer of the fruit cup was covered with tin foil to prevent temperature fluctuations from affecting the test results and extending the “berry cup’s” longevity. Sampling times were designated as 9:30–10:00, 11:30–12:00, 13:30–14:00, 15:30–16:00, and 17:30–18:00, during which the replaced buffer was discharged from the cup. Each time, the buffer liquid from the fruit cup was collected, and it was filled using a syringe. The collected samples were immediately frozen using liquid nitrogen and stored at −80 °C for subsequent analysis of glucose, fructose, and sucrose.

2.8. Statistical analysis

Data analysis was performed using Microsoft Excel 2010 and SPSS 26.0. Graphical representations were created using Prism 9. Pearson correlation analysis was conducted using R statistical software.

3. Results

3.1. Effects of GA on physicochemical characteristics of grape berries

To investigate the effects of GA on Shine Muscat grapes, various key attributes of grape berries during development were evaluated in both GA-treated and CK groups. These parameters included berry weight, berry size, LD, TD, vertical/horizontal ratio (a measure of fruit shape index), TA (g/L), and TSS (Brix) of the berries. As illustrated in Figure 1A, the trend of changes in berry weight was analogous between GA-treated groups and the CK group. Notably, at each developmental stage, significantly higher berry weights were observed under T1, T2, and T3 treatments than those under the CK treatment. Furthermore, berry weights were substantially higher under T1 and T2 treatments than those under the T3 treatment. At 20–100 DAT, the berry size exhibited a comparable trend to the berry weight for both GA-treated groups and the CK group. The berry weight and volume were more rapidly increased under the T2 treatment than those under any other treatment (as shown in Figure 1A, B). Specifically, the berry weights were 11.3 g and 12.05 g in 2022 and 2021, respectively, marking a remarkable 52.2% (2022) and 42.8% (2021) increase, respectively, over those in the CK group at 100 DAT (Figure 1A). Similarly, the berry volumes were 10.59 cm3 and 10.63 cm3 in 2022 and 2021, respectively, representing a substantial 21.8% (2021) and 59.2% (2022) increase, respectively, than those in the CK group at 100 DAT (Figure 1B). Moreover, the berry size was significantly higher in T1, T2, and T3 treatment groups than in the CK group (Figure 1C). Moreover, no significant increase in TD was observed between GA-treated groups and the CK group (Figure 1D). In contrast, LD was substantially higher in the three GA-treated groups than in the CK group at 20–100 DAT (Figure 1E). However, similar LD values were observed for T1 and T2 treatment groups. The fruit shape index was evidently higher in T1 and T2 treatment groups than in the other groups, with that in the CK group being the lowest. At 21–70 DAT, a slight increase in fruit shape index in the Tl group was evident compared to that in the T2 group (Figure 1F).
The TSS content in grape berries is illustrated in Figure 2A. It displayed an increasing trend as the fruits ripened. In 2021, at 20 DAT, remarkably higher TSS contents were observed in GA-treated groups than those in the CK group. At 40 DAT, the TSS content demonstrated the following trend: T2 > T3 > T1 > CK. At 60, 80, and 100 DAT, significantly higher TSS contents were observed under T1 and T2 treatments than those under T3 and CK treatments. Moreover, in 2022, TSS contents increased across all treatments as the fruits matured. At 60 DAT, a significant difference in TSS content was observed between the T3 and CK groups. However, at 20, 40, 80, and 100 DAT, no significant difference was observed between GA-treated groups and the CK group.
An evident declining trend in TA during fruit development was observed, particularly at 20–60 DAT (Figure 2B). In 2021, no significant differences were observed among the treatments at 20, 80, and 100 DAT. However, at 40 DAT, CK and T3 groups exhibited higher TA compared to T1 and T2 groups. At 60 DAT, the CK group displayed the highest TA among the treatments. In 2022, no significant differences in TA were recorded among the treatments at 20 DAT. At 40 DAT, the overall trend observed was as follows: CK > T3 > T1 > T2, with TA in the CK group being significantly higher than that in GA-treated groups. At 100 DAT, a significant difference in TA was observed under the T3 and T1 (or T2) treatments.
The TSS/TA ratio demonstrated a consistent increasing trend during the course of fruit development. In 2021, varying GA treatments yielded minimal influence on the TSS/TA ratio across different time periods. However, at 100 DAT, the TSS/TA ratio in CK, T1, and T2 groups was significantly higher than that in the T3 group; however, no significant differences were observed among CK, T1, and T2 groups (Figure 2C).

3.2. Effect of GA on soluble sugar content in grape berries

As illustrated in Figure 2D, the soluble sugar composition of grape berries primarily consists of glucose and fructose. Notably, the sucrose content was negligible. The glucose and fructose contents within all four groups exhibited a consistent increasing trend; however, throughout the shelf life, the fructose content consistently surpassed that of glucose. The increase in total sugar content was significantly higher at 40–60 DAT than during the maturation stage (80–100 DAT). Prior to reaching 60 DAT, the total sugar content in T1 and T2 treatment groups was significantly higher than that in the CK group. Moreover, the total sugar content in the T3 group was the lowest at 80 DAT among all four treatment groups. During the ripening stage, the T2 group exhibited the highest total sugar content among all four treatment groups.
The fructose and glucose contents within each treatment group gradually increased throughout the shelf life of the fruit. Notably, the changes in glucose and fructose contents in fruits at 80 DAT closely paralleled those observed in the total sugar content. Upon attaining complete maturation, the contents of fructose and glucose in the T2 treatment group exceeded those observed in other treatment groups (glucose: 78.66 ± 1.30 mg/g; fructose: 63.49 ± 1.03 mg/g). Notably, among all groups, the lowest glucose (65.36 ± 1.40 mg/g) and fructose (57.58 ± 1.01 mg/g) contents were observed in the T3 treatment group.

3.3. Effect of GA on organic acid content in grape berries

As illustrated in Figure 2D, the primary organic acids present in grape berries include tartaric acid, malic acid, and citric acid, with citric acid content being relatively minimal. During the period of shelf life, we noted a distinctive V-shaped trend in the tartaric acid content, which peaked at 20 DAT and plummeted at 80 DAT across all treatments. Prior to 60 DAT, significant differences were observed in the tartaric content among all four treatments; however, these differences were not evident at 100 DAT. At 100 DAT, we observed the lowest and highest tartaric acid contents under T1 (4.31 ± 0.08 mg/g) and T2 (4.38 ± 0.02 mg/g) treatments, respectively.
The malic acid content in grape berries initially increased and then decreased over the course of shelf life. Notably, a significant difference among the four treatments was evident at 40 DAT alone. During the maturation phase, the lowest malic acid content was observed under the T1 treatment (1.83 mg/g), which was 8.74% lower than that observed under the CK treatment.
The change trends in citric acid content paralleled those of malic acid content. At 40 DAT, the citric acid content in all groups reached its peak during the shelf life. Notably, significant differences in citric acid content were observed across all groups at 40 and 60 DAT. At 100 DAT, the lowest citric acid content was noted under the T2 treatment (0.012 mg/g), which was 8.33% lower than those observed under the T1 and T3 treatments.

3.4. Analysis of sugar unloading in phloem

Our results demonstrated that the primary sugars in grape berries were fructose and glucose, with relatively minimal sucrose content at the veraison stage. Notably, the fructose content exceeded that of glucose. As depicted in Figure 2E, the changes in fructose and glucose unloading in the phloem exhibited a double-peak curve at five distinct time points following GA treatment. In the morning, sugar unloading increased rapidly, reaching its peak at 12:00, closely aligned with the increase in photosynthetic rate. Subsequently, a sharp decline was observed at 14:00. Furthermore, a resurgence in sugar unloading was observed with an increase in photosynthetic rate, reaching a minor peak at 16:00, followed by a gradual decline at sunset. Notably, the morning unloading volume was approximately 1.1 times that of the afternoon unloading. However, the changes in fructose and glucose unloading in the phloem exhibited a unimodal curve at the five time points under the CK treatment. In the morning, sugar unloading increased rapidly, reaching its peak at 14:00, followed by a rapid decrement.
The maximum fructose unloading (3.89 mg/g) was observed between 9:30 and 10:00 in the T3 treatment fruits, significantly surpassing that observed in the CK treatment fruits (3.51 mg/g). For the T1 treatment, the peak fructose unloading (4.21 mg/g) occurred between 11:30 and 12:00. Similarly, the peak fructose unloading for the CK treatment reached 4.00 mg/g at 13:30–14:00 and that for the T2 treatment reached 3.89 mg/g and 3.30 mg/g at 15:30–16:00 and 17:30–18:00, respectively.
Furthermore, the peak glucose unloading for the T3 treatment occurred at 9:30–10:00 and 11:30–12:00, measuring 2.98 mg/g and 3.11 mg/g, respectively. For the CK treatment, the highest glucose unloading occurred at 13:30–14:00, following a pattern similar to that of fructose. Notably, the T3 treatment exhibited the highest glucose unloading (3.11 mg/g) at 15:30–16:00. In addition, at 17:30–18:00, glucose unloading was significantly higher in the T2 group than that observed in other groups. Overall, these findings indicate that the highest sugar (fructose+glucose) unloading (32.25 mg/g) was observed in the T2 treatment group during the veraison stage.
Similarly, in vivo analysis of sugar unloading at the maturation stage resembled that at the veraison stage. The primary sugars in grape berries were fructose and glucose, with relatively minimal amounts of sucrose owing to incomplete transformation. As illustrated in Figure 2E, the changes in fructose and glucose unloading in the phloem exhibited a unimodal curve at each time point during the maturation stage. Sugar unloading reached its peak at 10:00 in the morning, gradually declining to a minimum at 14:00. Subsequently, sugar unloading increased again with an increase in photosynthetic rate, achieving a minor peak at 16:00 and finally decreasing at sunset. The unloading value observed in the morning was 1 fold higher than that observed in the afternoon. Among the treatments, the T2 treatment demonstrated the highest fructose unloading at the first four time points, while the highest unloading under the T3 treatment was observed at 17:30–18:00. Moreover, the maximum glucose unloading under the T2 treatment was evident at 9:30–10:00 (3.09 mg/g), 13:30–14:00 (2.92 mg/g), and 17:30–18:00 (2.72 mg/g). The T3 treatment group exhibited the highest glucose unloading at 11:30–12:00 (2.91 mg/g) and 15:30–16:00 (3.03 mg/g).
Based on the sugar unloading data at various time points for each treatment, the T2 treatment displayed the highest sugar unloading (fructose+glucose) at the maturation stage, with a value of 31.26 mg/g.

3.5. Effect of GA on endogenous hormone contents in grape berries

Four types of endogenous hormones—GA, ABA, IAA, and CTK—were analyzed. The GA content demonstrated a fluctuating trend at 0–48 h after GA3 expansion treatment. It reached a peak at 0–4 h, followed by a sharp decline and a rapid increase at 24 h. In contrast, the GA content under the CK treatment exhibited a gradual increase at 48 h. Notably, the overall trend in GA content at 1–24 h was as follows: T3 > T2 > T1 > CK (Figure 3A). The ABA contents were significantly higher in CK treatment fruits than those in fruits subjected to the GA treatment. The lowest and highest ABA contents in berries were observed at 4 h and 24 h after GA treatment, respectively. At 48 h, the overall trend in ABA content was as follows: CK > T1 > T3 > T2 (Figure 3B). The IAA content in fruits displayed a trend similar to that observed for GA. The highest IAA content in berries under the T1 treatment was observed at 8 h, and at 48 h, the following overall trend was noted: T3 > T2 > T1 > CK (Figure 3C). The CTK content in berries increased rapidly at 0–2 h under CK, T2, and T3 treatments, followed by a slower upward trend that gradually decreased at 4 h. Moreover, the CTK content showed an upward trend under the T1 treatment at 0–8 h. At 48 h, the following trend was observed: T2 > T3 > T1 > CK (Figure 3D).

3.6. Correlation analysis

The correlation analysis depicted in Figure 4 examines the relationships among the contents of endogenous GA, soluble sugars, organic acids, and key parameters, namely, berry weight, berry shape index, TD, LD, volume, TSS, and TA. Through correlation analysis, we observed a significant positive correlation between the GA content and the key parameters, including berry weight, berry shape index, TD, LD, and volume, particularly with TD. Furthermore, a significant negative correlation was noted between GA and TSS contents at 60, 80, and 100 DAT, as well as between the GA content and TA at 20 and 40 DAT. In contrast, a significant positive correlation was observed between the GA content and TA at 60 and 80 DAT (Figure 4A).
The contents of reducing sugars, glucose, and fructose in fruits displayed a significant positive correlation with GA contents at 40 and 60 DAT (Figure. 4B); however, a significant negative correlation between them was observed at 80 and 100 DAT. Notably, the GA content was significantly negatively correlated with tartaric acid content at 20 DAT, malic acid content at 40, 60, 80, and 100 DAT, and citric acid content at 80 and 100 DAT (Figure. 4B). In contrast, the GA content was significantly positively correlated with tartaric acid content at 60, 80, and 100 DAT, malic acid content at 20 DAT, and citric acid content at 20, 40 and 60 DAT (Figure 4B).
The heatmap presented in Figure 4C indicates that IAA and CTK contents were significantly positively correlated with GA3 contents at 1, 2, 4, 8, 24, and 48 h. Furthermore, the correlation analysis highlighted that ABA contents were significantly negatively correlated with GA3 contents and IAA and CTK contents at 1, 2, 4, 8, 24, and 48 h. Notably, a significant positive correlation was observed between IAA and CTK contents.

4. Discussion

4.1. GA3 application by cluster dipping at different periods results in berry size enlargement in grapes

Our study conclusively showed that GA3 effectively enhanced the size of Shine Muscat berries. Intriguingly, we found that the effect of GA3 on LD was more significant than that on TD. Subsequently, we conducted additional analyses to measure berry weight and volume in order to investigate the enlargement effect. As anticipated, GA3 significantly enhanced both berry weight and volume, particularly under the seedlessness + expansion treatment (T2 treatment). This observation potentially indicates that the concentrations of compounds, such as sugars and acids, might have been regulated through expansion dilution.

4.2. GA application alters the concentrations of sugars and acids

To verify the hypothesis mentioned earlier, we conducted analyses of soluble sugar and organic acid contents. The application of GA led to a significant increase in TSS content in berries during 2021, which is consistent with the findings of Tyagi et al. [14]in Sangiovese grapes at 79 d. However, this effect was not observed during 2022, consistent with previous findings for Sangiovese grapes at 49 d[14]. We further observed lower TA in GA-treated berries than that in untreated control at 40 and 60 DAT in 2021 and 40 DAT in 2022. This finding aligns with that of Gao et al.[2]in Cabernet Sauvignon grapes. However, no significant difference in TA was observed at 80 and 100 DAT in 2021 and 80 DAT in 2022, aligning with the findings of previous research on Sangiovese grapes[14] . In commercial settings, the TSS/TA ratio is considered the most reliable indicator of fruit flavor[15]. This study demonstrates that GA-treated berries [Rachis elongation (2022.4.28)+Seedlessness (2022.5.24)+Expansion treatment (2022.6.7), or Seedlessness (2022.5.24)+Expansion treatment (2022.6.7)] exhibited a higher TSS/TA ratio at 80 DAT in both 2021 and 2022, suggesting improved berry flavor. Subsequently, we investigated the effect of GA application on soluble sugar and organic acid contents in fruits to thoroughly explore the potential existence of the expansion dilution effect.
It is widely accepted that soluble sugars in grape berries primarily consist of glucose and fructose. Although some studies have reported the presence of sucrose in certain table grape varieties[16], sucrose has not been detected in Cabernet Sauvignon grapes[7]. In the present study, traces of sucrose were detectable, albeit in limited quantities. Monosaccharides serve as effective osmotica in plants. Herein, GA treatment [Seedlessness (2022.5.24)+Expansion treatment (2022.6.7)] increased the accumulation of monosaccharides (fructose and glucose) in grape berries at 100 DAT. However, berries treated with expansion alone exhibited the lowest values, suggesting an effect of expansion dilution. This phenomenon may be influenced by other phytohormones mediated by GA or from the competition for photosynthates due to seed development. Our findings indicate that appropriate GA3 application (T2 treatment) can enhance sugar accumulation, overcoming the potential impact of expansion dilution.

4.3. Interactions between signaling pathways of GA and other phytohormones contribute to the coordinated development of berries

GA, in conjunction with IAA and CTK, plays a role in promoting cell division and expansion, leading to increased cell growth and protein synthesis[17]. Therefore, these hormones coordinate the development and enlargement of fruits following fertilization. Notably, these three plant hormones, along with ABA, have significant regulatory functions in the establishment and maintenance of fruit sink strength [11]. The application of exogenous hormones could potentially trigger changes in endogenous hormone levels in plants [18]. Therefore, in this study, we assessed alterations in the levels of four endogenous hormones (GA3, CTK, IAA, and ABA). Our findings indicated that GA3 application led to increased levels of IAA and CTK compared with those in the CK group, while also exerting an inhibitory effect on ABA production. The correlation analysis unveiled a significant positive correlation between GA3 and IAA/CTK contents and a significant negative correlation between GA3 and ABA contents. These results, in conjunction with prior research, suggest a potential crosstalk between GA and other plant hormones such as auxin and CTK[19]. Interestingly, the growth of grape berries induced by auxin and CK treatment is also partially dependent on GA biosynthesis [20], underscoring the significance of the synergy between these plant hormones. GA3 has been substantiated to enhance fruit size by promoting cell division and expansion, contributing to the mechanisms of cell wall formation and relaxation[21]. Previous studies have revealed that auxin modulates grape berry ripening, which could be associated to cell expansion[22,23]. Numerous studies have also provided support for the notion that CTK could govern critical rate-limiting steps in nutrient distribution and utilization[24,25]. These findings collectively suggest the significant contributions of GA, IAA, and CTK in modulating fruit sink strength.

4.4. GA application increases the sink strength in grapes

The application of gibberellin can lead to an increase in auxin content. Gibberellin, along with auxin and CTK, plays a role in nutrient transport and attraction, ensuring that the fruit becomes a robust “reservoir” and is positioned favorably in the nutritional competition. Sink organs in plants serve as net importers of assimilates. Throughout various stages of plant development, all plant organs can function as sinks, thereby receiving assimilates. In terms of the transport of assimilates, the capacity of a sink organ to import assimilates is referred to as its sink strength. However, a substantial proportion of the imported assimilates can be utilized for respiratory processes within the sink organs. Consequently, measuring the sink strength of a sink organ solely based on parameters such as the absolute growth rate or net accumulation rate of dry matter falls short of accurately gauging the actual capacity of a sink organ to receive assimilates. This conventional approach instead represents an apparent sink strength. Measuring the import rate of assimilates, which takes into account the sum of the net carbon gain and the carbon loss through respiration within a sink organ, would provide a more accurate estimate of the true sink strength. In addition, CTK has the potential to enhance sink capacity by promoting cell proliferation or sustaining sink activity through the regulation of sucrolytic enzymes, thereby allowing the acquisition of more photoassimilates.
Sink strength is generally defined as the competitive ability of an organ to attract assimilates, represented by the product of sink size and sink activity [26]. In various fruit crops, including grapes, GA has been shown to increase fruit size by stimulating cell division and elongation [27]. Our study demonstrated that GA3 application significantly increased berry size, berry weight, LD, and TD. This suggests that GA3 can bolster sink strength by actively contributing to the enhancement of sink size, encompassing both cell number and size. This augmentation, in turn, ensures the accrual of a greater amount of photoassimilates and dry matter within the sink structure. Similar findings were observed in pear (Pyrus pyrifolia), where GA treatment led to a significant increase in sink demand and fruit size[28]. The transport of photosynthates occurs through a permeation-driven loading or unloading process within the phloem, which is governed by the pressure gradient existed between the source and the phloem reservoir. Our findings demonstrate that an appropriate application of GA3 can expedite phloem unloading, resulting in enhanced sugar unloading in berry fruits during both the veraison and maturation stages. This enhancement implies an increase in sink activity. Apart from sugar translocation, the accumulation of sugars stands as a pivotal determinant of sink strength. Our study shows that GA3 application leads to the accumulation of substantial concentrations of glucose and fructose in berry fruits. This aligns with the findings of previous studies [29,30], which collectively support the notion that GA3 can affect sink strength by expediting the accumulation of hexose sugars. These findings imply that the phytohormone GA3 actively participates in the regulation of sink size and activity, encompassing an increase in berry size, sugar phloem unloading, and sugar accumulation in the sink cells. Altogether, these mechanisms significantly contribute to the modulation of sink strength in grapes (as depicted in Figure 5).

5. Conclusions

Collectively, this study introduced two significant and novel findings that contribute to our understanding of grape physiology and growth regulation: (i) GA3 exerted significant effects on the contents of soluble sugars, organic acids, and endogenous hormones (IAA, CTK, and ABA). The application of GA3 led to enhanced sugar unloading during the softening and ripening stages. (ii) The appropriate application of GA3 plays a crucial role in orchestrating the modulation of sink size and activity, including the enhancement of berry size, the facilitation of sugar phloem unloading, and the accumulation of sugars within sink cells. These insights collectively exert a robust effect on the overall sink strength in grape development. Moreover, these novel findings significantly enhance our comprehension of the complex interplay between GA3 and soluble sugar, organic acid, and endogenous hormone contents. This study presents compelling empirical evidence that contributes to the broader body of knowledge supporting effective strategies for table grape cultivation and the proficient utilization of GA3, ultimately benefiting the grape industry in China.

Author Contributions

Conceptualization, X.J.L., Z.S.X., L.L. and B.L.; methodology, X.J.L., Z.H.C. and Z.S.X.; sample preparation and analysis, Z.H.C.; investigation, X.L.L. and G.W.Y; data curation, Z.H.C. and S.X.L.; writing—original draft preparation, L.L.; writing—review and editing, B.L.; Supervision, X.L.L., Y.S.W. and Z.H.; project administration, B.L.; funding acquisition, Z.H., L.L. and B.L. All authors approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of Shandong Province (Grant No. ZR2023MC101, ZR2021QC165 and ZR2022QC185), Agricultural Science and Technology Innovation Project of Shandong Academy of Agricultural Sciences (Grant No. CXGC2023A47), Agricultural Variety Project of Shandong Province-New Varieties Grape Cultivation of High Quality with Characteristic and Industrialization Technology Innovation & Promotion (Grant No. 2022LZGCQY019), Collaborative Extension Project of Major Agricultural Technology of Shandong Province (Grant No. SDNYXTTG-2023-18), Key R&D Program of Shandong Province (Grant No. 2022TZXD0011) and National Natural Science Foundation of China (32202462).

Data Availability Statement

Data will be made available on request.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Cheng, C.; Jiao, C.; Singer, S. D.; Gao, M.; Xu, X.; Zhou, Y.; Li, Z.; Fei, Z.; Wang, Y.; Wang, X. Gibberellin-induced changes in the transcriptome of grapevine (Vitis labrusca × V. vinifera) cv. Kyoho flowers. BMC Genomics. 2015, 16, 128. [Google Scholar] [CrossRef]
  2. Gao, X.T.; Wu, M.H.; Sun, D.; Li, H.Q.; Chen, W.K.; Yang, H.Y.; Liu, F.Q.; Wang, Q.C.; Wang, Y.Y.; Wang, J.; He, F. Effects of gibberellic acid (GA3) application before anthesis on rachis elongation and berry quality and aroma and flavour compounds in Vitis vinifera L. 'Cabernet Franc' and 'Cabernet Sauvignon' grapes. J. Sci. Food Agric. 2020, 100, 3729–3740. [Google Scholar] [CrossRef] [PubMed]
  3. Grimplet, J.; Tello, J.; Laguna, N.; Ibáñez, J. Differences in flower transcriptome between grapevine clones are related to their cluster compactness, fruitfulness, and berry size. Front. Plant Sci. 2017, 8, 632. [Google Scholar] [CrossRef] [PubMed]
  4. Korkutal, I. ; Bahar, E; Gokhan, O. The characteristics of substances regulating growth and development of plants and the utilization of gibberellic acid (GA3) in viticulture. World Journal of Agricultural Sciences,2008, 4: 321–325.
  5. Sun, D. The effects of peduncle and berry quality with gibberellin treatments in the Cabernet Franc grape (Vitis vinifera L. cv.). Thesis for M.S., Northeast Forestry University. 2017.
  6. Wang, M.; Huang, L.P.; Liu, X.T.; Zhao, Q.F.; Ma, X.H. Effects of gibberellin on grape inflorescence elongation and berry growth treatment before flowering, Anhui Nongxue Tongbao, 2018, 24, 35-36.
  7. Yang, L.L.; Niu, Z.Z.; Wei, J.G.; Zhao, Y.Z.; Chen, Z.; Xuan, L.F.; Chu, F. J. Effects of exogenous GA3 on inflorescence elongation and gibberellin anabolism related gene expression in grapevine. Molecular Plant Breeding, 2020, 18,5600-5606.
  8. Han, D.H.; Lee, C.H. The effects of GA3, CPPU and ABA applications on the quality of kyoho (Vitis vinifera L. X V. Labrusca L.) grape. Acta Horticulturae, 2004, 653, 193-197. [CrossRef]
  9. Zhao, J. The influence of GA3, CPPU and streptomycin on the denuclearization and quality of grape. Master thesis of Shenyang Agricultural University. 2018. [Google Scholar]
  10. Ho, L.C. Metabolism and compartmentation of imported sugars in sink organs in relation to sink strength. Annu. rev. 1988, 39, 355–378. [Google Scholar] [CrossRef]
  11. Li, Y.M.; Forney, C.; Bondada, B.; Leng, F.; Xie, Z.S. The molecular regulation of carbon sink strength in grapevine (Vitis vinifera L. ). Front. Plant Sci. 2021, 11, 606918. [Google Scholar] [CrossRef] [PubMed]
  12. Ma, Q.J.; Sun, M.H.; Lu, J.; Liu, Y.J.; Hu, D.G.; Hao, Y.J. Transcription factor AREB2 is involved in soluble sugar accumulation by activating sugar transporter and amylase genes. Plant Physiol. 2017, 174, 2348–2362. [Google Scholar] [CrossRef] [PubMed]
  13. Wu, Y.; Li, X.; Zhang, W.; Wang, L.; Li, B.; Wang, S. Aroma profiling of Shine Muscat grape provides detailed insights into the regulatory effect of gibberellic acid and N-(2-chloro-4-pyridinyl)-N-phenylurea applications on aroma quality. Food Res. Int. 2023, 170, 112950. [Google Scholar] [CrossRef] [PubMed]
  14. Tyagi, K.; Maoz, I.; Lapidot, O.; Kochanek, B.; Butnaro, Y.; Shlisel, M.; Lerno, L.; Ebeler, S.; Lichter, A. Effects of gibberellin and cytokinin on phenolic and volatile composition of Sangiovese grapes. Sci. Hortic. 2022, 295, 110860. [Google Scholar] [CrossRef]
  15. Parson, S.; Paull, R. E. Pineapple organic acid metabolism and accumulation during fruit development. Sci. Hortic. 2007, 112, 297–303. [Google Scholar] [CrossRef]
  16. Shiraishi, M.; Shinomiya, R.; Chijiwa, H. Preliminary genetic analysis of sucrose accumulation in berries of table grapes. Sci. Hortic. 2012, 137, 0–113. [Google Scholar] [CrossRef]
  17. Spartz, A.K.; Lee, S.H.; Wenger, J.P.; Gonzalez, N.; Itoh, H.; Inzé, D.; Peer, W.A.; Murphy, A.S.; Overvoorde, P.J.; Gray, W.M. The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J. 2012, 70, 978–990. [Google Scholar] [CrossRef]
  18. El-Mergawi, R.A.; Abd El-Wahed, M.S.A. Effect of exogenous salicylic acid or indole acetic acid on their endogenous levels, germination, and growth in maize. Bull Natl Res Cent. 2020, 44, 1–8. [Google Scholar] [CrossRef]
  19. Wang, W.; Bai, Y.; Koilkonda, P.; Guan, L.; Zhuge, Y.; Wang, X.; Liu, Z.; Jia, H.; Wang, C.; Fang, J. Genome-wide identification and characterization of gibberellin metabolic and signal transduction (GA MST) pathway mediating seed and berry development (SBD) in grape (Vitis vinifera L. ). BMC Plant Biol. 2020, 20, 384. [Google Scholar] [CrossRef] [PubMed]
  20. Lu, L.; Liang, J.; Zhu, X.K.; Li, T.Z.; Hu, J.F. Auxin- and cytokinin-induced berries set in grapevine partly rely on enhanced gibberellin biosynthesis. Tree Genet. 2016, 12, 41. [Google Scholar] [CrossRef]
  21. Upadhyay, A.; Maske, S.; Jogaiah, S.; Kadoo, N.Y.; Gupta, V.S. GA3 application in grapes (Vitis vinifera L.) modulates different sets of genes at cluster emergence, full bloom, and berry stage as revealed by RNA sequence-based transcriptome analysis. Funct. Integr. Genomics. 2018, 18, 439–455. [CrossRef]
  22. Wong, D.C.; Lopez Gutierrez, R.; Dimopoulos, N.; Gambetta, G.A.; Castellarin, S.D. Combined physiological, transcriptome, and cis-regulatory element analyses indicate that key aspects of ripening, metabolism, and transcriptional program in grapes (Vitis vinifera L.) are differentially modulated accordingly to fruit size. BMC Genomics. 2016, 17, 416. [CrossRef]
  23. Dal Santo, S.; Tucker, M.R.; Tan, H.T.; Burbidge, C.A.; Fasoli, M.; Böttcher, C.; Boss, P.K.; Pezzotti, M.; Davies, C. Auxin treatment of grapevine (Vitis vinifera L.) berries delays ripening onset by inhibiting cell expansion. Plant Mol. Biol. 2020, 103, 91–111. [CrossRef]
  24. Body, M.; Kaiser, W.; Dubreuil, G.; Casas, J.; Giron, D. Leaf-miners co-opt microorganisms to enhance their nutritional environment. J. Chem. Ecol. 2013, 39, 969–977. [Google Scholar] [CrossRef] [PubMed]
  25. Werner, T.; Holst, K.; Pörs, Y.; Guivarc'h, A.; Mustroph, A.; Chriqui, D.; Grimm, B.; Schmülling, T. Cytokinin deficiency causes distinct changes of sink and source parameters in tobacco shoots and roots. J. Exp. Bot. 2008, 59, 2659–2672. [CrossRef]
  26. Farrar, J. Sink strength: what is it and how do we measure it? Introduction. Plant Cell Environ.1993, 16 ,1015-1016. [CrossRef]
  27. Zhang, C.; Whiting, M. Pre-harvest foliar application of Prohexadione-Ca and gibberellins modify canopy source-sink relations and improve quality and shelf-life of ‘Bing’sweet cherry. Plant Growth Regul. 2011, 65, 145–156. [Google Scholar] [CrossRef]
  28. Li, J.; Yu, X.; Lou, Y.; Wang, L.; Slovin, J.P.; Xu, W.; Wang, S.; Zhang, C. Proteomic analysis of the effects of gibberellin on increased fruit sink strength in asian pear (pyrus pyrifolia). Sci. Hortic. 2015, 195, 25–36. [Google Scholar] [CrossRef]
  29. Lecourieux, F.; Kappel, C.; Lecourieux, D.; Serrano, A.; Torres, E.; Arce-Johnson, P.; Delrot, S. An update on sugar transport and signalling in grapevine. J. Exp. Bot. 2014, 65, 821–832. [Google Scholar] [CrossRef] [PubMed]
  30. Zhang, Y.; Zhen, L.; Tan, X.; Li, L.; Wang, X. The involvement of hexokinase in the coordinated regulation of glucose and gibberellin on cell wall invertase and sucrose synthesis in grape berry. Mol. Biol. Rep. 2014, 41, 7899–7910. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Berry weight(A) and volume(B) in 2021and 2022, berry growth at three stages: hardcore stage, softening stage, and harvest stage(C), transverse diameter(D), longitudinal diameter(E), and berry shape index(F) in 2021and 2022 of “Shine Muscat” under various gibberellin treatments.Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Figure 1. Berry weight(A) and volume(B) in 2021and 2022, berry growth at three stages: hardcore stage, softening stage, and harvest stage(C), transverse diameter(D), longitudinal diameter(E), and berry shape index(F) in 2021and 2022 of “Shine Muscat” under various gibberellin treatments.Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Preprints 99360 g001
Figure 2. Total soluble solid (TSS) content(A), titratable acidity (TA)(B), and TSS/TA ratio(C) in 2021 and 2022. Contents of soluble sugars and organic acids (D), and sugar phloem unloading in “Shine Muscat” grapes during veraison and harvest stages under various gibberellin treatments in “Shine Muscat” grapes. Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Figure 2. Total soluble solid (TSS) content(A), titratable acidity (TA)(B), and TSS/TA ratio(C) in 2021 and 2022. Contents of soluble sugars and organic acids (D), and sugar phloem unloading in “Shine Muscat” grapes during veraison and harvest stages under various gibberellin treatments in “Shine Muscat” grapes. Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Preprints 99360 g002
Figure 3. Contents of endogenous hormones—gibberellic acid (GA3) (A), abscisic acid (ABA) (B), indole-3-acetic acid (IAA) (C), and cytokinin (CTK) (D)—in “Shine Muscat” grapes under various gibberellin treatments. Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Figure 3. Contents of endogenous hormones—gibberellic acid (GA3) (A), abscisic acid (ABA) (B), indole-3-acetic acid (IAA) (C), and cytokinin (CTK) (D)—in “Shine Muscat” grapes under various gibberellin treatments. Different letters denote statistically significant differences among treatments at the same period, as determined by Duncan’s test (p < 0.05).
Preprints 99360 g003
Figure 4. Correlation analysis of GA3 content with key parameters (A), soluble sugar and organic acid contents (B), and endogenous hormone contents (C) from grapes under various GA3 treatments.
Figure 4. Correlation analysis of GA3 content with key parameters (A), soluble sugar and organic acid contents (B), and endogenous hormone contents (C) from grapes under various GA3 treatments.
Preprints 99360 g004
Figure 5. Mechanism through which GA3 regulates sink strength via modulation of sink size and activity, encompassing an increase in berry size, sugar phloem unloading, and sugar accumulation in the sink cells in grapes.
Figure 5. Mechanism through which GA3 regulates sink strength via modulation of sink size and activity, encompassing an increase in berry size, sugar phloem unloading, and sugar accumulation in the sink cells in grapes.
Preprints 99360 g005
Table 1. Experimental design of different gibberellin treatments.
Table 1. Experimental design of different gibberellin treatments.
NO. Treatment
CK Water
T1 Rachis elongation treatment + Seedlessness treatment + Expansion treatment
T2 Seedlessness treatment + Expansion treatment
T3 Expansion treatment
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated