The Effect of Biotic Stress in Plant Species Induced by ‘ Candidatus Phytoplasma solani’—An Artificial Neural Network Approach

: Infections with phytoplasma present one of the most significant biotic stresses influencing plant health, growth, and production. The phytoplasma ‘ Candidatus Phytoplasma solani’ infects a variety of plant species. This pathogen impacts the physiological and morphological characteristics of plants causing stunting, yellowing, leaf curling, and other symptoms that can lead to significant economic losses. The aim of this study was to determine biochemical changes in peony ( Paeonia tenuifolia L.), mint ( Mentha × piperita L.), and dill ( Anethum graveolens L.) induced by ‘ Ca . Phytoplasma solani’ in Serbia as well as to predict the impact of the biotic stress using artificial neural network (ANN) modeling. The phylogenetic position of the Serbian ‘ Ca . Phytoplasma solani’ strains originated from the tested hosts using 16S rRNA (peony and carrot strains) and plsC (mint and dill strains) sequences indicated by their genetic homogeneity despite the host of origin. Biochemical parameters significantly differed in asymptomatic and symptomatic plants, except for total antho-cyanidins contents in dill and the capacity of peony and mint extracts to neutralize superoxide anions and hydroxyl radicals, respectively. Principal Component Analysis (PCA) showed a correlation between different chemical parameters and revealed a clear separation among the samples. Based on the ANN performance, the optimal number of hidden neurons for the calculation of TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, and Car was nine (using MLP 8-9-13), as it produced high r 2 values (1.000 during the training period) and low SOS values. Developing an effective early warning system for the detection of plant diseases in different plant species is critical for improving crop yield and quality.


Introduction
Plants undergo a wide range of biotic and abiotic stresses that limit their proper growth [1].Biotic stress in plants is caused by living organisms, including fungi, bacteria, phytoplasma, viruses, nematodes, insects, and weeds, causing diseases or damage [2].Specific ion channels and kinase cascades and reactive oxygen species (ROS), as well as phytohormones (e.g., abscisic acid, salicylic acid, jasmonic acid, and ethylene), are activated in plants following exposure to stress [3].Phytoplasmas are plant pathogenic bacteria belonging to the class Mollicutes.They lack cell walls and exclusively inhabit the nutrient-rich phloem tissue of plants [4].Diseases caused by phytoplasmas pose a threat to agricultural production by leading to substantial losses in yield and the quality of plants.More than 600 diverse diseases have been reported worldwide, mainly on vegetable crops belonging to the families Apiaceae, Asteraceae, Cucurbitaceae, Fabaceae, and Solanaceae [5][6][7][8].The symptoms of diseases caused by phytoplasmas include yellowing or reddening of leaves, phyllody, stunting, virescence, shortened internodes, big bud, little leaf, witches' broom or proliferation, leaf curl, giant calyx, floral malformation, and vascular discoloration [7,9].Phytoplasmas are transmitted from plant to plant by phloemfeeding insect vectors, and they propagate within the cytoplasm of both insects and plants.Phytoplasma attack causes various changes in terms of physiological, biochemical, and metabolic pathways [4,5,7,10,11].For example, phytoplasma infections have been known to decrease chlorophyll synthesis [5,12] and protein content or to increase peroxidase activity [11].Our previous studies [5,7] indicated biochemical changes in common eveningprimrose (Oenothera biennis L.) and carrot (Daucus carota L.) resulting from infection by 'Candidatus Phytoplasma solani'.Changes in O. biennis resulted in a significant increase in the peroxidation of lipids, phenylalanine ammonia-lyase activity, total sugars, polyphenols, and anthocyanins content, along with a decrease in photosynthetic pigments and total flavonoids [5].The oxidative damage of membranes in carrot cells was accompanied by a decrease in the content of photosynthetic pigments, as well as a pronounced reduction in the level of glutathione (GSH) content.Conceivably, anthocyanidins were responsible for the enhanced antioxidative capacity [7].
Artificial neural networks (ANNs) have found applications in diverse domains for forecasting outcomes of intricate systems by processing input data.ANNs belong to a class of machine learning algorithms capable of learning from data and predicting results with precision and dependability, such as identifying plant diseases.ANNs have exhibited remarkable potential in predicting biotic stress in plants resulting from diverse pathogens, enabling real-time plant health monitoring and offering a pre-emptive alert mechanism for farmers, agronomists, and researchers to forestall disease outbreaks.
The present study aimed to determine biochemical changes induced by 'Ca.Phytoplasma solani' in the three latest host plants identified in Serbia, including peony (Paeonia tenuifolia L.), mint (Mentha × piperita L.), and dill (Anethum graveolens L.), and to utilize these findings to construct an artificial neural network (ANN) model in order to predict the impact of the biotic stress in these particular plant species.

Plant Material
Experiments performed in this study were conceived for three plant species, i.e., P. tenuifolia, Mentha × piperita, and A. graveolens with externally visible disease symptoms caused by 'Ca.Phytoplasma solani' (Figure 1) [13][14][15].Plants of each species with symptoms (S-symptomatic) and without symptoms (A-asymptomatic) were collected from Bački Petrovac (45 • 21 ′ 38 ′′ N 19 • 35 ′ 30 ′′ E) in the Bačka region of Vojvodina, Serbia, at the "full-bloom" phenophase.The study field belongs to the Institute of Field and Vegetable Crops, National Institute of the Republic of Serbia (Novi Sad, Serbia), for the purpose of reproduction and maintenance of seed collection.
A minimum of thirty plants of each tested species (per fifteen randomly selected asymptomatic and symptomatic plants) were collected at the growing sites.One part of the collected leaves (two per taken plant) was immersed in liquid nitrogen and further stored under −20 • C, while the other part was dried at ambient temperature in a well-ventilated place.Previously published results of biochemical changes in D. carota plants infected with 'Ca.Phytoplasma solani' (Figure 1) in Serbia were used for the comparison in all assays [7].plants infected with 'Ca.Phytoplasma solani' (Figure 1) in Serbia were used for the comparison in all assays [7].

Phylogenetic Position of the Tested Serbian 'Ca. Phytoplasma solani' Strains
Phylogenetic analysis was performed to check the position of the five Serbian 'Ca.Phytoplasma solani' strains originating from the same A. graveolens, P. tenuifolia, Mentha × piperita, and D. carota plants as those used for biochemical analysis, in relation to Phytoplasma spp.strains from different hosts and countries.Two separate phylogenetic trees were constructed based on the previously published sequences of the Serbian strains of interest.The first Neighbor-Joining (NJ) phylogenetic tree was constructed based on the 16S rRNA sequences for the three Serbian 'Ca.Phytoplasma solani' strains originating from P. tenuifolia (B15 and B18) [14] and D. carota (ML_NS-2016) [7], as well as with 23 comparative Phytoplasma spp.strains belonging to different 16Sr groups (Table 1) retrieved from the GenBank database.The tree was rooted with Acholeplasma laidlawii strain NCTC10116 (Table 1).

Phylogenetic Position of the Tested Serbian 'Ca. Phytoplasma solani' Strains
Phylogenetic analysis was performed to check the position of the five Serbian 'Ca.Phytoplasma solani' strains originating from the same A. graveolens, P. tenuifolia, Mentha × piperita, and D. carota plants as those used for biochemical analysis, in relation to Phytoplasma spp.strains from different hosts and countries.Two separate phylogenetic trees were constructed based on the previously published sequences of the Serbian strains of interest.The first Neighbor-Joining (NJ) phylogenetic tree was constructed based on the 16S rRNA sequences for the three Serbian 'Ca.Phytoplasma solani' strains originating from P. tenuifolia (B15 and B18) [14] and D. carota (ML_NS-2016) [7], as well as with 23 comparative Phytoplasma spp.strains belonging to different 16Sr groups (Table 1) retrieved from the GenBank database.The tree was rooted with Acholeplasma laidlawii strain NCTC10116 (Table 1).The second NJ phylogenetic tree was constructed based on the plsC gene (1-acylsn-glycerol-3-phosphate acyltransferase) sequences with the remaining two Serbian 'Ca.Phytoplasma solani' strains of interest, originating from A. graveolens (STOL2) [13] and Mentha × piperita (STOL3) [15], and 13 comparative Phytoplasma spp.strains from different hosts and countries (Table 2), which were also from GenBank.This tree was rooted with Acholeplasma laidlawii strain DSM 23060 (Table 2).Before the construction of each phylogenetic tree, sequences of all strains were aligned to the same size (710 nt for 16S rRNA; 399 nt for plsC) in the BioEdit v.7.0 program using the ClustalW Multiple alignment function.Trees were constructed in Mega7 software, using the bootstrap value of 1000.The Kimura two-parameter nucleotide substitution model was used to compute genetic distances [14].

Principal Component Analysis (PCA)
Principal Component Analysis (PCA) is a statistical method commonly used in exploratory data analysis [27] to identify patterns and relationships within a dataset.The transformation is based on the eigenvalue decomposition of a data correlation matrix [28], which ensures that the first principal component captures the highest possible variance in the data.This technique can reveal spatial relationships among processing parameters and facilitate the identification of underlying patterns within complex datasets.

Artificial Neural Network (ANN) Modeling
During the training cycle for artificial neural network (ANN) modeling, all input and output data were normalized to improve the network's behavior.The optimal number of hidden layers and neurons in each layer was determined through trial and error to achieve good performance.A multi-layer perceptron (MLP) model with three layers (input, hidden, and output) was chosen for this study, as it is known to be effective at approximating nonlinear functions [29].In this study, the optimal number of hidden neurons was found to be ten.The ANN scheme is presented in Figure 2. The Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm was used for ANN modeling, and the training process was repeated multiple times to achieve the best performance considering parameter variability.Successful training was determined when the learning and cross-validation curves (Sum of Squares vs. training cycles) approached zero.The ANN was tested using the best weights obtained during the training step, and the coefficient of determination (r 2 ) and Sum of Squares (SOS) were used to evaluate the performance (i.e., accuracy) of the ANN.Artificial neural networks show high performance in predicting TS, RG, PAL, LP, TP, TT, Tflav, Tpro, Tant, NBT, • OH, DPPH, Tcha, Tchb, and Car by calculating weight coefficients and biases with respect to characteristic input parameters (logical variables including plant type and treatment) [30,31].The basic equation (Equation ( 1)) for calculating the output data of artificial neural networks is as follows: where Y represents the output value, f 1 and f 2 represent the transfer function in the hidden and output layer, and X represents the matrix of the input layer [30].
Horticulturae 2024, 10, x FOR PEER REVIEW 6 of 15 training process was repeated multiple times to achieve the best performance considering parameter variability.Successful training was determined when the learning and cross-validation curves (Sum of Squares vs. training cycles) approached zero.The ANN was tested using the best weights obtained during the training step, and the coefficient of determination (r 2 ) and Sum of Squares (SOS) were used to evaluate the performance (i.e., accuracy) of the ANN.Artificial neural networks show high performance in predicting TS, RG, PAL, LP, TP, TT, Tflav, Tpro, Tant, NBT, •OH, DPPH, Tcha, Tchb, and Car by calculating weight coefficients and biases with respect to characteristic input parameters (logical variables including plant type and treatment) [30,31].The basic equation (Equation ( 1)) for calculating the output data of artificial neural networks is as follows: where Y represents the output value, f1 and f2 represent the transfer function in the hidden and output layer, and X represents the matrix of the input layer [30].

Global Sensitivity Analysis
Yoon's global sensitivity equation was used to calculate the relative impact of the input parameters on the output variables, according to the weight coefficients of the developed ANN models [32,33]: where w-weight coefficient in the ANN model, i-input variable, j-output variable, k-hidden neuron, n-number of hidden neurons, and m-number of inputs.

Global Sensitivity Analysis
Yoon's global sensitivity equation was used to calculate the relative impact of the input parameters on the output variables, according to the weight coefficients of the developed ANN models [32,33]: where w-weight coefficient in the ANN model, i-input variable, j-output variable, k-hidden neuron, n-number of hidden neurons, and m-number of inputs.

Phylogenetic Position of 'Ca. Phytoplasma solani' Strains
The NJ phylogenetic tree based on the 16S rRNA sequences of three tested Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia (B15 and B18) and D. carota (ML_NS-2016) and the comparative Phytoplasma spp.strains from different hosts/countries is presented in Figure 3. Based on the constructed tree, all tested/comparative 'Ca.Phytoplasma solani' strains were placed in the same tree cluster (the XII-A 16Sr Group) showing genetic homogeneity among themselves regardless of the host/country of origin.The remaining Phytoplasma spp.strains belonging to 16Sr Groups I (A, B, and D), II (A and D), VIII-A, X-A, XI-A, and XII-B were clearly separated in other tree branches/clusters, each corresponding to a specific 16Sr Group.Acholeplasma laidlawii strain NCTC10116 was placed on a monophyletic tree branch, as an outgroup (Figure 3).tested/comparative 'Ca.Phytoplasma solani' strains were placed in the same tree cluster (the XII-A 16Sr Group) showing genetic homogeneity among themselves regardless of the host/country of origin.The remaining Phytoplasma spp.strains belonging to 16Sr Groups I (A, B, and D), II (A and D), VIII-A, X-A, XI-A, and XII-B were clearly separated in other tree branches/clusters, each corresponding to a specific 16Sr Group.Acholeplasma laidlawii strain NCTC10116 was placed on a monophyletic tree branch, as an outgroup (Figure 3).The NJ phylogenetic tree that was constructed based on the plsC gene sequences for the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and comparative Phytoplasma spp.strains from different hosts/countries is presented in Figure 4.According to the constructed tree, all Serbian tested 'Ca.Phytoplasma solani' as well as comparative 'Ca.Phytoplasma solani' strains and stolbur phytoplasmas were genetically homogenous and placed in the same cluster.Based on the same gene, sequences of Strawberry lethal yellows phytoplasma (CPA) strain NZSb11 and 'Ca.Phytoplasma australiense' strains were identical; therefore, these two strains were grouped together within the second cluster of the tree.An outgroup Acholeplasma laidlawii strain DSM23060 was placed on a monophyletic tree branch.Groups I (A, B, and D), II (A and D), VIII-A, X-A, XI-A, and XII-B were clearly separated in other tree branches/clusters, each corresponding to a specific 16Sr Group.Acholeplasma laidlawii strain NCTC10116 was placed on a monophyletic tree branch, as an outgroup (Figure 3).The NJ phylogenetic tree that was constructed based on the plsC gene sequences for the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and comparative Phytoplasma spp.strains from different hosts/countries is presented in Figure 4.According to the constructed tree, all Serbian tested 'Ca.Phytoplasma solani' as well as comparative 'Ca.Phytoplasma solani' strains and stolbur phytoplasmas were genetically homogenous and placed in the same cluster.Based on the same gene, sequences of Strawberry lethal yellows phytoplasma (CPA) strain NZSb11 and 'Ca.Phytoplasma australiense' strains were identical; therefore, these two strains were grouped together within the second cluster of the tree.An outgroup Acholeplasma laidlawii strain DSM23060 was placed on a monophyletic tree branch.

and D. carota
Groups I (A, B, and D), II (A and D), VIII-A, X-A, XI-A, and XII-B were clearly separated in other tree branches/clusters, each corresponding to a specific 16Sr Group.Acholeplasma laidlawii strain NCTC10116 was placed on a monophyletic tree branch, as an outgroup (Figure 3).The NJ phylogenetic tree that was constructed based on the plsC gene sequences for the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and comparative Phytoplasma spp.strains from different hosts/countries is presented in Figure 4.According to the constructed tree, all Serbian tested 'Ca.Phytoplasma solani' as well as comparative 'Ca.Phytoplasma solani' strains and stolbur phytoplasmas were genetically homogenous and placed in the same cluster.Based on the same gene, sequences of Strawberry lethal yellows phytoplasma (CPA) strain NZSb11 and 'Ca.Phytoplasma australiense' strains were identical; therefore, these two strains were grouped together within the second cluster of the tree.An outgroup Acholeplasma laidlawii strain DSM23060 was placed on a monophyletic tree branch.The NJ phylogenetic tree that was constructed based on the plsC gene sequences for the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and comparative Phytoplasma spp.strains from different hosts/countries is presented in Figure 4.According to the constructed tree, all Serbian tested 'Ca.Phytoplasma solani' as well as comparative 'Ca.Phytoplasma solani' strains and stolbur phytoplasmas were genetically homogenous and placed in the same cluster.Based on the same gene, sequences of Strawberry lethal yellows phytoplasma (CPA) strain NZSb11 and 'Ca.Phytoplasma australiense' strains were identical; therefore, these two strains were grouped together within the second cluster of the tree.An outgroup Acholeplasma laidlawii strain DSM23060 was placed on a monophyletic tree branch.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and 13 comparative Phytoplasma spp.strains according to sequences of the plsC gene.Acholeplasma laidlawii strain DSM23060 from the GenBank database served as an outgroup.

Biochemical Analyses of Peony, Mint, Dill, and Carrot Infected by 'Ca. Phytoplasma solani'
There was a significant difference for most of the tested biochemical parameters when asymptomatic and symptomatic plants were compared, except for total anthocyanidins contents in dill and the capacity of peony and mint extracts to neutralize superoxide anions and hydroxyl radicals, respectively (Table 3).
The total sugars content accumulated in symptomatic leaves was twice the amount of the sugars in healthy leaves.Lipid peroxidation intensity was clearly and dramatically

Biochemical Analyses of Peony, Mint, Dill, and Carrot Infected by 'Ca. Phytoplasma solani'
There was a significant difference for most of the tested biochemical parameters when asymptomatic and symptomatic plants were compared, except for total anthocyanidins contents in dill and the capacity of peony and mint extracts to neutralize superoxide anions and hydroxyl radicals, respectively (Table 3).

Biochemical Analyses of Peony, Mint, Dill, and Carrot Infected by 'Ca. Phytoplasma solani'
There was a significant difference for most of the tested biochemical parameters when asymptomatic and symptomatic plants were compared, except for total anthocyanidins contents in dill and the capacity of peony and mint extracts to neutralize superoxide anions and hydroxyl radicals, respectively (Table 3).and 13 comparative Phytoplasma spp.strains according to sequences of the plsC gene.Acholeplasma laidlawii strain DSM23060 from the GenBank database served as an outgroup.

Biochemical Analyses of Peony, Mint, Dill, and Carrot Infected by 'Ca. Phytoplasma solani'
There was a significant difference for most of the tested biochemical parameters when asymptomatic and symptomatic plants were compared, except for total anthocyanidins contents in dill and the capacity of peony and mint extracts to neutralize superoxide anions and hydroxyl radicals, respectively (Table 3).
The total sugars content accumulated in symptomatic leaves was twice the amount of the sugars in healthy leaves.Lipid peroxidation intensity was clearly and dramatically higher in the leaves of symptomatic plants, reaching twice the values of LP for the asymptomatic dill and peony plants.Reduced glutathione content was enhanced slightly but significantly in all tested species.Although the activity of PAL was markedly higher in the symptomatic leaves (1.25-3.5 times higher), only in mint and dill were the amounts of phenolic compounds positively correlated with this trend (possibly because of higher tannin and proanthocyanidins contents); however, the amount of phenolics in peony plants was lower in symptomatic leaves (except for anthocyanidins content) than in those of the other two species.Photosynthetic pigment contents (chlorophyll a and b and carotenoids) significantly decreased (2.45-5.08 times) in the leaves of all tested species.Although there was a difference in antioxidant capacity between asymptomatic and symptomatic plants within species, the mint and dill plants had enhanced antioxidant activity against superoxide anions (NBT test higher than 80% of neutralized radicals), while peony plants had pronounced antioxidant activity against DPPH radicals (DPPH test higher than 90% of neutralized radicals).

Principal Component Analysis (PCA)
In this study, PCA was used to investigate the correlation between different chemical parameters tested, such as TS, RG, PAL, LP, TP, TT, Tflav, Tpro, Tant, NBT, • OH, DPPH, Tcha, Tchb, and Car.The plant samples, including mint AL, mint SL, dill AL, dill SL, peony AL, and peony SL, were observed, and the content of these variables was utilized in the PCA.The first three principal components generated by the PCA were utilized to display the results visually in a PCA graph, which effectively differentiated the six plant samples.By applying PCA to the provided dataset, it was possible to differentiate between the samples based on the processing parameters used.This technique was utilized as a tool in exploratory data analysis to characterize and distinguish the input parameters for the neural network (as shown in Figure 5).
According to the PCA, the carrot SR and carrot AR samples were characterized by the highest Car content, while peony SL and peony AL were noted to obtain the most prominent value of DPPH, TS, TT TP, PAL, and RG content.The sample mint AL was characterized by augmented Tflav, Tcha, and Tchb contents.was utilized in the PCA.The first three principal components generated by the PCA were utilized to display the results visually in a PCA graph, which effectively differentiated the six plant samples.By applying PCA to the provided dataset, it was possible to differentiate between the samples based on the processing parameters used.This technique was utilized as a tool in exploratory data analysis to characterize and distinguish the input parameters for the neural network (as shown in Figure 5).The results of the PCA (shown in Figure 5) reveal a clear separation between the samples, indicating good discrimination.The first three principal components (having eigenvalues of 5.15, 3.36, and 2.17), which account for 71.22% (the total variance

Neurons in the ANN Hidden Layer
One of the crucial steps in designing an ANN is to determine the appropriate number of hidden layers and neurons in each layer, which depend on the complexity of the inputoutput relationship.As this relationship becomes more complex, more neurons need to be added to the hidden layer(s) [32].The optimal number of hidden neurons was selected by minimizing the difference between the predicted and desired output values, using the Sum of Squares (SOS) as the performance indicator during testing.The used MLP model was identified according to StatSoft Statistica's notation, which indicates the number of inputs, neurons in the hidden layer(s), and outputs.Based on the ANN performance, it was observed that the optimal number of hidden neurons for the calculation of TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, and Car was nine (using MLP 8-9-13), as it produced high r 2 values (1.000 during the training period) and low SOS values.The optimal networks used for predicting TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, and Car were able to produce reasonably accurate output values across a wide range of process variables.In most cases, the predicted values were very close to the experimental or target values, as evidenced by the high r 2 values obtained from the ANN models.The complexity of the ANN models was quite high, with 211 weight biases required for TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, and Car calculations.However, the models were able to fit the experimental data quite well, owing to the high nonlinearity of the system under investigation [31,32].During the training period, the r 2 values between the experimental measurements and the ANN model outputs for TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, and Car were equal to 1.000.

Sensitivity Analysis
A sensitivity analysis was conducted to determine the relative influence (RI) of inputs (such as plant type and treatment) on TS, RG, PAL, LP, NBT, • OH, TP, TT, Tflav, Tpro, Tant, DPPH, Tcha, Tchb, and Car, The results, shown in Figure 6, indicate how outputs change with respect to infinitesimal changes in inputs, reflecting both experimental errors and the inputs' influence on the outputs.

Discussion
Phytoplasmas belonging to the "stolbur" group 16SrXII-A ('Ca.Phytoplasma solani') are known to affect a wide range of wild and cultivated plants worldwide, causing severe

Discussion
Phytoplasmas belonging to the "stolbur" group 16SrXII-A ('Ca.Phytoplasma solani') are known to affect a wide range of wild and cultivated plants worldwide, causing severe economic damage [34].Since its first report on pepper in Serbia in 1949, 'Ca.Phytoplasma solani' has been described on diverse host plants (carrot, corn, grapevine, peony, peppermint, potato, etc.) in this country [12,[35][36][37][38][39][40].Species designation within phytoplasmas was primarily based on conserved 16S rRNA gene; however, since it does not allow distinctions among closely related Ca.Phytoplasma species, recent studies are encouraging the additional use of less conserved, housekeeping genes to enhance the resolving power [41].However, the use of 16S rRNA sequences in this study allowed us to determine phylogenetic relatedness among the tested Serbian 'Ca.Phytoplasma solani' strains.The performed phylogenetic analysis with 16S rRNA sequences of the three Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia (B15 and B18) and D. carota (ML_NS-2016) and various reference Phytoplasma spp.strains belonging to different 16Sr Groups indicated an affiliation between the tested strains and the 16SrXII-A Group and genetic homogeneity with the other (tested/reference) 'Ca.Phytoplasma solani' strains placed within this group, regardless of the host or country of origin.Additionally, earlier, D. carota strain ML_NS-2016 was confirmed to be genetically homogenous with different 'Ca.Phytoplasma solani' strains originating from various countries (Bulgaria, Canada, China, France, Greece, Italy, Poland, Russia, Serbia, and Turkey) and hosts (corn, grapevine, parsley, parsnip, pea, periwinkle, potato, red sage, tobacco, tomato, and valerian), all belonging to the "stolbur" group (16SrXII-A) [7].Similar to phylogenetic analysis based on 16S rRNA, a phylogenetic tree constructed based on the sequences of plsC gene, encoding the 1-acyl-sn-glycerol-3-phosphate acyltransferase enzyme involved in glycerophospholipid metabolism [42], placed the tested Serbian (A. graveolens strain STOL2 and Mentha × piperita strain STOL3) and comparative 'Ca.Phytoplasma solani' strains in the same "stolbur" tree cluster.This gene was discriminatory enough to distinguish 'Ca.Phytoplasma australiense' (first added as Strawberry lethal yellows phytoplasma) from the "stolbur" group.Some future research on Serbian 'Ca.phytoplasma solani' strains should be directed toward the sequencing and analysis of multiple housekeeping genes such as tuf, secY, vmp1, stamp, etc. [8,43], which could eventually reveal some genetic differences.
Artificial neural networks (ANNs) have found application in diverse domains for forecasting the outcomes of intricate systems by processing input data.ANNs belong to a class of machine learning algorithms capable of learning from data and predicting results with precision and dependability, such as identifying plant diseases.ANNs have exhibited remarkable potential in predicting biotic stress in plants resulting from diverse pathogens, enabling real-time plant health monitoring and offering a pre-emptive alert mechanism for farmers, agronomists, and researchers to forestall disease outbreaks.Developing an effective early warning system for the detection of plant diseases in different plant species is critical for improving crop yield and quality.The identification and prevention of diseases of crops are essential for improving crop prediction.Based on our results, the ANN model can be considered as an efficient computational methodology for modeling and predicting the effect of the biotic stress in different plants induced by 'Ca.Phytoplasma solani'.

Figure 2 .
Figure 2. ANN topology with three layers (input, output, and hidden) with weights, biases, and transfer functions.

Figure 2 .
Figure 2. ANN topology with three layers (input, output, and hidden) with weights, biases, and transfer functions.

Figure 3 .
Figure 3. Neighbor-joining phylogenetic tree showing the position of the three tested Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia and D. carota and 23 comparative Phytoplasma spp.strains belonging to different 16Sr Groups.Acholeplasma laidlawii strain NCTC10116 from the GenBank database served as an outgroup.

Figure 3 .
Figure 3. Neighbor-joining phylogenetic tree showing the position of the three tested Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia

Figure 3 .
Figure 3. Neighbor-joining phylogenetic tree showing the position of the three tested Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia and D. carota and 23 comparative Phytoplasma spp.strains belonging to different 16Sr Groups.Acholeplasma laidlawii strain NCTC10116 from the GenBank database served as an outgroup.

Figure 3 .
Figure 3. Neighbor-joining phylogenetic tree showing the position of the three tested Serbian 'Ca.Phytoplasma solani' strains from P. tenuifolia and D. carota and 23 comparative Phytoplasma spp.strains belonging to different 16Sr Groups.Acholeplasma laidlawii strain NCTC10116 from the GenBank database served as an outgroup.
and 23 comparative Phytoplasma spp.strains belonging to different 16Sr Groups.Acholeplasma laidlawii strain NCTC10116 from the GenBank database served as an outgroup.

Figure 4 .
Figure 4. Neighbor-joining phylogenetic tree showing the position of the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and 13 comparative Phytoplasma spp.strains according to sequences of the plsC gene.Acholeplasma laidlawii strain DSM23060 from the GenBank database served as an outgroup.

Figure 4 .
Figure 4. Neighbor-joining phylogenetic tree showing the position of the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens

Figure 4 .
Figure 4. Neighbor-joining phylogenetic tree showing the position of the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and 13 comparative Phytoplasma spp.strains according to sequences of the plsC gene.Acholeplasma laidlawii strain DSM23060 from the GenBank database served as an outgroup.

Figure 4 .
Figure 4. Neighbor-joining phylogenetic tree showing the position of the two tested Serbian 'Ca.Phytoplasma solani' strains from A. graveolens and Mentha × piperita and 13 comparative Phytoplasma spp.strains according to sequences of the plsC gene.Acholeplasma laidlawii strain DSM23060 from the GenBank database served as an outgroup.

Figure 5 .
Figure 5. Biplot graph of chemical analysis for asymptomatic (A) and symptomatic (S) samples of the leaves of peony, mint, dill, and carrot.

Figure 5 .
Figure 5. Biplot graph of chemical analysis for asymptomatic (A) and symptomatic (S) samples of the leaves of peony, mint, dill, and carrot.

Table 1 .
List of Phytoplasma spp.strains from GenBank used to construct the phylogenetic tree based on 16S rRNA sequences.

Table 1 .
List of Phytoplasma spp.strains from GenBank used to construct the phylogenetic tree based on 16S rRNA sequences.

Table 1 .
Cont.Phytoplasma solani strains of interest in this study; b outgroup.
a Serbian Ca.

Table 2 .
List of Phytoplasma spp.strains from GenBank used to construct the phylogenetic tree based on plsC gene sequences.
a Serbian 'Ca.Phytoplasma solani' strains of interest in this study; b outgroup.

Table 3 .
[7]t.-difference between A and S according to the t-test (p < 0.05); nd-no difference.acarrotdata retrieved from a previous study by Mitrović et al.[7]and not the primary result of this work. *