HRMAS NMR Spectroscopy
Semi-solid state HRMAS NMR spectroscopy represents a very innovative and powerful technique which, despite its undisputable potential, it is still underutilized in the fields of food and agricultural chemistry. HRMAS permits to obtain NMR spectra by the direct exam of fresh and intact tissues, relatively rapidly and with a resolution very similar to that attainable via NMR in liquid-state and high resolution [24–25]. On this basis, such a technique was exploited to identify the primary metaboloma of Fiano and Pallagrello Nero grapes and research possible correlations with both soil electrical conductivity and the treatment based on the biodynamic preparation 500 biostimulant.
Initially, the most appropriate rotor spin rate was established. In fact, if, on one hand, the resolution of an NMR spectrum of a sample analysed in the semi-solid state can be improved by increasing the spin-rate, on the other hand, an excessively high rotational speed can alter the integrity of the sample, penalizing reliability, and spectral quality and reproducibility [24–25]. Therefore, through a preliminary set of tests, we found that the best compromise in consisted in a spin rate of 5 kHz.
The full 1H HRMAS NMR spectrum of a representative F berry is shown in Fig. 2a. It is evident the spectral region in which was performed the suppression of the residual solvent signal (at around 4.8 ppm) and which was digitally removed to eliminate the corresponding and non informative artifacts. The identification of the most intense proton signals found in both grapes types permitted to define their primary metaboloma [4, 22,44–46]. Figure 2b reports the assignment of the main proton signals detected, almost exclusively, in the two regions of alkyl (0.5–2.5 ppm) and hydroxyalkyl (2.5–5.5 ppm) protons. Most of predominant signals were attributed to multiplets in carbohydrates such as glucose, fructose and sucrose, followed by xylose but with a content significantly lower. Several amino acids were found in P and F berries, and consisted in proline, asparagine, GABA, glutamate, arginine, lysine, alanine, threonine, and valine. In addition, the metabolomic profile was also composed by alcohols, like methanol and ethanol, organic acids, such as tartaric, succinic and malic acids, and other compounds such as choline and lipid molecules (broadened signals in the alkyl region). Except for very weak and broadened signals, which were attributed to tryptophan and phenylalanine, there were no relevant peaks and multiplets at low field frequencies, and, in particular, in the aromatic region. This was explained by the fact that the examined tissue consisted in grape mesocarp, which notoriously lacks of phenolic and polyphenolic compounds (i.e., tannins and anthocyanins). In fact, the latter ones are expected to be much more concentrated in the grape skin and pips. The profile of the primary metabolome of F and G grapes was thus defined and, from a qualitative point of view, the same compounds were detected in both grape varieties. Subsequently, a semiquantitative comparison among the integrations of F and G spectra was conducted aiming to find significant, objective, and diagnostic characteristics related to both the soil spatial variability (based on ECa values detected via EMI technique) and the action exerted by the p500 application. It is important to underline that the studied grapes were harvested in the same vintage (September 2021), from the same vineyards and the same type of soil, managed by the same producer, and grown, practically, in the same microclimatic conditions.
Soil Spatial microvariability
In order to correlate soil vineyard microvariability with the primary metabolome of studied grapes, we elaborated the integrations of their HRMAS NMR spectra via Principal Component Analysis (PCA). The latter is an unsupervised multivariate statistical analysis enabling the rapid identification of statistical groups of observations characterized by specific peculiarities. It allows an effective exploration of very large data matrices, providing, in a single output (namely score plot), clear information on systematic variations induced by specific conditions and/or treatments, and revealing, at the same time, the variables responsible for these discriminations [43]. Therefore, HRMAS 1H NMR spectra of F samples collected from plants located in the soil sites characterized by a low (L, 9 mS/m) or high (H, 13 mS/m) ECa responses were thoroughly compared via PCA. Fig.s 3a and 3b show both the PCA score-plots and the most significant spectral comparisons between Fiano samples, for both control and p500 treated plants. The score-plot displayed in Fig. 3a refers to FC samples and results from the linear combination of the line PC1 and PC2, explaining 44.18 and 24.03% of total system variability, respectively. The fact that both L and H groups were placed in markedly different regions of the score-plot, indicated neat and peculiar differences in their primary metabolomas. Firstly, the metabolites mostly representing the PC1 and PC2 were individuated by examining the respective PCA loading-plot. Then, the metabolites responsible to significantly discriminate L from H were selected via ANOVA test. FC collected in L sites exhibited a lower content in amino acids, such as proline, glutamate, GABA, alanine, valine, and threonine, along with a lower content in ethanol and lipids (differentiation expressed along the PC1). Also the lower content in xylose permitted a collocation of L samples at higher PC2 values, thus contributing to the discrimination. On the right-side of Fig. 3a it is exhibited a comparison among the spectral regions (3 replicates per grape type) including the NMR signals mostly involved in the discrimination between the two types of grapes.
Figure 3b reports a comparison between L and H Fiano samples resulting from treatment with biodynamic p500. Again, the combination of PC1 and PC2 (respectively 40.93 and 23.93% of the total explained variability) implied a clear differentiation between FT_L and FT_H. The discrimination occurred only along the PC1 and was due to the fact that FT_H samples exhibited a higher content in fructose, lysine, lipids, and valine, accompanied by a lower amount of GABA, threonine, and alanine. Again, the spectral comparison in the right-side of the figure highlights the neat semiquantitative differences between the two types of grapes. These results proved that, in all cases, the metabolome of F samples changed significantly, as a function of soil electrical conductivity, thus proving that the factors which determined a difference of 4 mS/m between the ECa values of L and H sites are sufficient to affect the vine metabolism.
The score-plot exhibited in Fig. 3c informs on the metabolome of Pallagrello grapes (untreated with p500) as a function of different ECa responses. Only the combination PC4 and PC7 (explaining the 14,68% of total variability) allowed a discrimination between H and L, and was exclusively due to a significantly higher content of malic acid in PC_L samples.
In case of P samples treated with the biostimulant p500, it was detected a neat discrimination as a function of soil ECa responses (Fig. 3d). The discrimination occurred by combining the PC1 (43.47%) with the PC2 (27.41%) and was ascribed to the higher content in GABA, proline, asparagine, malate, glutamate, alanine, ethanol, lipids, and valine in PT_H samples. Interestingly and in agreement with the results concerning F berries treated with p500 (Fig. 3b), also in PT samples it was recorded an over-expression of several amino acids. Another common aspect is that valine and lipid compounds resulted higher than in L, for both FT_H and PT_H. Also for P grapes, the different ECa values revealed soil sites with different characteristics and affecting differently the grape metabolism.
The soil spatial variability identified by EMI technique is mainly attributable to the overall distribution in active clays, water, and organic matter which contribute to increase the apparent electrical conductivity [16]. Our results, based on sites with ECa values at least of 8 mS/m and characterized by an average difference of 4 mS/m, proved that investigated soil sites experimented specific soil conditions which, not only were indirectly revealed by EMI, but also significantly impacted on both vine physiology and grape metaboloma. In fact, in all of studied cases, it was verified that at different ECa values corresponded systematically a different metaboloma. In addition, such a discrimination resulted even more pronounced for samples subjected to biostimulant treatment. Interestingly, this outcoming wel l agrees with a recent paper focusing on Aglianicone black grapes [22] and proving the existence of a tight and systematic correlations between the composition of Aglianicone grapes, revealed via HRMAS NMR, and ECa soil data.
Effects induced by p500 treatment on grape metabolome
PCA was also exploited to investigate on possible effects induced by the action of the biostimulant p500 on the composition of the two varieties of grapes F (Fig. 4a) and P (Fig. 4b). Interestingly, both varieties of grapes showed an enhanced p500-dependent change in the primary metaboloma. The score-plot differentiating FC from FT grapes was obtained by combining PC1 and PC2 (56.67% of the total system variability explained; Fig. 4a). The discrimination occurred along the bisector direction between the second and the fourth quadrant and was determined by a higher content in proline, glutamate, GABA, alanine, ethanol, and lipids, in samples resulting from p500 treatment, accompanied by a lower content in sucrose and fructose. For P samples, the neat differentiation occurred along the PC1 (48% of the total variability of the system; Fig. 4b). Analysis of the loading plot, combined with the ANOVA test, revealed that the PC samples exhibited a significantly lower content in metabolites such as glucose, fructose, proline, glutamate, lysine, alanine, threonine, and lipids.
These outcomes suggest that, for both P and F grapes, the p500 treatment may promote the expression of several amino acids as well as can impact on carbohydrates contents. As compared to the control, the p500 induced an increase of proline, glutamate and lipids for both FT and PT. However, the influence induced by p500 on vine metabolism seems to vary as a function of the variety. A certain influence of biodynamic products on the expression of amino acids have been already observed previously. For example, Picone and coworkers [14], in agreement with our results on FT, have demonstrated that the concentration of GABA, in treated grapes, was significantly higher than the grapes subjected to traditional organic treatment, as well as it was proved that biodynamic treatments positively elicited the production of isoleucine and valine [47].
Heatmap (Fig. 5) not only confirmed the discriminability of explored samples, by clusterizing control and p500 samples into two distinct groups for both P and F samples, but also highlighted, in a direct and simple way, the semi-quantitative response (red to yellow scale) of the variables significantly involved in the grape type discrimination. In case of F samples, the map underlined a lower content in sucrose and fructose in FT samples, accompained by a larger amount of glutamate, proline, lipids, alanine ethanol and GABA (Fig. 5a). In case of Pallagrello grapes, the heat map confirmed the pronounced difference in the primary metabolome induced by p500 treatment. In fact, PT exhibited a lower abundance of fructose and glucose, while PC showed a pronounced expression of several compounds, including amino acids and lipids. Finally, it was applied the cross-validation test of PLS-DA to further prove the capacity of the used analytical approach to discriminate control from p500-treated Fiano and Pallagrello grapes in a reliable way. The test was successful since in each of five cross-validations, the "unknown" replicates were classified correctly at the 100%. In Supporting Fig. S3 they are reported the outputs of one representative cross-validation test for both Fiano (Supporting Fig. S3 a-c) and Pallagrello grapes (Supporting Fig. S3 d-f). The figure includes the PLS-DA Cross-validation score-plot, the PLS-DA Receiver Operator Characteristics (ROC) curve in which the Area Under the Curve corresponded to 1 and the PLS-DA cross-validation results. In conclusion, these outcomes showed unequivocally that biostimulant preparation 500, when applied with appropriate times, methods, and amounts, can influence the primary metabolism of vines, with effects which reflect on the berry composition and with possible implications on grape and wine quality. This finding is very important because permits to candidate the proposed HRMAS-based metabolomic approach as a reliable tool to prove objectively the biodynamic nature of the product, thus permitting to protect biodynamic producers as well as permit them to valorize and promote their products.
Titratable acidity, contents in total phenols and antioxidant agents
The titratable acidity of grape is an important parameter because is related to the pH, accounts for content in weak acids, impacts on organoleptic properties and influences the vinemaking process. It is mostly influenced by predominant organic acids which are represented by tartaric acid (0.2–1.0%, on fresh weight), followed by malic, succinic and citric acid [48]. It is reported that the content of tartaric acid, generally ranging between 3 and 9 g/l, depends on the variety, production area, vintage, and terroir [14, 48–51].
F grapes exhibited a significantly higher titratable acidity than P one (Table 1). This was expected, due to the fact that the pH of white grape berries is typically more acid than red ones. Moreover, samples FT_L, PT_H and PT_L showed the highest values for F and P variety, respectively. Importantly, all of these samples resulted from the p500 treatment, thus suggesting a possible role of this biostimulant in eliciting an accumulation of organic acids in grapes. The titratable acidity did not vary among P samples, with an average value of 3.86 g/l. Conversely, this parameter varied markedly among F samples, as a function of both p500 treatment and soil ECa values. In detail, FT_L was characterized by the highest titratable acidity, which corresponded to 5.62 g/l. Interestingly, the slightly but significantly higher titratable acidity of PC_L than PC_H is in line with the highest content in malic acid for PC_L already observed by HRMAS NMR.
The content in total phenols was evaluated in all investigated grapes through the Folin Ciocalteau assay. As expected, the black grape berries P exhibited significantly higher contents in total phenols than the white berries F, with average values of 1.455 and 0.655 mg GAE/g fw, respectively. These values are in line with the content in total phenols for red grapes, which generally ranges within 1.2 and 3.35 mg GAE/g fw, while, for white berries, the range is typically included between 0.32 and 2 mg GAE/g fw [52–55]. The highest contents in total phenols were observed for F and P grapes treated with p500, as compared to the respective control (Fig. 6). The assumption of food rich in total phenols is important due to the beneficial properties of these compounds. In fact, they may serve as powerful antioxidants, prevent inflammations, stimulate positively the nervous system, prevent cardiovascular diseases and seem to be related to anti-tumor activity [56]. On this basis, it is possible to conclude that the application of biodynamic treatment not only influenced the primary metabolome of F and P grapes, but also implied an increase in their nutraceutical and healthy value.
In addition, it was found also a neat correlation between the phenols content and site-specific ECa responses. In fact, except for FT_H and FT_L which did not vary at all among each other, FC_L was significantly higher than FC_H. Conversely, in case of Pallagrello samples, the total phenols responses changed significantly, as a function of ECa values, by exhibiting the highest contents in case of PC_H and PT_L (Fig. 6a).
The DPPH assay was carried out to quantify the free radical scavenging activity in studied grapes. It was detected a significant amount of antioxidants in both varieties (Fig. 6b) with higher average values for P grapes. This well agrees with the responses observed for the total phenols, since most of the antioxidant activity deriving from the assumption of grapes, is actually ascribable to phenolic and polyphenolic compounds. The antioxidant activity in red grape berries ranges between 2244 and 3154 µg of AEE/g fw [56–57], while the average content of antioxidant agents in white berried grapes typically is included within 1233 and 1408 µg of AEE/g fw. Again, in agreement with the total phenols, the samples FT_H and PT_L, both resulting from the treatment with the biodynamic product, gave, in absolute, the highest responses in antioxidant activity, with values of approximately 1718 and 3419 µg AAE/g fw, respectively. These results were even higher than those typical for red and white berries. About the correlation with ECa vineyard responses, in case of Pallagrello the values varied slightly, as a function of ECa values, with L samples giving the relatively higher values. In case of Fiano grapes, the changes ascribable to site-specific ECa values, were much more pronounced, with the highest values observed for FC_L and FT_H.
MRI analysis of Fiano and Pallagrello grape berries
MRI spectroscopy was exploited to examine the inner morphology of intact berries of Fiano and Pallagrello cultivars. In Fig. 7 they are shown the MRI images related to the central slices of representative Fiano (a) and Pallagrello nero (b) grape berries. MRI images clearly permitted to distinguish the three main morphological tissues composing this fruit. In fact, the thin outer darker layer is composed by the exocarp and the epidermis, whereas the mesocarp corresponds to the large area occupied by intense and whitish signals due to the high water content. The seeds are neatly identifiable in the endocarp by the presence of dark semi-ellipsoidal areas, characterized by a relatively darker and external area (the tegument) and relatively greyish and lighter area (the endosperm contained in the seeds). The latter ones appear significantly darker than the surrounding mesocarp tissue being them poorer in water and richer in oils and polyphenols. Furthermore, it was also possible to appreciate the vascular system and, in particular, the peripheral vascular bundles, branching out in proximity of the exocarp, the central fibrovascular bundles, developing below the cercine, and the bundles directly connected to the seeds (Fig. 7; supporting video 1 and 2). These images suggest that MRI can reveal anatomical details of intact in-vivo or ex-vivo samples, with a resolution of tens of microns. This is important because paves the way to its use to investigate on the early detection of treatment-related effects on inner grape morphological traits. It is important to underline that, to date, only one research group [58] has successfully attempted to use MRI to investigate on grape berries, by detecting interesting effects exerted on Nero di Troia grapes as a function of different fertilizers. Although this further proves the potential role of MRI in studying grapes, on the other hand, it certifies how the MRI potential to examine grape is still largely unexpressed and unexplored, thus requiring further validations.
Firstly, MRI images of F and P were accurately compared. They were only appreciated morphometric differences between the two cultivars, with the diameter and volume of P berries and seeds (detected in integer berry) appearing significantly larger than F grapes. However, from a qualitative and morphological point of view, no relevant and reproducible differences were appreciated as a function of both p500 treatment and soil electrical conductivity.
Therefore, the next step consisted in examining hydrogen relaxation times (T1 and T2) and self-diffusion coefficients to obtain information on grape microstructure and on behaviour, mobility and nature of the water molecules in the grapes mesocarp. Since these MRI parameters are notoriously very stable and constant, their possible and pronounced variation is diagnostic of an alteration in the tissue under examination [30]. These changes can be attributed, directly or indirectly, to factors such as: a sharp variation in the content of free water; change in microstructure and therefore in the confinement of water molecules; alteration in the consistency of the material; thickening of vascular bundles and so on [29–31].
In Fig. 8 and Sup. Table S2 are reported the results concerning T1, T2 and self-diffusion coefficients of F grapes. Differently than for MR images, MRI parametric data revealed singular and interesting responses depending on both p500 treatment and soil microvariability (Sup. Table S3). T2 values ranged between 39.56 to 53.22 ms and the differences among the F types resulted statistically significant in all cases (Fig. 8a, Sup. Table S2). Even though this is only a preliminary study, it suggests the potential use of this single parameter as an innovative and diagnostic element of quality, useful to serve as supporting data to discriminate, trace and identify grapes treated or not with preparation 500 and/or sampled from micro-areas characterized by different soil properties. Interestingly, while p500 led to the highest T2 value, in case of high ECa values, it was detected the lowest value in case of grapes treated with p500 but sampled at soil sites with the lowest electrical conductivity. For T1 relaxation times (Fig. 8b, Sup. Table S2), the values ranged between 1220 and 1489 ms. The T1 values for FT_H and FT_L did not differ from each other, while a neat difference was observed between FC_H and FC_L samples, which exhibited a gap of 269 ms, with FT_H characterized by a longitudinal relaxation time of 1489 ms. Finally, the diffusion coefficients varied from 1.213 to 1.413 mm2/s * 10− 3 (Fig. 8c, Sup. Table S2). Interestingly, the results followed a trend similar to that of T2 results with the longest relaxation time recorded for samples collected at high ECa values and resulting from p500 treatment. For samples collected at points characterized by low ECa values, no differences were found between treated (FT_L) and untreated plants (FC_L), with values averaging around 1.3 mm2/s * 10− 3.
Our results suggested that MRI proton relaxation times and diffusion coefficients are capable to appreciate grape structural properties and characteristics, which depend on conditions and factors affecting the vine phenological development during the fruit production. In fact, the application of p500 as well as different soil electrical conductivity were capable to induce pronounced changes in Fiano grapes, leading to relatively small but statistically significant inter-class differences (Sup. Table S3). In addition, the fact that all replicates (per berry and per ROI) exhibited a very low intra-class variance certified an excellent and promising reproducibility for such innovative analytical method. A variation in relaxation times and in water diffusivity may be explained by factors such as a change in the content and confinement of free water in the berries; alteration in the grape microstructure; development of a denser and thicker vascular system, with consequent impact on mesocarpic water mobility; modification in fruit consistency. The systematic change observed for these parameters permitted to conclude that they correlated with both p500 treatment and soil ECa, and in particular to conditions such as: (1) different bioaccessible water in soil, which directly influenced soil ECa values and indirectly the water content of the fruit, thus implying a higher dry weight; (2) a different and site-specific content of organic matter (excluding the role of active clayes content which were assumed constant within the examined F vineyard) (3) the impact of the biodynamic treatment on the bioavailability of water and/or nutrients, inducing an alteration that varied depending on spatial microvariability (4) the action of biodynamic treatment serving as a primer to shift and enrich the vine rhizospheric microbial community, with an influence on nutrient and water uptake as well as their respective physiology (5) the greater bioavailability of iron (presumably related to the induced presence of siderophore microorganisms) capable to imply a more efficient uptake of this ferromagnetic element, known to promote an artificial shortening (i.e. not related to the tissue itself) of relaxation times. Although, at the moment, we cannot discriminate among these factors, we can assume that most of these concurred to determine the observed results and hypothize to performe, in the future, futher experiments aimed to go deeper inside these important aspects. It is important to emphasize that the content in antioxidant activity and total phenols for FT_H sample resulted, absolutely, the highest. This outcome permits to infer that Fiano grapes treated with preparation 500 and developed in points characterized by relatively higher ECa values yielded grapes with a neatly superior nutraceutical quality. Consequently, we found that MRI T2 relaxation times and diffusion coefficients of FT_H showed a similar trend, exhibiting the highest values for both parameters. This original evidence suggests that these MRI parameters may serve as indirect markers of nutraceutical quality. Consistently, the absolutely lowest value found in T2 was that of FC_H samples, which was also the sample with the markedly poorest content in antioxidants and total phenols. Again, this finding suggests a possible correlation between nutraceutical quality of Fiano and MRI T2.
The analytical approach conducted for Fiano grapes was attempted identically also for Pallagrello grapes. This notwithstanding, it was not possible to obtain reliable and useful MRI data with an unacceptable variance. In fact, differently than for the Fiano variety, all examined Pallagrello berries tended to collapse drastically in the MRI tube before the end of all MRI experiments (each series of MRI experiments required 8 hours of continuous acquisition time) (Supporting Fig S4). At the end of the analysis we incurred in these problems: (a) a general collapse of the entire berry, whereas in MRI analysis it is required to mantain the sample motionless throughout the whole experiments [29–31] (b) a relevant amount of juice exuded from the berry on the bottom of the tube (c) a greater susceptibility to pathogens attack and (d) a drastically reduced consistency of the berry. Except for just few and not statistically relevant samples, all these conditions concurred to make Pallagrello samples (and presumably all grape varieties with similar characteristics) unsuitable for this type of MRI experiments.