#### Abstract

Multivariate statistical techniques were used to compare the relationship between yield and its related traits under noninoculated and inoculated cultivars with mycorrhizal fungus (*Glomus intraradices*); each one consisted of three wheat cultivars and four water regimes. Results showed that, under inoculation conditions, spike weight per plant and total chlorophyll content of the flag leaf were the most important variables contributing to wheat grain yield variation, while, under noninoculated condition, in addition to two mentioned traits, grain weight per spike and leaf area were also important variables accounting for wheat grain yield variation. Therefore, spike weight per plant and chlorophyll content of flag leaf can be used as selection criteria in breeding programs for both inoculated and noninoculated wheat cultivars under different water regimes, and also grain weight per spike and leaf area can be considered for noninoculated condition. Furthermore, inoculation of wheat cultivars showed higher value in the most measured traits, and the results indicated that inoculation treatment could change the relationship among morphological traits of wheat cultivars under drought stress. Also, it seems that the results of stepwise regression as a selecting method together with principal component and factor analysis are stronger methods to be applied in breeding programs for screening important traits.

#### 1. Introduction

Bread wheat (*Triticum aestivum* L.) is one of the most widely cultivated and important food crop in the world. Development of high yielding wheat cultivars is the major objective of breeding programs [1]. On the other hand, targeted efforts to breed genotypes for improved mycorrhizal symbiosis result in increased yield in crops under a wide range of environmental conditions and contribute toward sustainability of agricultural ecosystems in which soil-plant-microbe interactions will be better exploited. Screening genotypes via molecular biology and traditional breeding techniques can increase productivity of symbioses and eventually result in increased economic yield of crop plants [2].

Arbuscular mycorrhizal (AM) fungi colonize the roots of most monocotyledons and dicotyledons, including important crops such as rice, maize, and wheat despite their different root architecture and cell patterning. In nature, mineral nutrient acquisition and water uptake by plant roots are often assisted by symbioses with beneficial AM fungi. During this intimate association, the extraradical hyphal mycelium acquires minerals from the soil beyond the zone accessible by the roots. Inside the root, a considerable proportion of the minerals is delivered to the host in exchange for carbohydrates [3]. However, a few breeding programs and physiological researches have examined the morpho-physiological effects of mycorrhizal inoculation on wheat cultivars under drought stress. The benefit arising from the mycorrhizal symbiosis was not proportional to the extent of the root colonization in durum wheat genotypes; because various genotypes showed various responses to mycorrhizal inoculation [2, 4].

Aiming to providing methods in order to attain higher production, breeders have used important yield components as selection criteria to perform higher yields via indirect selection [5, 6]. Kumbhar et al. [7] illustrated high efficiency of tiller and kernel production on wheat yield. Grain weight/spike has been reported as the most closely variable contributing to grain yield, and it has often been used in selecting high yielding wheat genotypes [7]. 1000-grain weight has been shown as the main yield component for 20% of variation in wheat grain yield. Moghaddam et al. [8] found a negative correlation between plant height and grain yield of wheat due to the lower number of grains/spike. The results of Leilah and Al-Khateeb [9] showed that number of spike per square meter, 100-grain weight, grain weight/spike, and biological yield were the most effective variables influencing grain yield. The study of O. Alizadeh and A. Alizadeh [10] showed that using mycorrhizal inoculation caused better nutrient elements uptake in plant and generated plant growth regulators for increasing yield and yield components in corn. When *Azospirillum* and *Mycorrhiza* were used as biofertilizers in study conducted by Ardakani et al. [11], significant effects on wheat yield and its related traits were found. Statistical techniques have important roles in detecting relationship among plants’ traits and between the yield and its related components. One of the basic statistical methods to study the relationships between traits is simple correlation coefficient analysis [12]. Estimating correlation coefficients of different variables with grain yield is a suitable technique to decisions about the relative importance of these characters and their values as selection criteria [13]. On the other hand, the correlation coefficient may not give enough information about the relationship between different variables as much as statistical multivariate methods give. Therefore, other statistical techniques such as multiple regressions, factor analysis, principal component analysis, cluster analysis, and path coefficient analysis have been defined as appropriate techniques for interpreting these relationships in crop plants [13, 14]. The multivariate statistical analyses can provide more insights on the deep structure of data and traits’ relationships.

Attempts to identifying an ideal model for producing high-yielding wheat plants under different water regimes using mycorrhizal symbiosis in wheat breeding have rarely been made. This study was conducted to clarify and interpreting the relationship between wheat grain yield and its related components under drought conditions and mycorrhizal symbiosis with aiming to provide theoretical foundations guiding wheat breeders for researching the association of the main yield components and their influences on wheat plant productivity.

#### 2. Materials and Methods

##### 2.1. Experimental Procedures

The experiment was carried out in the greenhouse of Crop Production and Plant Breeding Department, College of Agriculture, Shiraz University, Shiraz, Iran, in 2010. Maximum, minimum, and the mean temperatures of this area were 20.9, 4.4, and 11°C, respectively, with 47.3% relative humidity and 602 mm annual rainfall. A factorial experiment based on completely randomized design with three replications was used. The studied factors were four water regimes (100%, 75%, 50%, and 25% of field capacity), three wheat cultivars (Darab 2 as a semi-resistant cultivar; Shiraz and Falat as sensitive cultivars), and mycorrhizal inoculation, including inoculated and noninoculated (control) treatments. The fungus used in this experiment was *Glomus intraradices* Schenck and Smith, provided by the Department of Soil Science, Shiraz University, Shiraz, Iran. Mycorrhizal inoculums were prepared through the trap culture in maize (*Zea mays* L.) with spores of *G. intraradices*. The mixture of trap culture medium was obtained from autoclaved soil/quartz sand (<1 mm) (4 : 1, v/v).

The soil samples used for planting were collected from Bajgah, Fars, Iran (1810 m asl; longitude 29°50′ and latitude 52°46′). The soil samples were air-dried, passed through 2 mm sieve, and mixed uniformly. The physicochemical properties of the soil were sandy loam, fine, mixed, mesic, calcixerollic xerochrepts, field capacity 25.3%, pH 7.9 (soil: distilled water, 1 : 1), electrical conductivity 0.5 d S m^{−1}, carbonate calcium equivalent 11.6%, total organic matter 1.34%, total Kjeldahl nitrogen 0.06%, Olsen phosphorus 15 mg kg^{−1}, extractable potassium 240 mg kg^{−1}, and DTPA-extractable Fe, Cu, Mn, and Zn were 5, 2, 11.3, and 1.7 mg kg^{−1}, respectively [15]. The day/night air temperature of the greenhouse was 25/15.

The pots were filled with 5 kg washed and sieved soil (mentioned previously) without purification or sterilization to simulate the real soil field properties. 150 mg N kg^{−1} soil (urea 46%) and some micronutrient elements such as Zn, Fe, Ca, and K up to 5 mg kg^{−1} were applied in each pot. The seeds were treated with ethanol 98% for about 20 s, washed three times with distilled water, and kept at 20°C for a week. About 5 cm of surface soil of each pot was removed and in mycorrhizal treatments, 50 g inoculums (containing spore numbers of 8 g^{−1} substrate and root colonization of 85 percent) were placed and incorporated with the remaining soil. Then, 3 cm of the removed soils was added to each pot; afterwards eight seeds were planted at equal distances in each pot. Finally, the rest of the removed soils were added to the pots. After germination, seedlings were thinned to four plants per pot.

Pots were daily weighed and based on the decreasing weight of each pot; decalcified water was added to each pot up to desired field capacity (FC) until the date of applying water regimes’ treatments. The water regimes were started at the tillering stage. The temperature during the experiment ranged from 15 to 28°C, with a 16/8 h light/dark period.

##### 2.2. Leaf Area Measurements

Plants’ leaves diameter were measured with a ruler, and leaves areas (X11) were calculated using the following equation [16, 17]:

##### 2.3. Total Chlorophyll Content

Total chlorophyll contents of the flag leaves (X8) were determined according to Iqbal et al. [6] procedure. Total chlorophyll content was extracted in 80% cold acetone, and the absorbance of the extractions was measured spectrophotometrically at 645 and 663 nm. Total chlorophyll was determined based on the following standard formula [18]: where, is the spectrophotometer reading at 645 and 663 wavelengths (nm), and Chl T is total chlorophyll content.

##### 2.4. Agronomic Traits

Grain yield (Y), tiller numbers (X1), spike length (X2), spikelet/spike (X3), spike weight/plant (X4), grains/spike (X5), grain weight/spike (X6), 100-grain weight (X7), biological yield (X9), and root weight (X10) were determined after the plants were harvested.

##### 2.5. Statistical Analyses

###### 2.5.1. Simple Correlation Coefficients

Pearson simple correlation coefficients were calculated, and the matrix of these correlations was studied [19].

###### 2.5.2. Multiple Linear Regressions

A multiple linear regression model was used for determining relative contribution of related components to the grain yield (Y) variations by applying the following equation [19]: where, is the dependent variable (yield), the ’s are independent variables (measured traits) affecting dependent one, is the intercept coefficient, and the ’s are the related coefficients of independent variables in predicting the dependent variable.

###### 2.5.3. Stepwise Regression

Stepwise regression [20] was used in order to determine the most important variables significantly contributed to total yield variability.

###### 2.5.4. Factor Analysis

The factor analysis method is consisted of the reduction of a large number of correlated variables to a much smaller number of uncorrelated variables. After extracting main factors, the matrix of factor loading was used to a varimax orthogonal rotation, and the communality or variance of uncorrelated variables was estimated by the highest correlation coefficient in each array as suggested by Seiller and Stafford [21].

###### 2.5.5. Principal Components Analysis

Principal components analysis is a mathematical procedure used to classify a large number of variables (items) into major components and determine their contribution to the total variation. The first principal component is accounted for the highest variability in the data, and each succeeding component accounts for the highest remaining variability as possible [22].

###### 2.5.6. Path Analysis

Path coefficient analysis was performed using simple Pearson correlation coefficients using grain yield/plant as dependent variable and the other characters as influential variables. The direct and indirect effects of influential variables on grain yield were calculated according to proposed method of Dewey and Lu [23].

###### 2.5.7. Cluster Analysis

Cluster analysis was used for arranging variables into different clusters to find the clusters that their cases within are more similar and correlated to one another comparing to other clusters. This procedure was performed using a measure of similarity levels and Euclidean distance [24, 25].

All statistical analyses were performed using SAS-9.1 [26] and Minitab-14 packages.

#### 3. Results

##### 3.1. Simple Correlations

Table 1 shows the minimum and maximum values, arithmetic mean, standard deviation, and standard error of means for all the estimated variables of wheat, separately for the inoculated and noninoculated plants. Results of simple correlation analysis (Table 2) show that all variables for either inoculated or noninoculated cultivars have a significant positive correlation with grain yield. The highest correlation with yield under both inoculated and noninoculated conditions was recorded for spike weight/plant ( for inoculated and 0.993 for noninoculated conditions). Overall, correlation coefficients with grain yield were higher in the inoculated plants than that of the noninoculated ones.

##### 3.2. Path Coefficients

The correlation coefficients were partitioned into direct and indirect effects in both inoculated and noninoculated plants (Tables 3 and 4). The highest direct effect (1.023 for the inoculation and 1.013 for the noninoculation) on yields for both conditions belonged to the spike weight/plant, while direct effects of the other variables were relatively low. The indirect effects of the spike weight/plant for both conditions were negative. Except for the spike weight/plant, the other variables had high indirect effects on the grain yield. Spikelets/spike showed the lowest direct contribution to the grain yield variations for both conditions but the highest indirect contribution through other variables.

##### 3.3. Multiple Linear Regressions

Regression coefficients and the probability of the estimated variables in predicting the wheat grain yield separately for the inoculated and noninoculated conditions are presented in Table 5. Based on these results, the predicting model equations for the grain yield/plant (Y) are formulated as follows.

*Model for Inoculated Condition*

where, Y is the grain yield, X1 is the tiller numbers, X2 is the spike length, X3 is the spikelet/spike, X4 is the spike weight/plant, X5 is the grains/spike, X6 is the grain weight/spike, X7 is the 100-grain weight, X8 is the total chlorophyll contents of the flag leaves, X9 is the biological yield, X10 is the root weight, and X11 is the leaf area.

*Model for Noninoculated Condition*

where, Y is grain yield, X1 is the tiller numbers, X2 is the spike length, X3 is the spikelet/spike, X4 is the spike weight/plant, X5 is the grains/spike, X6 is the grain weight/spike, X7 is the 100-grain weight, X8 is the total chlorophyll contents of the flag leaves, X9 is the biological yield, X10 is the root weight, and X11 is the leaf area. The formulas explained 99.5% and 99.2% () of the total variations of the grain yields for the inoculated and noninoculated conditions, respectively, and the remaining 0.5% and 0.8% are probably due to residual effects. The *t*-test for the variables revealed that the spike weight/plant and grain weight/spike contributed significantly in grain yields of both conditions, while the grains/spike, leaf area and root weight were significant only for the noninoculated condition.

##### 3.4. Stepwise Regression

Tables 6 and 7 show the entered or removed variables from the established model by stepwise regression separately from inoculated and noninoculated conditions. The partial and cumulative determination coefficient (), the probability value of entered variables to the model or removed variables from models, and standard error of the variables are also presented in these Tables (Tables 6 and 7). Under the inoculated condition, spike weight/plant, grains/spike, grain weight/spike, and total chlorophyll content of the flag leaf were entered into the model, and none were removed; while, under the noninoculated condition, seven variables, including spike weight/plant, grains/spike, grain weight/spike, total chlorophyll content of the flag leaf, biological yield, root weight, and the leaf area were entered into the model, and there were no variables to be removed. According to the results, 99.75% and 99.91% of the total variations in the grain yields were explained by the selected variables under the inoculated and noninoculated conditions, respectively. Due to their low relative contributions, the other variables were not included in the models. Therefore, based on the final step of stepwise regression analyses, the equations for prediction of grain yield (Y) were computed as follows.

*Final Model under Inoculated Condition*

where, Y is grain yield, X4 is the spike weight/plant, X5 is the grains/spike, X6 is the grain weight/spike, and X8 is the total chlorophyll contents of the flag leaves.

*Final Model under Noninoculated Condition*

where, Y is grain yield, X4 is the spike weight/plant, X5 is the grains/spike, X6 is the grain weight/spike, X8 is the total chlorophyll contents of the flag leaves, X9 is the biological yield, X10 is the root weight, and X11 is the leaf area.

##### 3.5. Factor Analysis

The first two of the twelve factors accounted for 59.2% and 60.1% of the total variations of the inoculated and noninoculated conditions, respectively (Table 8). The first factor was included for yield and spike weight/plant for both conditions, and it could explain 47.6% and 49.5% of the total variations in the dependent structure for the inoculated and noninoculated conditions, respectively; therefore, it can be named as grain yield factor. The second factor included spike length and biological yield for the inoculated condition which accounted for 11.6% of the total variability in the dependent structure, and it was considered as the spike length (Table 9). Under the noninoculated condition, the second factor accounted for 10.6% of the total variability and was consisted of total chlorophyll content of the flag leaf; therefore it could be named as total chlorophyll factor. The first two factors of the inoculated and noninoculated conditions are graphically depicted in Figures 1(a) and 1(b).

**(a)**

**(b)**

##### 3.6. Principal Component

The data presented in Table 10 and Figures 2(a) and 2(b) demonstrated that the increase in the number of components was associated with a decrease in eigenvalues. According to the results, the estimated wheat variables grouped into two main components that accounted for 72.4% and 70.3% of the total variations of the grain yields for the inoculated and noninoculated conditions, respectively. PC1 was moderately correlated with spike length (), spike weight/plant (), grains/spike (), biological yield (), and grain yield () in the inoculated plants and with spike weight/plant (), biological yield (), and yield () under the noninoculated condition. PC2 for the inoculated condition was moderately correlated with spikelets/spike () and total chlorophyll content (), while it was highly correlated with the leaf area (). In the noninoculated plants PC2 was moderately correlated with the spike length (), grain weight/spike (), total chlorophyll (), leaf area (), and yield (). PC1 accounted for about 63.5% and 60% of the variations in the grain yields; while PC2 explained 8.9% and 10.3% of the variations in the grain yields under the inoculated and noninoculated conditions, respectively. The first two components and their contributions in the variables for both conditions are graphically presented in Figures 3(a) and 3(b).

**(a)**

**(b)**

**(a)**

**(b)**

##### 3.7. Cluster Analysis

Hierarchical cluster analysis showed the similarity distance ranged between 72%–86% under the inoculated condition and 65%–82% under the noninoculated condition. The examined variables of wheat cultivars could be agglomerated into four and three clusters, respectively (Figures 4(a) and 4(b)). As the important variables affecting wheat grain yield, spike weight/plant, and 100-grain yield under the inoculated condition and the spike weight/plant, grain weight/spike, and the leaf area under the noninoculated condition were grouped together with the wheat grain yield into the third cluster.

**(a)**

**(b)**

Tables 11(a) and 11(b) outline the overall results of determining the most important variables affecting wheat grain yield for separately inoculated and noninoculated conditions.

#### 4. Discussion

Grain yield of wheat is the integration of many variables that affect plant growth throughout the growing period. Great efforts have been made to develop proper models that can predict wheat grain yield and distinguish the ideal- and high-yielding crop plants (ideotype). The knowledge of association and relationship between grain yield and its components under water deficit conditions would improve the efficiency of breeding programs by identifying appropriate indices to select wheat cultivars [27]. Because of wide variations in grain yield under normal and water stress conditions, simulating performance of wheat under soil moisture deficit represents special challenges for wheat modelers [28]. The results of the present study showed that spike weight/plant had the highest positive correlation with grain yield of wheat genotypes under both inoculated and noninoculated conditions. As well as spike weight/plant, other variables showed high positive correlations with yield. Based on the correlation coefficient analysis, all the variables had a high contribution with the wheat grain yield, but spike weight/plant was the most effective variable for both conditions. Moghaddam et al. [8] showed a negative correlation between the plant height and the grain yield. These investigators [8] attributed that increase in stem height caused the lower number of grains/spike. Kumbhar et al. [7] and Mohamed [29] had shown that grain weight/spike, biological yield, and the number of spikes were closely related to grain yield. In the study of Leilah and Al-Khateeb [9] grain yield had a high positive correlation with the number of spikes, number of grains/spike, 100-grain weight, weight of grains/spike, and the biological yield. In the study of Heidari et al. [1], grain yield was correlated positively with each of the biological yield, spikes/m^{2}, harvest index, and the grain weight/spike, while there was no correlation between the grain yield and the heading date or maturity. The differential relations of yield components to the grain yield may be attributed to environmental effects on plant growth [9, 30].

Since simple correlation coefficient can only determine the linear relationship between two related variables, but it cannot show how multiple variables are related to one another contributing to dependent variable (yield); path analysis was used for dividing the correlation coefficient of variables with yield into their direct and indirect effects via other variables. Spike weight/plant had the highest direct effect on grain yield under either condition. Path analysis shows that direct effect of biological yield, under inoculation condition, is positive, while, under noninoculated condition, it is negative, which is probably due to the effect of AM symbiosis on higher uptake of nutrient elements. The results also show that biological yield of the inoculated plants is significantly higher than those in noninoculated ones. Higher biological yield, under noninoculated condition, causes more consumption of plant energy and nutrient elements for more vegetative growth, resulting in decreasing grain yield, but inoculation of the mycorrhizal fungus compensated this problem by higher nutrient elements uptake from the soil. Based on the path coefficient analysis of wheat grain yield, Kumbhar et al. [7], Mohamed [29], and Leilah and Al-Khateeb [9] reported that the biological yield, harvest index, grains weight/spike, number of spikes/m^{2}, and 100-grain weight have the greatest impact in relation to the wheat grain yield. In the study of Heidari et al. [1], number of grains/spike had the largest direct and positive effect on the grain yield.

Linear regression analysis revealed that spike weight/plant and grain weight/spike are variables that significantly contributed in grain yield for both conditions, while the grains/spike, leaves area, and root weight were significant only for the noninoculated condition, so, these are variables determined as the most effective variables contributing to the grain yield by this statistical method. These results show that inoculation can affect relationship between the traits, and it can be useful for the breeding programs. Results of this study for regression analysis, under the noninoculated condition, are relative similar to the other studies such as Kumbhar et al. [7] and Leilah and Al-Khateeb [9] which showed that spike length, number of spikes/m^{2}, grain weights/spike, and the biological yield have contributed significantly to the grain yield. Asseng et al. [30] reported that increased kernel number had improved the potential yield of wheat under certain environmental conditions limited by water supply.

Stepwise regression analysis is a multiple statistical method that can screen or select the most important variables through a dependent variable such as the grain yield. Based on this method, spike weight/plant, number of grains/spike, grain weight/spike, and total chlorophyll content of the flag leaf are the most important variables contributing to the yield under the inoculated condition. Under the noninoculated condition, in addition to the four mentioned variables for the inoculated condition, three other variables, including the biological yield, root weight, and the leaves area are important. Relatively differences between the selected models by stepwise regression for the inoculated and noninoculated conditions show effectiveness of the root colonization on interrelationship of the variables. Leilah and Al-Khateeb [9] reported that the weight of grains/spike, harvest index, biological yield, number of spikes/m^{2}, and lastly spike length can explain 98.1% of the grain yield variations of wheat. Also, Mohamed [29] found that spike length, spike number, grain numbers/spike, spike weight, and straw yield are associated significantly with wheat grain yield. Heidari et al. [1] showed that the most important components for grain yield based on this method are grain weight/spike, spikes/m^{2}, and spikelets/spike.

Factor analysis showed that spike weight/plant under both inoculated and noninoculated conditions, spike length, and biological yield for the inoculated condition, and total chlorophyll content of the flag leaf for the noninoculated condition had the highest relative contribution to wheat grain yield (Table 9). Such results can be recognized by means of Figures 2(a) and 2(b). Similar results were obtained by Mohamed [29] who stated that factor analysis had classified the ten wheat variables into two main groups which accounted for 80.79% of the total variability in the dependence structure. Leilah and Al-Khateeb [9] showed that biological yield, harvest index, weight of grains/spike, spike length, and number of spikes/m^{2} had a high relative contribution to wheat grain yield.

Results of the principal component analysis revealed that, under the inoculated condition, spike length, spike weight/plant, grains/spike, biological yield, number of spikelets/spike, and total chlorophyll are important variables affecting grain yield, and under the noninoculated condition, spike length, spike weight/plant, biological yield, grain weight/spike, total chlorophyll, and leaves area are the most important factors. Harvest index, biological yield, spike diameter, number of spikes/m^{2}, spike length, grain weights/spike, and 100-grain weight in the study of Leilah and Al-Khateeb [9] were the most important factors in contributing to the yield. Also, Yin et al. [31] stated that the grain yield was divided into three components, namely, number of spikes/m^{2}, number of kernels/spike, and 1000-kernel weight.

Hierarchical cluster analysis showed that, under the inoculated condition, grain yield, spike weight/plant, and 100-grain yield are highly important variables for contributing to yield, while, under the noninoculated condition grain yield, spike weight/plant, grain weight/spike, and leaves area are highly important variables. Hierarchical cluster analysis method starts with the calculation of the distance of each variable in relation to other variables. Groups are then formed by the process of agglomeration division. In this process, all variables start individually in groups of one. The close groups are then gradually merged until finally all variables come to a single group. Repeated splitting of the groups will result in all the evaluated variables being in groups of their own. In the study of Leilah and Al-Khateeb [9] results proved that 100-grain weight, weight of grains/spike, harvest index, and the biological yield were the variables most closely related to the grain yield.

It has been well established that AM symbiosis protects host plants against negative effects of drought stress due to nutritional, physical, and cellular improvements [32]. In addition, the AM symbiosis increases host plant growth due to improved plant nutrient and water uptakes via external hyphae in inoculated roots [33]. The beneficial effects of different mycorrhizal fungi on plant growth, under drought conditions, have been demonstrated in wheat [34] and other plant species [35, 36]. Mouchesi et al. [37] showed higher production of yield and its related traits under inoculated condition in compare to noninoculated one due to higher uptake of nutrient elements and water by mycorrhizal roots.

#### 5. Conclusions

The multiple statistical procedures which have been used in this study showed that, under water stress condition and mycorrhizal inoculation, spike weight/plant and total chlorophyll content of the flag leaf are the most important variables contributing to wheat grain yield, while for the noninoculated condition, grain weight/spike and chlorophyll content of the flag leaf, grain weight/plant, and leaves area are also important. Therefore, spike weight per plant and chlorophyll content of the flag leaf can be used as selection criteria in the breeding programs for both the inoculated and the noninoculated wheat cultivars under different water regimes, and also grain weight per spike and leaf area can be considered for the noninoculated condition. Furthermore, the results indicated that root inoculation with the mycorrhizal fungus changed the relationship among morphological traits of wheat cultivars under drought stress.

Overall, the results indicate that mycorrhizal symbiosis can relatively change the impact of important yield-related traits through production of wheat grain yield which is due to the effect of mycorrhizal symbiosis on water and nutrient elements uptakes by root. Also, plant breeders can consider fewer variables as selection criteria under mycorrhizal symbiosis than nonmycorrhizal condition to find higher yielding cultivars in their breeding programs.

Also, using different statistical techniques for determining important traits contributing to grain yield of wheat showed that simple correlation cannot distinguish important variables. On the other hand, since path analysis uses results of simple correlation, it is not suitable for using it in selecting important variables. Similar to the results of correlation and path analysis, cluster analysis and multiple regression analysis could not clearly distinguish important traits either. It seems that the results of stepwise regression as a selecting method together with principal component and factor analysis are stronger statistical methods to be applied in breeding programs for screening important traits.