Relationship between biochemical and photosynthetic traits with Asian soybean rust

in this to investigate the relationship between and photosyntetic parameters in soybean with Two experiments performed in a randomized complete block with three treatments (water, Tween 20, and The evaluated traits were: severity, chlorophyll pools, concentration of phenolic compounds content, enzyme activity, and photosyntetic parameters. Based on Pearson correlation, correlation network and path analysis it was verified that the severity had high correlation with almost all traits evaluated meanwhile photosynthesis is weakly related with almost all traits. Therefore, the occurrence of ASR affects directly the traits related to enzymatic activity and phenolic compounds content in soybean plants inoculated with ASR. The management of ASR is important to keep in normal levels the rates of photosynthesis carried out by the plant, and thus not affect the yield. Besides that, understanding the biochemical mechanisms and ecophysiological responses that occur during the soybean- P. pachyrhizi interaction has a great importance for breeding programs, as it will help for obtaining cultivars or efficient methods in fungus control.


INTRODUCTION
Soybean [Glycine max (L.) Merril] currently represents one of the main export products of Brazilian agribusiness, contributing significantly to the trade balance, being the culture that has grown most in the last three decades. Soybean production worldwide in the 2015/2016 harvest was 312 million tons in a cultivated area of 119 million hectares (USDA 2017). In Brazil, soybean production in the 2015/2016 harvest was 95 million tons in an area of 33 million hectares (CONAB 2017). Parallel to this growth, it is worth mentioning the increase in phytosanitary problems affecting this crop, with 3926 TÁSSIA B. DE OLIVEIRA et al. more than 100 species of pathogens being reported, of which about 35 are of great economic importance (Dorrance et al. 2003). Among these diseases, Asian soybean rust (ASR), caused by the biotrophic fungus Phakopsora pachyrhizi H. Sydow & Sydow, is considered a serious problem for agribusiness in soybean producing countries (Goellner et al. 2010). Its rapid expansion and high potential for causing losses in soybean have made this disease one of the main diseases of the crop, which can reach levels that make harvesting inviable (Yorinori et al. 2005). Although chemical control is the most used measure to manage ASR, the indiscriminate use of agrochemicals causes irreversible damage to the environment, causing environmental imbalance and favoring the selection of new physiologically resistant races. Another way of handling ASR is the use of phytoregulators such as methyl jasmonate that stands out as a signaling compound that activates plant defense mechanisms in response to stresses, such as symptoms caused by pests and pathogens (Browse 2009). Jasmonate still plays an important role in the regulation of the expression of genes encoding stress-related proteins (Ruiz-Garcia and Gômez-Plaza 2013). Besides the management forms mentioned above, genetic breeding is a suitable alternative for the management of ASR using resistant varieties.
For the plant breeder to identify tolerant soybean varieties for ASR, it is first necessary to know the association between characters, especially if the selection in one of them presents difficulty due to the low heritability and/or difficulty in the measurement. Thus, the estimates of the correlation coefficients allow to evaluate the magnitude and direction of the relation between two characters, and consequently obtain gains for one of them through the indirect selection using the other character. In some cases, indirect selection, based on the correlated response, can lead to progress faster than the direct selection of the desired variable (Cruz et al. 2012).
Despite the great utility of correlation coefficients for a better understanding of the relationships between variables, they only estimate the association between variables and do not allow conclusions about cause and effect relationships. Studies about the correlation coefficient split are done by path analysis, developed by Wright (1923). This analysis consists of the study of the direct and indirect effects of characters on a main variable, whose estimates are obtained by means of regression equations, in which the variables are previously standardized (Cruz et al. 2012, Vencovsky andBarriga 1992).
In the literature there are few studies that report the association between biochemical and physiological variables with components related to gas exchange and the occurrence of ASR. Therefore, the goal of this work was to investigate the relationship between biochemical and photosyntetic parameters in soybean with ASR.

EXPERIMENTAL DESIGN
Two experiments were performed in a greenhouse at Embrapa Soja in the Warta District, in Londrina, Parana State. Soybean plants of BRS 361 genotype, considered susceptible to ASR, were used. Seeds were obtained from Embrapa Soja Germplasm Bank. The plants were kept in greenhouse (relative humidity 65%, 28 °C) until the V5 stage (Fehr et al. 1971) and then submitted to the exogenous application of the different products using pressurized sprayer with CO 2 .
The experiments were set up in a randomized block design with three treatments (water, Tween 20 e methyl jasmonate) with six replications, and each experimental plot consisted of a pot containing five plants. It was not observed interaction between experiments, then it was used the mean between APPLICATION OF PATH ANALYSIS IN SOYBEAN 3927

EVALUATION OF ENZYMATIC ACTIVITY
To determine the activity of the enzymes chitinase (CHI) and β-1,3-glucanase (GLU), leaves in the V4 and V5 stages were collected, using two plants per sample. The samples were collected at 48, 96, 144 and 192 hours after inoculation (hai). Leaves of non-inoculated plants were also collected in the same period. The samples were stored individually in aluminum foil and immediately frozen in liquid nitrogen and stored in ultra-freezer at a temperature of -80 °C. From each leaf extract obtained, enzymatic trails were performed in triplicate. The activity of CHI and GLU was determined. The total protein concentration in each sample was determined following the method of Bradford (1976).
The foliar extract was obtained by maceration of 0.3 g of leaf tissue crushed in a liquid nitrogen mortar with addition of 2% (p/v) polyvinylpolypyrrolidone (PVPP). The obtained powder was homogenized in 2 mL of 50 mM sodium phosphate buffer (pH 6.5) and 1 mM phenylmethylsulfonyl fluoride (PMSF). The homogenate was centrifuged at 20,000 x g, at 4 °C for 25 min. The supernatant was used to determine the enzymatic activity.
CHI activity was determined according to the method proposed by Roberts andSelitrennikoff (1988) modified by Harman et al. (1993), using as substrate p-nitrophenyl-β-D-N,N'diacetylchitobiose. The incubation system consisted of adding a 20 μL aliquot of the supernatant obtained as described above, to a mixture containing 470 μL of 50 mM sodium acetate buffer (pH 5.0) and 10 μL of PVPP (2 mg mL -1 ) and maintained at 37 °C for 2 h. The reaction was paralyzed by adding 0.5 mL of 0.2 M sodium carbonate and the absorbances was determined at wavelength 410 nm. To calculate the enzymatic activity was used the molar extinction coefficient of 7 x mM -1 cm -1 with the activity of CHI expressed in mM of p-nitrophenyl produced by min -1 mg -1 protein.
experiments to perform the statistics analysis in this study.

CHEMICALS USED
The products used were: water as a control treatment; the surfactant (Tween 20 -Polyoxyethylene 20 sorbitan monolaurate) at the concentration of 0.02%; and methyl jasmonate (MJ) at the concentration of 1.25 mM prepared in 0.02% Tween 20 (MJ + T20).
Plants treated with MJ were kept in another greenhouse for 48 hours under the same temperature and humidity conditions to avoid induction in the control plants (Water and Tween 20), since MJ is known as a volatile regulatory phyto (Seo et al. 2001).

PLANTS INOCULATION WITH P. pachyrhizi
After 24 hours of spraying the products described above, the plants were inoculated. To prepare the inoculum of P. pachyrhizi the fungus was previously kept in the cultivar BRS 316, in a greenhouse, at Embrapa Soja. The spores were collected with light beats on the leaves with uredinospores with pustules, on white paper and packed in microtubes. Then, a suspension was prepared with distilled water containing 0.01% Tween 20, and adjusted to a final concentration of 1.4x10 5 urospores/ ml. After inoculation, the plants were kept under nebulization to stimulate the infectious process as it requires a high temperature, 23 °C to 25 °C and relative humidity greater than 95% in the absence of light for 12 hours.
EVALUATION OF ASR SEVERITY ASR severity was evaluated 16 days after plant inoculation in V2 leaves, using a scale proposed by Godoy et al. (2006). This scale takes into account the percentage of infected leaf area and ranges from 0.6% to 78.5%.
GLU activity was determined according to the method described by Lever (1972), with modifications: 3,5-dinitrosalicylic acid (DNS) to replace p-hydroxybenzoic acid hydrazide (Miller 1959). The reaction consisted of 230 μL of 100 mM sodium acetate buffer (pH 5.0), 250 μL of the laminarin substrate solution (4 mg mL -1 ) and 20 μL of the plant extract. The system was incubated at 45 °C for 30 min, and after this period 1 mL of DNS was added to the reaction and then this mixture was heated at 100 °C for 5 min. After cooling on ice to 30 °C the absorbance of the samples was determined at 540 nm and the results are expressed as absorbance units per min -1 mg -1 protein.

ANALYSIS OF PHENOLIC COMPOUNDS
The samples were composed of a V5 leaf collected 48, 96 and 144 hai. Leaves of not inoculated plants were also collected in the same periods. The samples were stored individually in aluminum foil packaging in liquid nitrogen and stored in a freezer (-20 °C) until extraction and analysis, which were performed according to the standard methodology of Laboratório de Química e Ecologia of EMBRAPA Soja, described below.
At the time of extraction, the leaves were crushed in a liquid nitrogen mortar and placed in a Falcon tube. Subsequently, the leaves were weighed (500 mg), and 5 mL of 90% methanol was added, and the solution was then subjected to the ultrasonic bath for 20 min. Then the samples were centrifuged at 9,880 x g at 4 °C for 12 min, vacuum dried, resolubilized in 1.5 mL of 80% methanol and homogenized manually. And then, the plant extracts were filtered through a 0.45 μm Millipore® membrane and transferred to tubes of the automated sampler and analyzed HPLC (High Performance Liquid Chromatography) -Shimadzu -Prominence model.
The methanolic extracts from the samples were analyzed on C18 column (250 mm long and 4.6 mm internal diameter, 5 μM particles). Aliquots of 10 μL were injected automatically into the equipment, with CBM-20ª controller, SPD-20A detector, DGU 20A5 degasser, LC-20AT pump, SIL-20A automatic sampler and CTO 20A furnace. The mobile phase was composed of two solvents: (A) 2% acetic acid (HOAc) and (B): A mixture of methanol, acetic acid and MilliQ® water (MeOH:HOAc:H 2 O; 18:1:1). The linear gradient system used in the analysis started from the initial condition with 75% of A and 25% of B reaching, and after 40 min, the inverse situation, in other words, 25% of A and 75% of B, where it was maintained during 5 min. At 45 min, it returned to the initial situation, remaining for 5 min, before the next injection. The solvent flow was 1 mL min -1 and ultraviolet region analysis was obtained at wavelength 260 and 280 nm.
The peak area of each identified substance was measured and multiplied by the correction factor of each standard compound, taking into account the initial tissue mass and volume of resolubilization. In order to obtain the correction factor of each compound, different concentrations of isoflavones and phenolic acids standards (6.25, 12.5, 25, 50, 100 μg mL -1 ) were injected into the HPLC. From the values of the areas of the standard compounds and their concentrations, a dispersion graphic was constructed and a linear regression equation and its correction factor was obtained.

EVALUATION OF PHOTOSYNTETIC PARAMETERS AND CHLOROPHYLL POOLS
Evaluations of photosyntetic parameters and chlorophyll pools were performed on central V4 leaflets, previously marked at 39, 87, 135 and 183 hai, between 9 and 10 o'clock in the morning. A portable photosynthesis analyzer was used, model LCpro-SD of ADC brand. Active photosynthetic radiation of 1044 μmol m -2 .s -1 was used. The following physiological parameters were evaluated: C i -internal CO 2 concentration; E -transpiration rate; G s -stomatal conductance to water vapor; Anet CO 2 assimilation rate. The chlorophyll pools was determined indirectly through the SPAD index, using the chlorophyllometer SPAD-502. From the SPAD index obtained in the chlorophyllometer, the total chlorophyll pools was calculated using the equation described below, which was generated in calibration and validation tests performed at Embrapa Soja (Unpublished data): CC=((SPAD * 0,0007)-0,0071) Where CC is the chlorophyll pool expressed in mg.cm -2 .

STATISTICAL ANALYSIS
Initially, analysis of variance (ANOVA) was performed for each character according to the statistical model described below: Where: Y ij is the observation in the j-th block, evaluated in the i-th treatment; µ is the overall mean of the experiment; B j is the fixed effect of the j-th block; T i is the fixed effect of the i-th treatment; e ij is the random error associated with observation Y ij .
Tukey test was performed to verify significant statistic difference among treatments. Subsequent, the phenotypic correlations (r P ) between the pairs of characters were estimated according in the equation below: Where: COV P(xy) is the phenotypic covariance between the variables X e Y; 2 Px σ is the phenotypic variance of the variable X; 2 Pŷ σ is the phenotypic variance of the variable Y. Phenotypic correlations were tested by the t-test with n-2 degrees of freedom.
In order to graphically express the functional relation of the estimates of the phenotype correlation coefficients between the variables (r ij ) the correlation network was used, in which the proximity between the nodes (traces) is proportional to the absolute value of the correlation between these nodes. The thickness of the edges was controlled by applying a cut-off value of 0.60, which means that only |r ij | ≥ 0.60 have their edges highlighted. Finally, positive correlations were highlighted in green, while negative correlations were represented in red.
Subsequently, a multicollinearity analysis of the X'X correlation matrix was performed, which revealed severe multicollinearity (condition number > 1000) according to Montgomery et al. (2015) classification. Therefore, to carry out the path analyzes a constant k = 0.10 was added to the diagonal of the X'X matrix for path analysis, which provided weak multicolinearity (condition number < 100).
Path analysis, considering the SEV and A as main dependent variables, was performed by the model described below: The coefficient of determination (R²) of the path analysis was obtained by the equation

ANOVA AND ESTIMATES OF GENETIC PARAMETERS
In order to verify if there was a significant difference between the treatments and to estimate the genetic and experimental parameters, the analysis of variance for the 15 characters was performed.
Significant difference was verified for seven traits (DZI, MD, MG, SEV, C i , GLU e CHI), while for the other traits there was no significant difference between the treatments (Table I). CV e higher than 20% was observed for all phenolic compounds content evaluated, except for MG (Table I). However, the CV e for the variables related to photosyntetic parameters and chlorophyll pools, enzymatic activity and severity were low (below 20%).
The coefficient of determination (H 2 ) was high for most of the evaluated variables, except for GI and R, while intraclass correlation was high only for variables DZI, MD, MG, SEV, C i , GLU e CHI (Table II).
The quadratic treatment component was high for most of the evaluated variables, while the coefficient of genetic variation was high only for the variables DZI, MD, MG, C, SEV e CHI (Table  II).
Based on Tukey test, it was observed that methyl jasmonate was able to reduce ASR effects on soybean plants, decreasing the severity (Table  III). In addition the application of methyl jasmonate increased enzymatic activity which GLU increased 25% and CHI increased 100% compared with other treatments. On the other hand, phenolic compunds content decreased with the application of methyl jasmonate such as DZI, GI, MD, MG and C.

PHENOTYPIC CORRELATIONS
In order to verify the linear relationship between severity, phenolic compounds content, enzymatic activity, photosyntetic parameters, and chlorophyll pools, the correlation between these traits was estimated and the correlation network was used to graphically express these results. It was observed that of the 105 possible correlations 82 (78%) were considered high (above 0.5), where 42 (40%) were negative and 40 (38%) were positive (Table IV).

3932
TÁSSIA B. DE OLIVEIRA et al. evaluated variables, except for R and C i (Figure 1). The SEV has significant and positive correlations with DZI and MD, and negative with phenolic compounds content (R, E and GI) and E. Chlorophyll pools has significant and positive correlations with phenolic compounds content (DZI, MD, MG and C) and SEV, and negative with all compounds related to photosyntetic parameters (E, C i and G s ) and enzymatic activity (GLU and CHI). Substances related to enzymatic activity had a significant and high correlation with photosyntetic parameters.

PATH ANALYSIS
Path analyzes were performed in order to verify how the variables related to phenolic compounds content, enzymatic activity and photosyntetic parameters affect directly or indirectly the main variables (SEV and A). It was observed that the estimation coefficient was of high magnitude (= 0.98 e 0.97) and the residual effect of low magnitude (0.12 and 0.18) for SEV and A respectively (Tables  V and VI).

ANOVA AND ESTIMATES OF GENETIC PARAMETERS
The existence of a significant effect for the treatments of variables DZI, MD, MG, SEV, GLU and CHI showed that inoculation with the P. pachyrhizi fungus causes changes in enzymatic activity, phenolic compounds content in soybean plants inoculated with ASR. Thus, understanding the biochemical mechanisms and ecophysiological responses that occur during the soybean P. pachyrhizi interaction has a great importance for breeding programs, as it will help for obtaining resistant cultivars or efficient methods in fungus control.
Some of the evaluated phenolic compounds (GI, R and C) presented CV e estimates higher than their respective estimates for CV g . In addition, these traits presented low intraclass correlation estimates. Other variables that presented low intraclass correlation estimates were E and A. These results indicate that the non-significant effect among treatments occurred due to the absence of correlation between the replicates (low r), which inflated the residual variance and even made it impossible to obtain the genetic parameter estimates for some traits (Cruz et al. 2012) such as GEI, CC and G s . From these results it is possible to recommend the use of a greater number of repetitions to evaluate GI, R, GEI, C, CHL, E, G s and A variables in future experiments to increase accuracy in their estimates.
The estimation of H² makes it possible to know the fraction of the phenotypic variation between treatments that is genetic nature. The high estimates of this parameter for the variables DZI, MD, MG, SEV, C I , GLU and CHI and moderate for E and A are encouraging for breeding programs.      From the results reported here, other researches can be conducted to obtain genetic gains in these traits and may guide breeders on the genetic and physiological mechanisms of soybean defense against the P. pachyrhizi pathogen. The application of methyl jasmonate reduced ASR drastically (Table III). Findings in this study are according to results reported in the literature which proved that methyl jasmonate can efficiently reduce diseases severity in many crops such as sweet cherry fruit (Yao and Tian 2005b), Alternaria brassicicola, Botrytis cinerea and Plectosphaerella cucumerina in Arabidopsis thaliana (Thomma et al. 2000), Monilinia fructicola in peach (Yao and Tian 2005a), Erysiphe necator in grapes (Belhadj et al. 2006) and Colletotrichum acutatum in plum (Cao et al. 2008).

Estimates of the direct (main diagonal) and indirect effects of the variables on the main variable
Methyl jasmonate can reduce disease severity affecting fungal growth suppression as well as spores germination suppression, which both were reported for Alternaria alternata (Chen et al. 2014), and A. brassicicola in Arabidopsis thaliana (Thomma et al. 1998).
Another factor that should be considered is the increase of enzymatic activity with the methyl jasmonate application. According to Cruz et al. (2013) chitinase and β-1,3-glucanase are important to decrease the P. pachyrhizi colonization in soybean leaves tissue. These enzymes get involved to the plant defense against pathogen by the chitin and β-1,3-glucan hydrolytic action present on the fungal cell wall, and consequently the oligosaccharides liberation that may act as elicitors of the plant defense mechanism against fungal (Mohammadi et al. 2002). The overexpression of genes affecting chitinase and β-1,3-glucanase in plants have increased the resistant of plant for several pathogens (Van Loon 1997). Perhaps the enzyme activities increase favored the ASR severity reduction by methyl jasmonate application (Table III).
Despite that, there is no significant difference among treatments for photosyntetic parameters, plants treated with water had a slight decrease on them. Berger et al. (2007) proposed that plants reduce the photosynthesis and assimilation for increasing the respiration and other processes related with defense.

PHENOTYPIC CORRELATIONS
In a breeding program, knowledge about the relationship between variables is essential to achieve the desired genetic progress in the selection stages reducing cycle time and resources (Borém and Miranda 2013). In this work, we evaluated phenolic compounds content, enzymatic activity and photosyntetic parameters, which are difficult to measure and require high cost equipment and great demand for labor and time. Thus, estimates of phenotypic correlations reported here may guide physiologists and breeders about which selection strategies should be used in soybean. Besides that, a trait that has low H² and/or greater difficulty of measurement can be evaluated indirectly with the measurement of another correlated variable (Cruz et al. 2012).
The large number of traits evaluated in this work provided 105 estimates of phenotypic correlations, and the interpretation of the correlation matrix is a difficult task. Thus, we sought to use the correlation network to graphically visualize the linear relationships between phenolic compounds content, enzymatic activity and photosyntetic parameters. The two-dimensional network allows to detect complex phenotypic patterns that are difficult to extract with other techniques (Figure 1). The efficiency of this innovative technique has already been reported in studies that evaluated a large number of metabolic compounds in tomato (DiLeo et al. 2011, Ursem et al. 2008.
In the correlation network, positive correlations were represented in green color and negative correlations in red. The greater was the thickness of the traces, the greater was the magnitude of the correlations. It was possible to observe a great association of the phenolic compounds content, enzymatic activity and photosyntetic parameters, except for R, C i and A, that are further away at the ends of the two-dimensional network.
A was the trait that smalest correlated with the others, presenting a positive linear relation with R and negative with C i . The negative correlation between the traits A and C i was also verified in other studies (Concenço et al. 2008, Ferreira et al. 2011 and indicates that when the photosynthetic rate is considered higher the CO 2 concentration tends to be lower. It is important to highlight the evaluation of the R, since this compound correlated negatively with SEV and positively with A. It has been reported by Hoffmann-Campo et al. (2006) that R has a negative effect on the intake of soybean leaves by Anticarsia gemmatalis Hübner (Lep.: Noctuidae). However, there are still no reports on the negative effect of this compound on soybean diseases.

PATH ANALYSIS
Although important, the phenotypic correlation coefficient can produce misunderstandings about the relationship between two traits, and may not be an actual measure of cause and effect. A high or low correlation coefficient between two traits may be the result of the effect that a third trait or a group of traits has on the pair, not giving the exact relative importance of the direct and indirect effects of these factors (Cruz et al. 2012). Therefore, we performed the path analysis, which investigates the cause and effect relationship. This analysis provides a detailed knowledge of the influences of the traits involved, and justify the existence of positive and negative correlations, of high and low magnitude, between traits.
However, to obtain the direct and indirect effects of path analysis, it is necessary a well conditioned matrix X'X. In the presence of multicollinearity, the variances associated to the path coefficient estimators can reach exceedingly high values, becoming unreliable. Besides that, the parameter estimates may assume no reasonable values or no coherence with the biological phenomenon (Cruz et al. 2012). One way to evaluate multicollinearity is the condition number, which results from the ratio of the largest to the smallest eigenvalue of the correlation matrix. According to the criterion presented by Montgomery et al. (2015), the matrix of estimates of phenotypic correlations showed a severe multicollinearity, since the condition number was higher than 1000. Thus, it was necessary to carry out the path analysis, through the inclusion of a constant K = 0.10, in the diagonal of the correlation matrix.
The high estimates of the determination coefficients obtained when SEV and A were considered as main dependent traits, besides , coumestrol (C), severity (Sev), chlorophyll (CC), C i -internal CO 2 concentration; Etranspiration rate; G s -stomatal conductance to water vapor; A -net CO 2 assimilation rate; β-1,3-glucanase (GLU), chitinase (CHI). low estimates of residual effect in both cases, show that almost all the variation in these traits can be explained by the other evaluated ones. These results are important because A presented a low estimate of H², besides being difficult to measure and to be subject to the variations of the environment. Although SEV has a high estimate of H², it is extremely difficult to measure because it is established from scores and can become biased when more than one person performs the evaluations, which are common situations in soybean breeding programs.
For evaluating the direct effects of the path analyzes performed it was verified that there was a high direct negative effect of SEV on A and vice versa. This association was not identified by the correlation estimates and reiterates the importance of path analysis to verify the association between variables. The CO 2 assimilation rate reduction in inoculated soybean plants may be due to the destruction of leaves tissues that caused yellowing and premature defoliation. However, we did not observed visual alteration during the evaluation period and no leaf blade abscission. According to Goellner et al. (2010) visual symptoms caused by P. pachyrhizi can be observed after five through eight days after inoculation, and consequently it may explain why it was not observed any symptoms since the evaluation period was from one to eight days after inoculation. Therefore, the most important hypothesis to explain the CO 2 assimilation rate reduction could be the occurrence of any damage in the photosynthetic mechanism due to metabolic changes happen before visual changes.
It was possible to verify that R has direct effects of high magnitude on both traits (SEV and A). This relationship is translated into cause and effect, because R also correlated positively with A and negatively with SEV. The negative relationship between SEV and R is due to since R can induce the hydrolysis of membrane phospholipids altering the barriers and permeability of cell membranes (Gazzoni et al. 1997).
E also presented high correlation as direct effects on SEV and high H², however, no correlation was observed between E and A. These results demonstrate that E should be evaluated in segregating populations where selection will be performed for resistance to P. pachyrhizi. The indirect selection aims for increasing the photosynthetic efficiency of the genotypes can be performed by means of the evaluation of C i , which has high correlation with A, and high H². Situations where environment stress happens excitation energy increases, and then it exceeds the required content for photosynthesis metabolism (Asada 1999). The increase of excitation energy is associated with the CO 2 assimilation rate decrease due to the reduction of stomatal conductance to water vapor and the increase of production of reactive oxygen species (Mateo et al. 2004). Reactive oxygen species induced stomatal closure, and consequently it can explain the stomatal conductance to water vapor reduced (Munemasa et al. 2007). According to Loehrer et al. (2008) that studied the interaction between P. pachyrhizi and Arabidopsis observed that inoculated plants accumulated hydrogen peroxide which may be influenced the stomatal closure. The stomatal closure influences the CO 2 influx decrease, and consequently in the lack of losses in the water vapor via transpiration. This fact may cause damage when plants are infected by fungus due to the need of assimilates for plant defense plus the extra drain which are the pathogenic agents (Berger et al. 2007). On the other hand, the stomatal closure when plants are infected by fungus is a strategy used by plants to avoid the pathogenic agents penetration through the stomata (Gudesblat et al. 2009). This strategy may be useful since the most etiologic agents that causes rust use the stomata to start the infectious process. Unfortunately, this strategy is not useful for P. pachyrhizi due to the penetration by this Based on results reported here, no trait was shown to be promising for the indirect selection of genotypes with high photosynthetic efficiency and resistant to P. pachyrhizi. However, this research will contribute with future work aimed at the establishment of selection indexes in order to obtain multiple gains in SEV and A in soybean plants.

CONCLUSIONS
Methyl jasmonate was efficient to reduce ASR effects on soybean plants.
The occurrence of ASR affects directly and with high magnitude traits related to enzymatic activity, phenolic compounds content and photosyntetic parameters in soybean, reducing the concentration of some traits and increasing in others.
The management of ASR is important to maintain the photosynthetic rates carried out by the plant at normal levels, and thus not affect the soybean yield at the end of the reproductive cycle.