Profile of Semiquinone Radicals, Phytohormones and Sugars in Pistacia vera L. cv. Kirmizi Development

The main aim of this study was to investigate metabolic responses of fruits, leaves, and shoots of pistachio trees (Pistacia vera L. cv. Kirmizi) during their development. Electron paramagnetic resonance spectroscopy revealed significant increase in generation of semiquinone radicals in fruits and leaves of pistachio, while the flower cluster thinning application was conducted in relation to the control; especially at the second term of the plant material collection. Moreover, flower abscission caused an increase in the level of phytohormones such as indole-3-acetic acid and abscisic acid in fruits at the first term of fruits sampling. In turn, high-performance liquid chromatography analysis revealed differences both in the profile as well as the contents of soluble sugars detected in pistachio organs. The highest total sugar content was found in fruits of pistachio where the flower cluster thinning application was made early. Total sugar levels were higher also in leaves and shoots of the above-mentioned pistachio plants in relation to the control until the third sampling time. The importance of high levels of fructose in the fruits following the flower abscission was observed. Additionally, analyses of mineral elements in organs showed that copper and phosphorus contents in fruits were higher after the flower abscission in relation to the control. To conclude, our findings signal on contribution of semiquinone radicals, paramagnetic manganese ions, phytohormones, nutrients, and sugars in pistachio organs development on the background of the flower cluster thinning which was applied before fructification.


Introduction
Pistacia L. belongs to the family the Anacardiaceae (cashew family) [1] or order Sapindales [2]. It contains nine species and five subspecies [3]. Pistacia vera L. is believed to be the most ancestral species, whereas other species probably derived from it [4,5]. Recent research results of Karcı et al. [6] provided evidence that simple sequence repeats markers can be used to the characterization and the phylogenetic analysis of Pistacia species and cultivars. Moreover, P. vera is a small, dioecious tree originating from Central Asia and the Middle East. The Anacardiaceae family contains also important other species such as mango, pepper tree and sumac [7]. The genus Pistachio is a xerophytic plant [8]. Interestingly, P. vera is also a glycophyte highly tolerant to salt because pistachios have the capacity to maintain relatively high photosynthetic rates during treatment with salts [9]. For these reasons, P. vera trees are remarkably interesting from the point of plant physiology and molecular biology. These trees require specific climatic conditions such as long, hot, dry summers and chilling in the winter, but do not tolerate ground that freezes. Growth requires a temperature above 38 • C during the day and needs also low temperatures under 7 • C during winter months for break bud dormancy [10]. It should also be mentioned that the environment of pistachio trees needs to be arid, because these trees do not do well in areas of high humidity. Besides, P. vera grows do well in all soil types but really thrives in deep, sandy loam. Additionally, there are reports of pistachio trees being resistant to cold and wind but also sensitive to very extreme climate conditions, such as extreme drought, prolonged frost and excessive dampness and high humidity [11]. Key abiotic factors such as high temperature [12] and water deficit were found to lead to increased shell spliting of pistachio nuts and allow pathogenic fungi to colonize the nuts and contaminate them with mycotoxins [13,14].
P. vera is one of the most popular tree nuts in the world, which grows in defined climatic conditions. Source literature points to two centers of diversity of cultivated pistachio. As reported by Kafkas [5] (2019), the first comprises the Mediterranean region of Europe, Northern Africa, and the Middle East countries. In turn, the second comprises the Eastern part of Zagros Mountains from Crimea to the Caspian Sea. Moreover, pistachio cultivation extended westward from its center of origin to Italy, Spain, and other Mediterranean regions of Southern Europe as well as to North Africa, the Middle East, China, more often the United States and Australia [15][16][17].
It should be noted that in the literature, there are few molecular and physiological studies regarding the development of pistachio trees, and they are fragmentary. Pistachio needs almost 10 years to start an early fruit period and 20 years to enter full production period [11]. In the case of these trees, alternate bearing phenomena occurs, which is characterized by the fact that in the bearing year the production is very good, while in the alternative year the production is limited [18]. Research by Gündeşli [19] revealed the important role of phytohormones on the alternate bearing of pistachio (P. vera cv. Uzun) during different growth periods. The above author reported that gibberellin and abscisic acid (ABA) metabolites occurring in different organs play a key role in the control of pistachios during embryo development and flower bud abscission. Additionally, the importance of free polyamines in this phenomenon was also raised by Gündeşli [19]. Moreover, Okay et al. [20] reported that externally application of plant growth regulators to pistachio leads to regular yields, increases flower bud formation and prevents flower bud abscission. The above cited authors carefully analyzed seasonal changes of plant growth regulators to determine some physiological mechanisms like formation of flower buds, flower and fruit pouring and alternate bearing. Transcriptome analysis of P. vera L. cv. Bianca inflorescence buds in bearing and non-bearing shoots was carried out by Benny et al. [21], revealing differentially expressed genes, most of which being involved in sugar metabolism, plant hormone pathways, secondary metabolism and oxidative stress pathways.
It is extremely important and suitable to provide additional evidence regarding the mechanisms of pistachio tree fruiting versus time. The present work thus highlights for the first-time changes in the level of semiquinone radicals (stable radicals of organic origin) generation and paramagnetic center manganese ions (Mn 2+ ) in P. vera L. cv. Kirmizi organs at different developmental stages, where the flower cluster thinning was made before fructification. The electron paramagnetic resonance (EPR) method is one of the most sensitive techniques for detecting particles with unpaired electrons at concentrations as low as a fraction of ppm. Particles with unpaired electrons are rare in nature, but they include some transition metal ions (for the example Mn 2+ ) and free radicals that have been studied in the present work. It should also be mentioned that EPR spectroscopy is a unique technique not only used to detect, identify, quantify free radicals, study molecular structures, geometry, and dynamics, but to observe labeled species in situ in biological system and understanding redox processes, reaction kinetics and catalysis.
Moreover, in the present study it was also important to understand the sequence of generation of these semiquinone radicals and the time-dependent aspect of these changes. It is well known that semiquinone and phenoxyl radicals are formed in plant cells as a result of the oxidation of the hydroxyl groups of phenols and polyphenols being their generation positively associated with other markers of oxidation, such as protein carbonyls and peroxides [22]. Our previous research has shown that these radicals are involved in the direct defense response of plants against biotic factors as well as being involved in strengthening the cell wall [23][24][25][26][27][28][29]. Additionally, we used EPR spectroscopy in order to highlight biochemical and physical properties of apple fruits in relation to differences in semiquinone radical concentration between organs of different apple tree varieties [30]. Moreover, one of the goals of the present work was also to examine the levels of signalling molecules such as phytohormones (i.e., indole-3-acetic acid, IAA; ABA and salicylic acid, SA), soluble sugars (sucrose, monosacharides and galactose) and nutrient elements in pistachio organs where the flower cluster thinning application was made before fructification. Changes in the level of signaling molecules, i.e., phytohormones and sugars that are involved in the plant signalling network affect the regulation of metabolism of pistachio trees. As reported by Kozlowski [31] in perennial woody plants, most carbohydrates are produced in leaves but some are also synthesized in buds, twigs, stems, flowers and fruits. Internal competition for carbohydrates from sources to sink organs occurs, due to changes in rates of carbohydrate movement and reversals in direction of carbohydrate transport, as the relative sink strengths of various organs change. Previous studies of Morkunas et al. [32] have provided an extensive body of information on the contribution of carbohydrates to metabolic reactions of plant cells. In plants, sugars are essential components as respiratory substrates for the generation of energy and metabolic intermediates that are then used for the synthesis of macromolecules [33,34]. Sugars have also important hormone-like functions as physiological signals, which cause activation or repression of many plant genes, and this in turn leads to specific metabolic effects.
In the published literature, there are few reports on molecular and metabolomic changes occurring during the development of pistachio trees. Identifying these changes both at the molecular and metabolic levels and understanding the interrelationship between the studied parameters will afford new knowledge in the context of the development of this plant. In this study, we focused mainly on pistachio fruits, but additional analyses were also performed on leaves and shoots to obtain comprehensive information on the relationships between all these organs.

Plant Material and Growth Conditions
This study was conducted in 2018 in pistachio orchard ( Figure 1) (geographical location: 36 • 58 35.6" northern latitude and 37 • 51 34.9" eastern longitude). It should be mentioned that 'Kirmizi' variety mainly grows in the high-altitude parts of the Gaziantep region and is preferred due to its early ripening (i.e., beginning of September) [35]. with a coffee grinder and stored at −80 °C until analyzed. The analyses were carried out in the Department of Plant Physiology, Faculty of Agronomy and Bioengineering of Poznań University of Life Sciences, Institute of Molecular Physics Polish Academy of Sciences in Poznań and Department of Plant Physiology and Biotechnology, Faculty of Biological and Veterinary Sciences of Nicolaus Copernicus University in Poland and in the Plant Nutrition Physiology and Chromatography Laboratories of the Cukurova University, Faculty of Agriculture, Department of Horticulture in Turkey.   Samples of leaves, fruits ( Figure 2) and shoots from 30-year-old 'Kirmizi' pistachio trees from Nizip region of Gaziantep, in 4 different periods from the beginning of May to the end of July from the thinning applied and non-applied branches, were used. Applied trees were thinned 1 month before first sampling date. Leaf, fruit ( Figure 2) and shoot samples of this variety were brought from Gaziantep to Adana province in the cold chain condition and after harvesting immediately treated with liquid nitrogen and dried by freeze drying method (Ilshın FD5508 Bondıro Vacuum Freeze Dryer). All samples were homogenized with a coffee grinder and stored at −80 • C until analyzed. Samples of leaves, fruits ( Figure 2). and shoots from 30-year-old 'Kirmizi' pistachio trees from Nizip region of Gaziantep, in 4 different periods from the beginning of May to the end of July from the thinning applied and non-applied branches, were used. Applied trees were thinned 1 month before first sampling date. Leaf, fruit and shoot samples of this variety were brought from Gaziantep to Adana province in the cold chain condition and after harvesting immediately treated with liquid nitrogen and dried by freeze drying method (Ilshın FD5508 Bondıro Vacuum Freeze Dryer). All samples were homogenized with a coffee grinder and stored at −80 °C until analyzed. The analyses were carried out in the

Climatic Data of Research Area
Climatic data for the research area were obtained from the General Directorate of Meteorology. Monthly and yearly average temperatures, maximum and minimum temperatures, the average hours of sunshine and mean rainfall values in the period of 2018-2019 years for Gaziantep province are given in Table 1.

Electron Paramagnetic Resonance (EPR)
Samples (leaves, shoots, and fruits) of P. vera were lyophilized and the lyophilized plant material was transferred to EPR-type quartz tubes of 4 mm in diameter. EPR measurements were performed at room temperature with a Bruker ELEXSYS X-band spectrometer (Bruker, Billerica, MA, USA). The EPR spectra were recorded as first derivatives of microwave absorption. A microwave power of 2 mW and a magnetic field modulation of about 2 G was used for all the experiments to avoid signal saturation and deformation. EPR spectra of free radicals and Mn 2+ were recorded in the magnetic field range of 1400-4200 G.
To determine the number of paramagnetic centres in the samples, the spectra were doubleintegrated and compared with the intensity of the standard Al 2 O 3 :Cr 3+ single crystal with a known spin concentration [24,29,[37][38][39][40][41]. Before and after the first integration, some background corrections of the spectra were made to obtain a reliable absorption signal before the second integration. Free radicals gave signals with one g-values of 2.0045 ± 0.0005 while signals of Mn 2+ were recorded with g-values of 2.00 ± 0.01 in leaves, shoots, and fruits.

Detection of Phytohormones
Samples (leaves, shoots, and fruits) of P. vera were homogenized in liquid nitrogen and extracted twice with 20 mL 80% (v/v) methanol containing 20 mg/L butylated hydroxytoluene. Subsequently deuterated internal standards of ABA, [ 2 H6]ABA, IAA, [ 2 H2]IAA and SA, [ 2 H4]SA (100 ng of each) were added as internal standards. The mixture was shaken overnight in a cold room. The extract was then centrifuged, and the supernatant was collected in evaporation flasks. The remaining pellet was resuspended two times in 10 mL of 80% (v/v) methanol, shaken for two hours and centrifuged again. The combined extracts were evaporated until all methanol was removed. The flasks were washed with additional H 2 O and the extract was collected in 50 mL Falcon tubes. The appropriate amount of HCl was added for acidification until pH 2 was reached. This was checked with colored pH test strips. Then the extract was centrifuged, and the supernatant was collected into Erlenmeyer for solvent partition. The supernatant was partitioned three times against ethyl acetate. The organic phase was collected in evaporation flasks and dried under vacuum. The dry residue was then resolved in 5 mL 1 M formic acid (FA) and loaded on a Discovery ® DSC-18 SPE, 500 mg SPE cartridge (Supelco Inc., Bellefonte, PA, USA). These cartridges (columns) were preconditioned with 4 mL methanol and allowed to equilibrate with 4 mL of 1 M FA. The columns were then washed with 4 mL 1 M FA, Agronomy 2021, 11, 2115 6 of 20 4 mL of 20% methanol in 1 M FA and finally phytohormones were eluted with 4 mL of 80% methanol in water. The eluate was evaporated to dryness, resolved in 200 µL of 20% acetonitrile in water and applied to high-performance liquid chromatography (HPLC) for further sample purification. HPLC was done with a SUPELCOSIL ABZ+ PLUS column (250 × 4.5 mm, 5 µm particle size, Supelco). The samples were analyzed with a linear gradient of 20-80% acetonitrile in 0.1 M FA in 20 min, flow rate 1.0 mL/min at temperature 22 • C. The fractions collected at 12 ± 0.5 min were evaporated to dryness and transferred to appropriate small glass tubes with 100% methanol. The methanol was then evaporated in N 2 stream. Next the residue was methylated with diazomethane dissolved in diethyl ether. After this, the diethyl ether was evaporated with N 2, and the residue was dissolved in 30 µL of 100% methanol. 1 µL of this sample was analyzed by gas chromatography (GC)/mass spectrometry (MS) (Auto-System XL coupled to a TurboMass, Perking-Elmer) using an MDN-5 column (30 m × 0.25 mm, 0.25 µm phase thickness; Supelco). The GC temperature program was set at 60 • C for 1 min, 60-250 • C at 10 • C/min (5 min hold), flow rate 1.5 mL/min, injection port was 280 • C, electron potential 70

Determination of Soluble Sugar Content
Sugar contents (i.e., glucose, fructose, sucrose, and total sugar concentration) in dried leaf, shoot and fruit samples of the Kırmızı variety were determined according to the extraction method developed by Miron and Scahffer [43] using HPLC (HP 1100 Series), RID (Refractive Index Detector) detector and Transgenomic CARBOSep COREGEL 87C (250 × 4.6 mm, 5 µ) column in 3 repetitions. Distilled water was used as mobile phase and HPLC conditions were 1.0 mL/min flow and 85 • C oven temperature. Sugar contents in the samples were determined qualitatively and quantitatively using a refractive index detector according to the calibration curves created using the external standard and the retention time of the standard.

Nitrogen Analysis
Nitrogen in leaf samples was determined by the Kjeldahl method, which is a wet burning method recommended by Jacobs [44]. Dried and homogenized samples were weighed 200 mg for each repetition and 1 Kjeldahl tablet and 12 mL of sulfuric acid were added to nitrogen combustion tubes and put into the incinerator. Samples were burned for a total of 1 h: 10 min at 30 • C, 10 min at 60 • C and 40 min at 100 • C. The glass tube that completes the burning process was placed in the distillation device and the distillation process was started by placing 25 mL of boric acid in the flask into the other compartment. 0.1 N HCl was added dropwise to the sample, which took the color of onion peel after distillation, and titrated until it turned into the color of boric acid before the reaction.

Phosphorus Analysis
The phosphorus content was determined by the method proposed by Barton [45]. 200 mg of dried and homogenized samples were taken and burned in as furnace at 550 • C for 8 h. The ash produced after burning was dissolved in 1/3 HCl and filtered in blue band. The mixture prepared with 0.5 µL Bartone' solution and 4 mL of distilled water was vortexed. Vortexed samples were placed on plate and absorbance were measured in a spectrophotometer at a wavelength of 430 nm. 50 ppm of phosphorus were used for standard preparation with following concentrations 0%, 8.33%, 16.6%, 33.3%, 50%, 66.6%, 83.3%, 100%.

Macro and Micronutrient Analysis
200 mg of dried and homogenized samples were taken and burned in an ash furnace at 550 • C for 8 h. The ash produced after the combustion was dissolved in 1/3 HCl and filtered in the blue band and the K, Ca, Mg, Fe, Zn, Mn, and Cu concentrations were determined in the Varian brand Atomic Absorption Spectrometer [46].

Statistical Analysis
Three biological replicates per experimental variant were performed for each analysis. The normality of the distribution of the 16 observed traits was tested with Shapiro-Wilk's normality test [47] to check whether the analysis of variance (ANOVA) met the assumption that the ANOVA model residuals followed a normal distribution. The homogeneity of variance was tested using Bartlett's test. Multivariate normality and homogeneity of variance-covariance matrices were tested by Box's M test. A two-way (part, flower) multivariate analysis of variance (MANOVA) was made. Following this, two-way analyses of variance (ANOVA) were performed in order to verify the null hypotheses of a lack of part and flower effects as well as part × flower interaction effect in terms of the values of the 16 observed traits, independently for each trait. The arithmetic means and standard deviations were calculated. Moreover, Fisher's least significant differences (LSDs) were estimated at a significance level of α = 0.05. The relationships between the observed traits were estimated using Pearson's correlation coefficients. The GenStat v. 18 statistical software package was used for all the analyses.

Effect of Flower Cluster Thinning Application on the Concentration of Semiquinone Radicals and Mn 2+ in Organs of Pistachio Trees
In the first two sampling dates (i.e., 3rd May and 1st June), EPR analyses showed significant increased generation of semiquinone radicals in fruits and leaves of pistachio trees, where the flower cluster thinning application was made in relation to the control (i.e., without any flower cluster thinning application) (Figure 3a). In turn, at the third period (26th June), a significant decrease in the concentration of semiquinone radicals is visible in fruits and shoots, especially in trees receiving the flower cluster thinning application. Moreover, the same trend was recorded in fruits and shoots in the recent period of collection of plant material from trees with the flower cluster thinning. Besides, the highest level of radicals was observed in the shoots of pistachio trees without any flower cluster thinning application.

Macro and Micronutrient Analysis
200 mg of dried and homogenized samples were taken and burned in an ash furnace at 550 °C for 8 h. The ash produced after the combustion was dissolved in 1/3 HCl and filtered in the blue band and the K, Ca, Mg, Fe, Zn, Mn, and Cu concentrations were determined in the Varian brand Atomic Absorption Spectrometer [46].

Statistical Analysis
Three biological replicates per experimental variant were performed for each analysis. The normality of the distribution of the 16 observed traits was tested with Shapiro-Wilk's normality test [47] to check whether the analysis of variance (ANOVA) met the assumption that the ANOVA model residuals followed a normal distribution. The homogeneity of variance was tested using Bartlett's test. Multivariate normality and homogeneity of variance-covariance matrices were tested by Box's M test. A two-way (part, flower) multivariate analysis of variance (MANOVA) was made. Following this, two-way analyses of variance (ANOVA) were performed in order to verify the null hypotheses of a lack of part and flower effects as well as part × flower interaction effect in terms of the values of the 16 observed traits, independently for each trait. The arithmetic means and standard deviations were calculated. Moreover, Fisher's least significant differences (LSDs) were estimated at a significance level of α = 0.05. The relationships between the observed traits were estimated using Pearson's correlation coefficients. The GenStat v. 18 statistical software package was used for all the analyses.

Effect of Flower Cluster Thinning Application on the Concentration of Semiquinone Radicals and Mn 2+ in Organs of Pistachio Trees
In the first two sampling dates (i.e., 3rd May and 1st June), EPR analyses showed significant increased generation of semiquinone radicals in fruits and leaves of pistachio trees, where the flower cluster thinning application was made in relation to the control (i.e., without any flower cluster thinning application) (Figure 3a). In turn, at the third period (26th June), a significant decrease in the concentration of semiquinone radicals is visible in fruits and shoots, especially in trees receiving the flower cluster thinning application. Moreover, the same trend was recorded in fruits and shoots in the recent period of collection of plant material from trees with the flower cluster thinning. Besides, the highest level of radicals was observed in the shoots of pistachio trees without any flower cluster thinning application. In turn, the concentration of Mn 2+ showed fluctuations. The flower cluster thinning application caused decrease in the level of Mn 2+ in leaves of pistachio trees in the first term (3rd May), while an increase of these radicals occurred in shoots of pistachio trees with the flower cluster thinning as compared to the control (Figure 3b). At the third period (26th June), increases in concentrations of Mn 2+ were recorded in fruits and leaves. The highest level of these ions (21 × 10 14 Mn 2+ g −1 DW) was found in leaves of pistachio because of the flower cluster thinning application.

IAA Content
At the first stage (in May), the IAA content in fruits of pistachio trees where the flower cluster thinning application was made, was 4-fold higher than in the control. A similar general tendency was noted in leaves and shoots at certain time points. The highest level of IAA was recorded in the shoots and fruits of pistachio in the control (Figure 4). In turn, fruits of pistachio trees at the third time point (where the flower cluster thinning application was done) displayed 3-fold lower concentrations of IAA than the control trees.

ABA Content
After the flower cluster thinning application, the highest level of ABA was found in pistachio fruits (575 ng ABA•g −1 DW) at the first stage (in May) when nuts were collected ( Figure 5). Also, a high ABA content was noted in pistachio shoots in the last period of collecting samples, i.e., in July. Attention should also be drawn to the lower ABA level in fruits of pistachio trees from 1st June and leaves from 3rd May, where the flower cluster thinning application to pistachio trees was made. In turn, the concentration of Mn 2+ showed fluctuations. The flower cluster thinning application caused decrease in the level of Mn 2+ in leaves of pistachio trees in the first term (3 May), while an increase of these radicals occurred in shoots of pistachio trees with the flower cluster thinning as compared to the control (Figure 3b). At the third period (26 June), increases in concentrations of Mn 2+ were recorded in fruits and leaves. The highest level of these ions (21 × 10 14 Mn 2+ g −1 DW) was found in leaves of pistachio because of the flower cluster thinning application.

IAA Content
At the first stage (in May), the IAA content in fruits of pistachio trees where the flower cluster thinning application was made, was 4-fold higher than in the control. A similar general tendency was noted in leaves and shoots at certain time points. The highest level of IAA was recorded in the shoots and fruits of pistachio in the control (Figure 4). In turn, fruits of pistachio trees at the third time point (where the flower cluster thinning application was done) displayed 3-fold lower concentrations of IAA than the control trees. In turn, the concentration of Mn 2+ showed fluctuations. The flower cluster thinning application caused decrease in the level of Mn 2+ in leaves of pistachio trees in the first term (3rd May), while an increase of these radicals occurred in shoots of pistachio trees with the flower cluster thinning as compared to the control (Figure 3b). At the third period (26th June), increases in concentrations of Mn 2+ were recorded in fruits and leaves. The highest level of these ions (21 × 10 14 Mn 2+ g −1 DW) was found in leaves of pistachio because of the flower cluster thinning application.

IAA Content
At the first stage (in May), the IAA content in fruits of pistachio trees where the flower cluster thinning application was made, was 4-fold higher than in the control. A similar general tendency was noted in leaves and shoots at certain time points. The highest level of IAA was recorded in the shoots and fruits of pistachio in the control (Figure 4). In turn, fruits of pistachio trees at the third time point (where the flower cluster thinning application was done) displayed 3-fold lower concentrations of IAA than the control trees.

ABA Content
After the flower cluster thinning application, the highest level of ABA was found in pistachio fruits (575 ng ABA•g −1 DW) at the first stage (in May) when nuts were collected ( Figure 5). Also, a high ABA content was noted in pistachio shoots in the last period of collecting samples, i.e., in July. Attention should also be drawn to the lower ABA level in fruits of pistachio trees from 1st June and leaves from 3rd May, where the flower cluster thinning application to pistachio trees was made.

ABA Content
After the flower cluster thinning application, the highest level of ABA was found in pistachio fruits (575 ng ABA·g −1 DW) at the first stage (in May) when nuts were collected ( Figure 5). Also, a high ABA content was noted in pistachio shoots in the last period of collecting samples, i.e., in July. Attention should also be drawn to the lower ABA level in fruits of pistachio trees from 1st June and leaves from 3rd May, where the flower cluster thinning application to pistachio trees was made.

SA Content
It should be emphasized that after the flower cluster thinning application, the level of SA in pistachio organs (fruits, leaves and shoots) was lower than in the control at certain time points (Figure 6). For example, in June, a significant reduction in the SA concentration was noted in the fruits of pistachio trees before applying the flower cluster thinning as compared to the control. In turn, a significant reduction in the SA level was also observed in leaves especially in samples from the first, third and fourth time of plant material collection. Besides, the highest level of SA was recorded in fruits.

Changes in Sugar Levels after the Flower Cluster Thinning Application in Pistachio Organs
The total sugar content in pistachio organs after the flower cluster thinning application was generally higher than in the control (Figure 7). The significantly higher level of total sugars in the fruits of the pistachio trees following flower cluster thinning than in the control at the first time of nut harvest should draw attention (Figure 7a). A similar tendency at the next time point (on 1st June) in fruits was observed. Statistical analysis showed highly significant differences in these results (Table S1). Besides, at the last sampling time, the total sugar concentration was about 1.8 fold higher in shoots of the pistachio trees with flower cluster thinning than in the control. It should be emphasized that an extremely high fructose content in the nuts of pistachio trees after flower cluster thinning application was observed (2.7-fold higher than in the control) (Figure 7c). Moreover, sucrose concentration in fruits of pistachio trees following flower cluster thinning application that were collected in May, was also 1.2-fold higher than in the control (Figure 7d).

SA Content
It should be emphasized that after the flower cluster thinning application, the level of SA in pistachio organs (fruits, leaves and shoots) was lower than in the control at certain time points (Figure 6). For example, in June, a significant reduction in the SA concentration was noted in the fruits of pistachio trees before applying the flower cluster thinning as compared to the control. In turn, a significant reduction in the SA level was also observed in leaves especially in samples from the first, third and fourth time of plant material collection. Besides, the highest level of SA was recorded in fruits.

SA Content
It should be emphasized that after the flower cluster thinning application, the level of SA in pistachio organs (fruits, leaves and shoots) was lower than in the control at certain time points (Figure 6). For example, in June, a significant reduction in the SA concentration was noted in the fruits of pistachio trees before applying the flower cluster thinning as compared to the control. In turn, a significant reduction in the SA level was also observed in leaves especially in samples from the first, third and fourth time of plant material collection. Besides, the highest level of SA was recorded in fruits.

Changes in Sugar Levels after the Flower Cluster Thinning Application in Pistachio Organs
The total sugar content in pistachio organs after the flower cluster thinning application was generally higher than in the control (Figure 7). The significantly higher level of total sugars in the fruits of the pistachio trees following flower cluster thinning than in the control at the first time of nut harvest should draw attention (Figure 7a). A similar tendency at the next time point (on 1st June) in fruits was observed. Statistical analysis showed highly significant differences in these results (Table S1). Besides, at the last sampling time, the total sugar concentration was about 1.8 fold higher in shoots of the pistachio trees with flower cluster thinning than in the control. It should be emphasized that an extremely high fructose content in the nuts of pistachio trees after flower cluster thinning application was observed (2.7-fold higher than in the control) (Figure 7c). Moreover, sucrose concentration in fruits of pistachio trees following flower cluster thinning application that were collected in May, was also 1.2-fold higher than in the control (Figure 7d).

Changes in Sugar Levels after the Flower Cluster Thinning Application in Pistachio Organs
The total sugar content in pistachio organs after the flower cluster thinning application was generally higher than in the control (Figure 7). The significantly higher level of total sugars in the fruits of the pistachio trees following flower cluster thinning than in the control at the first time of nut harvest should draw attention (Figure 7a). A similar tendency at the next time point (on 1st June) in fruits was observed. Statistical analysis showed highly significant differences in these results (Table S1). Besides, at the last sampling time, the total sugar concentration was about 1.8 fold higher in shoots of the pistachio trees with flower cluster thinning than in the control. It should be emphasized that an extremely high fructose content in the nuts of pistachio trees after flower cluster thinning application was observed (2.7-fold higher than in the control) (Figure 7c). Moreover, sucrose concentration in fruits of pistachio trees following flower cluster thinning application that were collected in May, was also 1.2-fold higher than in the control (Figure 7d). Agronomy 2021, 11, x FOR PEER REVIEW 10 of 21

Effect of Flower Cluster Thinning Application on the Concentration of Mineral Elements in Pistachio Organs
The Cu content in fruits was higher after flower cluster thinning application as compared to the control (i.e., without cluster thinning) (Figure 8). In leaves, the Cu content was higher at the second and third period in the pistachio trees after flower cluster thinning and then decreasing for the last sampling point (Figure 8a). In shoots, the content of Cu was higher on the second harvest point in plants after the flower cluster thinning, while it was lower for the third and fourth harvest dates compared to the appropriate controls. Moreover, the Mg content in fruits was always lower after the flower cluster thinning application in relation to the control (Figure 8b). Also, in leaves of the above pistachio trees, the content of the Mg was lower at the first and second harvest dates. A similar tendency was noted in shoots at the last time point. It shoud be stressed that the flower cluster thinning leads to a strong accumulation of the Mn in fruits of pistachio trees at the last stage (in July); the concentration of the Mn was 13 times higher in these fruits than in the control (Figure 8c). In leaves and shoots of pistachio trees with flower abscission, the Mn concentration was generally lower than in the control. In fruits, the N content was lower at the first time point in pistachio trees after the flower cluster thinning, while it was higher at all other time points in the treated trees compared to the untreated ones ( Figure 8d). In leaves, the content of N was generally higher at the first, third and fourth dates in treated trees in relation to the appropriate controls. As for the shoot content, it was noticeable that the N level was lower on the second harvest date in plants after the flower cluster thinning, while it was higher at the last time points in treated trees. Analyzes of the P content in organs of pistachio trees with flower cluster thinning revealed its levels in fruits were higher in relation to the control, except the last time point (Figure 8e). In leaves and shoots, the P level was also higher on the first and second harvest dates in the treated pistachio trees in relation to the appropriate controls. Additionally, higher levels of this element in leaves and shoots were also observed in 25th July and June 26th, respectively. Moreover, in the fruits of pistachio trees as a result of flower cluster thinning higher the Fe contents were found (Figure 8f). Statistical analysis showed highly significant differences in these results (Table S1).

Effect of Flower Cluster Thinning Application on the Concentration of Mineral Elements in Pistachio Organs
The Cu content in fruits was higher after flower cluster thinning application as compared to the control (i.e., without cluster thinning) (Figure 8). In leaves, the Cu content was higher at the second and third period in the pistachio trees after flower cluster thinning and then decreasing for the last sampling point (Figure 8a). In shoots, the content of Cu was higher on the second harvest point in plants after the flower cluster thinning, while it was lower for the third and fourth harvest dates compared to the appropriate controls. Moreover, the Mg content in fruits was always lower after the flower cluster thinning application in relation to the control (Figure 8b). Also, in leaves of the above pistachio trees, the content of the Mg was lower at the first and second harvest dates. A similar tendency was noted in shoots at the last time point. It shoud be stressed that the flower cluster thinning leads to a strong accumulation of the Mn in fruits of pistachio trees at the last stage (in July); the concentration of the Mn was 13 times higher in these fruits than in the control (Figure 8c). In leaves and shoots of pistachio trees with flower abscission, the Mn concentration was generally lower than in the control. In fruits, the N content was lower at the first time point in pistachio trees after the flower cluster thinning, while it was higher at all other time points in the treated trees compared to the untreated ones ( Figure 8d). In leaves, the content of N was generally higher at the first, third and fourth dates in treated trees in relation to the appropriate controls. As for the shoot content, it was noticeable that the N level was lower on the second harvest date in plants after the flower cluster thinning, while it was higher at the last time points in treated trees. Analyzes of the P content in organs of pistachio trees with flower cluster thinning revealed its levels in fruits were higher in relation to the control, except the last time point (Figure 8e). In leaves and shoots, the P level was also higher on the first and second harvest dates in the treated pistachio trees in relation to the appropriate controls. Additionally, higher levels of this element in leaves and shoots were also observed in 25th July and June 26th, respectively. Moreover, in the fruits of pistachio trees as a result of flower cluster thinning higher the Fe contents were found (Figure 8f). Statistical analysis showed highly significant differences in these results (Table S1).

Correlation Coefficients between Observed Traits in Relation to Pistachio Organs
All the observed traits had normal distribution. The results of the MANOVA test indicated that all the organs (Wilk's λ = 0.0056; F = 39.37; p < 0.0001), flower (Wilk's λ = 0.5460; F = 2.65; p = 0.004) and organ × flower interaction (Wilk's λ = 0.3595; F = 2.13; p = 0.002) were significantly different regarding all the 16 quantitative traits. ANOVA indicated that the main effects of organs were significant for all the studied traits (Table S1). The main effects of flower were significant for SA, total sugar content, sucrose, Mn, and P (Table S1). The effects of organ × flower interaction were statistically significant for Cu, Mn,  (Tables S1 and S2). Statistically significant (at 0.05 level) correlation coefficients were observed for 95 from 120 pairs of the observed traits (Table S3, Figure 9).

Correlation Coefficients between Observed Traits in Relation to Pistachio Organs
All the observed traits had normal distribution. The results of the MANOVA test indicated that all the organs (Wilk's λ = 0.0056; F = 39.37; p < 0.0001), flower (Wilk's λ = 0.5460; F = 2.65; p = 0.004) and organ × flower interaction (Wilk's λ = 0.3595; F = 2.13; p = 0.002) were significantly different regarding all the 16 quantitative traits. ANOVA indicated that the main effects of organs were significant for all the studied traits (Table S1). The main effects of flower were significant for SA, total sugar content, sucrose, Mn, and P (Table S1). The effects of organ × flower interaction were statistically significant for Cu, Mn, N and Fe (Table S1, Table S2). Statistically significant (at 0.05 level) correlation coefficients were observed for 95 from 120 pairs of the observed traits (Table S3, Figure 9).

Discussion
The results of this study show for the first time selected molecular and metabolic responses of P. vera L.cv. Kirmizi during different developmental stages. Firstly, within the first two terms collected fruit from pistachio trees with the flower cluster thinning, Figure 9. Heatmaps for Pearson's correlation coefficients between the observed traits (r 0.001 = 0.23). Each cell denotes correlation coefficient between pair of observed traits. The correlation coefficients ranged from −1 (blue) to 1 (red).

Discussion
The results of this study show for the first time selected metabolic responses of P. vera L. cv. Kirmizi during different developmental stages. Firstly, within the first two terms collected fruit from pistachio trees with the flower cluster thinning, free radicals accumulation in relation to the control was demonstrated. We analyzed production of these radicals because it was important to know the responses of pistachio trees to hand-thinning of flowers during their development. Our study showed that free radicals gave signals with one g-values of 2.0045 ± 0.0005. The g-value of 2.0045 ± 0.0005, indicated these radicals are semiquinone-derived radicals [22,48]. Besides, apart from semiquinone radicals, EPR spectroscopy revealed also presence of paramagnetic Mn 2+ ions and differences in concentration of these ions in the organs of the pistachio trees with the flower cluster thinning application and the control pistachio trees. EPR spectroscopy is a method for studying materials with unpaired electrons. It is particularly useful for studying metal complexes and organic radicals. We reported here that the concentration of paramagnetic Mn 2+ ions demonstrated fluctuations. This result may indicate that Mn 2+ ions may change oxidation degree. However, Mn 3+ in contrast to Mn 2+ ions is extremely difficult to detection by EPR method applied by us (9 GHz). For Mn 3+ ions there are observable only forbidden transitions [49]. Besides, statistical analysis showed a strict correlation between observed parameters, i.e., the level of semiquinones and paramagnetic centers (Mn 2+ ) detected by EPR spectroscopy, and the concentration of Cu, Mg, Mn, and Fe measured by atomic absorption spectrometer. Increase in the level of semiquinone radicals in fruits of pistachio trees may be a consequence of a manual flower cluster thinning. Free radicals are molecular species with unpaired electrons; these radicals are highly reactive and can either extract an electron from molecules or donate an electron to other molecules thus acting as a reductant or an oxidant [50]. The higher levels of semiquinone radicals recorded in pistachio fruit trees may be related to the phenolic compound content and changes in the plant cell redox status as a consequence of manual abscission of the flowers; although phenolics were not analyzed in this study. For example, Yamasaki and Grace [51] using EPR spectroscopy detected phytophenoxyl radicals and have provided evidence for the redox coupling of plant phenolics with ascorbate in the H 2 O 2 -peroxidase system. An important implication of the results described by the above authors is that metal ions may influence the nature of plant phenolics in vivo by altering the lifetime of phenoxyl radicals as the oxidized products. The results in this study clearly indicate that zinc was a spin stabilization agent. Additionally, the above EPR studies revealed that ascorbate can act as a secondary reductant in the CGA peroxidase reaction, but not as a primary electron donor to the peroxidase. Additionally, results of our previous study and other researchers have demonstrated the dependency between to the presence of semiquinone radicals and phenolics using various experimental model systems [22,37,52]. It should be noted that research of Tomaino et al. [53] reported the occurrence of a broad spectrum of phenolic compounds in pistachio, e.g., gallic acid, catechin, eriodictyol-7-O-glucoside, naringenin-7-O-neohesperidoside, quercetin-3-O-rutinoside, eriodictyol, genistein-7-O-glucoside, genistein, daidzein, apigenin, epicatechin, quercetin, naringenin, uteolin, kaempferol, cyanidin-3-O-galactoside and cyanidin-3-O-glucoside. It has been shown that semiquinone radical formation is positively associated with markers of oxidation such as protein carbonyl and total peroxides [54]. Earlier research from Pearce et al. [23] suggested that semiquinone radicals intervene as free radical scavengers, quenching free radical chain reactions involved in the biodegradation of lignocelluloses, and thus protecting cell walls from decay. On the other hand, semiquinone radicals were reported to play an important role in the defense mechanism of plants [24,25,29,55,56]. It is believed that these stable free radicals may be associated with the deposition of lignin-like polyphenolic materials, where they may perform a key function in restricting pathogen damage and spread [57][58][59]. Thus, strong reduction of semiquinone radicals recorded by us in the shoots of pistachio trees at the second and third time point may indicate the involvement of these radicals in the lignification process. Additionally, an interesting issue to consider at this point is the participation of ascorbic acid (AsA) as a multifunc-tional metabolite with strong reducing properties that allows the neutralization of free radicals, including ROS and the reduction of molecules oxidized by ROS in cooperation with glutathione in the Foyer-Halliwell-Asada cycle. Studies have demonstrated that AsA is important for the growth and development of woody plants [60].
From a physiological point of view, it can be assumed that higher levels of semiquinone radicals in fruits of pistachio trees after flower removal may be also associated with higher metabolism of these fruits. In parallel, high concentrations of soluble sugars were also noted in those fruits being significantly higher than in the controls. This trend was clearly visible at the first and second time points of nut harvest. Therefore, probably larger amounts of sucrose and monosaccharides were redirected to these organs, intensifying their metabolism. The regulatory function of sugars affects all phases of plant life cycle interacting within phytohormones, controlling the processes of growth and development of plants [61,62]. Grossman and DeJong [63] reported that removal of flowers or subsequent thinning of fruit buds leads to even fruiting (increases its regularity) and increases the growth potential of tree fruits. Regulation of fruiting is one of the keys for obtaining fruits with appropriate and desired quality characteristics [64][65][66][67]. In the present study, the content of sucrose as well as glucose and fructose in fruits of pistachio trees following flower cluster thinning application, harvested in May and the first date in June were higher than in the control. It should be emphasized that the immense fructose peak is intriguing in fruits of pistachio trees as a result flower cluster thinning especially at the first term of fruits sampling. Fructose may act both as a signaling molecule as well as an oxidant. This high level of fructose may be related to the hydrolysis of sucrose. Thus, our results prove that both fructose as well sucrose play an important role in the early ripening of pistachio fruits. Stein and Granot [68] reported that most of the fructose found in plants, including also wood plants originates in carbon assimilated during photosynthesis. In photosynthesis, CO 2 is fixed in the chloroplasts via the Calvin cycle to yield triose phosphates (triose-P). Triose-P may then be exported to the cytosol, where two triose-P molecules are combined to create one molecule of fructose 1,6-biphosphate (F1,6BP). In turn, as reported Dennis and Blakeley [69], F1,6BP can be dephosphorylated to form F6P, which is isomerized to yield glucose 6-phosphate (G6P). Moreover, sucrose is the primary sugar transported from photosynthetic tissues through the phloem to non-photosynthetic tissues (sink tissues), where it serves as a main carbon source for metabolic pathways. In turn, Marra et al. [70] revealed that the deepest decrement of starch concentration after bud-break until bloom was in the trunk of in pistachio (P. vera L.), whereas, starch concentration started increasing after the initial growth flush.
A positive interdependence was found in this study between the fruit saccharide content and ABA concentration. Evidence was also provided here for the involvement of IAA in the fruit development of pistachio trees after flower cluster thinning. Our data have shown extremely high accumulation of IAA in fruits of pistachio trees with flower removal at the first term collection of plant materials. A similar upward trend in IAA content during fruit set was also reported by Lee et al. [71]. These latter authors also observed an increase in the activity of the soluble acid invertase which was accompanied by an increase in endogenous IAA content. Auxin affects any stage of fruit development. As reported by Nitsch [72], in order to produce flowers and fruits, plants have first to be shifted from a vegetative to a reproductive condition and this change is believed to be under hormonal control. It is well known that IAA, ABA and SA are phytohormones involved in the signalling and regulation of many crucial processes in plants, including fruit growth [73]. Endogenous concentrations of IAA are particularly high at fruit set and during initial growth developmental stages ( [73] and references cited in this work). Our analyses using GC/MS-SIM revealed that SA levels were generally lower in fruits, leaves and shoots of pistachio trees after flower cluster thinning as compared to the control organs. Earlier research results documented that some physiological events such as flower bud and embryo formation, development, bud abscission, fruit set and growth are regulated via phytohormones [74][75][76]. Our data confirmed that, mainly, IAA and ABA are involved in the development of pistachio fruits, because their levels were much higher than in the control fruits at the first harvest date of these organs.
Moreover, as reported by Zhang [77], thinning a cluster prior to nut growth will result in a higher percentage of filled nuts on the thinned cluster. Therefore, it is suspected that the intrinsic capacity of trees to store carbohydrates initially determines the percentage of filled nuts. Production of blank nuts is strongly affected by alternate bearing. Probably, the sugar status of a tree entering a crop year, sets the limits of the crop load a tree can set up and carry through maturation to splitting. Kazankaya et al. [78] revealed that kernel sugar composition of nut crops changed based on species, variety, genotype, and accession was affected by different ecological conditions. These authors also showed that the main sugar occurring in many pistachio varieties was sucrose, followed by glucose, fructose, and maltose. In turn, in walnut genotypes, the main sugar was glucose followed by sucrose, fructose and maltose. As for hazelnut, sucrose was the main sugar followed by glucose, maltose, and fructose. This observation indicates the key role of sucrose in the ripening process of pistachio fruits. As mentioned above, pistachio trees are a biennial-bearing crop. Trees that produce too much fruits due to limited carbohydrate and nutrient intake, produce small fruits with low firmness [79]. Fruit thinning contributes to an increase in the leaf to fruit ratio, which has a significant impact on increasing fruit weight, extract content as well as increase the sugar to organic acid ratio [80]. However, it should be noted that too strong thinning may lead to intensified vegetative growth, which may significantly weaken the setting of flower buds in the next year [63]. In turn, too late thinning of the fruits may lead to a weakening of fruit growth, for example its small size [81].
In the present study, we also found that among the analyzed mineral elements (Cu, P, Mg, Mn, N and Fe) in pistachio organs, Cu, P and Fe contents in fruits of pistachio trees after flower cluster thinning, were always higher than in the control fruits. Accumulation of these elements may indicate that they were redirected from other organs, e.g., leaves and shoots, to the fruits. These mineral elements play a key role in the physiological processes of plants. Mineral nutrition indeed strongly affects plant metabolism and ultimately affects plant growth, fruit development, fruit quality, nutrition, and postharvest storage [82]. Besides, mineral nutrition, such as N, P, and Fe, not only directly affects fruit quality but also influences the absorption and the content of other nutrient elements. Moreover, these trace elements mediate biochemical reactions by acting as cofactors for many enzymes, as well as acting as centers for stabilizing structures of enzymes and proteins.
Taken together, our results indicate on participation of semiquinone radicals, paramagnetic manganese ions, IAA, ABA, SA, nutrients, and sugars in P. vera organs development on the background of the flower cluster thinning which was done before fructification.

Conclusions
The obtained results for the first time revealed significant increase in generation of semiquinone radicals and fluctuations in paramagnetic metal ions (Mn 2+ ) in fruits of pistachio trees (Pistacia vera L. cv. Kirmizi) as a result flower cluster thinning during development. Additionally, our research shows the important role of signalling molecules such as phytohormones IAA and ABA in the fruit development of pistachios, especially in the early stages. As a result of flower abscission, an increase in the levels of IAA and ABA in the fruits at the first term of fruits sampling was observed. In parallel, the significantly higher level of total sugars, including sucrose and its monosaccharides in the fruits of the pistachio trees with flower cluster thinning than in the control at the first time of harvest nut draws attention. Changes in the level of mineral elements in pistachio organs indicated that these trace elements are very essential for cell functions at molecular and metabolic levels and can be transported from organ to organ depending on the demand given pistachio organ.

Supplementary Materials:
The following are available online at https://www.mdpi.com/article/ 10.3390/agronomy11112115/s1, Table S1: Mean squares from two-way analysis of variance for observed traits of pistachio, Table S2: Mean values and standard deviations of 16 observed traits, Table S3: Correlation coefficients between all pairs of observed traits.
Author Contributions: I.M. created the concept and participation in research design, directed the implementation of molecular and biochemical research, analyzed and interpreted the data, wrote and prepared the manuscript, M.Z.D., a MA student of Erasmus+ program of I.M. and N.E.K. collected plant material and took part in all molecular and biochemical analyses of plant material, A.W. took part in the preparation of samples for analysis of semiquinone radical, manganese ions and phytohormone, preparing all the results for statistical analysis, prepared the figures and the graphical abstract, was involved in compiling the entire manuscript, prepared references, W.B. contribution to the research methodology and measurements of semiquinone radicals and paramagnetic centers (Mn 2+ ) concentrations by EPR, calculation of the results from EPR, J.K. contribution to the research methodology and measurements of phytohormones, J.B. contributed to performing statistical analysis and prepared Tables S1-S3, Figure 9