Prerpints.org logo
Preprint
Article

First Attempt to Find Spectral Patterns in Cork Oak Leaves Revealing the Interaction of the Plant with Phytophthora cinnamomi

Altmetrics

Downloads

122

Views

51

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

19 July 2023

Posted:

20 July 2023

You are already at the latest version

Alerts
Abstract
Cork oak tree (Quercus suber L.) symbolises Montado landscape in Portugal and is a central element in the country's social and economic history. In recent decades, the loss of thousands of cork oaks has been reported, revealing the ongoing decline of these agro-forestry ecosystems. This emblematic tree of the Mediterranean Basin is host to the soil-born root pathogen Phytophthora cinnamomi, an active cork oak decline driver. In this framework, the early diagnosis of trees infected by the oomycete by non-invasive methods should contribute to the sustainable management of cork oak ecosystems and motivated this work. Using leaf reflectance spectroscopy it was possible to distinguish plants inoculated with P. cinnamomi from non-inoculated plants in a long course assay and to argue about the factors contributing to this differentiation. Spectral acquisitions in the visible-near infrared (400–1100 nm) were performed at 63, 78, 91, 126 and 248 days after inoculation (DAI), on the abaxial and adaxial faces of the leaves, and the spectral data was arranged in a three-way array of samples × wavelengths × days. The separation of the two plant groups were attributed to variations in water content (4v (O−H)) and/or a shift in the red edge and to structural modifications in the epidermal layer and/or the mesophyll of the leaves observed in the visible range. These spectral signatures can assist in the field identification of cork oaks that are interacting with P. cinnamomi.
Keywords: 
Subject: Biology and Life Sciences  -   Plant Sciences

1. Introduction

On the Portuguese mainland, cork oaks (Quercus suber L.) belong to the group of native evergreen woodlands known as montados, cork oak forests and holm oak forests, with a territorial occupation of 1 million hectares, representing 1/3 of the country’s forest [1]. Between 1990 and 2015, the loss of canopy cover associated with these plant species was quantified indicating forest degradation [2]. Bento et al. [3] emphasises that determining the percentage of forest cover is important for assessing the sensitivity of forests to drought highlighting the vulnerability of forests located in central, SW and NE Portugal. These reports are consistent with the deteriorating trend documented over the last few decades, threatening the survival of these highly productive agroforestry systems, which are critical to the Iberian Peninsula’s rural economy [4,5]. In the Mediterranean woodlands, the oomycete Phytophthora cinnamomi infects the roots of oaks [6] and shrubs [7,8] with the abundance of the parasite contributing to tree death along with other factors [9]. The ecological dynamics of the oomycete are adapted to variations in temperature and humidity across seasons and their abundance and spread will therefore be driven by climate change [10] and originate from recognised disease loci and through sporadic transmission of disease over longer distances [11].
The significance of remote sensing for Earth’s surface has grown considerably since the introduction of early weather satellites like Landsat and Nimbus. These satellites have revolutionized the observation of forest and crop canopies through the utilization of visible and infrared imaging, enabling the analysis of light reflection, absorption, and transmission by plant leaves. The technique has been widely used in agriculture to monitor crop growth, nutrient status, and stress responses [12]. In their comprehensive survey, Ustin and Jacquemoud [13] delve into the intricate relationship between leaf structure, biochemistry, and the optical properties of plants. Within the visible to near infrared band (400-1100 nm), the spectral reflectance characteristics primarily rely on the absorption of chlorophyll in the blue and red regions, as well as the absorption bands of water in the infrared [13,14]. Moreover, reflectance plateaus in the infrared region are associated with cellular structure, leading to higher reflectivity values compared to the lower reflectivity in the visible range caused by pigment absorption [15]. These differences in reflectance between the visible and infrared spectra serve as indicators of vegetation health and condition, particularly at the red edge transition occurring between 680 and 740 nm [16]. Indeed, many applications for crop disease detection rely on quantifying changes in these characteristics [17].
In the context of plant disease monitoring, visible-near infrared (Vis/NIR) spectroscopy serves as a valuable tool. While both remote sensing and Vis/NIR spectroscopy rely on the fundamental principle of reflectance, there are notable differences between the two approaches. Vis/NIR spectroscopy is typically performed under controlled conditions or by directly contacting the leaf, whereas remote sensing captures spectra or multi-spectral images from an airborne platform.
Thus, Vis/NIR spectroscopy acts as a testing ground for remote sensing, enabling measurements in controlled conditions and facilitating the examination of the correlation between spectral bands and vegetation traits, either structural or related with pathologies. It can also be employed as a standalone technique for proximal or contact measurements.
For instance, Slaton et al. [14] conducted a study to explore the relationship between near-infrared reflectance (NIRR) at 800 nm and leaf structural characteristics in alpine angiosperms. They discovered a strong correlation between NIRR and leaf bicoloration, leaf cuticle thickness, and the exposed mesophyll cell surface area per unit leaf surface area. This finding established a clear link between reflectance and plant structure and function. In the same line of research, Buitrago at al. [18] estimated leaf traits from infrared measurements (1.4-16.0 μm, outside the visible-near infrared range used in this work) taken from various plant species, and partial least squares regression models. Results showed that leaf thickness, water content, cellulose, lignin, and stomata density could be estimated with moderate to good accuracy.
The changes induced in the spectra by the presence of pathogens are surprisingly fast. This has been shown in a investigation of Kuska et al. [19], who used a microscope with a hyperspectral camera to analyze the interactions between Blumeria graminis f. sp. hordei and barley genotypes with different susceptibility levels. Changes in hyperspectral signatures, particularly in the green and red range, were observed within 48 hours after inoculation. Specific reflectance patterns indicated resistance, and hyperspectral imaging allowed for early detection of hypersensitive response (HR) spots.
Therefore, there exists a correlation between the emergence of infections and concurrent alterations in spectral reflectance. Numerous researchers have extensively utilized this relationship, typically employing various statistical multivariate analyses on spectral data to differentiate between healthy and infected plants or even among different types of plant diseases. The primary goal of these investigations is to enable the early detection of diseases (prior to visible symptoms) and/or their widespread occurrence through remote sensing. This approach has been applied to virtually all major disease categories [20]. For instance, it has been employed to identify infestations such as soybean cyst nematode [21], aphids in winter wheat [22], and mealybug in tomato plants [23]. Furthermore, it has been utilized to diagnose fungal diseases like sudden death syndrome [21], yellow rust and powdery mildew in winter wheat [22], target spot in tomato [24], and laurel wilt [25]. Additionally, bacterial diseases such as bacterial spot in tomato [24] and viral infections like Huanglongbing in citrus [26] have also been investigated using this approach.
Detecting Phytophthora cinnamomi early and accurately in cork oak is crucial for effective management of decline. However, traditional methods for oomycete detection involve time-consuming laboratory analyses and can be costly. Therefore, it is particularly relevant to consider the investigations conducted on Phytophthora infestans, such as those related to tomato late blight [24], for the present study. Newby et al. [27] monitored the effects of Phytophthora cinnamomi infection in Australian native vegetation through a two-factor experiment involving infection and water stress. They observed contrasting spectral responses between the inoculated and water stress groups. Hardoim et al. [28] utilized Fourier Transform Infrared Spectroscopy and found that metabolic patterns in cork oak roots were altered as early as 2 hours post-inoculation (hpi). Aubard et al. [5] assessed the decline of oak stands in Portugal and developed an effective monitoring tool for forest health and productivity using the NDVI (Normalized Difference Vegetation Index). They successfully correlated long-term trends in NDVI productivity. Another aspect explored in this study is the analysis of spectral characteristics over time following pathogen inoculation. This approach was also applied by Afonso et al. [29] in their research on citrus trees infected with Citrus Tristeza Virus. They found that discrimination between control and infected groups was not possible on certain measurement days, but became evident when analyzing the temporal trend.
The cork industry and the implementation of efficient disease management strategies would greatly benefit from the development of non-destructive and rapid methods for detecting oomycete-infected plants. In this study, we investigated the feasibility of utilizing VIS/NIR reflectance spectroscopy on cork oak leaves to identify plants that have been inoculated with P. cinnamomi and are engaged in an interaction with the oomycete. To achieve this, we employed chemometric techniques to establish classification models for the cork oak-P. cinnamomi interaction. Spectrum data were collected from leaves of both inoculated and non-inoculated cork oak plants to support our analysis.

2. Materials and Methods

The methodological design defined for this project is outlined in Figure 1. The scheme presents the different experimental approaches followed, which included the spectroscopic acquisition and the measurement and quantification of physiological and biochemical parameters in non-inoculated cork oak plants and in cork oak plants inoculated with P. cinnamomi over a period of 248 days.

2.1. Biological Material

The sample originated from twelve seedlings germinated from acorns with provenance in six cork oak trees (references S1.1, S2.1, S4.1, S5.1, S7.1, and S8.1). These parental cork oak trees were from Cachopo in the Algarve region (Portugal) and showed signs of decline. After germination, seventy-seven day-hold cork oak plants were divided into two groups of 6 plants each. One group was kept as control and the plants in the second group were inoculated with PA 45 Phytophthora cinnamomi isolate. Both groups, control (S1.1.e; S2.1.3; S4.1.3; S5.1.4; S7.1.4; S8.1.3) and inoculated (S1.1.1; S2.1.1; S4.1.2; S5.1.1; S7.1.3; S8.1.1), had one plant each germinated from seeds of the same mother tree. Since this same set of plants had already been analysed in a study dedicated to the leaf proteome, we refer the characterisation of the experimental procedure to the corresponding reference [30]. Then, the experimental procedures carried out during the inoculation of the plants with P. cinnamomi was described in detail in Coelho et al. [30], a manuscript dedicated to the study of the proteins associated with the immune response of cork oaks to this oomycete. Forty-eight hours after P. cinnamomi inoculation all plants were potted into a moisture of planting soil (PFLANZ-ERDE) and sand (proportion 2/3 for 1/3) in free-draining plastic containers (Top O 16 cm; Base O 13 cm; H 33 cm), transferred outside and watered regularly to container capacity. The assay was run for 248 days under environmental conditions.

2.2. Leaf Optical Properties Measurements

To gain insight of the changes occurring at the physiological level of cork oak plants inoculated with P. cinnamomi, the reflectance spectra of leaves from 6 control and 6 inoculated cork oak plants were captured, twice a week, for 6 months, and compared. Visible-near infrared (345 – 1037 nm) spectral data were registered through a portable setup consisting of lamp LS-1-LL (Ocean Optics, USA), spectrometer (USB4000, Ocean Optics, USA) and bifurcated fiber optic with leaf clip for reflectance measurements (PP Systems, USA). Spectral acquisitions were performed at five time points from the date of inoculation of the plants with P. cinnamomi. Measurements were taken at 63, 78, 91, 126 and 248 days after inoculation (DAI), on the abaxial and adaxial faces of the leaves (two leaves per plant), ensuring an exhaustive evaluation of the phenomenon under study.The calculations were conducted using Matlab R2019a (MathWorks, USA) and the PLS Toolbox for Matlab [31]. The Wilcoxon Rank Sum Test was applied using the ranksum function in Matlab. Additionally, the PLS Toolbox provided the necessary functions for Multiplicative Scatter Correction (mscorr), Savitzky-Golay filtering (savgol), and PARAFAC analysis (parafac).

2.3. Leaf Gas Exchange Measurements

Net photosynthetic rate (PN) (µmol CO2 m-2s-1), transpiration rate (Tr) (mmol H2O m-2s-1), stomatal conductance (gs) (mol H2O m-2s-1), intercellular CO2 concentration (Ci) (µmol CO2 mol-1) and water use efficiency (WUE) were determined for the two leaves per plant used in the spectroscopy assays taken at 63, 78, 91, 126 and 248 days after inoculation (DAI), with an infra-red gas analyser (IRGA) (LI-6400 Portable Photosynthesis System, LI-COR, USA) equipped with a 2x3 cm clear-bottom chamber (reference 6400-08, LI-COR, USA). The measurements were performed throughout the whole study period in the outside, under environmental conditions.

2.4. Quantification of Biochemical Parameters

At the end of the assay, the leaves of cork oak plants, inoculated and control, were collected and immediately frozen in liquid nitrogen and stored at -80˚C until further use for quantification of biochemical parameters. These were determined in random leaves collected from the whole cork oak plants. The quantitative determination of leaf pigments was determined by UV/VIS spectroscopy expressed in mass (mg) per unit area (m2). A circle 0.5 cm in diameter was removed from each leaf of the 12 plants in the sample, using a cutting cylinder. Pigment extractions were done in triplicate, performed on 3 leaf samples from each plant (3 biological replicates). The leaf material was submersed in 5 mL of methanol (pure solvent) for 48 hours in the dark at 4 ºC. Two aliquots were taken from each of the samples (2 technical replicates) and assayed in a spectrophotometer (model UV-160 A, Shimadzu Corporation, Japan). The content of chlorophyll-a (Chl a); chlorophyll-b (Chl b); total chlorophyll (Chla+b), and total carotenoids (tcar) were calculated according to the equations of Lichtenthaler [32,33]. Aditionally the Chl a/b ratio and Chl a + b/tcar ratio were calculated.
The quantitative determination of sugars started from samples prepared as follows: 1. random collection of 0.5 cm diameter discs from the leaves of the plants; 2. selection of 6 discs from each of the 12 leaf samples, corresponding to 100-200 mg fresh weight. Sugar determination was done in triplicate for each sample (3 technical replicates).
The whole-extract of soluble sugars was obtained by submerging the cork oak leaf discs in 1 mL ethanol 80%, followed by incubation at 80 ºC for 15 minutes.
This procedure was repeated 2 times then an aliquot of activated charcoal was added to 200 µL extract, centrifuged at 13,000 rpm for 1-2 minutes. Quantification was performed on the supernatant. The solid phase (plant leaves) was kept in the original eppendorf tube and stored at -20 ºC for later use. The determination of leaf starch content was carried out using the hydrolysis of the starch and the subsequent quantification of the obtained glucose equivalents. Cork oak leaves were ground in a mortar and pestle with quartz sand and distilled water (1 mL) and the resulting mixture was centrifuged at 15,000 rpm for 2 minutes; the supernatant was discarded and the pellet was re-suspended in water (1 mL). Prior centrifugation and resuspension procedures were repeated three times. Then, the re-suspended pellet was autoclaved for 20 min at 100 ºC. After cooling, one hundred microliters of the resulting starch suspension were submitted to enzymatic hydrolysis by the addition of 0.5 mL of amyloglucosidase (1.4 U) and α-amylase (2U) solution and incubation overnight at 37 ºC. After centrifugation of the hydrolysis solution at 15,000 rpm for 1 minute, the resulting glucose present in the supernatant was used for the quantification of the leaf starch content. Both the soluble sugars and the resulting glucose from starch cleavage samples were assayed enzymatically with the K-SUFRG–Sucrose/D-Fructose/D-Glucose Assay Kit (Megazyme International, Ireland; www.megazyme.com) at 340 nm, following the supplier’s instructions, on a Bio Tek Synergie 4 Plate Reader (BioTek Instruments, Inc., USA). Each sample was analysed in triplicate.

3. Results and Discussion

3.1. Biological Material Characterisation

As described in Coelho et al. [30], for the same group of plants under analysis, cork oak plants were in contact with the P. cinnamomi mycelium at the beginning of the experiment for 48 hours, with no (re)inoculation over the next eight months until the end of the experimental assay. Two hundred and forty-eight days after inoculation, the vegetative development of the plants looked similar in both experimental conditions, inoculated and non-inoculated. The inoculated plants were visually asymptomatic for leaf fall or yellowing, plant height or number of leaves. However, as stated in Coelho et al. [30], at the molecular level, there was a significant difference between the type and amount of proteins present in the leaves of inoculated plants when compared to control plants. Using the Arabidopsis proteome database as a reference, 424 proteins were confidently quantified in cork oak leaves, of which 80 proteins showed a p-value below 0.05 or a fold change greater than 2 or less than 0.5 in their levels between inoculated and control samples being considered as altered [30]. In the leaf proteome additional analysis, the protein profile observed in the inoculated plants expresses the defence response induced in the host by the oomycete and supports the identification of a homeostatic state in which the plant is interacting with the pathogen. It follows from these findings that the metabolism of the inoculated and non-inoculated plant groups, targeted in this new analysis, is different. The differences detected at the molecular level are a result of the effective establishment of infection in inoculated plants, although no leaf symptoms or variations in plant growth were visually observed.

3.2. Global Remarks on the Spectra Acquisition and Preprocessing

  • Signal to noise ratio
Since the aim of the spectral measurements was to detect potential variations between the control and inoculated plant groups, it is crucial to ensure that the spectra exhibit a suitable signal-to-noise ratio (SNR) across the entire range. Given that these differences might be subtle, it is essential to assess the SNR for each wavelength. SNR is computed by dividing the average signal by its standard deviation, obtained from a series of measurements. In a previous study, the same spectrometer was employed for non-invasive prediction of the internal quality of ’Jintao’ kiwifruit [29]. It was determined that the spectral ranges suitable for calibration exhibited an SNR greater than 200. Following a similar procedure, the SNR was calculated in this study, and the results are illustrated in Figure 2.
In Figure 2A, the signal-to-noise ratio (SNR) is depicted for the spectra obtained from the adaxial side of plant leaves in the control group across the five evaluation days, at 63, 78, 91, 126 and 248 days after inoculation (DAI). Each point represents the SNR value for a specific leaf at each wavelength. These individual values are then averaged and smoothed using a Savitzky-Golay filter with 200 points and a polynomial order of 2, resulting in the smooth dark curve. The yellow band indicates one standard deviation around this average. The dashed red line represents the threshold SNR value of 200, which serves as a criterion for identifying the optimal spectral range in terms of SNR. Based on this criterion, the range from 530 to 920 nm is determined to be the optimal spectral range. It is evident that beyond this range, the SNR curve rapidly declines. This decline can be attributed to two factors: (i) the spectral response of the spectrometer’s detector and (ii) the absorption caused by the pigments present in the leaves. For instance, in the range of 400 to 500 nm, both the low spectral response and the significant absorption by chlorophyll and carotenoids contribute to decreased SNR. The absorption by these pigments substantially reduces the reflected signal, leading to lower SNR values. Similarly, the drop observed above 900 nm is associated with the low quantum efficiency of silicon (which falls below 10% beyond 1000 nm) and the absorption caused by water around 960-980 nm. In Figure 2B, the average SNR curves are presented for the four datasets: control abaxial, control adaxial, inoculated abaxial, and inoculated adaxial. Interestingly, the smaller spectral features in the abaxial SNR curves closely resemble each other, as do the two adaxial SNR curves. This suggests that the fine details of the SNR are primarily determined by the spectral response of the spectrometer itself, while the larger features are influenced by the characteristics of the samples being measured. The distinction between the abaxial and adaxial leave sides becomes apparent in the 500-700 nm range, which can be attributed to the presence of pigments. On the other hand, the differences between the inoculated and control groups are not clearly evident in this overall plot combining data from all days. However, further analysis conducted separately for each day reveals discernible differences, as described below.
  • Experimental motivation for the choice of pre-processing
The spectral acquisition was conducted under sunlight and challenging conditions for the equipment. For instance, on hot days, the cooling of the light source was less effective, leading to gradual fluctuations in light intensity. Similarly, the spectrometer’s response and dark noise counts exhibited slow drifts over time. Consequently, using absolute reflectance values for comparing the two groups is not recommended. Instead, a preprocessing technique is necessary to mitigate the additive and multiplicative effects induced by environmental changes. Therefore, the spectra underwent Multiplicative Scatter Correction (MSC) as a preprocessing step. MSC was calculated from the average of the spectra from all days and all plants (i.e., including both groups). Prior to applying MSC, the spectral data was limited to the range of 530-920 nm, as determined by the SNR analysis. It is important to note that MSC causes a loss of direct interpretability for the absolute values of reflectance. However, the overall shape of the spectrum is preserved, allowing for the identification of local variations within specific spectral wavebands through changes in the spectrum shape and its relation to other bands. Figure 3 displays the superimposed spectra after applying MSC. The inoculated samples are shown in red, while the control samples are depicted in green. Figure 3A illustrates the reflectance on the abaxial side of the leaves, Figure 3B represents the reflectance on the adaxial side, and Figure 3C displays the average of the four datasets presented in the top figures.
In Figure 3C, distinct reflectance patterns are observed for the abaxial and adaxial sides. However, when considering the average spectra, it becomes challenging to differentiate between the control and inoculated groups. For the purpose of the following discussion, the spectra can be divided into three main bands. The first band corresponds to the visible range, spanning from 400 to 700 nm, which is influenced by the presence of pigments such as chlorophyll and carotenoids. The second band encompasses the near-infrared (NIR) range, ranging from 750 to 1000 nm, characterized by a plateau. The height of this plateau is primarily determined by the leaf structure, specifically the epidermal layers and the mesophyll [34]. Finally, there is a transition zone between the visible and NIR ranges, often referred to as the red edge. This region exhibits a sharp slope, connecting the low reflectance values in the visible range to higher reflectance levels in the NIR.
The observed phenomenon where the abaxial side reflects more in the visible range and slightly less in the NIR range compared to the adaxial side can be attributed to the optical characteristics of leaf tissues, which have naturally evolved to optimize their interaction with light from above. This optical engineering of the leaf tissues ensures that light is effectively channeled to the chloroplasts for optimal absorption. The palisade parenchyma, with its channellike structure, plays a key role in facilitating the conduction of light to the chloroplasts. In contrast, the spongy parenchyma, located beneath the palisade parenchyma, induces light scattering, which is particularly important for protecting the leaf from excessive infrared heating. However, when light is incident from below, the order of interaction with the two parenchyma layers is reversed. As a result, the optimal conduction of light to the chloroplasts is compromised, leading to lower absorption in the visible range and, consequently, higher reflection values. Additionally, the reflection of NIR light is also less efficient on the abaxial side, possibly due to the abaxial epidermis being less reflective than the adaxial epidermis.
  • Example of the analysis on one of the days
Although the overall comparison of spectra collected across all measurement days does not reveal distinct differences between the control and inoculated groups, conducting a specific analysis for each measurement day uncovers noticeable distinctions. These differences have been verified using the Wilcoxon Rank Sum Test (also known as the Mann Whitney U-test), which assesses whether the data in the control and inoculated groups originate from continuous distributions with equal medians. The test is performed for each wavelength individually. The application of the classical t-test across the entire spectral range was not feasible due to a significant portion of wavelengths not meeting the assumptions of normality and homogeneity of variances. Surprisingly, even with these violated assumptions, running the t-test produced results that were highly similar to those obtained with the nonparametric test. In Figure 4, the results for the second measurement day (78 DAI) are presented. The top plots, A and B, depict the MSC-transformed spectra for the abaxial and adaxial sides, respectively. The figure insets provide a closer view of a restricted range within the NIR plateau, specifically between 820 and 920 nm, where the differences between the control and inoculated groups are prominently visible. The bottom plots, C and D, illustrate the p-values obtained from the Wilcoxon Rank Sum Test for the data in plots A and B, respectively, as a function of wavelength. The p-values indicate the statistical significance of the differences observed between the two groups. The horizontal lines represent the threshold values of p = 0.05 (95% confidence level) and p = 0.01 (99% confidence level). It is important to note that the scale of the p-value axis is logarithmic. Regions of the curves below the p = 0.05 threshold are highlighted in red, indicating the wavelengths where the medians of the MSC reflectance for the control and inoculated groups significantly differ.
The differences observed in the spectra can be interpreted as follows:
1. The visible band (up to 700 nm), which is primarily influenced by pigments, does not exhibit significant differences between the control and inoculated groups. This suggests that there were no significant variations in the content of pigments between the two groups.
2. The regions highlighted in red, which correspond to the following spectral bands previously shown in Figure 3, are where notable differences between the control and inoculated groups can be observed:
(a) Around 750 nm (only on the abaxial side), coinciding with the presence of 4v (O − H), which represents the third overtone of the O–H stretching vibration.
(b) Around 830 nm, close to 3v + δ (O − H), which corresponds to the O–H combination band of the second overtone of stretching with the fundamental bending vibration.
(c) Around 920 nm, aligning with the presence of 4v (C −H), which represents the third overtone of the C–H stretching vibration.
These specific spectral bands indicate distinct molecular vibrations or combinations that may vary between the control and inoculated groups, suggesting differences in the chemical composition or molecular structures of the samples. The first two bands are related with water, while the third may reflect changes in, basically, any organic molecule(s).
However, there is an alternative interpretation that is both more plausible and straightforward. The broad absorption peak of 3v (O−H) at 960 nm is outside of the considered spectral range, but its effects may still be felt above 900 nm. As a result, the observed differences above 900 nm could also be attributed to variations in water content. Hence, adhering to the principle of Occam’s razor, the simplest and most likely explanation is that the spectral differences are indeed influenced by changes in water content.
  • Temporal analysis of the differences between the two groups
The analysis presented in Figure 4 can be extended to all five days of measurement, as demonstrated in Figure 5. In this figure, the p-values resulting from the Wilcoxon Rank Sum Test, which compares the MSC reflectances of the control and inoculated groups for each measurement day, are depicted. Each line in the figure corresponds to one of the measurement days (63, 78, 91, 126, and 248 DAI), while the columns represent the abaxial and adaxial sides, respectively. The p-values are plotted as a function of wavelength. The interpretation of the horizontal lines and colour codes remains consistent with Figure 4.
On the first day of measurement (63 DAI), there is a noticeable difference primarily around 750 nm. This difference can be attributed to variations in water content (4v (O−H)) and/or a shift in the red edge. The red edge is a sensitive region characterized by a steep slope, marking the transition from low reflectance in the visible range to high reflectance in the NIR plateau. Any alteration in pigments or the leaf structure that affects NIR scattering can lead to changes in the slope of the red edge, making it particularly useful for detecting NIR variations. Examining the p-values curve, we observe an additional region marked in red on the abaxial side, above 900 nm. This can also be influenced by changes in the off-range water absorption peak at 960 nm. In summary, the most plausible explanation for the differences observed on the first day is as follows: there are initial changes occurring in water levels (indicated by the variations above 900 nm), which slightly impact the NIR plateau. This change is further amplified in the slope of the red edge, making it appear as the primary source of differentiation between the control and inoculated groups.
The patterns observed on days 78, 91, and 126 DAI exhibit similarities that have been previously examined in the example of the analysis on one of the days. These patterns align with variations in water content, reflecting distinctions in the 3v+δ (O−H) combination band and the O–H combination band at 840 nm and 960 nm, respectively. Additionally, the 4v (C − H) band within the 910-930 nm range may contribute to these observations, although distinguishing the individual contributions of these bands can be challenging.
By the fourth day (126 DAI), differences attributable to water remain, but the distinctions across the near-infrared (NIR) plateau become more pronounced. With the exception of a band centred around 870 nm, all regions within the NIR plateau exhibit variations between healthy and inoculated plants. This suggests a significant disruption in the physical structure of the leaves. As previously mentioned, the palisade and spongy parenchyma tissues play a crucial role in light scattering within the NIR plateau. The observed differences across the entire spectral range of the plateau may indicate alterations in the structure of the parenchyma tissues, likely stemming from disparities in the initial water content. The final day of observation reveals a reversal of the observed patterns. On the abaxial side of the leaves, no notable differences are present except for occasional spikes of noise. On the adaxial side, most of the differences within the NIR plateau diminish, with the exception of the band ranging from 750 nm to 820 nm. However, a new development arises in the range of pigments, specifically between 620 nm and 700 nm, where differences between the two groups are apparent. The similarity observed in the NIR plateau between the inoculated and control plants could stem from two possible causes:
  • the inoculated plants have regained their initial structural integrity, or
  • the control plants have also experienced structural change.
The latter scenario is plausible considering that a summer season elapsed between the fourth and final measurements. On the other hand, the differences observed in the pigment range may also be attributed to two potential factors:
  • variations in pigment content between the two groups, or
  • structural changes in leaf anatomy altering the light path to and from the chloroplasts, thereby impacting reflectance levels.
In the latter hypothesis, the pigment content may be similar in both groups, but the reflectance levels differ. An important clue for discerning which scenario occurred lies in the absence of visible range differences on the abaxial side. This suggests that the pigment content remains unchanged, as any divergence would affect abaxial reflectance. Therefore, the most plausible explanation for the observations on the final day points to structural changes. The leaf structure of the control plants underwent alterations that made their NIR response resemble that of the infected plants. However, these same structural modifications resulted in visible reflectance differences between the two groups. How is this possible?
To comprehend why structural changes can lead to distinct alterations in visible and NIR reflectance patterns, Figure 6 provides valuable insights. In this figure, it is evident that visible light primarily interacts with the upper portion of the leaf, including the chloroplasts and extending down to the palisade mesophyll. Conversely, NIR light primarily interacts with the spongy mesophyll, where it is scattered. Based on the observations from the final day, the following conclusions can be drawn:
  • The pigments’ content is similar in both the control and inoculated groups.
  • The spongy mesophyll of the control group underwent alterations during the summer and acquired a structure similar to that of the inoculated group.
  • The upper epidermal layer and/or the palisade mesophyll experienced differential changes between the control and inoculated groups, resulting in distinct responses in the visible range.
These findings suggest that while the pigments’ content remains consistent, structural modifications in the spongy mesophyll account for the similarity in NIR reflectance between the two groups. Conversely, differential changes in the upper epidermal layer and/or the palisade mesophyll contribute to the divergent responses observed in the visible range. Indeed, the explanation provided aligns with the observed similarity in abaxial reflectances as well. In this scenario, where reflected light predominantly interacts with the spongy mesophyll, the comparable structure of this mesophyll layer in both the control and inoculated groups explains the absence of discernible differences in abaxial reflectances.
  • Vegetation Indices
Vegetation Indices (VI) are commonly employed to assess vegetation health through remote sensing. Therefore, it is important to explore whether VI can effectively detect distinctions between the two groups. The indices considered in this study are presented in Table 1.
It is crucial to note that the vegetation indices (VI) were derived from the raw reflectance data rather than the MSC-transformed data. Calculating the VIs directly from the raw data ensures that they are not distorted. This approach poses no issue since any potential spectrum fluctuations, which necessitated the application of MSC as a compensation transformation, are naturally accounted for by the nature of the VIs themselves. VIs are primarily defined by ratios, enabling them to compensate for such fluctuations. Relatively to the contents of Table 1, Normalized Difference Vegetation Index (NDVI) is the most widely adopted VI and with the greenness (health) of vegetation [35]. NDVI values range from 1 to -1. Healthy NDVI values are greater than 0.3 and lower NDVI values represent barren or drought-tolerant plants. The Green Normalized Difference Vegetation Index (GNDVI) method is used for estimating photosynthetic activity and is a commonly used vegetation index to determine water and nitrogen uptake into the plant canopy [36]. The Water Index (WI) is sensitive to changes in canopy water status. As the water content of vegetation canopies increases, the strength of the absorption around 970 nm increases relative to that of 900 nm [37]. We have included a similar index, calculated from the reflectance around 950 nm instead. This is called here Water Band Index (WBI). The Soil-adjusted vegetation index (SAVI) accounts for the differential red and near-infrared extinction through the vegetation canopy. The index is a transformation technique that minimizes soil brightness influences from spectral vegetation indices involving red and NIR wavelengths [38]. The indices O1 – O4 are simple ratio and/or differences between two bands, less used in the literature and that do not have a uniformly adopted nomenclature. O1 is a simple version of the Weighted Difference Vegetation Index (WDVI) [39] with the a coefficient taken as 1. Since it measures the difference in height between the red and NIR plateaus, and the later depends on the scattering by the leaf structure, it may be related with changes in leaf structure. O2 is the simple red/NIR ratio [40], while O4 is simply its inverse and usually called Simple Ratio NIR/RED Difference Vegetation Index (DVI). Finally, O3 is the green/red ratio [41]. The Photochemical Reflectance Index (PRI) is sensitive to changes in carotenoid pigments (e.g. xanthophyll pigments) in live foliage [42]. The Rededge Vegetation Stress Index (RVSI) identifies spectral changes in the red edge, which may also be traced to changes in leaf structure [43]. The Modified Chlorophyll Absorption in Reflectance Index (MCARI) minimizes the effects of nonphotosynthetic materials on spectral estimates of absorbed photosynthetically active radiation [44]. The Visible Atmospherically Resistant Index (VARI) is designed to emphasize vegetation in the visible portion of the spectrum, while mitigating illumination differences and atmospheric effects [45]. The Browning Reflectance Index (BRI) was for estimation of superficial scald development in anthocyanin-free fruit upon a background of Chl and Car absorption [46], but we have included it in the attempt of finding browning in the leafs.
Table 1. Vegetation Indices (VI) used in this work.
Table 1. Vegetation Indices (VI) used in this work.
Name Formula* Reference % Detection
NDVI (N - R)/(N + R) [35] 10
GNDVI (G - R)/(G + R) [36] 20
WBI R950/R900 [37] 80
SAVI [(1 + L)(N - R)]/(N + R + L), L = 0.5 [38] 40
O1 N - aR, a = 1 [39] 70
O2 R/N [40] 10
O3 G/R [41] 20
O4 N/R [47] 10
PRI (R531 - R570)/(R531 + R570) [42] 20
RSVI (R714 + R752)/2 - R733 [43] 40
MCARI [(R700 - R670)-0.2(R700 - R500)] X R700/R670 [44] 10
VARI (G - R)/(G + R - B) [45] 30
WI R900/R970 [37] 100
BRI (1/R550 - 1/R700)/N [46] 30
*The symbols R, G, B, and N correspond to the red, green, blue, and near-infrared bands, respectively, and were calculated based on the average values of the respective wavelength ranges (600-700 nm for red, 500-600 nm for green, 400-500 nm for blue, and 800-900 nm for near-infrared). Additionally, Rxyz represents the monochromatic reflection at a wavelength of xyz nm. The explanations of the acronyms in the first column can be found in the main text. The references providing a more detailed introduction or explanation of the VIs are listed in the third column. Lastly, the fourth column presents the success rate in detecting differences among the groups, expressed as a percentage. A 100% success rate corresponds to 10 successful detections (5 days × 2 sides, i.e., abaxial and adaxial), which can be easily understood by referring to Figure 8.
Figure 7 illustrates the results of the Vegetation Indices (VI) obtained from the adaxial side on the fourth day of measurements, corresponding to 126 DAI. The figure also displays the outcomes of the Wilcoxon Rank Sum Test comparing the two groups. Among the 14 VIs examined, eight of them demonstrate the ability to differentiate between the two groups. Notably, the fourth day exhibits the highest number of significant differences, which aligns with the observations depicted in Figure 5, where the fourth day exhibits more bands shown in red. Figure 8 is analogous to Figure 5, but it focuses on the VIs. Each row in the figure corresponds to one of the measurement days (63, 78, 91, 126, and 248 DAI), while the columns represent the abaxial and adaxial sides. The p-values are plotted against the VIs, which are shown along the x-axis of the bottom row. The interpretation of the horizontal lines and color codes remains consistent with Figure 5, where p-values below 0.05 are indicated by a red point encircled by a circumference.
Preprints 80004 g007
Figure 7. Vegetation Indices (VI) acquired from the adaxial side on the fourth day of measurements, or 126 DAI. Each subplot represents a specific VI, where the control values are represented in green and the inoculated values in red. The name of the parameter is displayed at the top of each subplot, and below it is the p-value resulting from the Wilcoxon Rank Sum Test, which compares the medians of the two groups. Significant p-values are highlighted in red. The significance levels are denoted by a star code: (*) indicates p < 0.05 (95% confidence level), (**) indicates p < 0.01 (99% C.L.), and (***) indicates p < 0.001 (99.9% C.L.). The boxes and error bars are defined by five horizontal segments, which, from top to bottom, represent the mean plus standard deviation, the third quartile, the median, the first quartile, and the mean minus standard deviation, respectively.
Figure 7. Vegetation Indices (VI) acquired from the adaxial side on the fourth day of measurements, or 126 DAI. Each subplot represents a specific VI, where the control values are represented in green and the inoculated values in red. The name of the parameter is displayed at the top of each subplot, and below it is the p-value resulting from the Wilcoxon Rank Sum Test, which compares the medians of the two groups. Significant p-values are highlighted in red. The significance levels are denoted by a star code: (*) indicates p < 0.05 (95% confidence level), (**) indicates p < 0.01 (99% C.L.), and (***) indicates p < 0.001 (99.9% C.L.). The boxes and error bars are defined by five horizontal segments, which, from top to bottom, represent the mean plus standard deviation, the third quartile, the median, the first quartile, and the mean minus standard deviation, respectively.
Preprints 80004 g007
The distribution of the red points in Figure 8 does not exhibit an obvious pattern. It does not align precisely with the distribution of the red and blue components in Figure 5. This discrepancy arises because the vegetation indices (VIs) convey nonlinear information that may not directly mirror the spectral information. However, upon a 8 more detailed analysis, the two figures do agree in the conclusions. As indicated in Table 1, the VIs that are most successful in discriminating between the two groups are WI (100%) and WBI (80%), both of which are associated with water content. The third most successful parameter is O1 or WDVI, which relates to changes in leaf structure. The success rates of all the other VIs are below 50%. Notably, SAVI and RSVI achieve a 40% success rate. While SAVI is a modified version of NDVI, RSVI targets the red edge, which, as mentioned earlier, may also be influenced by structural changes. In conclusion, the analysis of the VIs confirms the findings of the spectral analysis, highlighting that water content and leaf structural changes appear to be the primary distinguishing factors between the two groups.
Preprints 80004 g008
Figure 8. The p-values obtained by the Wilcoxon Rank Sum Test, on the comparison between the Vegetation Indices (VI) of control and inoculated groups from the data obtained in each measurement day. Each row corresponds to one of these days (63, 78, 91, 126 and 248 DAI), the left and right column correspond to the abaxial and adaxial sides, respectively. The p-values are shown as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL), at p = 0.01 (99% CL) and at p = 0.001 (99.9% CL). Note the scale is logarithmic. The points below the threshold p = 0.05 are represented in encircled red points, meaning that control and inoculated groups have different medians for that specific VI.
Figure 8. The p-values obtained by the Wilcoxon Rank Sum Test, on the comparison between the Vegetation Indices (VI) of control and inoculated groups from the data obtained in each measurement day. Each row corresponds to one of these days (63, 78, 91, 126 and 248 DAI), the left and right column correspond to the abaxial and adaxial sides, respectively. The p-values are shown as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL), at p = 0.01 (99% CL) and at p = 0.001 (99.9% CL). Note the scale is logarithmic. The points below the threshold p = 0.05 are represented in encircled red points, meaning that control and inoculated groups have different medians for that specific VI.
Preprints 80004 g008
  • PARAFAC analysis
Before delving into this section, it is important to mention that a Principal Component Analysis (PCA) was initially conducted on the spectroscopic results for each day. The objective was to identify two distinct clusters in the scores plot, which would effectively separate the two groups. Additionally, the loadings plot was examined to identify the most significant bands contributing to group separation. While this approach yielded successful results for some days, it was not consistently effective across all days. Furthermore, the information obtained from the PCA analysis did not provide substantial clarity or offer substantially distinct insights compared to those derived from Figure 5. Consequently, an alternative feature extraction technique was employed. Exploiting the fact that the data is organized in the form of a three-dimensional array (rows = samples) × (columns = wavelengths) × (sheets = days), a method called PARAFAC (Parallel Factor Analysis) was adopted [48]. PARAFAC facilitates a trilinear decomposition directly on the three dimensions of the data, providing a comprehensive overview of the entire process across the five days, rather than isolated analyses for each individual day. Furthermore, each dimension undergoes compression, similar to PCA. For instance, in the wavelength dimension, the more than 2000 variables were condensed to just three factors, effectively explaining the majority of data variability. Ultimately, each dimension is characterized by a reduced number of factors through the application of the PARAFAC technique. This reduction in dimensionality allows for a more concise representation of the data while still capturing the essential information and variability present in each dimension. After applying a Savitzky-Golay filter to the MSC-transformed data, the following results were obtained. The filter had a width of 155 points, approximately equivalent to 30 nm, and a polynomial order of 2 with a derivative order of 1. This filtering technique emphasized the shape of the spectra, enabling improved discrimination between the groups. However, as a trade-off, the clarity regarding the specific wavelengths’ contributions was reduced. The adoption of a large filter width was necessary to effectively remove noise from the data. The model presented in this study was ultimately reduced to two factors. The primary criterion for determining the number of factors was the core consistency check, which was performed using the loadings obtained from the PARAFAC model. When three factors were considered, the core consistency check did not pass, leading to the decision to retain only two components. While including more factors could potentially account for additional variability in the data, it would primarily capture noise. The first two components explain 44% of the total data variability on the adaxial side and 43% on the abaxial side. Although there is still unexplained variability, these two components are sufficient to capture the differences between the two groups.
Figure 9 displays the PARAFAC results obtained from the adaxial side data and Figure 10 shows the corresponding results for the abaxial side.
In the upper right panel, the scores of the samples are plotted in the first two components, where control leaves are represented by green markers and inoculated leaves by red markers. The upper left panel shows the wavelength loadings in the first component, while the lower left panel illustrates the wavelength loadings in the second component. Finally, the lower right panel exhibits the loadings associated with the third dimension (time) in the first two components. The upper right plots in these figures reveal the most notable observation. A distinct separation between the two groups is evident along the second factor. Generally, the control samples exhibit values of F2 < 0, while the inoculated samples show values of F2 > 0, a trend that holds true for both adaxial and abaxial sides. Additionally, an interesting observation is that the first factor, F1, appears to capture variations within the same tree, while the second factor, F2, captures variations between the trees. This distinction arises because the segments connecting each pair of leaves from the same tree align predominantly with the F1 axis, indicating that tree-related variations are primarily associated with variations in F1. This pattern is more pronounced in the control leaves compared to the inoculated leaves, particularly on the adaxial side. The reduced parallelism observed in the inoculated group suggests internal changes resulting from the infection. The area plots (upper right and lower left) illustrate the loadings across the wavelength dimension. F1, which corresponds to intra-tree variations, exhibits significant contributions primarily below 750 nm. On the other hand, F2, associated with inter-tree variations, demonstrates substantial contributions primarily above 750 nm. Consequently, intra-variations are predominantly linked to changes in pigments, while inter-variations (which enable group differentiation) are primarily associated with water content and leaf structural changes. These findings align with the previous conclusions drawn from the direct spectral analysis as well as the examination of vegetation indices, reinforcing the overall consistency of the results. The final subplot in the lower left show cases the loadings along the time dimension, which corresponds to the order of the measurement days. Its interpretation follows the previously assigned meanings of F1 and F2. For instance, the first and last days exhibit the greatest separation in terms of the first factor, indicating that the leaves from these days are the most dissimilar in terms of pigment content. This can be attributed to the natural aging process of the leaves. Conversely, in terms of water and structural changes (F2), the first and fifth days appear to converge towards a similar state by the end of the experiment. This observation is consistent with the spectral analysis in Figure 5, where the first and fifth rows exhibit relatively similar patterns, with fewer red bands compared to the other three days. However, caution is required when interpreting other observations due to the derivative spectra used in the analysis performed by PARAFAC. For instance, in both Figure 9 and Figure 10, it is the third day that appears to be the furthest from the rest in terms of F2, suggesting it is the day with the most significant differences between the two groups. On the contrary, Figure 5 and Figure 8 indicate that the fourth day exhibits more discriminative bands. The resolution to this apparent contradiction lies in recognizing that PARAFAC results are based on derivative data, and it is plausible that the rate of changes is higher on the third day, leading to a more pronounced separation along the F2 axis.

3.3. Leaf Gas Exchange Measurements

Globally, both plant groups exhibited similar values for the net photosynthetic rate (PN), transpiration rate (Tr), stomatal conductance (gs), intercellular CO2 concentration (Ci) and water use efficiency (WUE) (Table 2). Indeed, no statistically significant differences (p<0.05) were found between the control and the inoculated plants at any of the sampling dates.

3.4. Variation Pattern of Physiological Parameters and Biochemical Parameters

The pigment content [chlorophyll-a (chl a); chlorophyll-b (chl b); total carotenoids (tcar)] was calculated according to Lichtenthaler (1987)[32,33] and the results obtained on average for the control and for the inoculated plants are shown in Table 3.
In our sample, chlorophyll a content was lower in inoculated plants than in control plants with no statistically significant differences. For the other pigments also no values were found that differed significantly between the inoculated and non-inoculated samples.
Although there were no statistically significant differences between the two plant samples, a behaviour trend was observed for chlorophylls similar to that observed by Daoudi et al. [49] in trials with cork oaks from humid and semi-arid regions subjected to water stress.
Soluble sugar concentration expressed as mg per cm2 of leaf area was slightly higher in inoculated plants for fructose and glucose, with no statistically significant differences between control and inoculated plants (Table 4).
The plants inoculated with P. cinnamomi did not accumulate soluble sugars in values significantly different from control plants. The results obtained for the leaf proteome in this same assay stated that “the dynamics of proteins associated with sugar metabolism and sugar signaling reveals a metabolic tendency in favour of energy production and reducing power as opposed to the accumulation of sucrose and carbohydrates as reserve substances” [30].

4. Conclusions

During a plant-pathogen interaction, the homeostatic state of the participants varies over time. The immune response of plants to the oomycete differs from plant to plant and the plasticity of spectroscopic profiles in inoculated plants should exceed that of control plants.
The provided description of the observed spectral changes between 63 and 255 DAI is consistent with the measurements. Here is a summary of the key points:
1. At 63 DAI, slight changes in the mesophyll, primarily influenced by water content, are observed. These subtle differences are predominantly detected at the red edge around 750 nm.
2. By 78 and 91 DAI, clear differentiation between the control and inoculated groups becomes apparent, particularly concerning the water bands. The variations in water content contribute significantly to the distinguishable patterns observed.
3. At 126 DAI, the impact of water content changes extends throughout the entire mesophyll structure. This is evidenced by a broad range in the NIR region where reflectance differences are observed.
4. At 248 DAI, following the summer season, the pigments’ content remains similar between the control and inoculated groups. However, distinct alterations in the leaf structure become evident, leading to the following observations:
(a) The spongy mesophyll of the control group undergoes changes that render it structurally similar to that of the inoculated group.
(b) Differential changes occur in the upper epidermal layer and/or the palisade mesophyll between the control and inoculated groups, resulting in differentiated responses in the visible range.
5. Several vegetation indices exhibit the capability to detect the infection, particularly those related to water content such as WI (100% success rate) and WBI (80% success rate). Additionally, O1 (or WDVI) achieves a 70% success rate, which may indicate sensitivity to structural changes in the leaves caused by the infection. Interestingly, the widely used NDVI parameter proves to be almost ineffective in detecting the infection, with only a 10% success rate. This emphasizes the importance of considering specific vegetation indices tailored to the target application, as different indices capture different aspects of plant health and physiological changes.
6. The PARAFAC analysis on derivative spectra provided further confirmation of the conclusions derived from the direct spectral analysis and the vegetation indices. It revealed two factors, F1 and F2, which represent different directions of variation in the data. F1 captures the natural physiological variations within the leaves and is primarily influenced by changes in pigment bands in the visible spectrum. On the other hand, F2 is associated with the differences between the control and inoculated plants and is predominantly determined by variations in water content and leaf structural changes, as reflected in the NIR plateau region of the spectra. These findings align with the previous observations and highlight the distinct contributions of pigments and water/structural changes in discriminating between the control and inoculated plants.
Overall, this description aligns well with the measured data, providing a comprehensive understanding of the observed changes over the course of time.
The results of spectroscopy confirm the data from the analysis of the leaf proteome carried out on the same plants [30]. It is confirmed that it is possible to differentiate cork oaks inoculated by P. cinnamomi from non-inoculated ones, even when the plants are asymptomatic, by the analysis of proteins present in the leaves and by spectroscopic methods.

Author Contributions

Conceptualization, A.C.C. and R.G.; methodology, A.C.C., R.G., A.M.C.; formal analysis, A.C.C., R.G., A.M.C., G.S.; investigation, A.C.C., R.P., R.G., A.M.C., A.B.; resources, A.C.C., R.G.; writing—original draft preparation, R.G., A.C.C., A.M.C., A.B.; writing—review and editing, A.C.C., R.G., A.M.C., G.S., A.B., R.P.; visualization, R.G., A.C.C.; supervision, R.G., A.C.C.; project administration, A.C.C.; funding acquisition, R.G., A.C.C., G.S., A.M.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Fundação para a Ciência e a Tecnologia (FCT), grant numbers UID/Multi/00631/2013, UID/Multi/00631/2019, UIDB/00631/2020 CEOT BASE and UIDB/00631/2020 CEOT PROGRAMÁTICO.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 6.º Inventário Florestal Nacional – IFN6, Relatório Final. ICNF (Instituto da Conservação da Natureza e das Florestas): Lisboa, Portugal, 2015.
  2. Acácio, V.; Dias, F.S.; Catry, F.X.; Bugalho, N.; Moreira, F. Canopy Cover Loss of Mediterranean Oak Woodlands: Long-term Effects of Management and Climate. Ecosystems 2021, 24, 1775–1791. [Google Scholar] [CrossRef]
  3. Bento, V.A.; Russo, A.; Vieira, I.; Gouveia, C.M. Identification of forest vulnerability to droughts in the Iberian Peninsula. Theor. Appl. Climatol. 2023, 152, 559–579. [Google Scholar] [CrossRef]
  4. Duque-Lazo, J.; Navarro-Cerrillo, R.M.; van Gils, H.; Groen, T.A. Forecasting oak decline caused by Phytophthora cinnamomi in Andalusia: Identification of priority areas for intervention. For. Ecol. Manag. 2018, 417, 122–136. [Google Scholar] [CrossRef]
  5. Aubard, V.; Paulo, J.A.; Silva, J.M.N. Long-Term Monitoring of Cork and Holm Oak Stands Productivity in Portugal with Landsat Imagery. Remote Sens. 2019, 11, 525. [Google Scholar] [CrossRef]
  6. Serrano, M.S.; Rios, P.; Gonzalez, M.; Sánchez, M.E. Experimental minimum threshold for Phytophthora cinnamomi root disease expression on Quercus suber. Phytopathol Mediterr. 2015, 54(3), 461–464. [Google Scholar] [CrossRef]
  7. Moreira, A.C.; Martins, J.M.S. Influence of site factors on the impact of Phytophthora cinnamomi in cork oak stands in Portugal. For. Pathol. 2005, 35(3), 145–162. [Google Scholar] [CrossRef]
  8. Gómez-Aparicio, L.; Ibáñez, B.; Serrano, M.S.; De Vita, P.; Ávila, J.M.; Pérez-Ramos, I.M.; García, L.V.; Sánchez, M.E.; Marañón, T. Spatial patterns of soil pathogens in declining Mediterranean forests: implications for tree species regeneration. New Phytol. 2012, 194, 1014–1024. [Google Scholar] [CrossRef]
  9. Matías, L.; Abdelaziz, M.; Godoy, O.; Gómez-Aparicio, L. Disentangling the climatic and biotic factors driving changes in the dynamics of Quercus suber populations across the species’ latitudinal range. Divers. Distrib. 2019, 25, 524–535. [Google Scholar] [CrossRef]
  10. Serrano, M.; Romero, M.A.; Homet, P.; Gómez-Aparicio, L. Climate change impact on the population dynamics of exotic pathogens: The case of the worldwide pathogen Phytophthora cinnamomi. Agric. For. Meteorol. 2022, 322, 109002. [Google Scholar] [CrossRef]
  11. Cardillo, E.; Abad, E.; Meyer, S. Iberian oak decline caused by Phytophthora cinnamomi: A spatiotemporal analysis incorporating the effect of host heterogeneities at landscape scale. For. Pathol. 2021, 51, e12667. [Google Scholar] [CrossRef]
  12. Bauriegel, E.; Giebel, A.; Geyer, M.; Schmidt, U.; Herppich, W.B. Early detection of Fusarium infection in wheat using hyper-spectral imaging. Comput. Electron. Agric. 2011, 75(2), 304–312. [Google Scholar] [CrossRef]
  13. Ustin, S. L.; Jacquemoud, S. How the optical properties of leaves modify the absorption and scattering of energy and enhance leaf functionality. In Remote sensing of plant biodiversity; Publisher: Springer International Publishing, 2020; pp. 349–384. [Google Scholar] [CrossRef]
  14. Slaton, M.R.; Raymond Hunt Jr, E.; Smith, W.K. Estimating near-infrared leaf reflectance from leaf structural characteristics. Am J Bot. 2001, 88(2), 278–84. [Google Scholar] [CrossRef]
  15. Gausman, H.W.; Allen, W.A. Optical parameters of leaves of 30 plant species. Plant Physiol., 1993, 52, 57–62. [Google Scholar] [CrossRef]
  16. Dawson, T.P.; Curran, P.J. Technical note: A new technique for interpolating the reflectance red edge position. Int. J. Remote Sens. 1998, 19(11), 2133–2139. [Google Scholar] [CrossRef]
  17. Franke, J.; Menz, G. Multi-temporal wheat disease detection by multi-spectral remote sensing. Precis. Agric., 2007, 8, 161–172. [Google Scholar] [CrossRef]
  18. Buitrago, M.F.; Groen, T.A.; Hecker, C.A.; Skidmore, A.K. Spectroscopic determination of leaf traits using infrared spectra. Int J Appl Earth Obs Geoinformation 2018, 69, 237–250. [Google Scholar] [CrossRef]
  19. Kuska, M.T.; Brugger, A.; Thomas, S.; Wahabzada, M.; Kersting, K.; Oerke, E.-C.; Steiner, U.; Mahlein, A.-K. Spectral Patterns Reveal Early Resistance Reactions of Barley Against Blumeria graminis f. sp. hordei. Phytopathology 2017, 107, 1388–1398. [Google Scholar] [CrossRef] [PubMed]
  20. Zahir, S.A.D.M.; Omar, A.F.; Jamlos, M.F.; Azmi, M.A.M.; Muncan, J. A review of visible and near-infrared (Vis-NIR) spectroscopy application in plant stress detection. Sens. Actuator A Phys. 2022, 338(1), 113468. [Google Scholar] [CrossRef]
  21. Bajwa, S.G.; Rupe, J.C.; Mason, J. Soybean Disease Monitoring with Leaf Reflectance. Remote Sens. 2017, 9, 127. [Google Scholar] [CrossRef]
  22. Shi, Y.; Huang, W.; Luo, J.; Huang, L.; Zhou, X. Detection and discrimination of pests and diseases in winter wheat based on spectral indices and kernel discriminant analysis. Comput. Electron. Agric. 2017, 141, 171–180. [Google Scholar] [CrossRef]
  23. Canário, D.; Figueiredo, E.; Franco, J.C.; Guerra, R.M.F.D.N. Detecting early mealybug infestation stages on tomato plants using optical spectroscopy. Eur. J. Hortic. Sci., 2017, 82(3), 341–348. [Google Scholar] [CrossRef]
  24. Lu, J.; Ehsani, R.; Shi, Y.; Castro, A.I.; Wang, S. Detection of multi-tomato leaf diseases (late blight, target and bacterial spots) in diferent stages by using a spectral-based sensor. Sci. Rep. 2018, 8, 2793. [Google Scholar] [CrossRef] [PubMed]
  25. Abdulridha, J.; Ehsani, R.; Castro, A. Detection and Differentiation between Laurel Wilt Disease, Phytophthora Disease, and Salinity Damage Using a Hyperspectral Sensing Technique. Agriculture 2016, 6, 56. [Google Scholar] [CrossRef]
  26. Sankaran, S.; Mishra, A.; Maja, J.M.; Ehsani, R. Visible-near infrared spectroscopy for detection of Huanglongbing in citrus orchards. Comput. Electron. Agric., 2011, 77(2), 127–134. [Google Scholar] [CrossRef]
  27. Newby, Z.; Murphy, R.J.; Guest, D.I.; Ramp, D.; Liew, E.C.Y. Detecting symptoms of Phytophthora cinnamomi infection in Australian native vegetation using reflectance spectrometry: complex effects of water stress and species susceptibility. Australas. 2019, 48, 409–424. [Google Scholar] [CrossRef]
  28. Hardoim, P.R.; Guerra, R.; Rosa da Costa, A.M.; Serrano, M.S.; Sánchez, M.E.; Coelho, A. C. Temporal metabolic profiling of the Quercus suber–Phytophthora cinnamomi system by middle-infrared spectroscopy. For. Path. 2016, 46(2), 122–133. [Google Scholar] [CrossRef]
  29. Afonso, A.M.; Guerra, R.; Cavaco, A.M.; Pinto, P.; Andrade, A.; Duarte, A.; Power, D.M.; Marques, N.T. Identification of asymptomatic plants infected with Citrus tristeza virus from a time series of leaf spectral characteristics. Comput. Electron. Agric., 2017, 141, 340–350. [Google Scholar] [CrossRef]
  30. Coelho, A.C.; Pires, R.; Schütz, G.; Santa, C.; Manadas, B.; Pinto, P. Disclosing proteins in the leaves of cork oak plants associated with the immune response to Phytophthora cinnamomi inoculation in the roots: A long-term proteomics approach. PLoS ONE 2021, 16(1), e0245148. [Google Scholar] [CrossRef]
  31. PLS Toolbox 9.0. Eigenvector Research, Inc., Manson, WA USA 98831; 2021, software available at http://www.eigenvector.com.
  32. Lichtenthaler, H. K. Chlorophylls and carotenoids: Pigments of photosynthetic biomembranes. Meth. Enzymol. 1987, 148(C), 350–382. [Google Scholar] [CrossRef]
  33. Lichtenthaler, H.K.; Buschmann, C. Chlorophylls and carotenoids: Measurement and characterization by UV-VIS spectroscopy. Curr. protoc. food anal. chem. 2001, Supplement 1, 1–8 F4.3.1-F4. [Google Scholar] [CrossRef]
  34. Karabourniotis, G.; Liakopoulos, G.; Bresta, P.; Nikolopoulos, D. The Optical Properties of Leaf Structural Elements and Their Contribution to Photosynthetic Performance and Photoprotection. Plants. 2021, 10(7), 1455. [Google Scholar] [CrossRef] [PubMed]
  35. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. Third ERTS-1 Symposium NASA, NASA SP-351, Washington DC, 309-317, 1974.
  36. Gitelson, A.A.; Kaufman, Y.J.; Merzlyak, M.N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sens Environ. 1996, 58(3), 289–298. [Google Scholar] [CrossRef]
  37. Peñuelas, J.; Fillela, I.; Biel, C.; Serrano, L.; Savé, R. ”The Reflectance at the 950-970 Region as an Indicator of Plant Water Status.”. Int. J. Remote Sens. 1993, 14, 1887–1905. [Google Scholar] [CrossRef]
  38. Huete, A.R. A soil-adjusted vegetation index (SAVI). Remote Sens Environ. 1998, 25(3), 295–309. [Google Scholar] [CrossRef]
  39. Baret, F.; Guyot, G. Potentials and limits of vegetation indices for LAI and APAR assessment. Remote Sens Environ. 1991, 35(2-3), 161–173. [Google Scholar] [CrossRef]
  40. Ceccato, P.; Flasse, S.; Tarantola, S.; Jacquemoud, S.; Grégoire, J.M. Detecting vegetation leaf water content using reflectance in the optical domain. Remote Sens Environ. 2001, 77(1), 22–33. [Google Scholar] [CrossRef]
  41. Gamon, J.A.; Surfus, J.S. Assessing leaf pigment content and activity with a reflectometer. New Phytologist, 1999, 143, 105–117. [Google Scholar] [CrossRef]
  42. Gamon, J.; Peñuelas, J.; Field, C. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens Environ. 1992, 41, 35–44. [Google Scholar] [CrossRef]
  43. Merton, R. Monitoring community hysteresis using spectral shift analysis and the red-edge vegetation stress index. In Proceedings of the Seventh Annual JPL Airborne Earth Science Workshop (pp. 12-16). JPL Pasadena, CA, USA, 1998.
  44. Daughtry, C.S.; Walthall, C.L.; Kim, M.S.; De Colstoun, E.B.; McMurtrey III, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens Environ. 2000, 74(2), 229–239. [Google Scholar] [CrossRef]
  45. Gitelson, A.A.; Kaufman, Y.J.; Merzlyak, M.N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sens Environ 1996, 58(3), 289–298. [Google Scholar] [CrossRef]
  46. Merzlyak, M.N.; Gitelson, A.A.; Chivkunova, O.B.; Solovchenko, A.E.; Pogosyan, S.I. Application of reflectance spectroscopy for analysis of higher plant pigments. Russ. J. Plant Physiol 2003, 50, 704–710. [Google Scholar] [CrossRef]
  47. Zheng, G.; Moskal, L.M. (2009). Retrieving leaf area index (LAI) using remote sensing: theories, methods and sensors. Sensors 2009, 9(4), 2719–2745. [Google Scholar] [CrossRef] [PubMed]
  48. Bro, R. PARAFAC. Tutorial and applications. Chemometr Intell Lab Syst. 1997, 38(2), 149–171. [Google Scholar] [CrossRef]
  49. Daoudi, H.; Derridj, A.; Hannachi, L.; Mévy, J.P. Comparative drought responses of quercus suber seedlings of three algerian provenances under greenhouse conditions. Revue d’Ecologie (Terre et Vie) 2018, 73(1), 57–70. [Google Scholar] [CrossRef]
Preprints 80004 g001
Figure 1. Representative scheme of the project’s methodological design. (C: six non-inoculated cork oak plants, control plants; I: cork oak plants inoculated with P. cinnamomi; chl-a: chlorophyll a; chl-b: chlorophyll b; tcar: total carotenoids; gs: Stomatal conductance; Ci: intercellular CO2 concentration; Tr: Transpiration rate; WUE: Water use efficiency; Tair: Air temperature; PARo: external Photosynthetically Active Radiation; PN: Net photosynthetic rate).
Figure 1. Representative scheme of the project’s methodological design. (C: six non-inoculated cork oak plants, control plants; I: cork oak plants inoculated with P. cinnamomi; chl-a: chlorophyll a; chl-b: chlorophyll b; tcar: total carotenoids; gs: Stomatal conductance; Ci: intercellular CO2 concentration; Tr: Transpiration rate; WUE: Water use efficiency; Tair: Air temperature; PARo: external Photosynthetically Active Radiation; PN: Net photosynthetic rate).
Preprints 80004 g001
Preprints 80004 g002
Figure 2. A: Signal-to-noise ratio calculated from adaxial side of the control group. The points represent the value of SNR for each leaf, the dark curve is an average over the all points and the yellow band represents one standard deviation around the mean. The dashed red line shows the threshold SNR > 200. B: Average SNR curves for the four datasets: control and inoculated, either abaxial or adaxial leaf sides.
Figure 2. A: Signal-to-noise ratio calculated from adaxial side of the control group. The points represent the value of SNR for each leaf, the dark curve is an average over the all points and the yellow band represents one standard deviation around the mean. The dashed red line shows the threshold SNR > 200. B: Average SNR curves for the four datasets: control and inoculated, either abaxial or adaxial leaf sides.
Preprints 80004 g002
Preprints 80004 g003
Figure 3. A: Superposition of all MSC treated spectra, abaxial side. The inoculated samples are represented in red and the control samples in green. B: Same as A, but for the adaxial side. C: Average of the four datasets represented in A and B. The nominal positions of the main absorption bands are marked with dashed lines.
Figure 3. A: Superposition of all MSC treated spectra, abaxial side. The inoculated samples are represented in red and the control samples in green. B: Same as A, but for the adaxial side. C: Average of the four datasets represented in A and B. The nominal positions of the main absorption bands are marked with dashed lines.
Preprints 80004 g003
Preprints 80004 g004
Figure 4. Spectra for the second day of measurements, or 78 DAI. A and B: MSC-transformed spectra for the abaxial and adaxial sides, respectively. Figure insets: zoom on the 820 – 920 nm range. C and D, show the p-values for the Wilcoxon Rank Sum Test on the data of A and B, respectively, as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL) and at p = 0.01 (99% CL). Note the scale is logarithmic. The curve values below the threshold p = 0.05 are represented in red, meaning that control and inoculated groups have different MSC reflectance medians.
Figure 4. Spectra for the second day of measurements, or 78 DAI. A and B: MSC-transformed spectra for the abaxial and adaxial sides, respectively. Figure insets: zoom on the 820 – 920 nm range. C and D, show the p-values for the Wilcoxon Rank Sum Test on the data of A and B, respectively, as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL) and at p = 0.01 (99% CL). Note the scale is logarithmic. The curve values below the threshold p = 0.05 are represented in red, meaning that control and inoculated groups have different MSC reflectance medians.
Preprints 80004 g004
Preprints 80004 g005
Figure 5. The p-values obtained by the Wilcoxon Rank Sum Test, on the comparison between the MSC reflectances of control and inoculated groups from the data obtained in each measurement day. Each line corresponds to one of these days (63, 78, 91, 126 and 248 DAI), the left and right column correspond to the abaxial and adaxial sides, respectively. The p-values are shown as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL) and at p = 0.01 (99% CL). Note the scale is logarithmic. The curve values below the threshold p = 0.05 are represented in red, meaning that control and inoculated groups have different MSC reflectance medians.
Figure 5. The p-values obtained by the Wilcoxon Rank Sum Test, on the comparison between the MSC reflectances of control and inoculated groups from the data obtained in each measurement day. Each line corresponds to one of these days (63, 78, 91, 126 and 248 DAI), the left and right column correspond to the abaxial and adaxial sides, respectively. The p-values are shown as a function of the wavelength. The horizontal lines are at p = 0.05 (95% CL) and at p = 0.01 (99% CL). Note the scale is logarithmic. The curve values below the threshold p = 0.05 are represented in red, meaning that control and inoculated groups have different MSC reflectance medians.
Preprints 80004 g005
Preprints 80004 g006
Figure 6. Interaction of light with a plant’s leaf. When light, particularly in the visible spectrum, enters the leaf from the upper side (adaxial side), it undergoes significant absorption by chloroplasts in the palisade mesophyll. However, a portion of this visible light remains unabsorbed and is transmitted to the spongy mesophyll, where additional chloroplasts contribute to its absorption. The presence of air channels within the spongy mesophyll plays a crucial role by inducing strong light scattering. As a result, some of the transmitted photons are redirected back towards the palisade mesophyll, where most of them are absorbed, while a small fraction escapes as reflected light. Notably, the absorption is more pronounced for blue and red light. On the other hand, near infrared (NIR) light interacts weakly with the palisade mesophyll and is not absorbed. Consequently, the majority of NIR photons reach the spongy mesophyll, where they are scattered back towards the palisade mesophyll and largely transmitted back into the atmosphere. This process leads to a significant net reflectance. The existence of the NIR plateau, therefore, depends on the leaf’s structural characteristics.
Figure 6. Interaction of light with a plant’s leaf. When light, particularly in the visible spectrum, enters the leaf from the upper side (adaxial side), it undergoes significant absorption by chloroplasts in the palisade mesophyll. However, a portion of this visible light remains unabsorbed and is transmitted to the spongy mesophyll, where additional chloroplasts contribute to its absorption. The presence of air channels within the spongy mesophyll plays a crucial role by inducing strong light scattering. As a result, some of the transmitted photons are redirected back towards the palisade mesophyll, where most of them are absorbed, while a small fraction escapes as reflected light. Notably, the absorption is more pronounced for blue and red light. On the other hand, near infrared (NIR) light interacts weakly with the palisade mesophyll and is not absorbed. Consequently, the majority of NIR photons reach the spongy mesophyll, where they are scattered back towards the palisade mesophyll and largely transmitted back into the atmosphere. This process leads to a significant net reflectance. The existence of the NIR plateau, therefore, depends on the leaf’s structural characteristics.
Preprints 80004 g006
Preprints 80004 g009
Figure 9. PARAFAC results obtained from the adaxial side data. Upper Right: plot of the samples’ scores in the first two components. Control leaves are depicted in green, while inoculated leaves are shown in red. The lines join the leaves from the same tree. Upper left: wavelength loadings in the first component. Lower left: wavelength loadings in the second component. Lower right: plot of the loadings associated with the third dimension (time) in the first two components.
Figure 9. PARAFAC results obtained from the adaxial side data. Upper Right: plot of the samples’ scores in the first two components. Control leaves are depicted in green, while inoculated leaves are shown in red. The lines join the leaves from the same tree. Upper left: wavelength loadings in the first component. Lower left: wavelength loadings in the second component. Lower right: plot of the loadings associated with the third dimension (time) in the first two components.
Preprints 80004 g009
Preprints 80004 g010
Figure 10. PARAFAC results obtained from the abaxial side data. Upper Right: plot of the samples’ scores in the first two components. Control leaves are depicted in green, while inoculated leaves are shown in red. The lines join the leaves from the same tree. Upper left: wavelength loadings in the first component. Lower left: wavelength loadings in the second component. Lower right: plot of the loadings associated with the third dimension (time) in the first two components.
Figure 10. PARAFAC results obtained from the abaxial side data. Upper Right: plot of the samples’ scores in the first two components. Control leaves are depicted in green, while inoculated leaves are shown in red. The lines join the leaves from the same tree. Upper left: wavelength loadings in the first component. Lower left: wavelength loadings in the second component. Lower right: plot of the loadings associated with the third dimension (time) in the first two components.
Preprints 80004 g010
Table 2. Leaf gas exchange in control and inoculated plants determined under environmental conditions. The temporal global mean (5 samplings) and standard deviation of 12 independent replicates (2 leaves per plant; 6 control/inoculated plants) are shown on the table.
Table 2. Leaf gas exchange in control and inoculated plants determined under environmental conditions. The temporal global mean (5 samplings) and standard deviation of 12 independent replicates (2 leaves per plant; 6 control/inoculated plants) are shown on the table.
Physiological parameters* Control Samples (N=6x2x5) Inoculated Samples (N=6x2x5)
Mean Std. Deviation Mean Std. Deviation
PN (µmol CO2 m-2s-1) 7.76 0.65 9.41 0.47
gs (µmol m-2s-1) 0.1008 0.0094 0.0895 0.0081
Ci (µmol CO2 mol-1) 240.59 14.62 227.24 14.05
Tr (mmol H2O m-2s-1) 2.8692 0.3566 2.7638 0.3180
WUE (µmol CO2/mmol H2O) 3.62 0.30 3.63 0.27
Environmental conditions** Control Samples Inoculated Samples
Tair (ºC) 27.83 0.37 30.11 0.34
PARo (µE m-2s-1) 1485.917 117.3678 1512.050 117.2944
RH (%) 35.24 1.78 30.39 1.66
VPDL (kPa) 2.60 0.13 3.19 0.14
*Net photosynthetic rate (PN); transpiration rate (Tr); stomatal conductance (gs); intercellular CO2 concentration (Ci) and water use efficiency (WUE). **Tair: Air temperature; PARo: external Photosynthetically Active Radiation; RH: Air relative humidity; VPDL: leaf-to-air water vapor pressure deficit.
Table 3. Results obtained on average for chlorophyll-a (Chl a), chlorophyll-b (Chl b), and total carotenoids (tcar) in the inoculated and control plants. The mean and standard deviation are shown as well as pigment ratios or combinations.
Table 3. Results obtained on average for chlorophyll-a (Chl a), chlorophyll-b (Chl b), and total carotenoids (tcar) in the inoculated and control plants. The mean and standard deviation are shown as well as pigment ratios or combinations.
Pigments Control samples Inoculated samples
Mean Std. Deviation Mean Std. Deviation
Chl a (mg/cm2) 0.0454 0.0131 0.0435 0.0047
Chl b (mg/cm2) 0.0082 0.0033 0.0091 0.0014
tcar (mg/cm2) 0.0120 0.0033 0.0122 0.0018
Pigment ratios/combinations Control samples Inoculated samples
Chl a+b (mg/cm2) 0.0536 0.0526
Chl a/b 5.5 4.8
Chl a+b/tcar 4.5 4.3
Table 4. Results obtained on average for glucose, fructose, sacarose and starch in the inoculated and control plants. The mean and standard deviation are shown.
Table 4. Results obtained on average for glucose, fructose, sacarose and starch in the inoculated and control plants. The mean and standard deviation are shown.
Sugars Control samples Inoculated samples
Mean Std. Deviation Mean Std. Deviation
Glucose (mg/cm2) 0.0265 0.0080 0.0299 0.0152
Fructose (mg/cm2) 0.0135 0.0074 0.0206 0.0191
Sacarose (mg/cm2) 0.3619 0.1010 0.3571 0.0623
Starch (mg/cm2) 0.1198 0.0558 0.1341 0.0862
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.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated